3
_f @ s d Z dd ZdS )a
**********************************************
*** Program for the computation of the ***
*** line profile ***
**********************************************
This code is part of the LIFELINE program.
Copyright (C) 2020-2021 University of Liege (Belgium)
Enmanuelle Mossoux (STAR Institute)
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public
License as published by the Free Software Foundation, either version 3 of the License, or any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied
warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
details.
You should have received a copy of the GNU General Public License along with this program.
If not, see .
Run: from compute_profile import profile
p=profile(direct, direct2, ipar, nbr_bin_profile, q1, T, cs_interp, T_cs, energy, M, add, add2)
Read the files created by [r/a]shock[_coriolis].py in direct
# Parameters
# ==========
direct - directory where the shock characteritics are saved
direct2 - directory where to save files
ipar - index of the binary system for which we read the ray-tracing files
nbr_bin_profile - number of bins in the line profile
q1 - interpolation function of the emissivity
T - tabulated temperatures used to define q1
cs_interp - interpolation function of the cross-section
T_cs - tabulated temperatures used to define cs_interp
energy - rest energy of the line
M - orbital phase
add[2] - additional name when reading the set of binary systems
# Output
# ======
p - numpy array containing the emission in each bin
# Versions
# ========
v1 - 21/10/19
c
, C sp dd l }
dd l}ddlm} dd lj}t| d t| | d d}xB|j j }|}|s\P yt
|d P W qF t k
r Y qFX qFW t| d t| | d d}xB|j j }|}|sP yt
|d P W q t k
r Y qX qW t| d t| | d d}xH|j j }|}|s"P yt
|d P W n t k
rJ Y nX q
W g }dt| | }t|| d d }|j d
t| d |j d t|
d |j d
g }xt
t|D ]} |jd | dkr|jt
|d |} |} n*|jt
|j |j j } |j j } dd | D } dd | D } t| d t| d t| | d d} xD| j j } | s~P yt
| d P W n t k
r Y nX qjW xt
tt| D ]v}!| |! dkr| |! d |d k s| |! d |d% krqn||j| |! d }"|!dkr4| j j } | dd d }#| dd d }$dd |#D }#dd |$D }$d}%xvt
tt|#D ]b}&|#|& t|k r|t|}'n,|#|& t|kr|t|}'n ||#|& }'|%|'|$|& }%qW |
j| |! |
j|% d|" r|| | |! |
j|% d|" || < qW | j |j t|d& d t|| d qW |j |j |j |j |jd d| }(|(|
j| d |d |( })t|dkrtd |)d }*d }+n$||
j| }*t||
j| d }+|jd|) |* |j| | gdd gdd! |j d|+ |j d" |j d# |j || d$ |j |S )'N ) constantez/bin_parz.datarz/emiss_part_shock_parz/temp_emiss_parz/line_profile_parwz# Energy (keV) =
z5# Phases (mean anomaly, radian, M=0 -> conjunction): z8# tangential velocity (km/s) | luminosity (10^27 erg/s)
c S s g | ] }t |qS )float).0xr r 8/local/network/htdocs/Lifeline/codes3/compute_profile.py
h s zprofile..c S s g | ] }t |qS r )r )r r r r r
r i s z/ray_tracing_par_binghh!w> c S s g | ] }t |qS r )r )r r r r r
r s c S s g | ] }t |qS r )r )r r r r r
r s g g $@ g0O(@g @@czNo X-ray emissiong ?g?)r z Energy (keV)zNormalized emissionz.pdfr ) numpymath
constantesr matplotlib.pyplotpyplotopenstrreadlinesplitr
ValueErrorwriterangeintappendlenlog10minmaxisfiniteexpclosefigurearrayprintnansumplotylimxlabelylabelsavefig),directdirect2Ziparuseriparnbr_bin_profileq1T cs_interpT_csenergyMaddadd2npr r pltfbinlineZbin1fhistoZhisto1ftemptemp1ZvtangZ nom_debutZfprofileZemissZibinZfacteurZtemperfRTZligneitempqZtemperRTZfactRTtauZitemp2cssZ
wavelengthwaveZ y_to_plotZy_limitr r r
profile1 s
* ,
&0,
rI N)__doc__rI r r r r
- s