Calculate N_H along LOS to NS over the period of a binary.
The binary contains optical star with a spherical wind and NS with accretion wake where density is increased by some factor compared to the wind.
The wake side at the current moment coincides with the locus of the wind particles that "met" NS in the past.
The extension of the wake is set by the fraction of the NS orbit, "remembered" by the wake. The width of the wake equals the accretion radius (variable along the orbit) times user-defined factor. The head of the wake is missing so far.
User parameters should be listed in XRB.ini. See default parameters at the bottom (in the future). XRB_Blondin.ini contains parameters for a binary from paper Blondin et al (1990)
See results for GX 301−2 in folder GX_dif_geo/.
How to run
- to plot ionization map, run:
mpiexec -n 8 python xi_map_mpi.py [phase]
where phase is a float number >=0 and <=1
- to calculate NH(phase) for the whole period, run:
mpiexec -n 8 python NH_period_inclined_mpi.py
It will produce PDF files and ASCII *dat files in folder data/
- to plot different models NH(phase), run:
python plot_many_NH_period.py [dirname] [parameter names]
where
-
dirname is the path to the folder with *dat files made beforehand by NH_period_inclined_mpi
-
parameter names is optional argument: parameters separated by spaces to indicate in the legends. Name of parameters may be found in the headers of dat-files.
