PyDataLondon2026-FastForwardingModels/cv_wind.pf at main · AustenWallis/PyDataLondon2026-FastForwardingModels · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
System_type(star,cv,bh,agn,previous)                   cv

### Parameters for the Central Object
Central_object.mass(msol)                       0.8
Central_object.radius(cm)               7.25182e+08
Central_object.radiation(yes,no)                  yes
Central_object.rad_type_to_make_wind(bb,models)                   bb
Central_object.temp                           40000
Binary.mass_sec(msol)                           0.4
Binary.period(hr)                               3.2

### Parameters for the Disk (if there is one)
Disk.type(none,flat,vertically.extended,rmin>central.obj.rad)                 flat
Disk.radiation(yes,no)                          yes
Disk.rad_type_to_make_wind(bb,models,mod_bb)                   bb
Disk.temperature.profile(standard,readin)             standard
Disk.mdot(msol/yr)                            1e-08
Disk.radmax(cm)                         2.17555e+10

### Parameters for Boundary Layer or the compact object in an X-ray Binary or AGN
Boundary_layer.radiation(yes,no)                   no

### Parameters describing the various winds or coronae in the system
Wind.number_of_components                         1
Wind.type(SV,star,hydro,corona,kwd,homologous,shell,imported)                   SV
Wind.coord_system(spherical,cylindrical,polar,cyl_var)          cylindrical
Wind.dim.in.x_or_r.direction                     30
Wind.dim.in.z_or_theta.direction                   30

### Parameters associated with photon number, cycles,ionization and radiative transfer options
Photons_per_cycle                            900000
Ionization_cycles                                 25
Spectrum_cycles                                   12
Wind.ionization(on.the.spot,ML93,LTE_tr,LTE_te,fixed,matrix_bb,matrix_pow,matrix_est)            matrix_bb
Line_transfer(pure_abs,pure_scat,sing_scat,escape_prob,thermal_trapping,macro_atoms_escape_prob,macro_atoms_thermal_trapping)     thermal_trapping
Wind.radiation(yes,no)                          yes
Surface.reflection.or.absorption(reflect,absorb,thermalized.rerad)              reflect
Wind_heating.extra_processes(none,adiabatic,nonthermal,both)            adiabatic
Atomic_data                     data/standard80.dat

### Parameters for Domain 0
Wind.mdot(msol/yr)                            1e-09
SV.diskmin(units_of_rstar)                        4
SV.diskmax(units_of_rstar)                       12
SV.thetamin(deg)                                 20
SV.thetamax(deg)                                 65
SV.mdot_r_exponent                                0
SV.v_infinity(in_units_of_vescape                    3
SV.acceleration_length(cm)                    7e+10
SV.acceleration_exponent                        1.5
SV.gamma(streamline_skew;1=usually)                    1
SV.v_zero_mode(fixed,sound_speed)                fixed
SV.v_zero(cm/s)                              600000
Wind.radmax(cm)                               1e+12
Wind.t.init                                   40000
Wind.filling_factor(1=smooth,<1=clumped)                    1

### Parameters defining the spectra seen by observers

Central_object.rad_type_in_final_spectrum(bb,models,uniform,mono)                   bb
Disk.rad_type_in_final_spectrum(bb,models,uniform,mono,mod_bb)                   bb

### The minimum and maximum wavelengths in the final spectra and the number of wavelength bins
Spectrum.nwave                                10000
Spectrum.wavemin(Angstroms)                     850
Spectrum.wavemax(Angstroms)                    1850

### The observers and their location relative to the system
Spectrum.no_observers                             4
Spectrum.angle(0=pole)                           10
Spectrum.angle(0=pole)                           30
Spectrum.angle(0=pole)                           60
Spectrum.angle(0=pole)                           80
Spectrum.orbit_phase(0=inferior_conjunction)                  0.5
Spectrum.orbit_phase(0=inferior_conjunction)                  0.5
Spectrum.orbit_phase(0=inferior_conjunction)                  0.5
Spectrum.orbit_phase(0=inferior_conjunction)                  0.5
Spectrum.live_or_die(live.or.die,extract)              extract
Spectrum.type(flambda,fnu,basic)              flambda

### Parameters for Reverberation Modeling (if needed)
Reverb.type(none,photon,wind,matom)                 none

### Other parameters
Photon_sampling.approach(T_star,cv,yso,AGN,tde_bb,min_max_freq,user_bands,cloudy_test,wide,logarithmic)                   cv