3
ëÚ_; ã @ sz d Z dd„ Zdd„ Zdd„ Zdd„ Zd d
„ Zd d „ Zd
d„ Zdd„ Zdd„ Z dd„ Z
dd„ Z d dd„Z dd„ Z
d d „ Zd S ) a¿
***************************************************
*** Program for the computation of the shock ***
*** in radiative colliding wind binaries ***
*** with Coriolis deflection ***
***************************************************
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 rshock_coriolis import rshock
rien=rshock(Mdot1, Mdot2, vinf1, vinf2, mass_ratio, a, per, R1, R2, beta1, beta2, Teff_1, Teff_2, nbr_points_shock_2D, nbr_points_shock_3D, nstep_width, nbr_bin_profile, direct, cut_lim, wind_prim, wind_sec, ex, omega, phase, incl, ipar, mu, T)
Use the Antokhin et al. (2004) formalism for the shock creation and the Parkin & Pittard (2008) formalism for the Coriolis deflection
# Parameters
# ==========
Mdoti - Mass loss rate [Msun/yr]
vinfi - terminal velocity of the winds [km/s]
mass_ratio - mass ratio
a - semi-major axis [Rsun]
per - orbital period [days]
Ri - stellar radius [Rsun]
betai - index of the beta velocity law
Teff_i - efective temperature [K]
nbr_points_shock_2D - number of points to create the 2D shape of the shock
nbr_points_shock_3D - number of points to create the 3D shape of the shock
nstep_width - number of points to discretize the inside of the shock
direct - directory where to save files
cut_lim - Limit of the tangential velocity to see the emission
wind_prim - h5 file containing the wind distribution of the primary
wind_sec - h5 file containing the wind distribution of the secondary
ex - eccentricity
omega - argument of the periapsis [radian]
phase - phase of the system
incl - inclination of the orbit
ipar - index on the set of stellar parameter
mu - Mean molecular weight
T - Temperatures of emissivities
# Versions
# ========
v1 - 16/03/18
v2 - 26/03/19 - Taking the radiative inhibition into account
c C sò dd l } | d |d |d k r|d} |jd|| || ggƒ\}
}|d }| j|
j| d d |
j| d d ƒ}|
j| d dkoš|
j| d dkr||d }|d }d}|
j| d }|
j| d }|d }xz|dkoæ|dkrP|dkr"|
j| d }|
j| d }|d }|}d }qØd}|
j| d }|
j| d }|d }|}qØW | j|
j| d d |
j| d d ƒ}|| d |d |d k röd}|j|| || ggƒ\}
}|d }| j| j| d d | j| d d ƒ} | j| d dkrö| j| d dkrö|d }|d }d}| j| d }| j| d }|d }x~|dkrÊ|dkrÊ|dkrš| j| d }| j| d }|d }|}d }n,d}| j| d }| j| d }|d }|}qNW | j| j| d d | j| d d ƒ} | d |d |d krâ|| d |d |d krâ|j| | || ggƒ\}
}|d }| j| j| d d | j| d d ƒ} | j| d dkr†| j| d dkr†|d }|d }d}| j| d }| j| d }|d }x~|dkrZ|dkrZ|dkr*| j| d }| j| d }|d }|}d }n,d}| j| d }| j| d }|d }|}qÞW | j| j| d d | j| d d ƒ} |jd| | || ggƒ\}
}|d }| j|
j| d d |
j| d d ƒ}|
j| d dkrâ|
j| d dkrâ|d }|d }d}|
j| d }|
j| d }|d }x~|dkr¶|dkr¶|dkr†|
j| d }|
j| d }|d }|}d }n,d}|
j| d }|
j| d }|d }|}q:W | j|
j| d d |
j| d d ƒ}|| || S )
Né é g ÐcAg ð?ÚuxÚuyé ÚupÚdown)ÚnumpyÚqueryÚsqrtÚiloc)ÚxÚyÚR1ÚR2ÚMdot1ÚMdot2ÚdÚtree_prim_windÚ
tree_sec_windÚ
dwind_primÚ dwind_secÚnpÚvv1ÚrienÚ ind_closeÚvv2Ú
ind_closepÚ
ind_closemÚstadeÚdw1Údw2© r! ú8/local/network/htdocs/Lifeline/codes3/rshock_coriolis.pyÚcomp_eta9 sÆ *&
* *(
*8 *(
* *(
*r# c C s: dd l } t| d||||||| |
| ƒ }
| |d| j|
ƒ S )Nr g g ð?)Úmathr# r
)r r r r Úv1Úv2r r r r r r r$ Úetar! r! r" Úfind_x0§ s r( c C sb dd l }dd l}y<||j|ƒ | d|j ||j|ƒ d |j|ƒ }W n d}Y nX |S )Nr g @r g )r r$ ÚsinÚpi) r ÚyyZddeltaZaangleZsslopeZMMdotr r$ Úretr! r! r" Ú ddzeta_dy® s <
r- c C sÌ ddl m} dd l}dd l}| d| kr~|j|ƒ|j|ƒ dkr\|t|
|| g||||fd}|| | |j|ƒ }d }
d
}nD|dkrd}| || d ||| d }d |jd|jd| ƒ ƒ |jd| ƒ|| || |jd| ƒ || d|jd| ƒ d
|| d |jd| ƒ |||jd| ƒ d|jd| ƒ }|dk rbd}|
|| d d||| } | d| | d dd || | d
|| d }| dkr¨|
| d d|| } | d| | d dd | | d
|d }n˜| |d || d }| || d ||| d }d |jd|jd| ƒ ƒ |jd| ƒ|| || |jd| ƒ || d|jd| ƒ d
|| d |jd| ƒ |||jd| ƒ d|jd| ƒ }|dk rþd}|
| d d|| } | d| | d dd || | d
|| d }| dkr¨|
|| d d||| } | d| | d dd | | d
|d }|| | |j|ƒ }
| }||
|gS )Nr )Úodeintgš™™™™™É?g )Úargsg ð?Úyesr g @g @g @g ð¿g ð¿)Úscipy.integrater. r r$ r) r- Úcosr
) Úy_newZz_newÚxstagÚy_vecZ anglesignZ deltasignZ slopesignZMdotZvsignZetasignZrhosignÚbetaZ cote_chocZ dzeta_oldÚy_oldr r r ÚradiusÚcrashingr. r r$ ZdzetaÚsigmaÚc1Úc2Zz0Zsigma_0r! r! r" Úsigmas¹ sB
¨
8
4 ¨
8
0r= c ! C sò dd l }dd l}ddlm} dd l}
ddlm} dd l} d| jd< t |d dƒ}
|
j
d| d ƒ |
j
d
ƒ |
j ƒ |
j j
|d }|j|| d kƒd }| d t|ƒk r®dg}| d t|ƒkrÐtt|ƒƒd
g}||d d }d}d}| |d}|j|ƒ}g }xr gQ²@³-(;é g )ÚrefgêWôTáßDr gfi´íÍ ð?g¾´dbÒùA?Fg-<¥Ò}#ð?g ð?é éÿÿÿÿrF rF ) r r$ Ú sun_abundr> ÚpyatomdbÚ
constantesr? ÚosÚenvironÚopenÚwriteÚcloseÚutilÚload_user_prefsÚwhereÚminÚmaxÚintÚlenÚarrayÚrangeÚappendÚatomicÚ Z_to_massÚsumÚargminÚabsÚtakeÚarangeÚisnanÚfloor)!ÚkTÚkT_tabÚcool_tabÚdirectÚkT_ionÚion_fracÚsunabundr r$ r> rH r? rJ ÚfuserÚsettingÚindÚ
cool_tab_hereÚnzÚsum_vecÚnumber_fractionÚ
mass_fractionÚzzÚZÚfrac_totÚ
ion_frac_hereÚZ_vecÚXhÚ
atomic_weightÚAzÚAznÚXzÚmuÚmubarr! r! r" Ú integrand_Tä sd
$
2& r} c ! C sî dd l }dd l}ddlm} dd l}
ddlm} dd l} d| jd< t |d dƒ}
|
j
d| d ƒ |
j
d
ƒ |
j ƒ |
j j
|d }|j|| d kƒd }| d t|ƒk r®dg}| d t|ƒkrÐtt|ƒƒd
g}||d d }d}d}| |d}|j|ƒ}g }x )r? z'https://hea-www.cfa.harvard.edu/AtomDB/r@ z /userdatarA zfilemap=z /filemap
z3atomdbroot=https://hea-www.cfa.harvard.edu/AtomDB/
)rB ghÉh÷Æ!w>r gQ²@³-(;rC g )rD gêWôTáßDr gfi´íÍ ð?g¾´dbÒùA?Fg-<¥Ò}#ð?g ð?rE rF rF rF ) r r$ rG r> rH rI r? rJ rK rL rM rN rO rP rQ rR rS rT rU rV rW rX rY rZ r[ r\ r] r^ r_ r` ra )!rb rc rd re rf rg rh r r$ r> rH r? rJ ri rj rk rl rm rn ro rp rq rr rs rt ru rv rw rx ry rz r{ r| r! r! r" Úvar_rho! sd
$
2& r~ c C s& dd l }|j|d || d ƒ| S )Nr r )r r
)ZcenterxZcenteryr8 r r r! r! r" Ú plot_star^ s r c C s: dd l }|}| }|j|d |d ƒ}|j|| || ƒS )Nr r )r r
Úarctan2)r
Úzr ÚoppÚadjÚhypor! r! r" Úangleyzb s
r… c C s dd l }| |j|ƒ S )Nr )r$ r2 )ÚvecÚangler$ r! r! r" Úprojxi s rˆ c C s dd l }| |j|ƒ |j|ƒ S )Nr )r$ r) r2 )r† Úangle1Úangle2r$ r! r! r" Úprojyl s r‹ c C s dd l }| |j|ƒ |j|ƒ S )Nr )r$ r) )r† r‰ rŠ r$ r! r! r" Úprojzo s rŒ çš™™™™™Ù?c C sx dd l }| j}|| |j|| ddd }|| |j|| ddd }|j|j|d |d ƒdd}| j||kƒ d S )Nr r )Úaxisr )r Ú trianglesÚrollrS r
Úset_mask) Útriangr r
Úalphar r ZxtriZytriZmaxir! r! r" Ú
apply_masks s r” c
C s† dd l }| \}}}|\}}} g }
g } xXttt|ƒƒƒD ]D} | j|j|| ||f||fƒƒ |
j|j|| ||f|| fƒƒ q6W | |
fS )Nr )r rW rT rU rX Úinterp)
Úp1Úp2r r Úx1Úy1Úz1Úx2Úy2Úz2ZreszZresyÚir! r! r" Úfind_xy~ s
$rŸ c 9 C sâ ddl m} ddlm} ddlm} dd l}dd l}ddlm } dd l
j }dd l j
}ddlm} dd l}
dd l} dd l} dd l}
ddlm} ddlm} | \ }}}}}}}}}}}}} } } } } }!}"}#}$}%}&}'}(})}*}+},}-}.}/|
jd ƒ |j}0|#d
t|.ƒ t|/ƒ d }1| j|%ƒj d d
}2| j|&ƒj d d
}3||j |j |2j d d … d |2j d d … d fƒd–t t!|2j d d … d ƒƒfƒj"dd}4||j |j |3j d d … d |3j d d … d fƒd—t t!|3j d d … d ƒƒfƒj"dd}5|| || }6d}7d}8|dƒ}9|dƒ}:|dƒ};|dƒd }<|dƒ}=||: }||: }||: }||9 d˜ }||9 d™ }d |0d |d |<|d d }>|>d | }?t#||gƒ}@|j$|"|j%d!d }A|A|j&t%|"ƒd" ƒd |@ |j&t%|"ƒd" ƒ }Bg }Cg }Dt'|1ƒT}ExL|ED ]D}F|Fj(ƒ d d#krä|Fj(ƒ }G|Dj)t%|Gd ƒƒ |Cj)t%|Gd ƒƒ qäW W d Q R X |j*|Dƒd$ }D|j+|D|j*|Cƒd% ƒ |j,d&ƒ |j-d&ƒ |j.d'ƒ |j/d(ƒ |j0|#d) t|.ƒ t|/ƒ d* ƒ |j1ƒ |0d" }H|6d krp|0d d |6 }I|0d" }Hd |0d" |I |j2|0d" ƒ }J|J|I|H }K|H|K }HxLt3|Kƒd+krTt3|Jƒd+krTd |H|I |j2|Hƒ }J|J|I|H }K|H|K }Hq
W |H|0d" k rp|0d" |H }Hg }Lg }Mg }Nt'|#d, ƒ˜}Ox|OD ]ˆ}F|Fj(ƒ }P|Pd d d-kr|Lj)t%|Pd ƒƒ |Mj)t%|Pd ƒƒ g }Qx4t4t |Pd ƒd ƒD ] }R|Qj)t%|P|Rd ƒƒ qìW |Nj)|Qƒ qW W d Q R X |j*|Lƒ}Lt#||gƒd.|,d d$ d/ d0 dš|+d/ d2 k rnt5d3ƒ d4}Sd›|0 }T|'}U|'}V|V}W|V|W ||j6|Wƒ d ||j7|Wƒ }XxDt3|Xƒd+krò|W|X }W|V|W ||j6|Wƒ d ||j7|Wƒ }Xq°W |W||j6|Wƒ | d"|0 }Y|j7|Wƒ| d ||j7|Wƒ }Z|j8d || ƒ|j6|Wƒ d ||j7|Wƒ }[|j9|[|Zƒ}\|\d"|0 kr~d6}\|d |d d ||j7|\ƒ }]|]}^|j:||] ƒ}_|j;| |^ ||^ d7ƒ}`td9„ t4d?ƒD ƒ\}Š}‹}Œ}d@d9„ t4dAƒD ƒ\}Ž}}}‘}’}“|U}”|j$| ƒdB |0 dC| }•dB|0 dC| }–|j$d"|! ƒd6 }—|j$d"|! ƒd6 }˜|V}™|'}šdD}›d6}œCxXt4|Sd ƒD C]F}|U|T| }V|V}W|V|W ||j6|Wƒ d ||j7|Wƒ }XxDt3|Xƒd+kr¼|W|X }W|V|W ||j6|Wƒ d ||j7|Wƒ }XqzW |W||j6|Wƒ | d"|0 }žt3|ž|Y ƒ}Ÿ|j7|Wƒ| d ||j7|Wƒ }Z|j8d || ƒ|j6|Wƒ d ||j7|Wƒ }[|j9|[|Zƒ}\|\d"|0 k rTd6}\|d |d d ||j7|\ƒ }]|j8|<|> d"|] d | ƒ} |]|? |> }¡|t=||] | d" |]|||||||4|5|2|3f dE}¢|¢d }¢|dk
r>t'|#dF t|)ƒ dGƒ}£|£j>t|¢ƒdH ƒ |£j>t|¢| ƒdI ƒ |£j>t|¢|] ƒdJ ƒ |£j1ƒ t#|¢|gƒ}¢|4j?|¢|] d6ggƒ\}¤}¥|¥d }¥|j8|2j |¥ dK d |2j |¥ dL d ƒ}¦|2j |¥ dK dk rž|2j |¥ dL dk rž|¥d }§|¥d }¨dM}©|2j |§ dK }ª|2j |§ dL }«|¥d }§x~|ªdk rr|«dk rr|©dMk rB|2j |¨ dK }ª|2j |¨ dL }«|¨d }¨|¨}¥dN}©n,dM}©|2j |§ dK }ª|2j |§ dL }«|§d }§|§}¥
qöW |j8|2j |¥ dK d |2j |¥ dL d ƒ}¦|5j?d |¢|] d6ggƒ\}¤}¥|¥d }¥|j8|3j |¥ dK d |3j |¥ dL d ƒ}¬|3j |¥ dK dk rö|3j |¥ dL dk rö|¥d }§|¥d }¨dM}©|3j |§ dK }ª|3j |§ dL }«|¥d }§x~|ªdk rÊ|«dk rÊ|©dMk rš|3j |¨ dK }ª|3j |¨ dL }«|¨d }¨|¨}¥dN}©n,dM}©|3j |§ dK }ª|3j |§ dL }«|§d }§|§}¥ qNW |j8|3j |¥ dK d |3j |¥ dL d ƒ}¬|j:| |¦ ƒ}|¢|k
s |¢|]| k
rH|dk
rBdO}›t5dPƒ t5dQƒ t5dPƒ d6}n||j|H k
r\d6}|j$| ƒ|H | dR|0 d" }®|j$| dS ƒd |0|H | dS d }¯|¯|j2|¯ƒ }°d d |6 |®|j2|®ƒ d }±g }²xH|±D ]@}³|j@t3|°|³ ƒ|jAt3|°|³ ƒƒkƒd }´|²j)|¯|´d ƒ
qÖW |j*|²ƒ}²|j6|²ƒ|j6|²|® ƒ }µ|µ|j6|®ƒ |] }¶|j@|3d d/kƒd }·|¢}¸d6}¹g }ºg }»xÈ|¶D ]À}¼|j8|¸d |¹d ƒ}½|j8|]|¸ d |¹d ƒ}¾tB|¸|¹|||||]|4|5|2|3ƒ }¿|¿dkr |¸|]|½d |j8|¿ƒ |½d |j8|¿ƒ |¾d |¼ }À|À|¼|¹ |¸ }¸|¼}¹|ºj)|¸|] ƒ |»j)|¿ƒ qzW |j*|»ƒ}»|j*|ºƒ}º|j8|ºd |¶|] d ƒ}Á|¶|] |Á }Â|º|Á }Ã|jC|Â|Ã}®|ºdd … |ºd dœ… }Ä|¶dd … |] |¶d d… |] }Å|j8|Åd |Äd ƒ}Æ|jC|Å|Æ |Ä|Æ ƒ}Ç|j*|Çd gƒ}È|¢|kr|j*|0d" gƒ}È|j |È|Çfƒ}Ç|Ç|® }É|0d" |Éd< |j@d d||] |j8|ºd |¶|] d ƒ | |j7|Ƀ |$k ƒd }Êd d6d6g}Ë|º|] }º|dk-r2 x’t4|Êdž d ƒD ]|}Ìt5t|̃dT t|ÊdŸ ƒ ƒ |º|Ì d6k r4|j:|¶|Ì |]|º|Ì ƒ|_k r4|j2|_ƒ|]|º|Ì }Ít3|̓t%|jD|¶|Ì ƒƒ |¶|Ì< |jE|jF|º|Ì |¶|Ì d6g|˃|jGjH|º|Ì |¶|Ì d6gddU|jGjH|ËddU ƒ}Î|0|jE|jF|º|Ì |] |¶|Ì d6g|˃|jGjH|º|Ì |] |¶|Ì d6gddU|jGjH|ËddU ƒ }Ï|j8|º|Ì d |¶|Ì d ƒ}Ð|j8|]|º|Ì d |¶|Ì d ƒ}Ñd6}Òd6}Ó|0d" }Ô|Ìdkr||jE|jF|º|Ìd |¶|Ìd d6g|˃|jGjH|º|Ìd |¶|Ìd d6gddU|jGjH|ËddU ƒ}Ò|0|jE|jF|º|Ìd |] |¶|Ìd d6g|˃|jGjH|º|Ìd |] |¶|Ìd d6gddU|jGjH|ËddU ƒ }Ó|jE|jF|º|Ì |º|Ìd |¶|Ì |¶|Ìd d6g|˃|jGjH|º|Ì |º|Ìd |¶|Ì |¶|Ìd d6gddU|jGjH|ËddU ƒ}Ôt3|Î|Ò ƒ}Õt3|Ï|Ó ƒ}Öd6}×|j$|!ƒd6 }Ø|Ð|krd|j;d6|Ð|] ||] d dVƒ}Ùg }×|j@|2d d |2d d d/ |Ðk ƒd }·|2j |· }Ú|j@|Úd d |Úd d d/ ||] kƒd }·|Új |· }Ú|j@|Úd |¶|Ì k ƒd }·|Új |· }Úx’t4dVƒD ]„}ÛtI|Ð|Ù|Û |] |΃}ÜtJ|Ð|Ù|Û |] |Îd6ƒ}Ý|4j?|Ü|] |Ý|] ggƒ\}¤}¥|¥d }¥t#||j8|2j |¥ dK d |2j |¥ dL d ƒƒ}Þ|2j |¥ dK dkr|2j |¥ dL dkr|¥d }§|¥d }¨dM}©|2j |§ dK }ª|2j |§ dL }«|¥d }§x~|ªdkrÔ|«dkrÔ|©dMkr¤|2j |¨ dK }ª|2j |¨ dL }«|¨d }¨|¨}¥dN}©n,dM}©|2j |§ dK }ª|2j |§ dL }«|§d }§|§}¥qXW t#||j8|2j |¥ dK d |2j |¥ dL d ƒƒ}Þd1|+d. |Þ|j6|Ô|Î ƒ d0 d }ß|j@|D|ßkƒd }´|ßtK|Dƒk rRdg}´|ßt#|Dƒkrrt t!|Dƒƒd g}´|C|´d }àd | |Þd |0 |Ð|Ù|Û |] d }á|×j)dW|8d |;d d t3|Þ|j6|Ô|Î ƒ ƒd |]dX |7 |á |à ƒ qnW |j@t3|×|Ù ƒtKt3|×|Ù ƒƒkƒd }´|×|´d |] }×|j$|!ƒd"|! d |× d"|! |!d |×d"|! }Ød6}â|j$|!ƒd6 }ãt3|j6|0|Ô |Ï ƒƒd+kr`|j;d6|Ñ|] ||] d dVƒ}Ùg }â|j@|3d d |3d d d/ |Ñk ƒd }·|3j |· }ä|j@|äd d |äd d d/ ||] kƒd }·|äj |· }ä|j@|äd |¶|Ì k ƒd }·|äj |· }äx¬t4dVƒD ]ž}ÛtI|Ñ|Ù|Û |] |0|Ï ƒ|] }ÜtJ|Ñ|Ù|Û |] |0|Ï d6ƒ}Ý|5j?|]|Ü |] |Ý|] ggƒ\}¤}¥|¥d }¥|j8|äj |¥ dK d |äj |¥ dL d ƒ}Þ|3j |¥ dK dkrì|3j |¥ dL dkrì|¥d }§|¥d }¨dM}©|3j |§ dK }ª|3j |§ dL }«|¥d }§x~|ªdkrÀ|«dkrÀ|©dMkr|3j |¨ dK }ª|3j |¨ dL }«|¨d }¨|¨}¥dN}©n,dM}©|3j |§ dK }ª|3j |§ dL }«|§d }§|§}¥qDW |j8|3j |¥ dK d |3j |¥ dL d ƒ}Þ|Þ|krú|}Þd1|+d. |Þ|j6|0|Ô |Ï ƒ d0 d }ß|j@|D|ßkƒd }´|ßtK|Dƒk rJdg}´|ßt#|Dƒkrjt t!|Dƒƒd g}´|C|´d¡ }àd | |Þd |0 |Ñ|Ù|Û |] d }á|âj)dW|8d |;d d t3|Þ|j6|0|Ô |Ï ƒ ƒd |]dX |7 |á |à ƒ qPW |j@t3|â|Ù ƒtKt3|â|Ù ƒƒkƒd }´|â|´d |] }â|j$|!ƒd"|! d |â d"|! |!d |âd"|! }ã|ÌdkrŽ|Ð|× |åk r||æ}×|Ñ|â |çk rŽ|è}ât#|œ|×|âƒ}œtI|Ð|Ù|Û |΃}ÜtJ|Ð|Ù|Û |Îd6ƒ}Ý|4j?|Ü|] |Ý|] ggƒ\}¤}¥|¥d }¥tK||j8|2j |¥ dK d |2j |¥ dL d ƒƒ}é|2j |¥ dK dk r$|2j |¥ dL dk r$|¥d }§|¥d }¨dM}©|2j |§ dK }ª|2j |§ dL }«|¥d }§x~|ªdk rò|«dk rò|©dMk rÂ|2j |¨ dK }ª|2j |¨ dL }«|¨d }¨|¨}¥dN}©n,dM}©|2j |§ dK }ª|2j |§ dL }«|§d }§|§}¥ qvW tK||j8|2j |¥ dK d |2j |¥ dL d ƒƒ}ét#dY|éƒ}étI|Ñ|Ù|Û |0|Ï ƒ|] }ÜtJ|Ñ|Ù|Û |0|Ï d6ƒ}Ý|5j?|]|Ü |] |Ý|] ggƒ\}¤}¥|¥d }¥tK||j8|3j |¥ dK d |3j |¥ dL d ƒƒ}ê|3j |¥ dK dk rÈ|3j |¥ dL dk rÈ|¥d }§|¥d }¨dM}©|3j |§ dK }ª|3j |§ dL }«|¥d }§x~|ªdk r–|«dk r–|©dMk rf|3j |¨ dK }ª|3j |¨ dL }«|¨d }¨|¨}¥dN}©n,dM}©|3j |§ dK }ª|3j |§ dL }«|§d }§|§}¥ qW tK||j8|3j |¥ dK d |3j |¥ dL d ƒƒ}êt#dY|êƒ}ê|j$|!ƒ}ë|j |Ž|×|Ð|×|ë |! d |j6|0d" |Î ƒ |Õ |– t3|j6|Ô|Î ƒƒ|! fƒ}Ž|j |Ž|â|Ð|â|ë |! d |j6|0d" |Ï ƒ |Ö |– t3|j6|Ô|Ï ƒƒ|! fƒ}ŽtI|Ð|× |΃}ìtJ|Ð|× |Îd6ƒ}ítI|Ñ|â |0|Ï ƒ|] }îtJ|Ñ|â |0|Ï d6ƒ}ïtI|Ð|Ø |΃}ðtJ|Ð|Ø |Îd6ƒ}ñtI|Ñ|ã |0|Ï ƒ|] }òtJ|Ñ|ã |0|Ï d6ƒ}ó|0d" }Ôd6}ôd6}õ|Ìdk!r|jE|jF|ì|ö |í|÷ d6g|˃|jGjH|ì|ö |í|÷ d6gddU|jGjH|ËddU ƒ}Ô|é|j7|Ô|Î ƒ }ô|Ì|Êd¢ k rœ|Ô}ø|jE|jF|î|ù |ï|ú d6g|˃|jGjH|î|ù |ï|ú d6gddU|jGjH|ËddU ƒ}Ô|ê|j7|0|Ô |Ï ƒ }õ|Ì|Êd£ k!r|Ô}û|ì}ö|í}÷|î}ù|ï}ú|Ð|× }å|Ñ|â }ç|×}æ|â}è|j8|éd |ôd ƒ}ü|j8|êd |õd ƒ}ý|j |‘|ôgfƒ}‘|j |’|õgfƒ}’d6}þ|j$|!ƒd6 }ÿ|j$|!ƒd6 } d | |éd |0 |Ðd |j$|!ƒd6 }|×d6k$râd1|+d. |üd0 d }þ|j;|×d"|! |×|!ƒ}ÿd | |éd |0 |Ð|× d }|j;d |þd$ |!ƒ}d d¤|; d |éd | d[|=d }g }x8|D ].}|j)d\tL||D|C|#|Ld$ |N|*ƒ ƒ "q|W ||jM|ƒ|jN|dd … |d d¥… ƒ |] }||× t#|ƒ }|||d$ d]d d^}E|ÿd6 |ÿ }|j@||] tK|ƒk ƒd }´t!|´ƒdk#rjtK|ƒgt!|´ƒ ||´< |j@||] t#|ƒkƒd }´t!|´ƒdk#r°t#|ƒgt!|´ƒ ||´< |E|ƒ} g }|g}x„t4t t!| ƒƒƒD ]n}|dk#ràd\tO| | d$ |D|C|#|Ld$ |N|*ƒ } |j)| t3|Ø|d |Ø| ƒ | ƒ #qàW |j*|ƒt#|ƒd_ }xrt4t t!| ƒƒƒD ]\}|dk$r‚||d |d¦ d ||d }
|j)|d§ |
ƒ $q‚W d6} |j$|!ƒd6 } |j$|!ƒd6 }
d | |êd |0 |Ñd |j$|!ƒd6 }|âd6k(rDd1|+d. |ýd0 d } |j;|âd"|! |â|!ƒ} d | |êd |0 |Ñ|â d }|j;d | d$ |!ƒ}d d¨|; d |êd | d[|=d }g }x8|D ].}|j)d\tL||D|C|#|Ld$ |N|*ƒ ƒ %qâW ||jM|ƒ|jN|dd … |d d©… ƒ |] }||â t#|ƒ }|||d$ d]d d^}E| d6 | }|j@|tK|ƒk ƒd }´t!|´ƒdk&rÐtK|ƒgt!|´ƒ ||´< |j@|t#|ƒkƒd }´t!|´ƒdk'rt#|ƒgt!|´ƒ ||´< |E|ƒ}
g }|g}x„t4t t!|
ƒƒƒD ]n}|dk'rBd\tO|
| d$ |D|C|#|Ld$ |N|*ƒ } |j)| t3|ã|d |ã| ƒ | ƒ 'qBW |j*|ƒt#|ƒd_ }xrt4t t!|
ƒƒƒD ]\}|dk'rä||d |dª d ||d }
|j)|d« |
ƒ 'qäW t#|»|Ì d |»|Ì gƒ|»|Ì< |þdk)r†x¦t4t t!|؃ƒƒD ] }tPtJ|Ð|Ø| |Îd6ƒd6|¢|¶|Î|Ô|Î |Ô||é|»|Ì || |d|—| |˜| |||]||›ƒ}|d d¬k)rJt!|ƒd k)r&|d |—|< |d |˜|< n$|d d |—|< |d |˜|< x tQ|t%ƒ)sj|d })qLW |j |“|gfƒ}“(q~W nŽd | |éd |0 |Ð|× d }xjt4t t!|؃ƒƒD ]V}|j |“|gfƒ}“|tJ||Îd6ƒ |é |j7|Ô|Î ƒ |—|< tJ||Îd6ƒ|˜|< )qºW | dk+rhxÞt4t t!|ッƒD ],}tPtJ|Ñ|ã| |Ïd6ƒd6|¢|¶|Ï|0|Ô |Ï |Ô||ê|»|Ì || | d |—|t t!|؃ƒ |˜|t t!|؃ƒ |||]||›ƒ}|d dk+r,t!|ƒd k+r |d |—||! < |d |˜||! < n,|d d |—||! < |d |˜||! < x tQ|t%ƒ+sL|d }+q.W |j |“|gfƒ}“*q4W nšd | |êd |0 |Ñ|â d }xvt4t t!|ッƒD ]b}|j |“|gfƒ}“|tJ||Ïd6ƒ |ê |j7|0|Ô |Ï ƒ |—||! < tJ||Ïd6ƒ|˜||! < +qœW |Ì|Êd® k,rt| }|}|
}|}|é}|ê}|×}|â}||j$|!ƒd6 }||j$|!ƒd6 }|j |Š|ìgfƒ}Š|j |†|ígfƒ}†|j |‹|îgfƒ}‹|j |‡|ïgfƒ}‡|j |Œ|ðfƒ}Œ|j |ˆ|ñfƒ}ˆ|j ||òfƒ}|j |‰|ófƒ}‰|j || fƒ}|j ||fƒ}|j ||
fƒ}|j ||fƒ}q°W |¡|k -rBd6}nr|t3|j:|¶|Êd¯ |º|Êd° |¡ ƒƒ|jR||¡ ƒ k-r´t3|j:|¶|Êd± |º|Êd² |¡ ƒƒ|jR||¡ ƒ d` } xÈt4| ƒD ]º}t%|d || ƒda d6k.r&dbt|d || ƒ dc t| |Sd ƒ dd } t5| ƒ |º} |¶|j7|•| ƒ } |¶|j6|•| ƒ } |†|j7|•| ƒ }í|†|j6|•| ƒ } |‡|j7|•| ƒ }ï|‡|j6|•| ƒ }!|ˆ|j7|•| ƒ }ñ|ˆ|j6|•| ƒ }"|‰|j7|•| ƒ }ó|‰|j6|•| ƒ }#|dk:r0 xt4|Êd³ d ƒD
]ð}Ì|j |z|j*| |Ì ƒ|¡ |j7|ƒ |j*| |Ì ƒ|j6|ƒ |¡ gfƒ}z|j |{| |Ì gfƒ}{|j |||j*| |Ì ƒ|¡ |j6|ƒ |j*| |Ì ƒ|j7|ƒ gfƒ}||jj)|Š|Ì |¡ |j7|ƒ | |Ì |j6|ƒ |¡ ƒ |kj)|í|Ì ƒ |lj)|Š|Ì |¡ |j6|ƒ | |Ì |j7|ƒ ƒ |mj)|‹|Ì |¡ |j7|ƒ |!|Ì |j6|ƒ |¡ ƒ |nj)|ï|Ì ƒ |oj)|‹|Ì |¡ |j6|ƒ |!|Ì |j7|ƒ ƒ |md´ d |ndµ d |od¶ d |d k 1rœ|jE|jF|md· |nd¸ |od¹ g|˃|jGjH|mdº |nd» |od¼ gddU|jGjH|ËddU ƒ}ÎtS|nd½ |od¾ ƒ}$t3|md¿ ƒ|k 1rxt3tI||΃ƒt%|jD|mdÀ ƒƒ |mdÁ< tJ||Î|$ƒ|ndÂ< tT||Î|$ƒ|odÃ< |j |}|Œt |Ì|! ƒt |Ìd |! ƒ… |¡ |j7|ƒ |"t |Ì|! ƒt |Ìd |! ƒ… |j6|ƒ |¡ fƒ}}|j |~|ñt |Ì|! ƒt |Ìd |! ƒ… fƒ}~|j ||Œt |Ì|! ƒt |Ìd |! ƒ… |¡ |j6|ƒ |"t |Ì|! ƒt |Ìd |! ƒ… |j7|ƒ fƒ}|zdÄ d |{dÅ d ||dÆ d |d k 3rò|jE|jF|zdÇ |{dÈ ||dÉ g|˃|jGjH|zdÊ |{dË ||dÌ gddU|jGjH|ËddU ƒ}ÎtS|{dÍ ||dÎ ƒ}$|›dDk3rLt5dPƒ t5dQƒ t5dPƒ dO}›t3|zdÏ ƒ|k 3r†t3tI||΃ƒt%|jD|zdÐ ƒƒ |zdÑ< tJ||Î|$ƒ|{dÒ< tT||Î|$ƒ||dÓ< |zdÔ |jdÕ< |{dÖ |kd×< ||dØ |ldÙ< |zdÚ |}dÛ< |{dÜ |~dÝ< ||dÞ |dß< |j |}|t |Ì|! ƒt |Ìd |! ƒ… |¡ |j7|ƒ |#t |Ì|! ƒt |Ìd |! ƒ… |j6|ƒ |¡ fƒ}}|j |~|ót |Ì|! ƒt |Ìd |! ƒ… fƒ}~|j ||t |Ì|! ƒt |Ìd |! ƒ… |¡ |j6|ƒ |#t |Ì|! ƒt |Ìd |! ƒ… |j7|ƒ fƒ}|}dà d |~dá d |dâ d |d k 5rÚ|jE|jF|}dã |~dä |då g|˃|jGjH|}dæ |~dç |dè gddU|jGjH|ËddU ƒ}ÎtS|~dé |dê ƒ}$t3|}dë ƒ|k 5r¶t3tI||΃ƒt%|jD|}dì ƒƒ |}dí< tJ||Î|$ƒ|~dî< tT||Î|$ƒ|dï< |jE|jF|zdð |{dñ ||dò g|˃|jGjH|zdó |{dô ||dõ gddU|jGjH|ËddU ƒ}Î|0|jE|jF|zdö |] |{d÷ ||dø g|˃|jGjH|zdù |] |{dú ||dû gddU|jGjH|ËddU ƒ }Ï|0d" }Ô|Ìdk7r:|jE|jF|jdü |jdý |kdþ |kdÿ |ld |ld g|˃|jGjH|jd |jd |kd |kd |ld |ld gddU|jGjH|ËddU ƒ}ÔtS|kd |ld ƒ}$|j |wtI|j$|!ƒd6 |‘|Ì |Ôƒfƒ}w|j |wtJ|j$|!ƒd6 |‘|Ì |Ô|$ƒfƒ}x|j |wtT|j$|!ƒd6 |‘|Ì |Ô|$ƒfƒ}y|0d" }Ô|Ìdk8r||jE|jF|md
|md |nd |nd
|od |od g|˃|jGjH|md |md |nd |nd |od |od gddU|jGjH|ËddU ƒ}ÔtS|nd |od ƒ}$|j |wtI|j$|!ƒd6 |’|Ì |Ôƒfƒ}w|j |xtJ|j$|!ƒd6 |’|Ì |Ô|$ƒfƒ}x|j |ytT|j$|!ƒd6 |’|Ì |Ô|$ƒfƒ}y|j |€dd6|j$|!ƒ fƒ}€|j |€d d6|j$|!ƒ fƒ}€|j |p|t |Ìd" |! ƒt |Ìd d" |! ƒ… fƒ}p|j |q|t |Ìd" |! ƒt |Ìd d" |! ƒ… fƒ}q|j |v|Žt |Ìd" |! ƒt |Ìd d" |! ƒ… fƒ}v|j |r|“t |Ìd" |! ƒt |Ìd d" |! ƒ… fƒ}r/q W |‚j)|‘|Ì ƒ |„j)|’|Ì ƒ | |Êd d }%| |Êd d }&| |Êd d }'| |Êd } | |Êd } | |Êd } |dk-rÀ| d6k ;rv|j:|j8| d | d ƒ|]| ƒ|_k ;rv|j2|_ƒ|]| }Í|jU| |j8| d | d ƒ ƒ}(t3|Í|j7|(ƒ ƒt%|jD| ƒƒ } t3|Í|j6|(ƒ ƒt%|jD| ƒƒ } |jE|jF| | | g|˃|jGjH| | | gddU|jGjH|ËddU ƒ}ÎtS|&|'ƒ}$| d | d | d |d k rè|ê|j7|Ô|Î ƒ })nP|âdk?r|é|j7|Ô|Î ƒ })n0|é|× |ê|â d |× d |â |j7|Ô|Î ƒ })|j ||j$|!ƒd6 |‚| fƒ}|j ||j$|!ƒd6 |„| fƒ}|}×|}â|} |}
|}|j |q|fƒ}q|}|j |q|fƒ}q|Ð|× |j7|΃ d }*|Ñ|â |j7|σ d }+| tI|)|Ôƒ|Ÿ d },| tJ|)|Ô|$ƒ|Ÿ d }-| tT|)|Ô|$ƒ|Ÿ d }.tI|j8| d | d | d ƒ|× |΃tI|‚| |øƒ|Ÿ d }/tT|j8| d | d | d ƒ|× |Î|$ƒtT|‚| |ø|$ƒ|Ÿ d }0tJ|j8| d | d | d ƒ|× |Î|$ƒtJ|‚| |ø|$ƒ|Ÿ d }1tI|j8|]| d | d | d ƒ|â |0|Ï ƒ|] tI|„| |ûƒ|Ÿ d }2tJ|j8|]| d | d | d ƒ|â |0|Ï |$ƒtJ|„| |û|$ƒ|Ÿ d }3tT|j8|]| d | d | d ƒ|â |0|Ï |$ƒtT|„| |û|$ƒ|Ÿ d }4|j8|,|/ d |-|1 d |.|0 d ƒ}×|j8|,|2 d |-|3 d |.|4 d ƒ}â|›dOkBrºd6}×|j8|/d |1d |0d ƒ|j7|jE|jF|,|-|.g|˃|jGjH|,|-|.gddU|jGjH|ËddU ƒƒ d }5|j8|2d |3d |4d ƒ|j7|jE|jF|,|] |-|.g|˃|jGjH|,|] |-|.gddU|jGjH|ËddU ƒƒ d }6|j |r||* |5 fƒ}r|j |r||+ |6 fƒ}r|j$|!ƒt3|/|, ƒ |! tK|,|/gƒ }7|,|-|.f}8|/|1|0f}9|/|,kDrb|8|9 }9}8tV|8|9|7ƒ}:|j$|!ƒt3|2|, ƒ |! tK|,|2gƒ };|,|-|.f}8|2|3|4f}9|2|,kDrê|8|9 }9}8tV|8|9|;ƒ}<|j$|!ƒ}ë|j |v|×|j8| d | d | d ƒ|×|ë |! d |j6|0d" |Î ƒ |Õ |– t3|j6|Ô|Î ƒƒd" |! d6|j$|!ƒ fƒ}v|j |v|â|j8|]| d | d | d ƒ|â|ë |! d |j6|0d" |Ï ƒ |Ö |– t3|j6|Ô|Ï ƒƒd" |! d6|j$|!ƒ fƒ}v|j |~|j*|:d d d d … ƒfƒ}~|jW|:d |/ d |:d |1 d |7|0 d ƒ}|jW|7|/ d |:d |1 d |:d |0 d ƒt!|:d ƒd kGrB|j |}|7|¡ |j7|ƒ |j*|:d ƒ|j6|ƒ |¡ fƒ}}|j ||7|¡ |j6|ƒ |j*|:d ƒ|j7|ƒ fƒ}n¤|j |}|7d d d … |¡ |j7|ƒ |j*|:d d d d … ƒ|j6|ƒ |¡ fƒ}}|j ||7d d d!… |¡ |j6|ƒ |j*|:d d d d"… ƒ|j7|ƒ fƒ}|j |~|j*||š|T|S }U|U}?xÊt4|Sd ƒD ]¸}|U|T|d }@|@}W|@|W ||j6|Wƒ d ||j7|Wƒ }XxFt3|Xƒd+kLrZ|W|X }W|@|W ||j6|Wƒ d ||j7|Wƒ }XLqW |j7|Wƒ| d ||j7|Wƒ }Z|j8d || ƒ|j6|Wƒ d ||j7|Wƒ }[|j9|[|Zƒ}A|Ad"|0 kLrÎd6}A|d |d d ||j7|Aƒ }]t3|?|@ ƒ}B|j*|st |d |= ƒ d … ƒ|¡ |j7|Bƒ |j*|ut |d |= ƒ d … ƒ|j6|Bƒ |¡ }C|j*|st |d |= ƒ d … ƒ|¡ |j6|Bƒ |j*|ut |d |= ƒ d … ƒ|j7|Bƒ }D|C|st |d |= ƒ d …< |D|ut |d |= ƒ d …< |dkPrÖ|j*|}t ||> ƒ d … ƒ|¡ |j7|Bƒ |j*|t ||> ƒ d … ƒ|j6|Bƒ |¡ }E|j*|}t ||> ƒ d … ƒ|¡ |j6|Bƒ |j*|t ||> ƒ d … ƒ|j7|Bƒ }F|E|}t ||> ƒ d …< |F|t ||> ƒ d …< |j*|dt ||= ƒ d … ƒ|¡ |j7|Bƒ |j*|ft ||= ƒ d … ƒ|j6|Bƒ |¡ }G|j*|dt ||= ƒ d … ƒ|¡ |j6|Bƒ |j*|ft ||= ƒ d … ƒ|j7|Bƒ }H|G|dt ||= ƒ d …< |H|ft ||= ƒ d …< |j*|gt ||= ƒ d … ƒ|¡ |j7|Bƒ |j*|it ||= ƒ d … ƒ|j6|Bƒ |¡ }I|j*|gt ||= ƒ d … ƒ|¡ |j6|Bƒ |j*|it ||= ƒ d … ƒ|j7|Bƒ }J|I|gt ||= ƒ d …< |J|it ||= ƒ d …< xŽt4t | ƒƒD ]|}K|K|j&| d" ƒkZr†de|] }L|}t |d |> d"|K d |! ƒ d }M|t |d |> d"|K d |! ƒ d }Nxt4dfƒD ]‚}Ot%|Oƒ|B dg }P|j*|}t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j7|Pƒ |j*|t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|j6|Pƒ |¡ }E|j*|}t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j6|Pƒ |j*|t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|j7|Pƒ }F|j8|Ed' |M d |Fd( |N d ƒ|LkSrÈ|E|M |M|C|K |M|Ed |M }E|F|N |N|D|K |N|Fd |N }FP |j8|Ed) |M d |Fd* |N d ƒ}LQqzW |E|}t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ …< |F|t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ …< |j*|dt |d |= ƒ d … ƒ|¡ |j7|Pƒ |j*|ft |d |= ƒ d … ƒ|j6|Pƒ |¡ }G|j*|dt |d |= ƒ d … ƒ|¡ |j6|Pƒ |j*|ft |d |= ƒ d … ƒ|j7|Pƒ }H|G|dt |d |= ƒ d …< |H|ft |d |= ƒ d …< t |d |> d"|K d |! ƒ dk Xr|j*|}t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j7|Bƒ |j*|t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|j6|Bƒ |¡ }E|j*|}t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j6|Bƒ |j*|t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|j7|Bƒ }F|E|}t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ …< |F|t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ …< np|j*|}t |d |> d"|K d |! ƒ d … ƒ|¡ |j7|Bƒ |j*|t |d |> d"|K d |! ƒ d … ƒ|j6|Bƒ |¡ }E|j*|}t |d |> d"|K d |! ƒ d … ƒ|¡ |j6|Bƒ |j*|t |d |> d"|K d |! ƒ d … ƒ|j7|Bƒ }F|E|}t |d |> d"|K d |! ƒ d …< |F|t |d |> d"|K d |! ƒ d …< |j*|gt |d |= ƒ d … ƒ|¡ |j7|Bƒ |j*|it |d |= ƒ d … ƒ|j6|Bƒ |¡ }I|j*|gt |d |= ƒ d … ƒ|¡ |j6|Bƒ |j*|it |d |= ƒ d … ƒ|j7|Bƒ }J|I|gt |d |= ƒ d …< |J|it |d |= ƒ d …< Pqæde|] }L|}t |d |> d"|K d |! ƒ d }M|t |d |> d"|K d |! ƒ d }Nxt4dfƒD ]}Ot%|Oƒ|B dg }Pt |d |> d"|K d |! ƒ dk \râ|j*|}t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j7|Pƒ |j*|t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|j6|Pƒ |¡ }E|j*|}t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j6|Pƒ |j*|t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ … ƒ|j7|Pƒ }Fn|j*|}t |d |> d"|K d |! ƒ d … ƒ|¡ |j7|Pƒ |j*|t |d |> d"|K d |! ƒ d … ƒ|j6|Pƒ |¡ }E|j*|}t |d |> d"|K d |! ƒ d … ƒ|¡ |j6|Pƒ |j*|t |d |> d"|K d |! ƒ d … ƒ|j7|Pƒ }F|j8|Ed+ |M d |Fd, |N d ƒ|Lk^rÊ|E|Ed- |Ed. |C|K |Ed/ |Ed |Ed0 }E|F|Fd1 |Fd2 |D|K |Fd3 |Fd |Fd4 }FP |j8|Ed5 |M d |Fd6 |N d ƒ}LZqüW t |d |> d"|K d |! ƒ dk _rÖ|E|}t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ …< |F|t |d |> d"|K d |! ƒ t |d |> d"|K d |! ƒ …< n`|E|}t |d |> d"|K d |! ƒ d …< |F|t |d |> d"|K d |! ƒ d …< |j*|gt |d |= ƒ d … ƒ|¡ |j7|Pƒ |j*|it |d |= ƒ d … ƒ|j6|Pƒ |¡ }I|j*|gt |d |= ƒ d … ƒ|¡ |j6|Pƒ |j*|it |d |= ƒ d … ƒ|j7|Pƒ }J|I|gt |d |= ƒ d …< |J|it |d |= ƒ d …< |j*|}t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j7|Bƒ |j*|t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|j6|Bƒ |¡ }E|j*|}t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|¡ |j6|Bƒ |j*|t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ … ƒ|j7|Bƒ }F|E|}t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ …< |F|t |d |> |K|! d ƒ t |d |> d"|K d |! ƒ …< |j*|dt |d |= ƒ d … ƒ|¡ |j7|Bƒ |j*|ft |d |= ƒ d … ƒ|j6|Bƒ |¡ }G|j*|dt |d |= ƒ d … ƒ|¡ |j6|Bƒ |j*|ft |d |= ƒ d … ƒ|j7|Bƒ }H|G|dt |d |= ƒ d …< |H|ft |d |= ƒ d …< PqæW |@}?t |d |> ƒ}QKqÌW |™|- }'|'dk dr¨d"|0 |' }'x~t4|QƒD ]n}tS|~|d ||d ƒ}$|jE|jF|}|d |}|d | |~|d |~|d | ||d ||d | g|˃|jGjH|}|d |}|d | |~|d |~|d | ||d ||d | gddU|jGjH|ËddU ƒ}ÔtI||d |Ôƒ|w|d < tT||d |Ô|$ƒ|y|d < tJ||d |Ô|$ƒ|x|d < dq¶W ~t t!|sƒƒ}S|j |z|sfƒ}º|j |{|tfƒ}¶|j |||ufƒ}R|j |j|dfƒ}d|j |k|efƒ}e|j |l|ffƒ}f|j |m|gfƒ}g|j |n|hfƒ}h|j |o|ifƒ}i|^}]|j*|}t |Sd" |! | ƒ d … ƒ}S|j*|}d t |Sd" |! | ƒ … ƒ}T|j*|~t |Sd" |! | ƒ d … ƒ}U|j*|~d t |Sd" |! | ƒ … ƒ}V|j*|t |Sd" |! | ƒ d … ƒ}W|j*|d t |Sd" |! | ƒ … ƒ}X|j*|pt |Sd" |! | ƒ d … ƒ}Y|j*|pd t |Sd" |! | ƒ … ƒ}Z|j*|qt |Sd" |! | ƒ d … ƒ}[|j*|qd t |Sd" |! | ƒ … ƒ}\t t!|Tƒ| ƒ}]|j |Tt tX| |] d ƒƒt tX| |] d ƒ|] ƒ… |Tt tX| |] dh d ƒƒt tX| |] dh d ƒ|] ƒ… fƒ}^|j |Xt tX| |] d ƒƒt tX| |] d ƒ|] ƒ… |Xt tX| |] dh d ƒƒt tX| |] dh d ƒ|] ƒ… fƒ}_|j |Zt tX| |] d ƒƒt tX| |] d ƒ|] ƒ… |Zt tX| |] dh d ƒƒt tX| |] dh d ƒ|] ƒ… fƒ}`t t!|Sƒ|S ƒ}]g }ag }bg }cx†t4|SƒD ]x}|j |a|St tX| d" |! d ||] ƒƒt tX| d" |! d ||] ƒd"|! ƒ… fƒ}a|j |a|St tX| d" |! dh d ||] ƒƒt tX| d" |! dh d ||] ƒd"|! ƒ… fƒ}a|j |b|Wt tX| d" |! d ||] ƒƒt tX| d" |! d ||] ƒd"|! ƒ… fƒ}b|j |b|Wt tX| d" |! dh d ||] ƒƒt tX| d" |! dh d ||] ƒd"|! ƒ… fƒ}b|j |c|Yt tX| d" |! d ||] ƒƒt tX| d" |! d ||] ƒd"|! ƒ… fƒ}c|j |c|Yt tX| d" |! dh d ||] ƒƒt tX| d" |! dh d ||] ƒd"|! ƒ… fƒ}ciqþW |j |^|afƒ}d|j |_|bfƒ}e|j |`|cfƒ}f|jYd7||] didj}g|jYd8||] didj}h|j$deƒ|jZt#|fƒƒ dZ }i|j[|d|] |e|] ƒ}jt\|j|d|] |e|] dkdl |j]dƒ |j^ƒ }k|kj_ƒ }l|j`|j|f|idm |lja|gƒ |lja|hƒ |j+|`|adndo |j+|b|cdndo |j+|`|a dndo |j+|b|c dndo |jbdpdq}m|mjcdrƒ |mjddt#|pƒƒ |j.dsƒ |j/dtƒ |jed9d ƒ |jfd:dƒ |j0|#du t|)ƒ d* ƒ |j1ƒ t t!|Tƒ| ƒ}]|j |Td|]… |Tt tX| |] d/ ƒƒt tX| |] d/ ƒ|] ƒ… fƒ}^|j |Vd|]… |Vt tX| |] d/ ƒƒt tX| |] d/ ƒ|] ƒ… fƒ}n|j |Zd|]… |Zt tX| |] d/ ƒƒt tX| |] d/ ƒ|] ƒ… fƒ}`|j |\d|]… |\t tX| |] d/ ƒƒt tX| |] d/ ƒ|] ƒ… fƒ}ot t!|Sƒ|S ƒ}]g }ag }pg }cg }qx¶t4|SƒD ]¨}|j |a|St tX||] ƒƒt tX||] ƒd"|! ƒ… fƒ}a|j |a|St tX| d" |! d/ ||] ƒƒt tX| d" |! d/ ||] ƒd"|! ƒ… fƒ}a|j |p|Ut tX||] ƒƒt tX||] ƒd"|! ƒ… fƒ}p|j |p|Ut tX| d" |! d/ ||] ƒƒt tX| d" |! d/ ||] ƒd"|! ƒ… fƒ}p|j |c|Yt tX||] ƒƒt tX||] ƒd"|! ƒ… fƒ}c|j |c|Yt tX| d" |! d/ ||] ƒƒt tX| d" |! d/ ||] ƒd"|! ƒ… fƒ}c|j |q|[t tX||] ƒƒt tX||] ƒd"|! ƒ… fƒ}q|j |q|[t tX| d" |! d/ ||] ƒƒt tX| d" |! d/ ||] ƒd"|! ƒ… fƒ}qoqàW |j |^|afƒ}d|j |n|pfƒ}r|j |`|cfƒ}f|j |o|qfƒ}sdv}tx(t#|sƒ|t dk sr |tde }trqæW |jYd;||] didj}g|jYd<||] didj}h|jg|j$dwƒ|jZt#|sƒ|t ƒ|j&tK|sƒ|t ƒ dx |j&tK|sƒ|t ƒ ƒ}i|j[|d|] |r|] ƒ}jt\|j|d|] |r|] dkdl |j]dƒ |j^ƒ }u|uj_ƒ }v|vja|hƒ |j`|j|s|t |idm |vja|gƒ |j+|`|adndo |j+|b|cdndo |j+|`|a dndo |j+|b|c dndo |jbdpdq}m|mjcdyt|tƒ dz ƒ |mjdtK|q|t ƒt#|q|t ƒƒ |j.d{ƒ |j/dtƒ |jed=d ƒ |jfd>dƒ |j0|#d| t|)ƒ d* ƒ |j1ƒ t t!|Sƒ|S ƒ}]g }ag }pg }cxt4|SƒD ] }|j |a|St tX||] ƒƒt tX||] ƒd"|! ƒ… fƒ}a|j |a|St tX| d" |! d/ ||] ƒƒt tX| d" |! d/ ||] ƒd"|! ƒ… fƒ}a|j |p|Ut tX||] ƒƒt tX||] ƒd"|! ƒ… fƒ}p|j |p|Ut tX| d" |! d/ ||] ƒƒt tX| d" |! d/ ||] ƒd"|! ƒ… fƒ}p|j |c|Yt tX||] ƒƒt tX||] ƒd"|! ƒ… fƒ}c|j |c|Yt tX| d" |! d/ ||] ƒƒt tX| d" |! d/ ||] ƒd"|! ƒ… fƒ}cuq8W |j |^|afƒ}d|j |n|pfƒ}r|j |`|cfƒ}f|jYd?||] didj}g|jYd@||] didj}h|j$deƒ|jZt#|fƒƒ dZ }i|j[|d|] |r|] ƒ}jt\|j|d|] |r|] dkdl |j]dƒ |j^ƒ }w|wj_ƒ }x|xja|hƒ |j`|j|f|idm |xja|gƒ |j+|`|adndo |j+|b|cdndo |j+|`|a dndo |j+|b|c dndo |jbdpdq}m|mjcdrƒ |mjddt#|pƒƒ |j.d{ƒ |j/dtƒ |jedAd ƒ |jfdBdƒ |j0|#d} t|)ƒ d* ƒ |j1ƒ |}|~||p|qd~œ}y|j>|y|#d t|)ƒ d€ ƒ t t!|}ƒƒ}zd|œ }{tKd‚|jZ|j8t#|ºƒtK|ºƒ d t#|¶ƒtK|¶ƒ d t#|RƒtK|Rƒ d ƒ|{ ƒgƒ}|tK|dƒ}}d_t#|gƒ }~tK|hƒ}d_t#|hƒ }€|j |j |}|] |~|] ||] fƒdCt t!|}ƒƒfƒj"}~}~~~||dd}‚||j |j |d|] |e|] |f|] fƒdDt t!|dƒƒfƒj"dd}ƒ||j |j |g|] |h|] |i|] fƒdEt t!|gƒƒfƒj"dd}„||j |j |º|] |¶|] |R|] fƒdFt t!|ºƒƒfƒj"dd}…~d~e~f~g~h~i|{|j7|'ƒ |j6|(ƒ }†|{|j7|(ƒ }‡|{ |j6|'ƒ |j6|(ƒ }ˆ|j*|†|‡|ˆgƒ|jGjH|j*|†|‡|ˆgƒƒ }‰g }Šg }‹xlt4|"ƒD ]`}|Šj)g ƒ |‹j)g ƒ t'|#dƒ t|)ƒ d„ t|ƒ d… dGƒ}Œ|Œj>d†ƒ |Œj1ƒ {q–W t5dPƒ t5d‡ƒ t5dˆƒ t5d‡ƒ xtt4|zƒD ]d}t%|d ƒda d6k|rxdbt|d ƒ dc t|zƒ dd } t5| ƒ ||df |] }||df |] }Ž||d f |] }|p| }|q| }át3|v| ƒ}‘|w| |j7|'ƒ |j6|(ƒ |x| |j7|(ƒ |y| |j6|(ƒ |j6|'ƒ }’|j6|'ƒ| |j7|'ƒ| }“|j7|'ƒ|j6|(ƒ | |j7|(ƒ|Ž |j6|'ƒ|j6|(ƒ | }”|j7|(ƒ |j7|'ƒ | |j6|(ƒ|Ž |j7|(ƒ|j6|'ƒ | }•|j6|'ƒ|d |j7|'ƒ| }–|j7|'ƒ|j6|(ƒ |d |j7|(ƒ|Ž |j6|'ƒ|j6|(ƒ | }—|j7|(ƒ |j7|'ƒ |d |j6|(ƒ|Ž |j7|(ƒ|j6|'ƒ | }˜d }™|”d6k ~rÖ|j8|•d |“d ƒ|k ~rÖd6}™d }š|—d6k r|j8|˜d |–d ƒ|k rd6}š|™|š d6k|r(|j@|B|’kƒd }›t!|›ƒdk|r(|‘d6k|r(|Š|›dG j)|ƒ |‹|›dH j)d%d‰|ád dŠ|;d d |‘ d‹ ƒ |||||Ž||{|†|‡|ˆ|º|¶|R|]|€|p|r||||| | |4|5|2|3|‰|;|”|•|“|—|˜|–|ƒ|„|‚|…|}|~||€ƒ)}œt'|#dƒ t|)ƒ d„ t|›dI ƒ d… dŒƒ}Œ|Œj>dTjhtit|œƒƒd ƒ |Œj1ƒ |q(W |›dOk€rºt'|#dŽ t|)ƒ d dGƒ}n t'|#dŽ t|)ƒ d… dGƒ}|j>dƒ t'|#d‘ t|)ƒ d… dGƒ}ž|žj>d’ƒ t'|#d“ t|)ƒ d… dGƒ}Ÿ|Ÿj>d”ƒ |Ÿj>d•ƒ x|t4|"ƒD ]p}|Ÿj>dTjhtit|Š| ƒƒd ƒ |žj>dTjhtit|‹| ƒƒd ƒ |j>t|B| d‚ ƒd ƒ qHW |Ÿj1ƒ |žj1ƒ |j1ƒ ||›gS (J Nr )ÚKDTree)Úfsolve)Úinterp1d)r? )ÚAxes3D)Úascii)ÚRTi' z$/cooling_functions/cooling_function_z.tabT)Zdropr r
r iè )Úleafsizeg®Gáz®ï?gÍÌÌÌÌÌô?ÚMsunÚRsunÚmpÚGg @@Úkbg¤p=
×Óv@g õ@g @r rE g ð?)Údtypeg @z#depthghÉh÷Æ!w>gQ²@³-(;Úlogz lambda(T) (10^-23 erg cm^3/s)zTemperature (keV)z/cooling_func_z.pdfg:Œ0âŽyE>z/ionisation_fraction.tabú#g×£p=
×ã?g à?g „×—Ag\Âõ(\ó?g333333Ã?z3X-ray emission can be created only in the shock capé gú~j¼t“h?g é( c s s | ]
}g V qd S )Nr! )Ú.0Ú_r! r! r" ú s zrshock..é c s s | ]
}g V qd S )Nr! )r± r² r! r! r" r³ s c s s | ]
}g V qd S )Nr! )r± r² r! r! r" r³ s é c s s | ]
}g V qd S )Nr! )r± r² r! r! r" r³ s é c s s | ]
}g V qd S )Nr! )r± r² r! r! r" r³ s é g €v@g €f@Úno)r/ z/histograms/xstag_parrA z cm
z R1
z d
r r r r r0 Ú z
Crashing windg-C ëâ6?g @@ú )rŽ éÈ g .@g €@g ÐcAg "@g D@göJáÇ-µDÚcubic)ÚkindÚ
assume_sortedgš™™™™™ñ?gffffffæ?g Y@z* ú/z *g $@éd g ÀX@g @Úwhite)Úcolorg333333Ó?)r“ )ÚlevelsÚblack)Úcz%.2f)ÚformatzTemperature [keV]zz/Dzx/Dz/plots/temperature_xz_parg‚vIhÂ%<=g 4@g 3@z Density [z g/cm^3]zy/Dz/plots/density_parz/plots/temperature_xy_par)r r
r rb Úrhoz/plots/char_shock_parz.datgš™™™™™©?g jø@z/histograms/ray_tracing_parZ_binz.dataz# Temperature / Factor
z ***********z Ray tracingg¸ …ëQ¸ò?gq=
×£põ?g„þäqÙ‰EÚaÚ
z/histograms/bin_parz_crashing.dataz # Tangential velocity (km/s)
z /histograms/emiss_part_shock_parz:# Part of the emissivity which depends only on the shock.
z/histograms/temp_emiss_parz9# Temperature of the cell where the emission is created.
z*# Each line of a bin of the line profile.
rF rF g H~Ag H~Agš™™™™™ñ?gú~j¼t“x?rF rF rF rF rF rF rF rF gR¸ …ëÑ!@rF rF rF gR¸ …ëÑ!@rF rF rF g ð¿g ð¿rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF éþÿÿÿrF rÊ rF rÊ rF rÊ rF rÊ rF rÊ rF rF rF rÊ rF rÊ rF rÊ rF rÊ rF rÊ rF rÊ rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF rF )r r )r r rF rÊ )r r )r r rF rÊ )r r )r r rF rÊ rF rF rF rF rF rF rF )jZ
scipy.spatialr Zscipy.optimizer¡ Úscipy.interpolater¢ r r$ rI r? Úmatplotlib.pyplotÚpyplotÚmatplotlib.triÚtriZmpl_toolkits.mplot3dr£ Ú
matplotlibrJ ZpandasÚsysZ
astropy.ior¤ Z
raytracingr¥ Úsetrecursionlimitr* ÚstrZread_hdfZ reset_indexÚreshapeÚ concatenater rT rU ÚTrS r_ Úfloatra rL ÚsplitrX rV ÚplotÚyscaleÚxscaleÚylabelÚxlabelÚsavefigrN Útanr] rW Úprintr) r2 r
Úatan2ÚatanÚlinspacer r( rM r rQ Únanminr# r€ ÚsignÚarccosÚdotÚlinalgÚnormrˆ r‹ rR r} ÚcumsumÚmeanr~ r= Ú
isinstanceÚasinr… rŒ ÚacosrŸ r\ ÚroundÚCircleÚceilÚ
Triangulationr” ÚfigureÚgcfÚgcaÚ tricontourfÚ
add_artistÚcolorbarÚ set_labelÚset_climÚxlimÚylimÚsortÚjoinÚmap( Ú liste_paramr r¡ r¢ r r$ r? ÚpltrÏ r£ rÐ rJ ÚpdrÑ r¤ r¥ r r Úvinf1Úvinf2Ú
mass_ratioÚexÚomegarÈ Úperr r Úbeta1Úbeta2ZTeff_1ZTeff_2Únbr_points_shock_2DÚnbr_points_shock_3DZ nstep_widthÚnbr_bin_profilere Úcut_limÚ wind_primÚwind_secÚphaseÚinclÚiparrh r{ rÖ ÚM_conjÚatomÚionr* Zfcoolr r r r r6 ZSumXzZXhr| r§ r¨ r© Úgravr« Úmass_sumÚM2Úvmaxr† Zvvvrd rc ÚfÚlineÚline2ZthetainfÚasympZdeltafZ
deltathetarf ZZ_ionrg ÚfionZlignert ZjionZnbr_balZstep_balZ debut_balÚMÚEÚdEZT2ÚcoiÚsiiZphase2r Zd_vouluZ
ouverture_minZxstar1Zystar1Zxstar2Zystar2Zx_cdw1Zy_cdw1Zz_cdw1Zx_cdw2Zy_cdw2Zz_cdw2Zx_cdw1gZy_cdw1gZz_cdw1gZx_cdw2gZy_cdw2gZz_cdw2gZkT_allZrho_allZ sigma_allZx_cdZy_cdZz_cdZvol_cdZvx_cdZvy_cdZvz_cdZx_skewZy_skewZz_skewZxx_allZyy_allZzz_allZcote_cdZvt_vecZvt1_vecZvp1_vecZvt2_vecZvp2_vecZ y_cdw1g1_allZ y_cdw2g1_allZ
y_cdw1g1a_allZ
y_cdw2g1a_allZ x_cdw1g1_allZ x_cdw2g1_allZ
x_cdw1g1a_allZ
x_cdw2g1a_allZvol_cdoZkT_alloZrho_alloZvt1_vecoZvt2_vecoZ
sigma_alloZM_oldZphi_vecÚdphiZ dzeta_o_vecZy_o_vecZM_vouluZ phase_voulur9 ÚlmaxZibalZT1ZTelapsÚvorbZCMr4 Zfxstagr r r r r r r r r ÚskewZ theta_vecZ
theta1_vecZ
ThetaCotThetaÚpointZthe1_vecZ pointhererk Zr_testr5 Zind_halfZx_oldr7 Zx_vecZeta_vecZ y_vechereÚr1Úr2r' ZpenteZr_vecÚ sin_thetaÚ cos_thetarƒ r‚ r„ Z slope_vecÚfirstZ
delta1_vecZind_cutÚxaxisZcutZRprimeZthetahZ the1_vechZr1_vechZr2_vechZthetahm1Z the1_vechm1Z
slope_vechÚdthetaZdthe1Zl0_1Zl1_widthZl0Zdwind_prim_chocZil0ÚxhZyhÚvZkT0rl rÇ Zl0_2Zl2_widthZdwind_sec_chocZrl1Zl011Zrl2Zl022r% r& Z increasevolZx_cdw1g1Zy_cdw1g1Zx_cdw2g1Zy_cdw2g1Z x_cdw1g1aZ y_cdw1g1aZ x_cdw2g1aZ y_cdw2g1aZvt1Zvt2Z x_cdw1g11Z y_cdw1g11Zslope_vec1_endZ x_cdw2g11Z y_cdw2g11Zslope_vec2_endZvp1Zvp2ZkT0_1Z l_1_widthZkT1Zrho1_vecZrho0ZkT_widthÚconstÚtestrž ZwheretoZconst3Zconst2ZdrhoZkT0_2Z l_2_widthZkT2Zrho2_vecZisigmar: ZkT1_endZ rho1_vec_endZkT2_endZ rho2_vec_endZv1_endZv2_endZl0_1_endZl0_2_endZ
sigma1_endZ
sigma2_endZicutÚsentenceZxcutZycutZzcutZz_cdw1g1Zz_cdw2g1Z z_cdw1g1aZ z_cdw2g1aÚlatZx_cutm1Zy_cutm1Zz_cutm1r‡ Zvt1cdZsurf11Zsurf21Zx_cd1Zy_cd1Zz_cd1Zx_cdw11Zz_cdw11Zy_cdw11Zx_cdw21Zy_cdw21Zz_cdw21Zsurf12Zsurf22Zzint1Zparam1Zparam2Zxyint1Zzint2Zxyint2ZstepcdZstepallZMoldZMhZphasehZ
angle_to_turnZx_cdnZz_cdnZxx_allnZzz_allnZx_cdw1nZz_cdw1nZx_cdw2nZz_cdw2nZiligneÚdistanceZ xx_alln_oldZ zz_alln_oldZiangleZ anglehereZfinZz_vecZx_balZx_capZy_balZy_capZz_balZz_capZkT_balZkT_capZrho_balZrho_capZ nbr_layerZ x_cap_layerZ z_cap_layerZ kT_cap_layerZ x_bal_layerZ z_bal_layerZ kT_bal_layerZx_layerZz_layerZkT_layerZcircle1Zcircle2rà Ztriang2ÚfigÚaxÚcbarZ y_cap_layerZ
rho_cap_layerZ y_bal_layerZ
rho_bal_layerZy_layerZ rho_layerÚfactorZfig2Úax2Zfig3Zax3ÚtableZ
nbr_pt_totZdeltaZZZndeltaZZÚxmin1Zxmax2Úymin1Zymax2Zensemble_pointsÚtreeZtreecdw1Ztreecdw2ZtreevecZdeltaXZdeltaYZdeltaZZdir_obsZ
temp_emissZ emiss_part1ZfRTÚxxr+ rq rb ZvolZvradÚyt1Zzt1Zpt1Úyt2Zzt2Zpt2Zvis1Zvis2Zind_vradZ ray_tracingZfbinZfhistoZftempr! r! r" ÚrshockŒ sô
D
VV $ 0
"
"
6(. . "00" (. . " 0
*(
* *(
*
&
(
8
F
. NZ "
^j~
*
.
0(
0& (R&6 *
.
*(
*
* (V&6
0(
0
0(
0
RR
NN,
",,:
*6 0 ,
",,:
*6 0
^"* z"2
:8 4
HD4040.V$h*d.V
$ h*d.V$Vb
–&**
ž&** 2228
<(&&J( $( $nVJV $(
0"" HPPTXX66
z‚ 00
„ˆ&@R<:TPR<:TP404040 ( *0 .$fb
^Z^Z^Z
00ÆÂ:66:NNfb ,ÎÊRVŠ†00fb "
00,ÎΊ†:NN:,RT00fb ÆÂNNfb $
æ "( &&&&&&&&&&€€€dldldr
PPPPDdDdDdDjX
"
DdDdDj
\
6>>@ 6 * &
R"FH&JL,, >Š0
$$&
rD N)r )Ú__doc__r# r( r- r= r} r~ r r… rˆ r‹ rŒ r” rŸ rD r! r! r! r" Ú6 s n +==