Skip to content

Commit

Permalink
This is Python 81.
Browse files Browse the repository at this point in the history
  • Loading branch information
jhmatthews committed May 18, 2016
2 parents 87481a3 + a24345a commit c41bf0f
Show file tree
Hide file tree
Showing 37 changed files with 817 additions and 143 deletions.
5 changes: 4 additions & 1 deletion .travis.yml
Original file line number Diff line number Diff line change
Expand Up @@ -49,11 +49,14 @@ script:
- cd ../examples/travis
- Setup_Py_Dir
- py balmer_test > balmer_test.out
- cd ..
- cd ../core/
- Setup_Py_Dir
- py -i cv_macro_benchmark
- py -i cv_standard
- py -i -d sv_detailedmode
- py -i fiducial_agn
- py -i 1d_sn
- cd ../beta/
- py -i agn_ss_2010_modela.pf
- py -i ulx1
- py -i ngc5548
26 changes: 26 additions & 0 deletions examples/README.md
Original file line number Diff line number Diff line change
@@ -0,0 +1,26 @@
### examples

This folder contains a number of standard parameter files, used for testing and porting to the user.

### core

Core files which are standard tests

* cv_standard.pf -- Standard CV model (Shlosman & Vitello 1993; Long & Knigge 2002).
* fiducial_agn.pf -- QSO model from Higginbottom et al. (2013).
* 1d_sn.pf -- simple Supernovae model as presented in Kerzendorf & Sim (2013).
* star.pf -- simple spherical stellar wind model

More performance intensive core matom models:

* cv_macro_benchmark.pf -- Macro-atom CV run from Matthews et al. (2015).
* m16_agn.pf -- Macro-atom AGN run from Matthews et al. (2016).


### beta

Files of astrophysical interest not currently fully tested.

* ngc5548.pf -- designed to represent rough params for NGC 5548
* ulx1.pf -- First attempt at a ULX model
* agn_ss_2010_modela.pf -- AGN model A from Sim et al. (2010).
File renamed without changes.
59 changes: 59 additions & 0 deletions examples/beta/ngc5548.pf
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
Wind_type(0=SV,1=Sphere,2=Previous,3=Proga,4=Corona,5=knigge,6=thierry,7=yso,8=elvis,9=shell) 0
Atomic_data data/standard78
photons_per_cycle 20000000
Ionization_cycles 30
spectrum_cycles 20
Coord.system(0=spherical,1=cylindrical,2=spherical_polar,3=cyl_var) 1
Wind.dim.in.x_or_r.direction 100
Wind.dim.in.z_or_theta.direction 100
Wind_ionization(0=on.the.spot,1=LTE,2=fixed,3=recalc_bb,5=recalc_pow,6=pairwise_bb,7=pairwise_pow,8=matrix_bb,9=matrix_pow) 7
Line_transfer(0=pure.abs,1=pure.scat,2=sing.scat,3=escape.prob,6=macro_atoms,7=macro_atoms+aniso.scattering) 3
Thermal_balance_options(0=everything.on,1=no.adiabatic) 0
System_type(0=star,1=binary,2=agn) 2
disk.type(0=no.disk,1=standard.flat.disk,2=vertically.extended.disk) 1
Disk_radiation(y=1) 1
Wind_radiation(y=1) 1
QSO_BH_radiation(y=1) 1
Rad_type_for_disk(0=bb,1=models)_to_make_wind 0
Rad_type_for_agn(0=bb,1=models,3=power_law,4=cloudy_table)_to_make_wind 3
mstar(msol) 1e7
rstar(cm) 8.85667e+12
disk.mdot(msol/yr) 0.02
Disk.illumination.treatment(0=no.rerad,1=high.albedo,2=thermalized.rerad,3=analytic) 0
Disk.temperature.profile(0=standard;1=readin) 0
disk.radmax(cm) 1e15
lum_agn(ergs/s) 1e43
agn_power_law_index -0.9
wind.radmax(cm) 1e17
wind.t.init 1e5
wind.mdot(msol/yr) 5
sv.diskmin(wd_rad) 50
sv.diskmax(wd_rad) 100
sv.thetamin(deg) 70
sv.thetamax(deg) 82
sv.mdot_r_exponent 0
sv.v_infinity(in_units_of_vescape 1
sv.acceleration_length(cm) 1e16
sv.acceleration_exponent 1.0
wind.filling_factor(1=smooth,<1=clumpted) 1
Rad_type_for_disk(0=bb,1=models,2=uniform)_in_final_spectrum 0
Rad_type_for_agn(0=bb,1=models,3=power_law,4=cloudy_table)_in_final_spectrum 3
spectrum_wavemin 200
spectrum_wavemax 2600
no_observers 8
angle(0=pole) 20
angle(0=pole) 40
angle(0=pole) 60
angle(0=pole) 70
angle(0=pole) 75
angle(0=pole) 80
angle(0=pole) 85
angle(0=pole) 89
live.or.die(0).or.extract(anything_else) 1
Select_specific_no_of_scatters_in_spectra(y/n) n
Select_photons_by_position(y/n) n
spec.type(flambda(1),fnu(2),basic(other) 1
Extra.diagnostics(0=no) 0
Use.standard.care.factors(1=yes) 1
reverb.type(0=None,1=Photon,2=Wind) 0
Photon.sampling.approach(0=T,1=(f1,f2),2=cv,3=yso,4=user_defined,5=cloudy_test,6=wide,7=AGN,8=logarithmic) 7
59 changes: 59 additions & 0 deletions examples/beta/ulx1.pf
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
Wind_type(0=SV,1=Sphere,2=Previous,3=Proga,4=Corona,5=knigge,6=thierry,7=yso,8=elvis,9=shell) 0
Atomic_data std_79
photons_per_cycle 20000000
Ionization_cycles 30
spectrum_cycles 20
Coord.system(0=spherical,1=cylindrical,2=spherical_polar,3=cyl_var) 1
Wind.dim.in.x_or_r.direction 100
Wind.dim.in.z_or_theta.direction 100
Wind_ionization(0=on.the.spot,1=LTE,2=fixed,3=recalc_bb,5=recalc_pow,6=pairwise_bb,7=pairwise_pow,8=matrix_bb,9=matrix_pow) 9
Line_transfer(0=pure.abs,1=pure.scat,2=sing.scat,3=escape.prob,6=macro_atoms,7=macro_atoms+aniso.scattering) 5
Thermal_balance_options(0=everything.on,1=no.adiabatic) 0
System_type(0=star,1=binary,2=agn) 2
disk.type(0=no.disk,1=standard.flat.disk,2=vertically.extended.disk) 0
Disk_radiation(y=1) 0
Wind_radiation(y=1) 1
QSO_BH_radiation(y=1) 1
Rad_type_for_disk(0=bb,1=models)_to_make_wind 0
Rad_type_for_agn(0=bb,1=models,3=power_law,4=cloudy_table)_to_make_wind 3
mstar(msol) 10.0
rstar(cm) 8.85667e+5
disk.mdot(msol/yr) 1.25e-8
Disk.illumination.treatment(0=no.rerad,1=high.albedo,2=thermalized.rerad,3=analytic) 0
Disk.temperature.profile(0=standard;1=readin) 0
disk.radmax(cm) 1e7
lum_agn(ergs/s) 1e39
agn_power_law_index -0.9
wind.radmax(cm) 1e10
wind.t.init 1e6
wind.mdot(msol/yr) 5e-6
sv.diskmin(wd_rad) 100
sv.diskmax(wd_rad) 200
sv.thetamin(deg) 70
sv.thetamax(deg) 82
sv.mdot_r_exponent 0
sv.v_infinity(in_units_of_vescape 2
sv.acceleration_length(cm) 1e8
sv.acceleration_exponent 0.2
wind.filling_factor(1=smooth,<1=clumpted) 1
Rad_type_for_disk(0=bb,1=models,2=uniform)_in_final_spectrum 0
Rad_type_for_agn(0=bb,1=models,3=power_law,4=cloudy_table)_in_final_spectrum 3
spectrum_wavemin 2
spectrum_wavemax 50
no_observers 8
angle(0=pole) 20
angle(0=pole) 40
angle(0=pole) 60
angle(0=pole) 70
angle(0=pole) 75
angle(0=pole) 80
angle(0=pole) 85
angle(0=pole) 89
live.or.die(0).or.extract(anything_else) 1
Select_specific_no_of_scatters_in_spectra(y/n) n
Select_photons_by_position(y/n) n
spec.type(flambda(1),fnu(2),basic(other) 2
Extra.diagnostics(0=no) 0
Use.standard.care.factors(1=yes) 1
reverb.type(0=None,1=Photon,2=Wind) 0
Photon.sampling.approach(0=T,1=(f1,f2),2=cv,3=yso,4=user_defined,5=cloudy_test,6=wide,7=AGN,8=logarithmic) 7
2 changes: 1 addition & 1 deletion examples/1d_sn.pf → examples/core/1d_sn.pf
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,7 @@ Wind_type(0=SV,1=Sphere,2=Previous,3=Proga,4=Corona,5=knigge,6=homologous,7=yso,
Atomic_data data/standard_sn_kurucz
photons_per_cycle 80000
Ionization_cycles 20
spectrum_cycles 1
spectrum_cycles 10
Coord.system(0=spherical,1=cylindrical,2=spherical_polar,3=cyl_var) 0
Wind.dim.in.x_or_r.direction 30
Wind_ionization(0=on.the.spot,1=LTE,2=fixed,3=recalc_bb,5=recalc_pow,6=pairwise_bb,7=pairwise_pow,8=matrix_bb,9=matrix_pow) 3
Expand Down
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
Wind_type(0=SV,1=Sphere,2=Previous,3=Proga,4=Corona,5=knigge,6=thierry,7=yso,8=elvis,9=shell) 0
Atomic_data data/h20_hetop_standard78
photons_per_cycle 100000
Ionization_cycles 1
spectrum_cycles 0
Ionization_cycles 30
spectrum_cycles 10
Coord.system(0=spherical,1=cylindrical,2=spherical_polar,3=cyl_var) 1
Wind.dim.in.x_or_r.direction 30
Wind.dim.in.z_or_theta.direction 30
Expand Down
File renamed without changes.
2 changes: 1 addition & 1 deletion examples/fiducial_agn.pf → examples/core/fiducial_agn.pf
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
Wind_type(0=SV,1=Sphere,2=Previous,3=Proga,4=Corona,5=knigge,6=thierry,7=yso,8=elvis,9=shell) 0
Atomic_data data/standard78
photons_per_cycle 200000
photons_per_cycle 2000000
Ionization_cycles 30
spectrum_cycles 20
Coord.system(0=spherical,1=cylindrical,2=spherical_polar,3=cyl_var) 1
Expand Down
72 changes: 72 additions & 0 deletions examples/core/m16_agn.pf
Original file line number Diff line number Diff line change
@@ -0,0 +1,72 @@
Wind_type(0=SV,1=Sphere,2=Previous,3=Proga,4=Corona,5=knigge,6=thierry,7=yso,8=elvis,9=shell) 0.0
Atomic_data data/h10_hetop_lohe1_standard78
photons_per_cycle 20000000.0
Ionization_cycles 20.0
spectrum_cycles 5.0
Coord.system(0=spherical,1=cylindrical,2=spherical_polar,3=cyl_var) 1.0
Wind.dim.in.x_or_r.direction 100.0
Wind.dim.in.z_or_theta.direction 200.0
Wind_ionization(0=on.the.spot,1=LTE,2=fixed,3=recalc_bb,5=recalc_pow,6=pairwise_bb,7=pairwise_pow) 9.0
Line_transfer(0=pure.abs,1=pure.scat,2=sing.scat,3=escape.prob,6=macro_atoms,7=macro_atoms+aniso.scattering) 7.0
Thermal_balance_options(0=everything.on,1=no.adiabatic) 0.0
System_type(0=star,1=binary,2=agn) 2.0
disk.type(0=no.disk,1=standard.flat.disk,2=vertically.extended.disk) 1.0
Disk_radiation(y=1) 1.0
Wind_radiation(y=1) 1.0
QSO_BH_radiation(y=1) 1.0
Rad_type_for_disk(0=bb,1=models)_to_make_wind 0.0
Rad_type_for_agn(0=bb,1=models,3=power_law,4=cloudy_table)_to_make_wind 3.0
mstar(msol) 1000000000.0
rstar(cm) 8.85667e+14
disk.mdot(msol/yr) 5.0
Disk.illumination.treatment(0=no.rerad,1=high.albedo,2=thermalized.rerad,3=analytic) 0.0
Disk.temperature.profile(0=standard;1=readin) 0.0
disk.radmax(cm) 1e+17
lum_agn(ergs/s) 1e+45
agn_power_law_index -0.9
Torus(0=no,1=yes) 0.0
wind.radmax(cm) 1e+20
wind.t.init 100000.0
wind.mdot(msol/yr) 5.0
sv.diskmin(wd_rad) 50
sv.diskmax(wd_rad) 100.0
sv.thetamin(deg) 70
sv.thetamax(deg) 82.0
sv.mdot_r_exponent 0.0
sv.v_infinity(in_units_of_vescape 1.0
sv.acceleration_length(cm) 1e19
sv.acceleration_exponent 0.5
wind.filling_factor(1=smooth,<1=clumped) 0.01
Rad_type_for_disk(0=bb,1=models,2=uniform)_in_final_spectrum 0.0
Rad_type_for_agn(0=bb,1=models,3=power_law,4=cloudy_table)_in_final_spectrum 3.0
spectrum_wavemin 500.0
spectrum_wavemax 7500.0
no_observers 20.0
angle(0=pole) 20.0
angle(0=pole) 40.0
angle(0=pole) 60.0
angle(0=pole) 70.0
angle(0=pole) 71.0
angle(0=pole) 72.0
angle(0=pole) 73.0
angle(0=pole) 74.0
angle(0=pole) 75.0
angle(0=pole) 76.0
angle(0=pole) 77.0
angle(0=pole) 78.0
angle(0=pole) 79.0
angle(0=pole) 80.0
angle(0=pole) 81.0
angle(0=pole) 82.0
angle(0=pole) 83.0
angle(0=pole) 84.0
angle(0=pole) 85.0
angle(0=pole) 89.0
live.or.die(0).or.extract(anything_else) 1.0
Select_specific_no_of_scatters_in_spectra(y/n) n
Select_photons_by_position(y/n) n
spec.type(flambda(1),fnu(2),basic(other) 1.0
reverb.type 0.0
Extra.diagnostics(0=no) 0.0
Use.standard.care.factors(1=yes) 1.0
Photon.sampling.approach(0=T,1=(f1,f2),2=cv,3=yso,4=user_defined), 7.0
46 changes: 46 additions & 0 deletions examples/core/star.pf
Original file line number Diff line number Diff line change
@@ -0,0 +1,46 @@
Wind_type(0=SV,1=Sphere,2=Previous,3=Proga,4=Corona,5=knigge,6=homologous,7=yso,8=elvis,9=shell) 1
Atomic_data data/standard78
photons_per_cycle 100000
Ionization_cycles 20
spectrum_cycles 20
Coord.system(0=spherical,1=cylindrical,2=spherical_polar,3=cyl_var) 0
Wind.dim.in.x_or_r.direction 30
Wind_ionization(0=on.the.spot,1=LTE,2=fixed,3=recalc_bb,6=pairwise_bb,7=pairwise_pow) 3
Line_transfer(0=pure.abs,1=pure.scat,2=sing.scat,3=escape.prob,6=macro_atoms,7=macro_atoms+aniso.scattering) 3
Thermal_balance_options(0=everything.on,1=no.adiabatic) 0
System_type(0=star,1=binary,2=agn) 0
disk.type(0=no.disk,1=standard.flat.disk,2=vertically.extended.disk) 0
Star_radiation(y=1) 1
Disk_radiation(y=1) 0
Boundary_layer_radiation(y=1) 0
Wind_radiation(y=1) 1
Rad_type_for_star(0=bb,1=models)_to_make_wind 0
Model_file kurucz91/kurucz91.ls
mstar(msol) 52.5
rstar(cm) 1.32e12
tstar 42000
disk.type(0=no.disk,1=standard.flat.disk,2=vertically.extended.disk) 0
Torus(0=no,1=yes) 0
wind.radmax(cm) 1e13
wind.t.init 40000
stellar_wind_mdot(msol/yr) 5.1e-6
stellar.wind.radmin(cm) 1.32e12
stellar.wind_vbase(cm) 2e+07
stellar.wind.v_infinity(cm) 2.25e+08
stellar.wind.acceleration_exponent 1
wind.filling_factor(1=smooth,<1=clumped) 1
Rad_type_for_star(0=bb,1=models,2=uniform)_in_final_spectrum 0
Model_file kurucz91/kurucz91.ls
spectrum_wavemin 850
spectrum_wavemax 1850
no_observers 1
angle(0=pole) 45
phase(0=inferior_conjunction) 0.5
live.or.die(0).or.extract(anything_else) 1
Select_specific_no_of_scatters_in_spectra(y/n) n
Select_photons_by_position(y/n) n
spec.type(flambda(1),fnu(2),basic(other) 1
Extra.diagnostics(0=no) 0
Use.standard.care.factors(1=yes) 1
reverb.type 0
Photon.sampling.approach(0=T,1=(f1,f2),2=cv,3=yso,4=user_defined,5=cloudy_test,6=wide,7=AGN,8=logarithmic) 2
File renamed without changes.
47 changes: 47 additions & 0 deletions examples/core/template_ionloop.pf
Original file line number Diff line number Diff line change
@@ -0,0 +1,47 @@
Wind_type() 9
Atomic_data data/standard78
photons_per_cycle 100000
Ionization_cycles 1
spectrum_cycles 0
Coord.system() 0
Wind.dim.in.x_or_r.direction 4
Wind_ionization() 7
Line_transfer() 5
Thermal_balance_options 0
System_type() 2
disk.type() 0
Star_radiation(y=1) 0
Disk_radiation(y=1) 0
Wind_radiation(y=1) 0
QSO_BH_radiation(y=1) 1
Rad_type_for_star(0=bb,1=models)_to_make_wind 0
Rad_type_for_agn()_to_make_wind 4
mstar(msol) 0.8
rstar(cm) 1e10
tstar 1000000
lum_agn(ergs/s) 1e+21
agn_power_law_index -0.9
low_energy_break(ev) 0.136
high_energy_break(ev) 20000
Torus(0=no,1=yes) 0
wind.radmax(cm) 1.00000000001e11
wind.t.init 10000
shell_wind_mdot(msol/yr) 0.00472e-17
shell.wind.radmin(cm) 1e11
shell.wind_v_at_rmin(cm) 1.00000
shell.wind.v_at_rmax(cm) 1.000010
shell.wind.acceleration_exponent 1
wind.filling_factor(1=smooth,<1=clumped) (1) 1
spec.type(flambda(1),fnu(2),basic(other) 2
reverb.type 0
Use.standard.care.factors(1=yes) 1
Photon.sampling.approach() 5
Num.of.frequency.bands 3
Lowest_energy_to_be_considered(eV) 0.0001
Highest_energy_to_be_considered(eV) 100000000
Band.boundary(eV) 2000
Band.boundary(eV) 10000
Band.minimum_fraction) 0.3
Band.minimum_fraction) 0.4
Band.minimum_fraction) 0.3

Loading

0 comments on commit c41bf0f

Please sign in to comment.