This HTML file was created by the program doc_tsuite.pl. doc_tsuite.txt contains a tab delimited list of the files.
agn_blr_albedo.in
measure rayleigh scattering of Lya
title measure rayleigh scattering of Lya
c model from Korista, K., & Ferland, G. 1998, ApJ, 495, 672
c
c commands controlling continuum =========
agn 6.683 -1.20 -1.20 -0.90
ionization parameter 1.0
c
c commands for density & abundances =========
init file = "ism.ini"
abundances old solar 84
hden 11.0
c
c commands controlling geometry =========
stop total column density = 23.75
c
c other commands for details =========
iterate
c
c commands controlling output =========
normalize to "feka" 2
print line faint -1
*print last iteration
print diffuse continua
print lines inward
punch continuum last "agn_blr_albedo.con" units kev
punch dr "agn_blr_albedo.dr"
punch reflected continuum last "agn_blr_albedo.ref"
punch emitted continuum last "agn_blr_albedo.emt"
c
c line would appear to change the rest of the spectrum
c agn_blr_albedo.in
c class blr
c ========================================
This model computes the albedo of a fairly
standard BLR cloud. This is the type of
model that was presented in the BLR albedo
paper by Korista & Ferland, 1998, ApJ 495, 672.
The print diffuse continua command enters continuum
fluxs into the emission line stack. The asserts then
check that these continua have the expected brightness.
agn_lex00_u0.in
intermediate-ionization x-ray ionized cloud from Lexington 2000
title intermediate-ionization x-ray ionized cloud from Lexington 2000
stop column density 16
print line sort wavelength
*constant temperature 720000
no induced processes
print lines column linear
normalise to "H 1" 1216
print line faint -1
print last iteration
hden 5
iterate
phi(h) 15.477 range 7.353 to 735.3
interpolate (-8 -3)
continue (-6 1.)
continue (-2 14.5)
continue (0.4771 12.7)
continue (0.8663 10.6)
continue (3.8663 7.6)
continue (6 1)
continue (9 -3)
element abundance helium -1
element abundance carbon -3.432
element abundance nitrogen -3.959
element abundance oxygen -3.097
element abundance neon -3.959
element abundance magnesium -4.4318
element abundance silicon -4.456
element abundance sulphur -4.7959
element abundance argon -5.4318
element abundance iron -4.398
element lithium off
element beryllium off
element boron off
element fluorine off
element sodium off
element aluminium off
element phosphrous off
element chlorine off
element potasium off
element calcium off
element scandium off
element titanium off
element vanadium off
element chromium off
element manganese off
element cobalt off
element nickel off
element copper off
element zinc off
c
c agn_lex00_u0.in
agn_lex00_u1.in
high-ionization x-ray ionized cloud from Lexington 2000
title high-ionization x-ray ionized cloud from Lexington 2000
stop column density 16
*constant temper 720000 K
no induced processes
print lines column linear
normalise to "H 1" 1216
print line faint -1
print last iteration
hden 5
iterate
phi(h) 16.477 range 7.353 to 735.3
interpolate (-8 -3)
continue (-6 1.)
continue (-2 14.5)
continue (0.4771 12.7)
continue (0.8663 10.6)
continue (3.8663 7.6)
continue (6 1)
continue (9 -3)
element abundance helium -1
element abundance carbon -3.432
element abundance nitrogen -3.959
element abundance oxygen -3.097
element abundance neon -3.959
element abundance magnesium -4.4318
element abundance silicon -4.456
element abundance sulphur -4.7959
element abundance argon -5.4318
element abundance iron -4.398
element lithium off
element beryllium off
element boron off
element fluorine off
element sodium off
element aluminium off
element phosphrous off
element chlorine off
element potasium off
element calcium off
element scandium off
element titanium off
element vanadium off
element chromium off
element manganese off
element cobalt off
element nickel off
element copper off
element zinc off
c
c agn_lex00_u1.in
agn_lex00_um1.in
low-ionization x-ray ionized cloud from Lexington 2000
title low-ionization x-ray ionized cloud from Lexington 2000
stop column density 16
no induced processes
print lines column linear
print lines sort intensity
normalise to "H 1" 1216
print line faint -2
print last iteration
hden 5
iterate
phi(h) 14.477 range 7.353 to 735.3
interpolate (-8 -3)
continue (-6 1.)
continue (-2 14.5)
continue (0.4771 12.7)
continue (0.8663 10.6)
continue (3.8663 7.6)
continue (6 1)
continue (9 -3)
element abundance helium -1
element abundance carbon -3.432
element abundance nitrogen -3.959
element abundance oxygen -3.097
element abundance neon -3.959
element abundance magnesium -4.4318
element abundance silicon -4.456
element abundance sulphur -4.7959
element abundance argon -5.4318
element abundance iron -4.398
element lithium off
element beryllium off
element boron off
element fluorine off
element sodium off
element aluminium off
element phosphrous off
element chlorine off
element potasium off
element calcium off
element scandium off
element titanium off
element vanadium off
element chromium off
element manganese off
element cobalt off
element nickel off
element copper off
element zinc off
c
c
c agn_lex00_um1.in
agn_reflector.in
model of Compton reflector
title model of Compton reflector
c
c commands controlling continuum =========
c this is a powerlaw that extends into the x-rays, but by default
c has exponential cutoff in FIR
power law -1. hi=1000000000
c the ionization parameter for this powerlaw
ionization parameter 1.
c add background so that FIR/radio does not have zero flux
background
iterate
c
c commands for density & abundances =========
hden 11
c
c commands controlling geometry =========
stop total column density 24
c
c other commands for details =========
constant temper 5
c
c commands controlling output =========
print line faint 1
print line inward
print diffuse continuum
c
c this is to produce lines with lots of contrast with the continuum
set PunchLWidth 100 km/sec
c a fig in part 3 of hazy uses data from this model. the fig shows lines
c with above and below linewidth
c set PunchLWidth c
punch continuum last "agn_reflector.con" units kev
punch continuum reflected last "agn_reflector.alb" units kev
punch overview last "agn_reflector.ovr"
punch dr last "agn_reflector.dr"
c
table lines "LineList_BLR.dat"
c agn_reflector.in
c class blr
c ========================================
c
This is a model of the Compton reflector in AGN. It is a
constant temperature since models of this region often make that assumption.
A plot in Part I of Hazy shows the incident and reflected portions of the
continuum. The code will complain that the cloud is Compton thick since it
is not really designed to simulate this situation.
agn_S_curve_grid.in
temperature across Spitzer thermal stability S curve
title temperature across Spitzer thermal stability S curve
c
c commands controlling continuum =========
c this is the generic AGN continuum
table agn
c the ionization parameter
ionization parameter 0.00 vary
grid from 0 to 3 in 0.25 dex steps
c add background so that FIR/radio does not have zero flux
background
iterate
c
c commands for density & abundances =========
c intended to make sim run as fast as possible
hden 8
init file "ism.ini"
element sulphur off
element chlorine off
element argon off
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c commands controlling output
c this produces a very large number of lines in the output mostly
c because Hbeta is very weak in hot gas. Choose Lya as the normalization
c line to make line list smaller
normalize to "totl" 1216
print line faint 10
punch grid "agn_S_curve_grid.grd" last
punch average "agn_S_curve_grid.avr" last
temperature hydrogen 1
end of temperatures
punch heating "agn_S_curve_grid.het" last
punch cooling "agn_S_curve_grid.col" last
c
asse hydro 2 temper 4.512 err0.05 grid 4.600 4.746 4.964 5.519 5.593 5.694 6.035 6.249 6.374 6.555 6.779 6.927
c
c agn_S_grid.in
c class BLR
c ========================================
c
This is part of a series of models that check the temperature
through the S curve in the Spitzer three-phase model of
ISM stability.
agn_warm_absorber.in
simple warm absorber model
title simple warm absorber model
c
c commands controlling continuum =========
table power law
ionization parameter 0
c
c commands for density & abundances =========
hden 6
abundances old solar 84
c
c commands controlling geometry =========
stop column density 22
c
c other commands for details =========
iterate to convergence
c
c commands controlling output =========
print line faint 10
plot continuum range 1 1000
punch continuum "agn_warm_absorber.con" last units kev
punch dr "agn_warm_absorber.dr" last
c this will be the default units, like the printout
punch line optical depths last "agn_warm_absorber_lab.lin"
c this will be the same, in rygbergs
punch line optical depths last "agn_warm_absorber_ryg.lin" units rydberg
c
c warrm_absorber
c class blr nlr
c ========================================
c
this is a simple warm absorber model. It makes a plot of the
transmitted continuum, and generates a list of lines with
significant optical depths
aperture_beam_int.in
test aperture beam command with intensity
title test aperture beam command with intensity
c
c commands for density & abundances =========
hden 0
init file "honly.ini"
c
c other commands for details =========
constant temper 4
c
c commands controlling continuum =========
laser 2
ionization -2
c
c commands controlling geometry =========
stop thickness 10
sphere
aperture beam
c
c commands controlling output =========
c
c aperture_beam_int.in
c class geometry
c ========================================
This is a homogeneous sphere that is especially simple. The model
is a test of the aperture command, a command that simulates observing
part of an extended object. In this case the aperture is a beam
contered on the center of the nebula, with a line of sight
extending through the object.
The code carries along a dummy emission line ("Unit 1") with a constant intensity
of 1e-10 erg cm-3 s-1. The line goes through all of the code's infrastructure,
and when the calculation is complete, the program confirms that the "luminosity"
of the line is the emitting volume times 1e-10. The aperture command
is verified by asserting that the emission line has the correct "luminosity".
In this case the inner radius is not specified so the returned value is unity.
aperture_beam_lum.in
test aperture beam command with luminosity
title test aperture beam command with luminosity
c
c commands controlling continuum =========
laser 2
q(h) 31.5
c
c commands for density & abundances =========
hden 0
init file "honly.ini"
c
c commands controlling geometry =========
stop thickness 10
radius 9
sphere
aperture beam
c
c other commands for details =========
constant temper 4
c
c commands controlling output =========
c
c aperture_beam_lum.in
c class geometry
c ========================================
This is a homogeneous sphere that is especially simple. The model
is a test of the aperture command, a command that simulates observing
part of an extended object. In this case the aperture is a beam
contered on the center of the nebula, with a line of sight
extending through the object.
The code carries along a dummy emission line ("Unit 1") with a constant intensity
of 1e-10 erg cm-3 s-1. The line goes through all of the code's infrastructure,
and when the calculation is complete, the program confirms that the "luminosity"
of the line is the emitting volume times 1e-10. The aperture command
is verified by asserting that the emission line has the correct "luminosity".
aperture_slit.in
test aperture slit command with luminosity
title test aperture slit command with luminosity
c
c commands controlling continuum =========
laser 2
q(h) 31.5
c
c commands for density & abundances =========
hden 0
init file "honly.ini"
c
c commands controlling geometry =========
radius 9 10
sphere
aperture slit
c
c other commands for details =========
constant temper 4
c
c commands controlling output =========
c
c aperture_slit.in
c class geometry
c ========================================
c
This is a homogeneous sphere that is especially simple. The model
is a test of the aperture command, a command that simulates observing
part of an extended object. In this case the aperture is a long
slit contered on the center of the nebula, extending beyond the outer
reaches of the matter.
The code carries along a dummy emission line ("Unit 1") with a constant intensity
of 1e-10 erg cm-3 s-1. The line goes through all of the code's infrastructure,
and when the calculation is complete, the program confirms that the "luminosity"
of the line is the emitting volume times 1e-10. The aperture command
is verified by asserting that the emission line has the correct "luminosity".
blr_f92.in
standard blr cloud in Ferland et al. 1992
title standard blr cloud in Ferland et al. 1992
c
c commands controlling continuum =========
table agn
ionization parameter -1
c
c commands for density & abundances =========
hden 11
abundances old solar 84
c
c commands controlling geometry =========
stop column density 25
c
c other commands for details =========
* stop zone 20
* element iron ionization -2 0 -2 -2 -2 - 2
iterate to convergence
c
c commands controlling output =========
plot continuum
normalize to "H 1" 1216 = 105.90
print line faint -2
punch overview last "blr_f92.ovr"
punch convergence reason "blr_f92.cvr"
punch element nitrogen last "blr_f92.nit"
punch element calcium last "blr_f92.ca"
punch dr "blr_f92.dr"
punch transmitted continuum "blr_f92.trn" last
c
c blr_f92.ini
c class blr
c ========================================
This is similar to one of the BLR models presented in Ferland et al. (1992)
for the well-studied Seyfert galaxy NGC 5548. It has a very large column
density and is marginally optically thick to electron scattering. The spectrum
is given relative to Lya, and the intensity of this line is reset to
produce a spectrum that is on the same intensity scale as that paper.
blr_fp89.in
final F+P 1989 BLR model table 3
title final F+P 1989 BLR model table 3
c
c commands controlling continuum =========
table agn
ionization parameter -0.5
c
c commands for density & abundances =========
hden 9.5
init file="c84.ini"
abundances old solar 84
c
c commands controlling geometry =========
stop column density 25.5
c
c other commands for details =========
iterate to convergence
failures 2
c
c commands controlling output =========
c set a relatively high faint level since lines relative to Lya=100
print line faint -1
print h-like departure coefficients
plot continuum
normalize to "H 1" 1216 = 85.60
punch heating "blr_fp89.het"
punch convergence reason "blr_fp89.cvr"
punch cooling "blr_fp89.col"
punch overview last "blr_fp89.ovr"
punch element nitrogen last "blr_fp89.nit"
punch results last "blr_fp89.rlt"
punch dr "blr_fp89.dr"
c
c blr_fp89.in
c class blr
c ========================================
Ferland and Persson (1989) presented this calculation of a BLR cloud.
The differences between the present predictions and
those given by FP are largely due to improved treatment
of Balmer line escape and destruction. The spectrum is given relative
to a Lya intensity of 100. The column density is VERY large, to
reproduce intensities of low-ionization lines, especially the Ca II lines.
blr_hizqso.in
high Z quasar cloud
title high Z quasar cloud
c
c commands controlling continuum =========
table agn
ionization parameter 1
c
c commands for density & abundances =========
c actual elec den will be about 100x larger
hden 8
abundances starburst 35
c
c commands controlling geometry =========
c this is more gramage than it appears because of high Z
stop column density 21
c
c other commands for details =========
c once elec den fail occurs due to tripping over He+ - He I-front
c code does recover so not a problem. The default for this is 10,
c and is set low because (hopefully) this would never occur
c hydrogen is a minor constituent in this high metallicity gas
c the electron density is well above the hydrogen density
failures 3
iterate 3
c
c commands controlling output =========
normalize "totl" 1909
punch dr "blr_hizqso.dr"
punch overview last "blr_hizqso.ovr"
punch element nitrogen last "blr_hizqso.nit"
punch element chlorine last "blr_hizqso.cl"
punch element iron last "blr_hizqso.fe"
punch element neon last "blr_hizqso.ne"
punch element calcium last "blr_hizqso.ca"
punch element zinc last "blr_hizqso.zn"
punch heating last "blr_hizqso.het"
punch cooling last "blr_hizqso.col"
c
c blr_hizqso.in
c class blr
c ========================================
This is a model of a very high metallicty BLR cloud. It checks the
intensities of some of the brigher lines, and is a check that the code
can converge a cloud with this high Z.
Secondary ionization is very important when H is highly ionized,
due to very high He abundance. Sec ionization becomes important
at the He+ - He ionization front, where H+/H is 1e-5.
blr_kk81.in
old blr
title old blr
c
c commands controlling continuum =========
interpolate (0 -5) (.05 -5) (.1 0) (1 -0.5) (7.353 -2.233)
continue (735 -3.233) (800 -15) (8.e6 -15)
f(nu) -7.32148
constant gas pressure
c
c commands for density & abundances =========
hden 9.60206
init file="c84.ini"
abundances he-1 c-3.699 n-4 o-3.1549 ne-4 na=-8 mg-4.5229
continue al-10 si-4.4229 s-10 ar-10 ca-10 fe-4.5229 ni=-8
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
failures 2
iterate to convergence
c
c commands controlling output =========
print line faint -1
print ages
normalise "H 1" 1216 100
plot continuum range -3
punch overview last "blr_kk81.ovr"
punch dr last "blr_kk81.dr"
punch continuum last "blr_kk81.con"
punch ionizing continuum "blr_kk81.ion"
c
c
c blr_kk81.in
c class blr
c ========================================
c
This is the 34standard34 BLR model presented by Kwan and Krolik (1981).
>>refer blr cloud Kwan, J., & Krolik, J. 1981, ApJ, 250, 478
Checks:
- Compare line intensities to previous versions of CLOUDY
by entering into table on page Error! Bookmark not defined..
- A code caution that the resulting total pressure was not constant
is to be expected. The KK calculation assumed constant gas pressure,
but radiation pressure was significant. Because of this the sum
of gas plus radiation pressure was not constant (but gas pressure was).
blr_level2.in
test dominant blr_level2 lines
title test dominant blr_level2 lines
c
c commands controlling continuum =========
table agn
ionization -2
c
c commands controlling geometry =========
stop column density 19
c
c commands for density & abundances =========
hden 10
element phosphorus abundance 0
c
c other commands for details =========
iterate
c
c commands controlling output =========
punch overview "blr_level2.ovr"
punch dr "blr_level2.dr"
c
c blr_level2.in
c class blr
c ========================================
This model checks predictions for the "blr_level2" lines. These are
lines that are normally very weak, have Opacity Project wavelengths,
and g-bar collision strengths. Phosphorus is given a large abundance
so that its blr_level2 lines are significant.
blr_n09_p18.in
BLR model, density 1e09 cm-3, flux of H-ion phots 1e18 cm2 s-1
title BLR model, density 1e09 cm-3, flux of H-ion phots 1e18 cm2 s-1
c
c commands controlling continuum =========
table agn
phi(h) 18
c
c commands for density & abundances =========
hden 9
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
init "c84.ini"
iterate convergence
c these are to try to speed things up
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n09_p18.ovr" last
punch dr "blr_n09_p18.dr" last
punch convergence reason "blr_n09_p18.cvr"
c
c blr_n09_p18.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n09_p18_Z20.in
BLR model, density 1e09 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20
title BLR model, density 1e09 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20
c
c commands controlling continuum =========
table agn
phi(h) 18
c
c commands for density & abundances =========
hden 9
abundances starburst 20
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate convergence
c these are to try to speed things up
init "c84.ini"
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n09_p18_Z20.ovr"
punch dr "blr_n09_p18_Z20.dr"
punch convergence reason "blr_n09_p18_Z20.cvr"
punch convergence error "blr_n09_p18_Z20.cve"
c
c blr_n09_p18_Z20.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n09_p20.in
BLR model, density 1e09 cm-3, flux of H-ion phots 1e20 cm2 s-1
title BLR model, density 1e09 cm-3, flux of H-ion phots 1e20 cm2 s-1
c
c commands controlling continuum =========
table agn
phi(h) 20
c
c commands for density & abundances =========
hden 9
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
init "c84.ini"
iterate convergence
c these are to try to speed things up
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n09_p18.ovr" last
punch dr "blr_n09_p18.dr" last
punch convergence reason "blr_n09_p18.cvr"
c
c blr_n09_p18.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n09_p20_Z20.in
BLR model, density 1e09 cm-3, flux of H-ion phots 1e20 cm2 s-1, Z=20
title BLR model, density 1e09 cm-3, flux of H-ion phots 1e20 cm2 s-1, Z=20
c
c commands controlling continuum =========
table agn
phi(h) 20
c
c commands for density & abundances =========
abundances starburst 20
hden 9
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
init "c84.ini"
iterate convergence
c these are to try to speed things up
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n09_p20_Z20.ovr" last
punch dr "blr_n09_p20_Z20.dr" last
punch convergence reason "blr_n09_p20_Z20.cvr"
c
c blr_n09_p20_Z20.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n09_p22.in
BLR model, density 1e09 cm-3, flux of H-ion phots 1e20 cm2 s-1
title BLR model, density 1e09 cm-3, flux of H-ion phots 1e20 cm2 s-1
c
c commands controlling continuum =========
table agn
phi(h) 22
c
c commands for density & abundances =========
hden 9
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
init "c84.ini"
iterate convergence
c these are to try to speed things up
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n09_p22.ovr" last
punch dr "blr_n09_p22.dr" last
punch convergence reason "blr_n09_p22.cvr"
c
c blr_n09_p22.in
c class blr
c ========================================
This is one of the models that sample the LOC plane.
blr_n09_p22_Z20.in
BLR model, density 1e09 cm-3, flux of H-ion phots 1e22 cm2 s-1, Z=20
title BLR model, density 1e09 cm-3, flux of H-ion phots 1e22 cm2 s-1, Z=20
c
c commands controlling continuum =========
table agn
phi(h) 22
c
c commands for density & abundances =========
abundances starburst 20
hden 9
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
init "c84.ini"
iterate convergence
c these are to try to speed things up
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n09_p22_Z20.ovr" last
punch dr "blr_n09_p22_Z20.dr" last
punch convergence reason "blr_n09_p22_Z20.cvr"
c
c blr_n09_p22_Z20.in
c class blr
c ========================================
This is one of the models that sample the LOC plane.
blr_n11_p20.in
BLR model, density 1e11 cm-3, flux of H-ion phots 1e20 cm2 s-1
title BLR model, density 1e11 cm-3, flux of H-ion phots 1e20 cm2 s-1
c
c commands controlling continuum =========
table agn
phi(h) 20
c
c commands for density & abundances =========
hden 11
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate convergence
c these are to try to speed things up
init "c84.ini"
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch convergence reason "blr_n11_p20.cvr"
punch overview "blr_n11_p20.ovr" last
punch dr "blr_n11_p20.dr" last
c
c blr_n11_p20.in
c class blr
c ========================================
c
This is one of the 5 models that sample the LOC plane.
blr_n11_p20_Z20.in
BLR model, density 1e11 cm-3, flux of H-ion phots 1e20 cm2 s-1, Z=20
title BLR model, density 1e11 cm-3, flux of H-ion phots 1e20 cm2 s-1, Z=20
c
c commands controlling continuum =========
table agn
phi(h) 20
c
c commands for density & abundances =========
hden 11
abundances starburst 20
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate convergence
c these are to try to speed things up
init "c84.ini"
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch cooling "blr_n11_p20_Z20.col"
punch heating "blr_n11_p20_Z20.het"
punch element calcium "blr_n11_p20_Z20.cal"
punch overview "blr_n11_p20_Z20.ovr"
punch dr "blr_n11_p20_Z20.dr"
punch temperature history "blr_n11_p20_Z20.tem"
c
c blr_n11_p20_Z20.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n12_p19.in
BLR model, density 1e12 cm-3, flux of H-ion phots 1e19 cm2 s-1
title BLR model, density 1e12 cm-3, flux of H-ion phots 1e19 cm2 s-1
c
c commands for density & abundances =========
c
c commands controlling continuum =========
table agn
phi(h) 19
c
c commands controlling geometry =========
hden 12
stop column density 23
c
c other commands for details =========
iterate convergence
c these are to try to speed things up
init "c84.ini"
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n12_p19.ovr" last
punch dr "blr_n12_p19.dr" last
c
c
c blr_n12_p19.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n12_p19_Z20.in
BLR model, density 1e12 cm-3, flux of H-ion phots 1e19 cm2 s-1, Z=20
title BLR model, density 1e12 cm-3, flux of H-ion phots 1e19 cm2 s-1, Z=20
c
c commands controlling continuum =========
table agn
phi(h) 19
c
c commands for density & abundances =========
hden 12
abundances starburst 20
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
failures 5
iterate convergence
c these are to try to speed things up
init "c84.ini"
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch convergence reason "blr_n12_p19_Z20.cvr"
punch overview "blr_n12_p19_Z20.ovr" last
punch dr "blr_n12_p19_Z20.dr" last
c
c blr_n12_p19_Z20.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n13_p18.in
BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1
title BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1
c
c this is a very low ionization cloud
c the conditions, and some lines, are surprisingly sensitive
c to the treatment of hydrogen molecules
c
c commands controlling continuum =========
table agn
phi(h) 18
c
c commands for density & abundances =========
hden 13
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate to convergence
c
c commands controlling output =========
normalize to "totl" 1216 100
init "c84.ini"
no level2 lines
punch overview "blr_n13_p18.ovr" last
punch dr "blr_n13_p18.dr" last
punch convergence reason "blr_n13_p18.cvr"
c
c blr_n13_p18.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n13_p18_Z20.in
BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20
title BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20
c
c this is a very low ionization cloud
c the conditions, and some lines, are surprisingly sensitive
c to the treatment of hydrogen molecules
c
c commands controlling continuum =========
phi(h) 18
table agn
c
c commands for density & abundances =========
hden 13
abundances starburst 20
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
c this sim needs roughly 10 iterations to converge
iterate convergence limit 15
c these are to try to speed things up
init "c84.ini"
no level2 lines
c abort if any failures occur - we want to pass this test
failures 1
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n13_p18_Z20.ovr" last
punch overview "blr_n13_p18_Z20.ovr1"
punch dr "blr_n13_p18_Z20.dr" last
punch convergence reason "blr_n13_p18_Z20.cvr"
c
c blr_n13_p18_Z20.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n13_p22.in
BLR model, density 1e13 cm-3, flux of H-ion phots 1e22 cm2 s-1
title BLR model, density 1e13 cm-3, flux of H-ion phots 1e22 cm2 s-1
c
c commands controlling continuum =========
table agn
phi(h) 22
c
c commands for density & abundances =========
hden 13
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate convergence
c these are to try to speed things up
init "c84.ini"
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n13_p22.ovr" last
punch dr "blr_n13_p22.dr" last
c
c blr_n13_p22.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_n13_p22_Z20.in
BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20
title BLR model, density 1e13 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20
c
c commands controlling continuum =========
table agn
phi(h) 22
c
c commands for density & abundances =========
hden 13
abundances starburst 20
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate convergence
c these are to try to speed things up
init "c84.ini"
no level2 lines
c
c commands controlling output =========
normalize to "totl" 1216 100
punch overview "blr_n13_p22_Z20.ovr"
punch dr "blr_n13_p22_Z20.dr"
punch heating "blr_n13_p22_Z20.het"
punch cooling "blr_n13_p22_Z20.col"
c
c blr_n13_p22_Z20.in
c class blr
c ========================================
This is one of the 5 models that sample the LOC plane.
blr_nf84.in
early model of blr
title early model of blr
c
c commands controlling continuum =========
ionization parameter -1.92
power law -1 100 0.01
c
c commands for density & abundances =========
hden 9.5
init, file="ism.ini"
abundances cameron
stop lyman optical 6
c
c commands controlling geometry =========
c
c other commands for details =========
constant gas pressure
iterate to convergence
c
c commands controlling output =========
print diffuse continuum
plot continuum
punch overview last "blr_nf84.ovr"
punch dr last "blr_nf84.dr"
punch results last "blr_nf84.rlt"
c
table lines "LineList_strong.dat"
c blr_nf84.in
c class blr
c ========================================
c
This is an example of a "conventional" BLR calculation. The parameters are
similar to those of Table 1 of Netzer and Ferland (1984). Notice that the
ratio of Lyalpha to Hbeta ratio is much larger than observed.
>>refer blr model Netzer, H., & Ferland, G. J. 1984, PASP, 96, 593
blr_nf84_45deg.in
early model of BLR, with illumination at 45 degree angle
title early model of BLR, with illumination at 45 degree angle
c
c commands controlling continuum =========
ionization parameter -1.92
power law -1 100 0.01
c this continuum is a beam coming in at 45 degrees
illuminate 45 degrees
c
c commands for density & abundances =========
hden 9.5
abundances cameron
stop lyman optical 6
c
c commands controlling geometry =========
c
c other commands for details =========
init file "ism.ini"
constant gas pressure
iterate to convergence
c
c commands controlling output =========
print diffuse continuum
plot continuum
punch overview "blr_nf84_45deg.ovr"
punch dr "blr_nf84_45deg.dr"
punch results last "blr_nf84_45deg.rlt"
c
c blr_nf84_45deg.in
c class blr
c ========================================
c
This is an example of a "conventional" BLR calculation. The parameters are
similar to those of Table 1 of Netzer and Ferland (1984). Notice that the
ratio of Lyalpha to Hbeta ratio is much larger than observed.
>>refer blr model Netzer, H., & Ferland, G. J. 1984, PASP, 96, 593
blr_rnfa.in
table 1 of Rees et al. ApJ 347, 648
title table 1 of Rees et al. ApJ 347, 648
c
c commands controlling continuum =========
table agn
ionization parameter -2
c
c commands for density & abundances =========
hden 10
init file="c84.ini"
abundances he -1 c -3.328 n -4.0088 o -3.0809 ne -4 na -20
continue mg -4.3768 al -5.5686 si -4.36653 s -4.76955
continue a -5.4202 ca -5.6383 fe -4.4815 ni -20
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate to convergence
c
c commands controlling output =========
normalise "H 1" 1216 = 100
print line faint -1
punch overview last "blr_rnfa.ovr"
punch continuum last "blr_rnfa.con" units microns
punch transmitted continuum last "blr_rnfa.trn" units rydbergs
punch element nitrogen last "blr_rnfa.nit"
punch dr last "blr_rnfa.dr"
punch results last "blr_rnfa.rlt"
c
c blr_rnfa.in
c class blr
c ========================================
c
This is the lower density cloud computed in Rees et al. (1989).
Table 1 of that paper lists the predictions, which were a mean of those
of Hagai Netzer's ION and roughly version 76 of CLOUDY. The lines
are generally still in good agreement with the predictions of that paper.
In particular the changes in the line fluxes shown in Figure 1
of that paper are reproduced quite well.
blr_rnfb.in
table 1 of Rees et al. ApJ 347, 648
title table 1 of Rees et al. ApJ 347, 648
c
c blr model from
c
c commands controlling continuum =========
table agn
ionization parameter -2
c
c commands for density & abundances =========
hden 12
init file="c84.ini"
abundances he -1 c -3.328 n -4.0088 o -3.0809 ne -4 na -20
continue mg -4.3768 al -5.5686 si -4.36653 s -4.76955
continue a -5.4202 ca -5.6383 fe -4.4815 ni -20
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
iterate to convergence
c
c commands controlling output =========
print line faint -1
punch overview last "blr_rnfb.ovr"
punch element nitrogen last "blr_rnfb.nit"
punch dr last "blr_rnfb.dr"
punch results last "blr_rnfb.rlt"
normalise "H 1" 1216 = 91.8
c
c blr_rnfb.in
c class blr
c ========================================
c
This is a very dense cloud, and was computed in Rees et al. (1989).
Table 1 of that paper lists the predictions, which were a mean of those of
Hagai Netzer's ION and roughly version 76 of CLOUDY. The lines are
generally still in good agreement with the predictions of that paper.
In particular the changes in the line fluxes shown in Figure 1 of that
paper are reproduced quite well. The fluxes of Lya and Hb are not
reproduced with great precision by this model because of changes in
collision rates for hydrogen and especially the form of the
escape probability function for subordinate lines. As Figure 1 of
RNF showed the line intensities are very sensitive to density for
these parameters.
c6_caseb.in
c6_caseb C VI case B
title c6_caseb C VI case B
c
c commands controlling continuum =========
table agn
ionization parameter 3
c
c commands for density & abundances =========
hden = 8
init file="ism.ini"
element carbon abundance 1
set eden 9
c
c commands controlling geometry =========
set dr -18
stop zone 1
c
c other commands for details =========
case b hummer no photoionization
no induced processes (OK)
constant temperature = 1.e5
atom h-like carbon levels large
c
c commands controlling output =========
punch results "c6_caseb.rlt"
c
c c6_caseb.in
c class limit
c ========================================
c
This test case compares the predictions of the multi-level H-like CVI atom
with the Storey and Hummer (1995) results. The set dr command sets the
zone thickness to 1 cm. The case b command sets Lyman line optical depths
to very large values.
Checks:
- Neutral fractions
- H* emissivity
- Relative line intensities
coll_coronal.in
model of active region of solar corona
title model of active region of solar corona
c
c commands controlling continuum =========
c temperature of corona
coronal equilibrium, 2.5e6 K
c assume the sun is an Atlas atmosphere at 5700K
blackbody 5770
luminosity solar linear 1
c
c commands for density & abundances =========
hden log=10
c
c commands controlling geometry =========
c its height above photosphere and thickness
radius 11 thickness = 10
c the gas fully covers the center of symmetry
sphere
c
c other commands for details =========
iterate
c
c commands controlling output =========
print last iteration
normalize to "o 8" 18.97A
punch dr "coll_coronal.dr" last
punch continuum "coll_coronal.con" last units angstroms
c
c coll_coronal.in
c class coll
c ========================================
c
This is a rough model of the solar corona.
The test checks that the table star atlas and coronal equilibrium commands work.
The gas is predominantly collisionally ionized.
coll_heat_only.in
test code in limit where ONLY mechanical heating is present
title test code in limit where ONLY mechanical heating is present
c
c commands controlling continuum =========
no photoionization
blackbody 5000 STE
c
c commands controlling geometry =========
stop column density 15
c
c commands for density & abundances =========
hden 10
c
c other commands for details =========
hextra -1.5
iterate
c
c commands controlling output =========
c output files
punch overview "coll_heat_only.ovr"
punch heating "coll_heat_only.het"
punch cooling "coll_heat_only.col"
punch continuum "coll_heat_only.con" units angstroms
c
c
c coll_heat_only.in
c class coronal
c ========================================
This test is an optically thin collisionally ionized gas with no
photoionization at all. Before c9601 the code would not do the ionization
stage trimming correctly. Original bug caught by Ronnie Hoogerweft.
coll_t4.in
coronal equilibrium at 10^4 K
title coronal equilibrium at 10^4 K
c
c commands controlling continuum =========
coronal 4
c
c other commands for details =========
iterate
c
c commands controlling geometry =========
stop column density 15
c
c commands for density & abundances =========
hden 10
c
c commands controlling output =========
punch overview "coll_t4.ovr"
punch continuum "coll_t4.con" units angstroms
c
c
c coll_t4.in
c class coronal
c ========================================
this tests conditions of collisional equilibrium at low densities
coll_t5.in
coronal equilibrium at 10^5 K
title coronal equilibrium at 10^5 K
c
c commands controlling continuum =========
coronal 5
c
c commands controlling geometry =========
stop column density 15
c
c commands for density & abundances =========
hden 10
c
c other commands for details =========
iterate
c
c
c commands controlling output =========
c normalize to strong HeII line
normalize to "He 2" 303.8A
punch overview "coll_t5.ovr"
punch continuum "coll_t5.con" units angstroms
c
table lines "LineList_He_like.dat"
c coll_t5.in
c class coronal
c ========================================
This is a test collisional ionization equilibrium
at 100000 K.
coll_t6.in
coronal equilibrium at 10^6 K
title coronal equilibrium at 10^6 K
c
c commands controlling continuum =========
coronal 6
c
c commands controlling geometry =========
stop column density 15
c
c commands for density & abundances =========
hden 10
c
c other commands for details =========
iterate
c
c commands controlling output =========
normalize to "c 5" 40.27A
c output files
punch overview "coll_t6.ovr"
punch continuum "coll_t6.con" units angstroms
c
c coll_t6.in
c class coronal
c ========================================
This test is an optically thin collisionally ionized gas.
coll_t7.in
coronal equilibrium at 10^7 K
title coronal equilibrium at 10^7 K
c
c commands controlling continuum =========
coronal 7
c
c commands controlling geometry =========
stop column density 15
c
c commands for density & abundances =========
hden 10
c
c other commands for details =========
c
normalize to "o 8" 18.97A
c commands controlling output =========
punch overview "coll_t7.ovr"
punch continuum "coll_t7.con" units angstroms
punch cooling "coll_t7.col"
c
c coll_t7.in
c class coronal
c ========================================
Test with only collisional ionization at a
high gas high temperature.
dynamics_orion_flow.in
Orion nebula blister with wind
title Orion nebula blister with wind
c
c commands controlling continuum =========
blackbody 40000
phi(h) 13.0
brems 6
phi(h) 10
c
c commands for density & abundances =========
hden 4
abundances hii region no grains
grains orion no qheat single
c
c commands controlling geometry =========
wind -7 km/s advection
c this is the account for continued matter beyond i front
double optical depths
sphere
stop AV 15
stop temperature linear 5
c
c other commands for details =========
iterate 3
magnetic field -4
cosmic rays background
c this has no effect on dynamics, but is to desaturate the lines as per flow
turbulence 8 km/sec no pressure
c try to speed up simulation without hurting dynamics
init "fast.ini"
c use Bakes & Tielens heating
set PAH Bakes
* init "c84.ini"
* failures 5
c
c commands controlling output =========
c want to print line intensities are surface brightness, per arcsec^2
print lines surface brightness arcsec
print diffuse continuum
print line faint -1.5
punch overview "dynamics_orion_flow.ovr" no hash
punch hydrogen 21 cm "dynamics_orion_flow.21cm" no hash
punch dr "dynamics_orion_flow.dr" no hash
punch continuum last "dynamics_orion_flow.con" units microns no hash
punch pressure "dynamics_orion_flow.pre" no hash
punch heating "dynamics_orion_flow.het" no hash
punch cooling "dynamics_orion_flow.col" no hash
punch molecules "dynamics_orion_flow.mol" no hash
punch wind "dynamics_orion_flow.wnd" last
c
c dynamics_orion_flow.in
c class wind hii
c ========================================
c
This is a model similar in spirit to the blister geometry HII region
model computed by Baldwin et al. (1991), but with a D-critical flow.
dynamics_veryfast.in
very fast dusty windy model
title very fast dusty windy model
c
c commands controlling continuum =========
table agn
ionization parameter -2
c
c commands for density & abundances =========
hden 9
c this is a set of commands to speed up calc
init file "fast.ini"
c
c commands controlling geometry =========
stop thickness 8.4
c this is a subsonic wind with advection
wind -5 advection
c
c other commands for details =========
c absolutely no convergence problems should occur
failures 1
c this degrades the continuum resolution by factor of 3, to
c speed up calc
set continuum resolution 3
iterate
c
c commands controlling output =========
punch pressure last "dynamics_veryfast.pre"
punch total opacity last "dynamics_veryfast.opc"
punch ionizing continuum last "dynamics_veryfast.ion"
punch continuum last "dynamics_veryfast.con"
c
c
c dynamics_veryfast.in
c class dynamics
c ========================================
c
This is meant to be a very fast calculation to
use when running extensive debug-enabled runtimes.
dynamics_wind.in
test of equations of motion in a wind
title test of equations of motion in a wind
c
c test of wind code
c radiative acceleration (e- only) is 9.54E-7 cm s^-2
c terminal velocity (e- only) is 7.6 km s^-1
c
c commands controlling continuum =========
table agn
luminosiy (total) 45
c
c commands for density & abundances =========
hden 4
c this only includes H and He
init file="hheonly.ini"
c
c commands controlling geometry =========
radius (parsecs) 1
stop thickness (parsecs) -1
wind 0.1
c
c other commands for details =========
c this will speed things up a bit
no level2
no radiation pressure
constant temperature 8
c
c commands controlling output =========
punch assets "dynamics_wind.asr"
punch pressure "dynamics_wind.pre"
punch dr "dynamics_wind.dr"
c
c dynamics_wind.in
c class limit dynamics
c ========================================
c
This tests the management of the radiative acceleration of an
electron scattering wind, and the resulting velocity. The parameters
were chosen so that electron scattering is the dominant opacity source,
so that the equations can be solved both numerically (in the example)
and analytically (the expected solution given above). In a realistic wind
the gas would be more neutral and line driving would dominate. The actual
acceleration is slightly below that due to Thomson scattering alone because
the gas opacity at high energies is slightly below Thomson.
Checks:
- The radiative acceleration is correct (e- 9.543910-7 cm s-2).
- The terminal velocity should be 7.57 km s-1.
- Force multiplier near unity (no line driving since so highly ionized).
- Thickness of cloud correct (R-Ro + dr/2 should be 3.086391017 cm).
feii_hin.in
test feii in high density limit
title test feii in high density limit
c
c commands controlling continuum =========
black 25000 K
ionization -10
c
c commands for density & abundances =========
hden 18
eden 18
c
c commands controlling geometry =========
constant tempera 4
stop zone 1
c
c other commands for details =========
atom feii
iterate
c want to test collision only limit
no induced processes
c
c commands controlling output =========
print last iteration
c
c
c feii_hin.in
c class limit
c ========================================
c
This checks that, at high particle densitites, in which
the gas should be in collisional
equilibrium, the level populations of the model Fe+ ion
go to the proper values, where the departure coefficients
are all equal to unity.
feii_hirad.in
feii in case of high radiation density limit
title feii in case of high radiation density limit
c
c commands controlling continuum =========
black 8000 K lte
c
c commands for density & abundances =========
hden 9
element iron abundance 2
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
atom feii
constant tempera 8000
iterate
c
c commands controlling output =========
print last iteration
punch feii relative populations range 0 200 "feii_hirad.lv1" last
punch feii relative populations range 201 371 "feii_hirad.lv2" last
c
c
c feii_hirad.in
c class limit
c ========================================
c
This checks that, at high radiation densitites, in which
the gas is irradiated by a blackbody in strict thermodynamic
equilibrium, the level populations of the model Fe+ ion
go to the proper values, where the departure coefficients
are all equal to unity.
feii_pump.in
test feii in continuum pumped limit
title test feii in continuum pumped limit
c
c commands controlling continuum =========
black 25000 K
ionization -5
c
c commands for density & abundances =========
hden 10
abundances old solar 84
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
constant tempera 4
atom feii
c
c commands controlling output =========
punch feii continuum "feii_pump.con"
c
c feii_pump.in
This is a constant temperature low ionization
cloud, with BLR-like densities, which includes the
large FeII atom. The tests check on the emission
predicted in the Fe II bands.
This model tests the large FeII model in the optically thin,
continuum pumped limit. The zone thickness is set to a
small value (1 cm) so that full continuum hits atom.
feii_ste.in
thermal equilibrium of FeII in STE limit
title thermal equilibrium of FeII in STE limit
c
c commands controlling continuum =========
black 8000 K lte
hden 9
c
c commands for density & abundances =========
element iron abundance 2
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
atom feii
iterate
c
c commands controlling output =========
punch feii continuum last "feii_ste.con"
punch feii column density last "feii_ste.col"
punch feii populations all last "feii_ste.pop"
print last iteration
c
c feii_ste.in
c class function
c ========================================
c
This model has a very high iron abundance, 100x H, and most Fe is
in the form of Fe+. It is irradiated by a blackbody in strict
thermodynamic equilibrium. The mode checks that the temperature
of the gas is equal to the radiation temperature, to confirm
that the thermal properties of the model FeII atom obey thermodynamics.
func_abund_fluc.in
this runs the standard, one command, test, which contains many asserts
c commands controlling continuum =========
blackbody 5
ionization -2
c
c commands for density & abundances =========
hden 3
c this varies the abundnaces between default and 1e-4 times default,
c over a 1e10 cm period - we will compute half a cycle
fluctuatsion abundances 10.3 0 -4
c
c commands controlling geometry =========
c stop at half a cycle so pick off lowest abundances
stop thickness 10
c
c other commands for details =========
c speed things up a bit
init "c84.ini"
c
c commands controlling output =========
print column density
punch overview "func_abund_fluc.ovr"
punch abundances "func_abund_fluc.abn"
c
c
c func_abund_fluc.in
c class functionality
c ========================================
c
This checks that the variable abundances option still works
func_distance.in
check that distance and "print flux earth" commands work
title check that distance and "print flux earth" commands work
c
c commands controlling continuum =========
c luminosity in ionizing radiation, this and distance
c were chosen so that flux at earth assocated with total
c luminoisity of object is 1 erg cm-2 s-1
luminosity 40
laser 2
c
c commands for density & abundances =========
hden 4
init file "honly.ini"
c
c commands controlling geometry =========
stop zone 1
radius 17
c
c other commands for details =========
distance 19.450395
constant temperature 4
c
c commands controlling output =========
c give flux at Earth
print lines, flux at Earth
c
c func_distance.in
c class function
c ========================================
c
Normally the code predicts the intensity or luminosity of
the emission lines. This test confirms that it can predict
the flux recieved at the Earth instead. The model is the simplest
and fastest that can be computed - a H-only constant temperature
single zone. The total luminosity is set to 1e40 erg/s, and the
ionization source is a laser at 2 ryd. With these set, the
total luminosity in ionizing radiation, the total luminoisity in
the incident continuum (the emission line labeled "Inci 0") will be 1e40.
The code will predict the flux at the Earth if both the distance to the
oject is specified with the distance command, and this is requested with
the print flux earth command. The distance was chosen so that the
total flux at the Earth will be 1 erg/s. This is asserted at the
end of the calculation.
func_globule.in
test of globule command
title test of globule command
c
c commands controlling continuum =========
black 40000 K
ionization parameter -1.5
c
c commands for density & abundances =========
globule initial density=5 scale depth=14
init file "ism.ini"
abundances hii region no qheat
c
c commands controlling geometry =========
c
c other commands for details =========
failures 2
no diffuse line pumping
c
c commands controlling output =========
print line faint -1.5
punch overview last "func_globule.ovr"
punch results last "func_globule.rlt"
punch dr "func_globule.dr" last
punch heating "func_globule.het" last
punch cooling "func_globule.col" last
c
c func_globule.in
c class function
c ========================================
c
This model uses the globule command, and tests that the zoning logic
works for this extreme case, and that the code is able to converge
the globule model.
func_grid.in
test grid command and fits output
title test grid command and fits output
c
c commands controlling continuum =========
c following gives range for variation of blackbody
c initial value ignored, will vary through above range
c value on blackbody command is ignored but needed to pass parser
blackbody 1.e5 K vary
grid, range from 4 to 6 with increment 1
luminosity 39
c
c commands for density & abundances =========
hden 0 vary
grid, range from 2 to 3 with increment 1
init "ism.ini"
c
c commands controlling geometry =========
stop zone 1
radius 4 parsecs
c
c other commands for details =========
constant temperature 4
c
c commands controlling output =========
punch xspec atable "func_grid.fit"
punch overview "func_grid.ovr" no clobber
punch linelist out to "func_grid.lin" read "LineList_HII.dat" no clobber
punch grid "func_grid.grd"
punch total opacity "func_grid.opc"
c
c
c func_grid.in
c class functionality
c ========================================
c
this model tests the grid command and fits output.
func_hotgas_coolstar.in
test very soft continuum, very hot gas
title test very soft continuum, very hot gas
c
c commands controlling continuum =========
c put in the cosmic background as the only continuum source
CMB
c gas has constant temperature of 1e6 K
constant temperature, t=6
c
c commands for density & abundances =========
hden 10.
c
c commands controlling geometry =========
stop zone 1
c specify a thin cell of gas - 1 cm thick
set dr 0
c
c other commands for details =========
c iterate since gas is optically thin
iterate
c
c commands controlling output =========
print last iteration
c
c func_hotgas_coolstar.in
c class limit
c ========================================
c
This is a test where the continuum source is very soft,
and does not extend into energies where the
code needs to work. There are special cases
used in this situation, for continuum addressing,
so this checks whether those still function.
func_ion_increase.in
test model where ionization increases with depth
title test model where ionization increases with depth
c
c commands controlling continuum =========
ionization parameter -4
table agn
c
c commands for density & abundances =========
hden 7 -4
init "hheonly.ini"
element oxygen on
element iron on
c
c commands controlling geometry =========
sphere
filling factor -5
radius 17 20
c
c other commands for details =========
constant temperature 4
c
c commands controlling output =========
punch dr "func_ion_increase.dr"
punch overveiw "func_ion_increase.ovr"
punch element oxygen "func_ion_increase.oxy"
punch element iron "func_ion_increase.fe"
c
c
c func_ion_increase.in
c class functionality
c ========================================
c
this model tests the case where the ionization increases with depth.
The density is a fast powerlaw to insure this.
func_map.in
map of heating vs cooling
title map of heating vs cooling
c
c commands controlling continuum =========
table agn
ionization parameter -2.5
c
c commands for density & abundances =========
hden 0
abundances old solar 84
c
c commands controlling geometry =========
set dr 0
stop lyman optical depth -7
stop zone 1
c
c other commands for details =========
c
c commands controlling output =========
c change plot to punch to generate plot for hazy
punch map file="func_map.dta"
set nmaps 100
c
c map.in
c class func limit
c ========================================
c
This is a test of the continuity of the code over a very large range
of temperature. It was used to produce one of the thermal maps shown in Hazy.
Checks:
- No breaks in the heating and cooling curves where
various approximations change.
func_set_ion.in
test impact of setting ionization
title test impact of setting ionization
c
c commands for density & abundances =========
init "ism.ini"
hden 0
element hydrogen ionization 0 1
element helium ionization 0 1 0
element carbon ionization 0 1
element oxygen ionization 0 0.1 0.2 0.3 0.2 0.1 0.1
c
c commands controlling continuum =========
blackbody 5
ionization parameter -2
c
c commands controlling geometry =========
c force model to be 1 cm thick
set dr 0
c
c other commands for details =========
constant temperature 4
stop zone 1
c
c commands controlling output =========
normalize to "O 3" 5007
punch dr "func_set_ion.dr"
c
c
c func_set_ion.in
c class function
c ========================================
c
this script exercies the option to specify the ionization of a species
func_stopline.in
test stop line command
title test stop line command
c
c commands controlling continuum =========
ionization parameter -2
blackbody 50000
c
c commands for density & abundances =========
hden 2
abundances old solar 84
init file "ism.ini"
c
c commands controlling geometry =========
c stop calc at line ratio, the point of this test
stop line "c 2" 157.6m reaches 0.002 relative to "o 3" 5007
c
c other commands for details =========
constant temperature 10000 K
c
c commands controlling output =========
normalise to "o 3" 5007
punch dr "func_stopline.dr" last
c
c func_stopline.in
c class functionality
c ========================================
c
this is an example of a simple calculation that stops when
a certain line ratio is reached
func_t10.in
test very soft continuum, very hot gas
title test very soft continuum, very hot gas
c
c commands controlling continuum =========
table agn
ionization parameters -2
c
c commands for density & abundances =========
hden 10.
c
c commands controlling geometry =========
stop zone 1
c specify a thin cell of gas - 1 cm thick
set dr 0
c
c other commands for details =========
c set constant temperature to high-T limit
constant temperature, t=10 log
c iterate since gas is optically thin
iterate
c
c commands controlling output =========
print last iteration
c
c
c func_t10.in
c class limit
c ========================================
c
This is a test of the highest temperature the code can do.
func_t3.in
test low temperature limit of code, 3K
title test low temperature limit of code, 3K
c
c commands controlling continuum =========
table agn
ionization parameters -2
c
c commands for density & abundances =========
hden 10.
c
c commands controlling geometry =========
stop zone 1
c specify a thin cell of gas - 1 cm thick
set dr 0
c
c other commands for details =========
c set constant temperature to low-T limit
constant temperature, t=3.0 linear
c iterate since gas is optically thin
iterate
c
c commands controlling output =========
print last iteration
c
c
c func_t3.in
c class limit
c ========================================
c
This is a test of the lowest temperature the code can do.
It runs a constant temperature of 3K
func_test.in
this runs the standard, one command, test, which contains many asserts
c commands controlling continuum =========
test
c commands for density & abundances =========
c commands controlling geometry =========
c other commands for details =========
c include this in one test to make sure it works
set continuum shield ferland
c check that LineList_BLR.dat still ok
table lines
c
c commands controlling output =========
c following will print physical constants used by the code
print constants
c print the past to the data directory
print path
c following prints column densities
print column density
punch overview "func_test.ovr"
c func_test.in
c class functionality
c ========================================
c
This runs the "test" command, which include several asserts.
func_trans_punch.in
first of func_trans_punch/transread pair, punch continuum
title first of func_trans_punch/transread pair, punch continuum
c
c commands controlling continuum =========
ionization parameter -1
table agn
c
c commands for density & abundances =========
hden 10
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
iterate
print last iteration
c
c commands controlling output =========
print incident continuum
punch ionization last "func_trans_punch.ion"
punch continuum last "func_trans_punch.con"
c this continuum will be used by transread.in
punch last transmitted continuum "func_trans_punch.trn"
c
init file "trans.dat"
c
c func_trans_punch.in
c class functionality
c ========================================
c
func_trans_punch.in and transread.in are a pair of
tests that check that the code can punch a transmitted
continuum then read it.
func_trans_read.in
second of transpunch/transread pair, used transmitted continuum
title second of transpunch/transread pair, used transmitted continuum
c the transmitted continuum produced by func_trans_punch"
c
c commands controlling continuum =========
ionization parameter -1
table read "func_trans_punch.trn"
brems 6
ionization parameter -6
c
c commands for density & abundances =========
hden 10
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
iterate
c
c commands controlling output =========
init file "trans.dat"
print last iteration
print incident continuum
punch ionization last "func_trans_read.ion"
c
c
c func_trans_read.in
c class functionality
c ========================================
c
func_trans_punch.in and func_trans_read.in are a pair of
tests that check that the code can punch a transmitted
continuum then read it.
grains_hot.in
test temperature of gas and dust in high energy density environment
title test temperature of gas and dust in high energy density environment
c
c commands controlling continuum =========
table power law
energy density temp 550
c
c commands for density & abundances =========
hden 5.0
init "ism.ini"
abundances orion no grains
grains orion no qheat
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
c
c commands controlling output =========
normalize to "O 8" 18.97A
punch overview "grains_hot.ovr" last
c
c
c grains_hot.in
c class limit
c ========================================
c
This tests the grains in an extreme condition - irradiation by an AGN
near the illuminated face of the molecular torus. The gas is
predominantly heated by the grains.
grains_lte.in
check that grains equilibriate at correct temp in ste limit
title check that grains equilibriate at correct temp in ste limit
c this also tests that the code works when H He ionization is ZERO
c
c commands controlling continuum =========
black 3 lte
c
c commands for density & abundances =========
hden 5
eden 5
init file "hheonly.ini"
element oxygen on
grain abund 1 type 1
grain abund 1 type 2
grain abund 1 type 3
grain abund 1 type 4
grain abund 1 type 5
grain abund 1 type 6
grain abund 1 type 7
grain abund 1 type 8
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate
c the equilibrium temp would not be correct since small H_2 molecule
c does not go to ste, being an equivalent two-level system
constant temperature 3
c
c commands controlling output =========
normalize to "O 1" 63.17m
c
c grains_lte.in
c class limit
c ========================================
c
This test radiations a set of grains with a true blackbody in
strict thermodynamic equilibrium. We expect the grains
(and everything else) to equilibriate at the blackbody temperature.
The gas temperature is forced to the radiation temperature because
the current molecule network (based on ISM approximations) does not
now go to LTE in the high radiation density limit. The calculation
asserts that all grain temperatures are very close to the radiation
temperature.
grains_qheat.in
cool atomic ISM with Si grain quantum heating
title cool atomic ISM with Si grain quantum heating
c
c commands controlling continuum =========
table ism
extinguish column=22 leak=0
c
c commands for density & abundances =========
hden 0
abundances old solar 84
constant density
c turns on grains with default abundance (log abundance ratio to ism = 0)
c and with quantum heating enabled by default (disabled with "no grains_qheat" option)
grains 0. "silicate_ism_10.opc"
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate
print last iteration
c this is included to that continuum pumping of H atom not important
case b
c
c commands controlling output =========
normalize to "C 2" 157.6m
print diffuse continuum
print line faint 1
punch grains temperature last "grains_qheat.qht"
punch continuum "grains_qheat.con" units microns
c
c
c grains_qheat.in
c class limit
c ========================================
c
the level 2 lines are left in, in this model. Many UV lines are
artifically too strong, for a physical reason. The calculated model
is optically thin, so continuum pumping is extremely important. In the real
ISM this process would occur but would be much weaker in the global
spectrum, because the lines would become optically thick, and
line shielding over a large column would limit the effect.
grains_temp.in
test all grain species temperature
title test all grain species temperature
c
c commands controlling continuum =========
black 50000
ionization parameter -2
c
c commands for density & abundances =========
hden 3
grains_temp "graphite_ism_01.opc"
grains_temp "silicate_ism_01.opc"
grains_temp "graphite_orion_01.opc"
grains_temp "silicate_orion_01.opc"
grains_temp "silicate_0m010.opc"
grains_temp "silicate_0m100.opc"
grains_temp "grey_ism_01.opc"
grains_temp "silicate_1m000.opc"
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
iterate
constant temper 4
c
c commands controlling output =========
print last iteration
c
c
c grains_temp.in
c class limit
c ========================================
c
This turns on all the grain species that the code knows about,
at least in its default state. A model of an ionized layer is
done and the asserts confirm the resulting grain temperatures.
grains_temp_all.in
test all grain species temperature
title test all grain species temperature
c
c commands controlling continuum =========
black 50000
ionization parameter -2
c
c commands for density & abundances =========
hden 3
c this turns on orion-style distributed grains
grains orion
c this turns on ism distributed grains
grains ism
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
constant temper 4
iterate
c
c commands controlling output =========
punch dr "grains_temp_all.dr" last
punch grain charge "grains_temp_all.grnchr" last
punch grain temperature "grains_temp_all.grntem" last
punch grain drift velocity "grains_temp_all.grndft" last
c
c grains_temp_all.in
c class limit
c ========================================
c
this is a relatively quick test of the new grains.
The Orion and ISM silicate and graphitic grains
are turned on and their equilibrium temperature
checked. This uses distributed grain sizes by default.
h2_cr.in
background cosmic ray ionization by suprathermal electrons only
title background cosmic ray ionization by suprathermal electrons only
c
c commands controlling continuum =========
c this is hot star continuum
black 30000
intensity -5. range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c background CR ionization rate should be 2.5e-17 s^-1,
c Tielens & Hollenbach 1985; McKee AstroPH 9901370
cosmic rays, background
c
c commands for density & abundances =========
hden 5
c
c commands controlling geometry =========
set dr 0
stop zone 5
stop temperature 400
c
c other commands for details =========
set h2 grain formation thermal pump
failures 5
atom h2
iterate
constant temperature 500K
c must turn off line photoexcitation
c since we want to test cr heating and ionization
no induced processes
c
c commands controlling output =========
print column density
print ages
punch overview "h2_cr.ovr" last
punch h2 lines "h2_cr.lin" last all
punch h2 rates "h2_cr.rat" last
punch h2 column density "h2_cr.col" last
c
//
c h2_cr.in
c class limit
c ========================================
c
This test conditions of cosmic ray ionization. Molecules and charge transfer
are disabled so that analytical estimates can be made.
h2_cr_grains.in
background cosmic ray ionization by suprathermal electrons only
title background cosmic ray ionization by suprathermal electrons only
c
c commands controlling continuum =========
c this is hot star continuum
black 30000
intensity -5. range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c background CR ionization rate should be 2.5e-17 s^-1,
c Tielens & Hollenbach 1985; McKee AstroPH 9901370
cosmic rays, background
c
c commands for density & abundances =========
abundances ism
hden 5
c
c commands controlling geometry =========
set dr 0
stop zone 5
stop temperature 400
c
c other commands for details =========
set h2 grain formation thermal pump
failures 5
atom h2
iterate
constant temperature 500K
c must turn off line photoexcitation
c since we want to test cr heating and ionization
no induced processes
c
c commands controlling output =========
normalize to "H2 " 2.121m
print column density
print ages
punch overview "h2_cr_grains.ovr" last
punch h2 lines "h2_cr_grains.lin" last all
punch h2 rates "h2_cr_grains.rat" last
punch h2 column density "h2_cr_grains.col" last
c
c h2_cr_grains.in
c class limit
c ========================================
c
This test conditions of cosmic ray ionization and grain formation pumping.
h2_hminus.in
H2 populations in H- dominated limit
title H2 populations in H- dominated limit
c this is one of a pair with h2_solomon.in
c
c commands controlling continuum =========
c this is hot star continuum
black 30000
intensity -5. range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c background CR ionization rate should be 2.5e-17 s^-1,
c Tielens & Hollenbach 1985; McKee AstroPH 9901370
cosmic rays, background
c
c commands for density & abundances =========
abundances ism no grains
hden 5
c
c commands controlling geometry =========
set dr 0
stop zone 5
stop temperature 400
c
c other commands for details =========
set h2 grain formation thermal pump
atom h2
iterate
constant temperature 500K
c
c commands controlling output =========
print column density
print ages
punch overview "h2_hminus.ovr" last
punch h2 lines "h2_hminus.lin" last all
punch h2 rates "h2_hminus.rat" last
punch h2 column density "h2_hminus.col" last
c
c h2_hminus.in
c class limit
c ========================================
c
This tests H2 in limit of H- formation and solomon destruction
h2_pdr_leiden_f1.in
low density flux model 1 as defined in e-mail
title low density flux model 1 as defined in e-mail
atom h2 matrix all
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
c
c commands for density & abundances =========
c hydrogen density
hden 3.
grains ism 1.16 no qheat
c
c commands controlling geometry =========
c
c other commands for details =========
c turn on the large H2 model
atom h2
c turn down the number of convergence failures - there should be none
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the gas kinetic temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
c default normalization line is Hbeta, which is not produced by this gas
normalize to "C 2" 157.6m
print line faint -4
punch overview "h2_pdr_leiden_f1.ovr"
punch line emissivity "h2_pdr_leiden_f1.ems"
H2 2.121m
12CO 647.2m
C 2 157.6m
O 1 63.17m
Si 2 34.81m
end of lines
punch leiden lines "h2_pdr_leiden_f1.lin"
punch leiden "h2_pdr_leiden_f1.lei"
punch dr "h2_pdr_leiden_f1.dr"
punch molecules "h2_pdr_leiden_f1.mol"
punch grain physical conditions "h2_pdr_leiden_f1.grn"
punch overview "h2_pdr_leiden_f1.ovr"
punch H2 lines "h2_pdr_leiden_f1.h2lin" all
punch H2 column density "h2_pdr_leiden_f1.h2col"
punch H2 populations matrix zone "h2_pdr_leiden_f1.pop"
punch H2 rates "h2_pdr_leiden_f1.rat"
c
c
c h2_pdr_leiden_f1.in
c class pdr
c ========================================
c
h2_solomon.in
H2 populations in solomon dominated limit
title H2 populations in solomon dominated limit
c this is a pair with h2_hminus
c
c commands controlling continuum =========
c this is hot star continuum
black 30000
intensity -5. range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c background CR ionization rate should be 2.5e-17 s^-1,
c Tielens & Hollenbach 1985; McKee AstroPH 9901370
cosmic rays, background
c
c commands for density & abundances =========
abundances ism
hden 5
c
c commands controlling geometry =========
set dr 0
stop zone 5
stop temperature 400
c
c other commands for details =========
set h2 grain formation thermal pump
atom h2
iterate
constant temperature 500K
c
c commands controlling output =========
c don't print so many lines
print line faint 10
print line column
print line sort wavelength
print column density
print ages
punch overview "h2_solomon.ovr" last
punch h2 lines "h2_solomon.lin" last all
punch h2 rates "h2_solomon.rat" last
punch h2 column density "h2_solomon.col" last
c
c h2_solomon.in
c class limit
c ========================================
c
This test H2 in case of grain formation and solomon destruction
h2_t2000.in
test large H2 molecule in shock-like conditions
title test large H2 molecule in shock-like conditions
c
c commands controlling continuum =========
c there is a continuum but it is unimportant
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
black 30000
intensity -5 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
* cosmic rays, background
c
c commands for density & abundances =========
hden 6
grains ism, abundance log 0.16 single, no qheat
init file="ism.ini"
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl=-7 ar-8 fe -6.60
c
c commands controlling geometry =========
stop zone 10
c force zone thickess to 1 cm
set dr 0
c need to lower stop temperature since gas is cooler than 4000K
stop temperature 3
c
c other commands for details =========
failures 2
atom h2
constant temperature 2000
iterate
c this stop Lymnan line excitation of H, which is then photoionized
c by Balmer continuum
case b
c
c commands controlling output =========
normalize to "H2 " 2.121m
c
punch pdr "h2_t2000.pdr"
punch h2 rates "h2_t2000.rat"
punch h2 column density "h2_t2000.col"
punch h2 lines "h2_t2000.lin" all last
punch h2 populations "h2_t2000.pop" all last
punch raw continuum "h2_t2000.raw"
punch continuum "h2_t2000.con"
punch overview "h2_t2000.ovr"
c
c h2_t2000.in
c class limit
c ========================================
This is a collisionally dominated H2 simulations. the temperature
has been fixed at 2000K and the large molecule turned on. The calculation
checks the returned value of the ortho to para densities.
cosmic rays and the incident continuum are small, density is high,
so populations should be close to LTE.
h2_t500.in
test large H2 molecule in PDR-like conditions
title test large H2 molecule in PDR-like conditions
c
c commands controlling continuum =========
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
blackbody, t = 75 K
intensity -2.7 (total)
c this is hot star continuum
black 30000
intensity -5 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c
c commands for density & abundances =========
hden 6
grains ism, abundance log 0.16 single, no qheat
init file="ism.ini"
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl=-7 ar-8 fe -6.60
turbulence 2.7 km/sec
c
c commands controlling geometry =========
stop zone 10
c force zone thickess to 1 cm
sphere
set dr 0
c stop when gas is fully neutral
stop efrac -10
c stop when gas is cold
stop temperature 10 linear
c
c other commands for details =========
failures 2
atom h2
constant temperature 500
iterate
case b
cosmic rays background
c
c commands controlling output =========
normalize to "12CO" 863.0m
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print line faint -2
print ages
c
punch pdr "h2_t500.pdr"
punch h2 rates "h2_t500.rat"
punch h2 lines "h2_t500.lin" all
punch h2 column density "h2_t500.col"
punch h2 populations "h2_t500.pop"
punch raw continuum "h2_t500.raw"
punch continuum "h2_t500.con"
punch overview "h2_t500.ovr"
c
table lines "LineList_PDR_H2.dat"
c h2_t500.in
c class limit
c ========================================
this is a dense molecular gas with no cosmic rays and the incident
continuum set to a small value. the lower levels are in LTE and
the populations check this.
h_casea.in
case A
title case A
c
c Seaton, M.J, 1959 MN 119, 90,
c 4pi j(beta) 5.56E-26
c for this model total H-beta=4.745
c b(2)=3.73E-3 (3)=3.69E-2 (4)=0.091 (5)=0.145 (6)=0.193
c n.b. very different results if not l-mixed
c
c commands controlling continuum =========
black body, T = 50000 K
ionization parameter -2
c
c commands for density & abundances =========
hden = 13;
init file="honly.ini"
c
c commands controlling geometry =========
stop lyman optical depth -6
stop zone 1
set dr 0
c
c other commands for details =========
no induced processes (OK) ;N.B. case A is a fiction; no incuded transitions
atom h-like collisions off ; must turn off all collisions
constant temperature = 10000 K
iterate
atom h-like levels large
c
c commands controlling output =========
c this should make the predicted totl Hbeta near unity for case a
normalize to "CaBo" 4861 = 2.23
print populations h-like
print departure coefficients
print last iteration
c
c h_casea.in
c class limit
c ========================================
c
Case A is a mathematical fiction; when the Lyman lines are
optically thin continuum pumping must be important if the gas
is ionized. Fluorescence is turned off with the
no induced processes command. The density is set to a
very high value (1015 cm-3) so that the 2s-2p states are well l-mixed,
in keeping with standard case A assumptions. As a result,
collisional excitation would dominate the level populations,
and hydrogen collisions must be turned off with the hydrogen collisions off
command. The Ly* optical depth is set to a small value.
The set dr command sets the zone thickness to 1 cm.
The abundances are set to a very small value so that the
electron density is equal to the hydrogen density.
Checks:
- Departure coefficients for H, He levels
- Neutral fractions
- H* emissivity
h_caseb_lon.in
low density case b
title low density case b
c
c commands controlling continuum =========
laser 2
ionization -2
c
c commands for density & abundances =========
hden 0
init "honly.ini"
c
c commands controlling geometry =========
stop thickness 0
c
c other commands for details =========
case b, tau 5
iterate
constant temper 4
c
c commands controlling output =========
print last iteration
punch dr "h_caseb_lon.dr"
c
c h_caseb_lon.in
c class limit
c ========================================
c
This is a hydrogen-only model, very thin thickness,
It is case B and checks thetwo photon luminosity.
h_caseb_lot.in
log density case B, T=5000, log n=2
title log density case B, T=5000, log n=2
c
c commands controlling continuum =========
black body, T = 2.e5 K
ionization parameter -1
c
c commands for density & abundances =========
hden = 2
init file="hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature = 500
c want H and he fully stripped so we can check H He spectra
case b hummer no photoionization
no induced processes (OK)
iterate
c
c commands controlling output =========
print diffuse continua
normalize to "Ca B" 4861
c
c does not agree with Hummer & Storey, as explained by Ferguson & Ferland
c until 03 jun 02, used small atom, which had smaller differences
c turn on large atom to get best answer
atom h-like levels large
c
punch results "h_caseb_lot.rlt" last
c
c h_caseb_lot.in
c class limit
c ========================================
c
This tests the ionization and emission line spectrum for
H case B at a low density and temperature.
h_caseb_n8.in
h_caseb_n8 high density case B
title h_caseb_n8 high density case B
c
c commands controlling continuum =========
black body, T = 2.e5 K
ionization parameter -1
c
c commands for density & abundances =========
hden = 8
init file="hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
case b hummer no photoionization
no induced processes (OK)
constant temperature = 20000
atom h-like levels large
c
c commands controlling output =========
normalize to "Ca B" 4861
punch results "h_caseb_n8.rlt"
c
c h_caseb_n8.in
c class limit
c ========================================
c
This test case compares the predictions of the multi-level hydrogen atom
with the Storey and Hummer (1995) results. The set dr command sets the
zone thickness to 1 cm. The case b command sets Lyman line optical depths
to very large values.
Checks:
- Neutral fractions
- H* emissivity
- Relative line intensities
h_casebn2.in
log density case B, T=5000, log n=2
title log density case B, T=5000, log n=2
c
c commands controlling continuum =========
black body, T = 2.e5 K
ionization parameter -1
c
c commands for density & abundances =========
hden = 2
init file="hheonly.ini"
c
c commands controlling geometry =========
c this forces unit emissivity when combined with density
set dr -4
stop zone 1
c
c other commands for details =========
atom h-like levels large
iterate
c want H and he fully stripped so we can check H He spectra
case b hummer tau 5 no photoionization
no induced processes (OK)
constant temperature = 5000
c
c commands controlling output =========
normalize to "Ca B" 4861
print populations h-like
print last iteration
punch results "h_casebn2.rlt" last column
c
c h_casebn2.in
c class limit
c ========================================
c
This test case compares the predictions of the multi-level hydrogen atom
with the Storey and Hummer (1995) results. The set dr command sets the
zone thickness to 1 cm. The case b command sets Lyman line optical depths
to very large values.
Checks:
- Neutral fractions
- H* emissivity
- Relative intensities of hydrogen lines
h_casec.in
case C
title case C
c
c commands controlling continuum =========
black body, T = 50000 K
ionization parameter -2
c
c commands for density & abundances =========
hden = 5
init file="hheonly.ini"
c
c commands controlling geometry =========
set dr -10
stop zone 1
c
c other commands for details =========
iterate
atom h-like levels large
stop lyman optical depth -6
constant temperature = 10000 K
c
c commands controlling output =========
print last iteration
normalize to "Ca B" 4861
print h-like populations to 50
print departure coefficients
c
c h_casec.in
c class limit
c ========================================
c
This is Case C, what really happens when optically thin gas is
irradiated by a continuum with Lyman line continuum fluorescence allowed.
Checks:
- Departure coefficients for H, He levels
- Neutral fractions
- H* emissivity
h_induc.in
constant temper black body limit from Ferland and Rees 1988
title constant temper black body limit from Ferland and Rees 1988
c tests whether departure coef are forced to unity by induced processes
c
c commands controlling continuum =========
black body, t = 50000 lte
c
c commands for density & abundances =========
hden 10
init file="hheonly.ini"
c
c commands controlling geometry =========
stop lyman continuum optical depth -6
set dr -10
stop zone 1
c
c other commands for details =========
iterate
atom h-like levels large
constant temper 50000
c
c commands controlling output =========
print populations h-like
print heating
print departure coef
print line faint -1
print last
c
c h_induc.in
c class limit
c ========================================
c
This example tests whether induced processes force level populations
of hydrogen to LTE when they are irradiated by a
blackbody in strict thermodynamic equilibrium. The density is
low enough value for radiation to dominate the rate equations
coupling levels with each other and the continuum. The expectation
is for all departure coefficients to equal unity. A small amount
of grains are included to check that the grain thermal balance
is handled properly in this radiation-dominated limit.
Checks:
- Departure coefficients exactly unity.
- Grain temperatures are exactly 5*104 K.
h_lrg_atom.in
h_lrg_atom case B with largest possible H atom
title h_lrg_atom case B with largest possible H atom
c
c commands controlling geometry =========
constant temperature = 10000
stop zone 1
set dr 0
c
c commands controlling continuum =========
black body, T = 2e5 K
ionization parameter -2
c
c commands for density & abundances =========
hden = 4
init file="honly.ini"
c
c other commands for details =========
c this is the largest possible hydrogen atom
atom h-like levels limit
c this forces case B, with the storey & hummer assumptions
case b hummer
c
c commands controlling output =========
normalize to "Ca B" 4861
punch results "h_lrg_atom.rlt"
c
c the balmer series, n->2
c
c the Paschen series, n->3
c
c the Pfund series, n->4
c
c the n->5 series
c
c the Balmer continuum
c h_lrg_atom.in
c class limit
c ========================================
This does a simple case B calculation of hydrogen line emission.
A large H atom is used, and the results are compared with the
Storey & Hummer predictions.
h_lym_thin.in
H only optically thin in Lyman continuum
title H only optically thin in Lyman continuum
c
c commands controlling continuum =========
blackbody 40000
ionization -2
c
c commands controlling geometry =========
stop lyman continuum optical depth -2 %Lyman lines thick but continuum thin
atom h-like levels large
c
c commands for density & abundances =========
hden 2
init file "honly.ini"
c
c other commands for details =========
iterate
constant temper 10000
no level2
c
c commands controlling output =========
print diffuse continuum
print line inward
print last iteration
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "h_lym_thin.ovr" last
punch fine optical depths last "h_lym_thin.fin" range 0.9 1.02 every 1
punch optical depths last "h_lym_thin.opt"
punch continuum "h_lym_thin.con" last units microns
punch dr "h_lym_thin.dr" last
c
c lots of continuum emission points turned on with print diffuse continuum
c energies dominated by two photon emission not included since very
c sensitive to Lya transport
c honlyoutsp.in
c class limit
c ========================================
c
This is a pure hydrogen cloud that is optically thin in the Lyman continuum.
The asserts check the emission in several H I lines and continua.
This should be close to what really happens in a low column density
cloud exposed to a continuum source that does not have strong Lyman lines.
(The continuum source used is a pure blackbody, and so has no lines).
So this is an example of "Case C" emission
>>refer H case C Ferland, G.J. 1999, PASP, 111, 1524
h_otsopen.in
test ots, inward fractions for pure hydrogen, open geo, filling factor
title test ots, inward fractions for pure hydrogen, open geo, filling factor
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 5
init file "honly.ini"
c
c commands controlling geometry =========
stop efrac 0.05 %stop when 5 percent ionized
constant temper 7500
filling factor 0.01
c
c other commands for details =========
iterate
diffuse ots
atom h-like levels large
c not iterating does not affect the solution
no level2
no induced processes (OK)
c
c commands controlling output =========
print diffuse continuum
print line inward
punch overview "h_otsopen.ovr" last
punch dr "h_otsopen.dr" last
punch continuum "h_otsopen.con" last units angstrom
normalize to "Ca B" 4861 = 1 % normalize to exact result
c
c this is an open geometry, so do not expect H atom to exactly be case B
c due to Lyman line leakage. we iterated, so have a good estimate of
c inward and outward fractions, which we test
* 99jun26, error had been just under 2%, increased above it when optical
* depths of 2s and 2p were brought together, so increased error
c h_otsopen.in
c class limit
c ========================================
c
This tests the total emission from a hydrogen Stromgren sphere
using the OTS approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked.
h_otspp.in
plane parallel conservation and hydrogenic emission for pure hydrogen
title plane parallel conservation and hydrogenic emission for pure hydrogen
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 5
init file "honly.ini"
c
c commands controlling geometry =========
sphere static (OK)
stop efrac 0.05 %stop when 5 percent ionized
c
c other commands for details =========
constant temper 7500
no level2
no induced processes (OK)
c needed to prevent 2s / Bal from becoming optically thick
turbulence 20 km/s
diffuse ots
atom h-like levels large
c
c commands controlling output =========
print diffuse continuum
normalize to "Ca B" 4861 = 1 % normalize to exact result for comp ion struc
punch overview "h_otspp.ovr"
punch dr "h_otspp.dr"
c
c lots of continuum emission points turned on with print diffuse continuum
c energies dominated by two photon emission not included since very
c sensitive to Lya transport
c h_otspp.in
c class limit
c ========================================
c
This tests the total emission from a plane parallel pure hydrogen
Stromgren sphere using the OTS approximation. The conservation
of the total number of ionizing photons, and the emitted spectrum,
are all checked.
h_otssp.in
spherical conservation and hydrogenic emission for pure hydrogen
title spherical conservation and hydrogenic emission for pure hydrogen
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 5
init file "honly.ini"
c
c commands controlling geometry =========
radius 13
c not iterating does not affect the solution
sphere static (OK)
stop efrac 0.05 %stop when 5 percent ionized
c
c other commands for details =========
constant temper 7500
diffuse ots
atom h-like levels large
no level2
no induced processes (OK)
c
c commands controlling output =========
print diffuse continuum
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "h_otssp.ovr"
punch dr "h_otssp.dr"
c
c lots of continuum emission points turned on with print diffuse continuum
c energies dominated by two photon emission not included since very
c sensitive to Lya transport
c h_otssp.in
c class limit
c ========================================
c
This tests the total emission from a spherical pure hydrogen Stromgren sphere
using the OTS approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked.
h_outopen.in
test inward fractions, open geo, filling factor, 2-photon emission
title test inward fractions, open geo, filling factor, 2-photon emission
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 0
init file "honly.ini"
c
c commands controlling geometry =========
stop efrac 0.05 %stop when 5 percent ionized
filling 0.01
c
c other commands for details =========
diffuse outward
iterate
atom h-like levels large
no level2
no induced processes (OK)
constant temper 4
c
c commands controlling output =========
print diffuse continuum
c normally will not print out inward parts of lines
print line inward
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "h_outopen.ovr" last
punch dr "h_outopen.dr" last
c
c this is an open geometry, so do not expect H atom to exactly be case B
c due to Lyman line leakage. we iterated, so have a good estimate of
c inward and outward fractions, which we test
c h_outopen.in
c class limit
c ========================================
c
This tests the total emission from m open geometry, a small filling factor,
hydrogen Stromgren sphere,
using the outward only approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked.
h_outpp.in
plane parallel H-only, close, test hydrogenic emission
title plane parallel H-only, close, test hydrogenic emission
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
c the hydrogen density
hden 5
c an ini file that will save time for H-only
init file "honly.ini"
c
c commands controlling geometry =========
stop efrac 0.05 %stop when 5 percent ionized
c need this to prevent lyman lines from escaping
sphere static (OK)
c
c other commands for details =========
c prevent continuum pumping of lyman lines
no induced processes (OK)
c must include this to keep Ha optically thin
turbulence 20 km/s
c this is the approximation we are testing
diffuse outward
c force a constant temperature, then give stopping criteria
c if stop were not present model would go to limiting number of zones
constant temper 7500
c increase number of levels to get full emission
atom h-like levels large
c
c commands controlling output =========
c need to turn this on so that continua at various energies will
c be included in the emission line list
print diffuse continuum
c normalize results to perfect case b - our predictions should
c be close to this
normalize to line "Ca B" 4861
punch overview "h_outpp.ovr"
punch dr "h_outpp.dr"
c
c lots of continuum emission points turned on with print diffuse continuum
c energies dominated by two photon emission not included since very
c sensitive to Lya transport
c h_outpp.in
c class limit
c ========================================
c
This tests the total emission from a plane parallel pure hydrogen Stromgren sphere
using the outward only approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked.
h_outsp.in
spherical conservation and hydrogenic emission for pure hydrogen
title spherical conservation and hydrogenic emission for pure hydrogen
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 5
init file "honly.ini"
c
c commands controlling geometry =========
radius 13
c not iterating does not affect the solution
sphere static (OK)
stop efrac 0.05 %stop when 5 percent ionized
c
c other commands for details =========
constant temper 7500
no level2
atom h-like levels large
no induced processes (OK)
diffuse outward
c
c commands controlling output =========
print diffuse continuum
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "h_outsp.ovr"
punch dr "h_outsp.dr"
c
c lots of continuum emission points turned on with print diffuse continuum
c energies dominated by two photon emission not included since very
c sensitive to Lya transport
c h_outsp.in
c class limit
c ========================================
c
This tests the total emission from a spherical pure hydrogen Stromgren sphere
using the outward only approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked.
h_t4_conemis.in
test continuous emission from model H atom
title test continuous emission from model H atom
c
c commands controlling continuum =========
laser 2
ionization -1
c
c commands for density & abundances =========
hden 8
init "honly.ini"
c
c commands controlling geometry =========
stop thickness -16
c
c other commands for details =========
atom h-like levels large
c increase default resolution stored in continuum_mesh.ini
set continuum resolution 0.1
iterate
constant temperature 4
c
c commands controlling output =========
print diffuse continuum
print last iteration
print lines column
punch continuum "h_t4_conemis.con" last no units microns
punch continuum "h_t4_conemis.con" last no title units microns
punch spectrum "h_t4_conemis.ncon" last no units microns
punch spectrum "h_t4_conemis.ncon" last no title units microns
c
c h_t4_conemis.in
c class limit
c ========================================
c
this tests the continuous emission from the model H atom. The
gas temperature is 10000 K and the continuous emissivity
is asserted for a range of wavelengths.
h_t4_conemis_lon.in
test low-den continuous emission from H atom, 2-nu is important
title test low-den continuous emission from H atom, 2-nu is important
c
c commands controlling continuum =========
laser 2
ionization -1
c
c commands for density & abundances =========
hden 0
init "honly.ini"
c
c commands controlling geometry =========
stop thickness 0
c
c other commands for details =========
atom h-like levels large
c increase default resolution stored in continuum_mesh.ini
set continuum resolution 0.1
case b
iterate
constant temperature 4
c
c commands controlling output =========
print last iteration
print diffuse continuum
print lines column
punch continuum "h_t4_conemis_lon.con" last units microns
punch two photon continuum "h_t4_conemis_lon.2nu" last rydbergs
punch spectrum "h_t4_conemis_lon.ncon" last units microns
c
c h_t4_conemis_lon.in
c class limit
c ========================================
c
This is a mate to hatomt10.in except that the density is low enough
for two-photon emission to be very imporant in the optical and uv.
h_t4_conemis_thick.in
test hydrogen atom continuous emissivity, used for plot in hazy
title test hydrogen atom continuous emissivity, used for plot in hazy
c
c commands controlling continuum =========
black 5
ionization parameter -2
c
c commands for density & abundances =========
c high density to suppress two-photon emission
hden 7
c turn off all elements, except for hydrogen, so that
c we really test hydrogen emission
init file= "honly.ini"
c
c commands controlling geometry =========
c this is necessary to stop model at H ionization front
c constant temperature models will go on to infinity without this
stop efrac -0.8
c
c other commands for details =========
turbulence 200 km/s ; needed to keep balmer lines thin
c increase continuum resolution by factor of ten
set continuum resolution 0.1
c line fluorescense would be very important with above turbulence,
c turn it off
no induced processes
c need large atom to get the many Rydgerg levels in the infrared
atom h-like levels large
c what full outward flux so use spherical geometry
sphere
constant temperature 4
c
c commands controlling output =========
c results will be relative to interplated Hummer&Storey Hbeta
c this tests hydrogen atom too
normalize "Ca B" 4861
c by default, the continuum flux at many wavelengths is not printed,
c this turns that information on
print diffuse continuum
c output the predicted continuum
punch emitted continuum last iteration "h_t4_conemis_thick.con" units micron
punch dr last "h_t4_conemis_thick.dr"
c
c h_t4_conemis_thick.in
c class limit
c ========================================
c
This checks that the predicted hydrogen continuum is in good agreement
with exact results in the optically thin nebular limit.
Checks:
- This output was used to generate figure h_t4_conemis_thick in Part I of HAZY.
- Continuum relative to Hbeta should agree with Ferland (1980)
filter averaged results.
- Hbeta should agree with Case B predictions, and Q(H) 4861.
he1n2t4.in
test hei atom vs Benjamin et al. 99
title test hei atom vs Benjamin et al. 99
c TODO departure coef for n=20 along triplets shows suspicious behavior
c is s-p-d progression
c
c >>refer hei intensity Benjamin, R.A., Skillman, E.D.,
c >>refercon & Smits, D.P. 1999, ApJ, 514, 307
c
c use Benjamin data, including rec coef and size of atom,
c and neglect collisions he does not include. If an attempt is made
c to set the size of the helium atom, and at the same time use the
c Benjamin command, cloudy will ignore the size command.
c
c commands controlling continuum =========
laser 2
ionization -2
c
c commands for density & abundances =========
hden 90 linear
set eden 2
init file "hheonly.ini"
element helium abundance -1
c
c commands controlling geometry =========
stop zone 1
set dr -2.95424
c
c other commands for details =========
atom he-like benjamin
atom he-like gbar off
case b no photoionzation
constant temper 4
iterate
c
c commands controlling output =========
c this increases the number of wavelength digits in the line output
c to make it easier to identify highly excited Rydberg lines
set line precision 6
print he-like departure helium
normalise to "He 1" 4471.47
print line faint 0.01
c
c
c he1n2t4.in
c class limit
c ========================================
c
This is one of the standard comparisons with the Benjamin et al.
(1999), ApJ, 514, 307, paper. The asserts give their answers.
the line wavelength precision is increased to differentiate between
various HeI lines
he1n2t4_best.in
test hei atom vs Benjamin et al. 99
title test hei atom vs Benjamin et al. 99
c >>refer hei intensity Benjamin, R.A., Skillman, E.D.,
c >>refercon & Smits, D.P. 1999, ApJ, 514, 307
c
c commands controlling continuum =========
laser 2
ionization -2
c
c commands for density & abundances =========
hden 90 linear
set eden 2
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr -2.95424
c
c other commands for details =========
c the no pdest option turns of dest probabilities, needed
c to prevent higher lyman lines from being destroyed by background
c opacity
case b no photoionzation no pdest
constant temper 4
iterate
c make he atom VERY large, n = 20 will be done with resolved levels
atom he-like levels 20
atom he-like collapsed levels 100
atom he-like gbar vriens
c
c commands controlling output =========
normalise to "He 1" 4471
print line faint 0.01
print he-like departure helium
c
c he1n2t4_best.in
c class limit
c ========================================
c
This is one of the standard comparisons with the Benjamin et al.
(1999), ApJ, 514, 307, paper. The asserts give their answers.
This turns on all processes and does not try to emulate their paper.
he1n2t4_Smits96.in
test hei atom vs Smits 96
title test hei atom vs Smits 96
c >>refer hei emission Smits, D.P., 1996, MNRAS, 278, 683
c
c commands controlling continuum =========
laser 2
ionization -3
c
c commands for density & abundances =========
c his h density was 0.9 * elec den
hden 1.954
set eden 2
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
c this is 1 / n_e n_He^+
set dr -2.9542
stop zone 1
c
c other commands for details =========
c make he atom VERY large, n = 20 will be done
c this is largest atom now possible without computing new rec coef
atom he-like levels 20
atom he-like collapsed levels 100
atom he-like collisional excitation off
atom he-like collisions fake off
case b no photoionization
constant temper 4
iterate
c
c commands controlling output =========
print he-like departure helium
normalise to "He 1" 4471
print line faint -1
c
c
c he1n2t4_Smits96.in
c class limit
c ========================================
c
This compares the He I model atom with the Smits (1996)
predictions. The thickness of the zone is set to force an
emissivity of unity.
he1n4t4.in
test hei atom vs Benjamin et al. 99
title test hei atom vs Benjamin et al. 99
c
c commands controlling continuum =========
laser 2
ionization -1
c
c commands for density & abundances =========
hden 9000 linear
set eden 4
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr -6.95424
c
c other commands for details =========
atom he-like gbar off
atom he-like benjamin
constant temper 4
case b no photoionzation
iterate
normalise to "He 1" 4471
print departure he-like helium
print line faint 0.01
c
c commands controlling output =========
punch fits "he1n4t4.fit"
punch diffuse continuum "he1n4t4.dif" last no units microns
punch diffuse continuum "he1n4t4.dif" last no title units microns
punch continuum "he1n4t4.con" last no units microns
punch continuum "he1n4t4.con" last no title units microns
c
c
c he1n4t4.in
c class limit
c ========================================
c
The Benjamin et al. 1999 model at 1e4 density,
turning off processes to simulate their calculations.
he1n4t4_best.in
the best we can do to predict the HeI emission spectrum
title the best we can do to predict the HeI emission spectrum
c
c commands controlling continuum =========
laser 2
ionization -1
c
c commands for density & abundances =========
hden 9000 linear
set eden 4
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr -6.95424
c
c other commands for details =========
atom he-like levels 30
atom he-like collapsed levels 70
atom he-like gbar vriens
constant temper 4
case b no photoionzation no pdest
iterate
c
c commands controlling output =========
normalise to "He 1" 4471
print departure he-like helium
print line faint 0.01
punch diffuse continuum "he1n4t4_best.dif" last no units microns
punch diffuse continuum "he1n4t4_best.dif" last no title units microns
punch continuum "he1n4t4_best.con" last no units microns
punch continuum "he1n4t4_best.con" last no title units microns
c
c
c he1n4t4_best.in
c class limit
c ========================================
c
This is the best and most complete model of He I that the code can do.
The results are currently the Benjamin et al. 1999 numbers.
he1n6t4.in
test hei atom vs Benjamin et al. 99
title test hei atom vs Benjamin et al. 99
c >>refer hei intensity Benjamin, R.A., Skillman, E.D.,
c >>refercon & Smits, D.P., 1999, ApJ, 514, 307
c
c commands controlling continuum =========
laser 2
ionization 0
c
c commands for density & abundances =========
hden 9.e5 linear
set eden 6
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr -10.95424
c
c other commands for details =========
atom he-like gbar off
c use benjamin rec coef, ignore collision he did not include
atom he-like benjamin
print he-like departure helium
case b no photoionization
constant temper 4
iterate
c
c commands controlling output =========
normalise to "He 1" 4471
print line faint -2
c
c
c he1n6t4.in
c class limit
c ========================================
c
This compares the He I model atom with the Benjamin et al (1999)
predictions. The thickness of the zone is set to force an
emissivity of unity.
he1n6t4_best.in
test hei atom vs Benjamin et al. 99
title test hei atom vs Benjamin et al. 99
c >>refer hei intensity Benjamin, R.A., Skillman, E.D.,
c >>refercon & Smits, D.P., 1999, ApJ, 514, 307
c
c commands controlling continuum =========
laser 2
ionization 0
c
c commands for density & abundances =========
hden 9.e5 linear
set eden 6
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr -10.95424
c
c other commands for details =========
c make he atom VERY large, n = 20 will be done
atom he-like levels 20
atom he-like collapsed levels 280
atom he-like gbar vriens
case b no photoionization no pdest
constant temper 4
iterate
c
c commands controlling output =========
normalise to "He 1" 4471
print he-like departure helium
print line faint 0.01
c
c check that lines still recognized
table lines "LineList_HeH.dat"
c
c he1n6r4_best.in
c class limit
c ========================================
c
This compares the He I model atom with the Benjamin et al (1999)
predictions. The thickness of the zone is set to force an
emissivity of unity.
he1n6t4_Smits96.in
test hei atom vs Smits 96
title test hei atom vs Smits 96
c >>refer hei emission Smits, D.P., 1996, MNRAS, 278, 683
c
c commands controlling continuum =========
laser 2
ionization 0
c
c commands for density & abundances =========
hden 5.95424
set eden 6
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr -10.95424
c
c other commands for details =========
c make he atom VERY large, n = 20 will be done
atom he-like levels 20
atom he-like collapsed levels 100
atom he-like collisional excitation off
case b no photoionization
constant temper 4
iterate
c
c commands controlling output =========
print he-like departure helium
normalise to "He 1" 4471
print line faint -2
print he-like departure helium
c
c
c he1n6t4_Smits96.in
c class limit
c ========================================
c
This compares the model He I line emission with the Smits (1996) results.
The thickness of the zone is set so that we can predict the precise
intensity of the normalization line, He I 4471.
he2_caseb.in
he2_caseb He II case B
title he2_caseb He II case B
c
c commands controlling continuum =========
table agn
ionization parameter 0
c
c commands for density & abundances =========
hden = 8
element helium abundance 1
init file="ism.ini"
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
atom h-like helium levels large
case b hummer no photoionization
no induced processes (OK)
constant temperature = 50000
c
c commands controlling output =========
punch results "he2_caseb.rlt"
c
c caseb_c6.in
c class limit
c ========================================
c
This test case compares the predictions of the multi-level hydrogen atom
with the Storey and Hummer (1995) results. The set dr command sets the
zone thickness to 1 cm. The case b command sets Lyman line optical depths
to very large values.
Checks:
- Neutral fractions
- H* emissivity
- Relative line intensities
heatomt10.in
continuous emission from HeI
title continuous emission from HeI
c
c commands controlling continuum =========
laser 1.9 ryd
ionization -1
c
c commands for density & abundances =========
c need high density to supress two-photon emission
hden 10
c force high density for electrons, He+ density will be the same
set eden 12
c set He/H ratio to 100, He and elec will have same density
element abundance helium 2
init "hheonly.ini"
c
c commands controlling geometry =========
c this thickness will get unit emission when combined with above densities
stop thickness -24
c
c other commands for details =========
c must make atom larger so get all edges in IR
atom he-like levels 15
c add large number of collapsed levels
atom he-like collapsed levels 50
iterate
constant temperature 4
c
c commands controlling output =========
print continuum
print last iteration
c
set continuum resolution 0.1
set nFnu diffuse outward, diffuse inward
print continuum
print last iteration
print line faint 1
print line column
punch spectrum "heatomt10.spc" last no units microns
punch spectrum "heatomt10.spc" last no title units microns
punch continuum "heatomt10.con" last no units microns
punch continuum "heatomt10.con" last no title units microns
punch diffuse continuum "heatomt10.dif" last no units microns
punch diffuse continuum "heatomt10.dif" last no title units microns
c
c heatomt10.in
c class limit
c ========================================
c
This tests continuous emission from the He I atom.
The laser is used so that the incident continuum
is not included in the total emission.
heatomt10lon.in
test low-den continuous emission from H atom, 2-nu is important
title test low-den continuous emission from H atom, 2-nu is important
c
c commands controlling continuum =========
laser 2
ionization -1
c
c commands for density & abundances =========
hden -8
element helium abundance log 3
set eden -5
init "hheonly.ini"
c
c commands controlling geometry =========
stop thickness log 10
c
c other commands for details =========
c force all helium to be in He+
atom he-like levels 20
atom he-like collapsed levels 20
no scattering opacity
element helium ionization -4 0 -4
case b
c increase default resolution stored in continuum_mesh.ini
set continuum resolution 0.1
iterate
constant temperature 4
c
c commands controlling output =========
normalize to "he 1" 4471
print lines column
print last iteration
print diffuse continuum
print lines column
punch continuum "heatomt10lon.con" last units microns
punch diffuse continuum "heatomt10lon.dif" last units microns
punch two photon continuum "heatomt10lon.2nu" last
punch spectrum "heatomt10lon.ncon" last units microns
c
c heatomt10lon.in
c class limit
c ========================================
c
This is a mate to hatomt10.in except that everything is He at low density
heiont10.in
continuous emission from HeII
title continuous emission from HeII
c
c commands controlling continuum =========
laser 5 ryd
ionization 1
c
c commands for density & abundances =========
hden 8
set eden 10
element abundance helium 2
init "hheonly.ini"
c
c commands controlling geometry =========
stop thickness -20
c
c other commands for details =========
constant temperature 4
iterate
c increase continuum resolution by factor of ten
set continuum resolution 0.1
set nFnu diffuse outward, diffuse inward
c
c commands controlling output =========
print continuum
print last iteration
print line faint 1
print line column
punch continuum "heiont10.con" last no units microns
punch continuum "heiont10.con" last no title units microns
punch spectrum "heiont10.spc" last no units microns
punch spectrum "heiont10.spc" last no title units microns
c
c heiont10.in
c class limit
c ========================================
c
This tests the He II continuous emission. The helium abundance
is very large so that He II overwhelms other emission sources.
The resolution of the continuum mesh is increased so that we get
a better representation of the continuous emission.
helike_ar.in
test he-like argon ion
title test he-like argon ion
c
c commands controlling continuum =========
laser 600
ionization -1
c
c commands for density & abundances =========
hden 7
element argon abundance -2
element argon ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element silicon off
element sulphur off
element oxygen off
element iron off
c
c commands controlling geometry =========
set dr -12
stop zone 1
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print he-like departure argon
print line faint -2
c This is Lyman alpha
normalize to "Ar17" 3.949A
c
c 3.994A is 2^3S to ground
c Totl is 2^3P to ground
c 3.365A is 3^1P to ground
c 21.54A is 3^3P to 2^3S
c 22.24A is 3^3D to 2^3P
c
c helike_ar.in
c class limit
c ========================================
c
test whether he-like ion populations and emission for argon
agree with published calcs.
helike_c.in
test he-like carbon ion
title test he-like carbon ion
c
c commands controlling continuum =========
laser 50
ionization -1
c
c commands for density & abundances =========
hden 7
element carbon abundance -2
element carbon ionization -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element nitrogen off
element oxygen off
element neon off
element magnesium off
element silicon off
element sulphur off
element argon off
element iron off
c
c commands controlling geometry =========
set dr -12
stop zone 1
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print he-like departure carbon
print line faint -2
c This is Lyman alpha
normalize to "C 5" 40.27A
c
c
c helike_c.in
c class limit
c ========================================
c
test whether he-like ion populations go to lte in high density limit.
The level populations should be in LTE, and the departure coefficients
should be unity.
helike_co.in
test he-like cobalt ion
title test he-like cobalt ion
atom he-like levels 10
atom he-like collapsed levels 20
print he-like departure cobalt
constant temper 5
iterate
hden 7
set dr -12
init file "ism.ini"
c must turn on since ism.ini turned it off
element cobalt on
c then set abundance and ionization after turning it on
element cobalt abundance -2
element cobalt ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
c now turn lots of elements off to save time
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element silicon off
element sulphur off
element oxygen off
element argon off
element iron off
laser 1300
ionization -1
stop zone 1
print line faint -2
helike_cu.in
test he-like copper ion
title test he-like copper ion
atom he-like levels 10
atom he-like collapsed levels 20
print he-like departure copper
constant temper 5
iterate
hden 7
set dr -12
init file "ism.ini"
c must turn on since ism.ini turned it off
element copper on
c then set abundance and ionization after turning it on
element copper abundance -2
element copper ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
c now turn lots of elements off to save time
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element silicon off
element sulphur off
element oxygen off
element argon off
element iron off
laser 1300
ionization -1
stop zone 1
print line faint -2
helike_fe.in
test he-like iron ion
title test he-like iron ion
c
c commands controlling continuum =========
laser 1300
ionization -1
c
c commands for density & abundances =========
hden 7
set dr -12
element iron abundance -2
element iron ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element silicon off
element sulphur off
element oxygen off
element argon off
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print line faint -2
print he-like departure iron
c This is Lyman alpha
normalize to "Fe25" 1.850A
c
c
c helike_fe.in
c class limit
c ========================================
c
test whether he-like ion populations and emission for argon
agree with published calcs.
helike_mg.in
test he-like magnesium ion
title test he-like magnesium ion
c
c commands controlling continuum =========
laser 250
ionization -1
c
c commands for density & abundances =========
hden 7
set dr -12
element magnesium abundance -2
element magnesium ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element nitrogen off
element carbon off
element neon off
element silicon off
element oxygen off
element sulphur off
element argon off
element iron off
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print he-like departure magnesium
print line faint -2
c This is Lyman alpha
normalize to "Mg11" 9.169A
c
c
c helike_mg.in
c class limit
c ========================================
c
test whether he-like ion populations and emission for magnesium
agree with published calcs.
helike_n.in
test he-like nitrogen ion
title test he-like nitrogen ion
c
c commands controlling continuum =========
laser 75
ionization -1
c
c commands for density & abundances =========
hden 7
element nitrogen abundance -2
element nitrogen ionization -5 -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element carbon off
element oxygen off
element neon off
element magnesium off
element silicon off
element sulphur off
element argon off
element iron off
c
c commands controlling geometry =========
stop zone 1
set dr -12
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print he-like departure nitrogen
print line faint -2
c This is Lyman alpha
normalize to "N 6" 28.79A
c
c helike_n.in
c class limit
c ========================================
c
test whether he-like ion populations go to lte in high density limit.
The level populations should be in LTE, and the departure coefficients
should be unity.
helike_ne.in
test he-like neon ion
title test he-like neon ion
c
c commands controlling continuum =========
laser 170
ionization -1
c
c commands for density & abundances =========
hden 7
element neon abundance -2
element neon ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element nitrogen off
element carbon off
element oxygen off
element magnesium off
element silicon off
element sulphur off
element argon off
element iron off
c
c commands controlling geometry =========
stop zone 1
set dr -12
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print he-like departure neon
print line faint -2
c This is Lyman alpha
normalize to "Ne 9" 13.45A
c
c
c helike_ne.in
c class limit
c ========================================
c
test whether he-like ion populations and emission for oxygen
agree with published calcs.
helike_ni.in
test he-like nickel ion
title test he-like nickel ion
atom he-like levels 10
atom he-like collapsed levels 20
print he-like departure nickel
constant temper 5
iterate
hden 7
set dr -12
init file "ism.ini"
c must turn on since ism.ini turned it off
element nickel on
c then set abundance and ionization after turning it on
element nickel abundance -2
element nickel ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
c now turn lots of elements off to save time
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element silicon off
element sulphur off
element oxygen off
element argon off
element iron off
laser 1300
ionization -1
stop zone 1
print line faint -2
helike_o.in
test he-like oxygen ion vs. Bautista & Kallman 2000 Table 1, column 3
title test he-like oxygen ion vs. Bautista & Kallman 2000 Table 1, column 3
c
c commands controlling continuum =========
laser 100
ionization -1
c
c commands for density & abundances =========
hden 7
element oxygen abundance -2
element oxygen ionization -5 -5 -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element silicon off
element sulphur off
element argon off
element iron off
c
c commands controlling geometry =========
stop zone 1
set dr -12
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print he-like departure oxygen
print line sort wavelength
print line faint -5
normalize to "o 7" 21.60A
c
c helike_o.in
c class limit
c ========================================
c
test whether he-like ion populations and emission for oxygen
agree with published calcs.
helike_si.in
test he-like silicon ion
title test he-like silicon ion
c
c commands controlling continuum =========
laser 350
ionization -1
c
c commands for density & abundances =========
hden 7
element silicon abundance -2
element silicon ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
init file "ism.ini"
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element oxygen off
element sulphur off
element argon off
element iron off
c
c commands controlling geometry =========
stop zone 1
set dr -12
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 5
iterate
c
c commands controlling output =========
print he-like departure silicon
print line faint -2
c This is Lyman alpha
normalize to "Si13" 6.648A
c
c
c helike_si.in
c class limit
c
test whether he-like ion populations and emission for silicon
agree with published calcs.
helike_zn.in
test he-like zinc ion
title test he-like zinc ion
atom he-like levels 10
atom he-like collapsed levels 20
print he-like departure zinc
constant temper 5
iterate
hden 7
set dr -12
init file "ism.ini"
c must turn on since ism.ini turned it off
element zinc on
c then set abundance and ionization after turning it on
element zinc abundance -2
element zinc ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
c now turn lots of elements off to save time
element helium off
element nitrogen off
element carbon off
element neon off
element magnesium off
element silicon off
element sulphur off
element oxygen off
element argon off
element iron off
laser 1300
ionization -1
stop zone 1
print line faint -2
hhe_otspp.in
plane parallel conservation and hydrogenic emission for pure H, He
title plane parallel conservation and hydrogenic emission for pure H, He
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 5
init file "hheonly.ini"
c
c commands controlling geometry =========
c expanding sphere to stop Lya from destroying HeI 23S
sphere
stop efrac 0.05 %stop when 5 percent ionized, so just beyond H I-front
c
c other commands for details =========
diffuse ots
atom h-like levels large
c that depend on having stable solution
iterate
constant temper 7500
no induced processes (OK)
c
c commands controlling output =========
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "hhe_otspp.ovr" last
punch dr "hhe_otspp.dr" last
c
c helium lines
c hhe_otspp.in
c class limit
c ========================================
c
This tests the total emission from a spherical pure H + He Stromgren sphere
using the outward only approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked.
hhe_otssp.in
pherical conservation and hydrogenic emission for hydrogen and helium
title pherical conservation and hydrogenic emission for hydrogen and helium
c
c commands controlling continuum =========
blackbody 40000
ionization -2
c
c commands for density & abundances =========
hden 5
init file "hheonly.ini"
c
c commands controlling geometry =========
diffuse ots
atom h-like levels large
radius 13
stop efrac 0.05 %stop when 5 percent ionized
sphere static
c
c other commands for details =========
constant temper 7500
iterate
no level2
no induced processes (OK)
c
c commands controlling output =========
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "hhe_otssp.ovr" last
punch dr "hhe_otssp.dr" last
c
c helium lines
c hheonlyoutpp.in
c class limit
c ========================================
c
This tests the total emission from a spherical pure H + He-like Stromgren sphere
using the outward only approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked. The geometry is plane paralel.
hhe_outpp.in
plane parallel conservation and H, He emission for pure H, He gas
title plane parallel conservation and H, He emission for pure H, He gas
c
c commands controlling continuum =========
blackbody 40000
ionization -2
c
c commands for density & abundances =========
hden 5
init file "hheonly.ini"
c
c commands controlling geometry =========
c open sphere to stop Lya from destroying HeI 23S
sphere
stop efrac 0.05 %stop when 5 percent ionized
c
c other commands for details =========
iterate
constant temper 7500
atom h-like levels large
no level2
no induced processes (OK)
c
c commands controlling output =========
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "hhe_outpp.ovr" last
punch dr "hhe_outpp.dr" last
c
c hhe_outpp.in
c class limit
c ========================================
c
This tests the total emission from a spherical pure hydrogen Stromgren sphere
using the outward only approximation. The conservation of the total number
of ionizing photons, and the emitted spectrum, are all checked.
hhe_outppff.in
plane parallel filling factor for pure H, He gas
title plane parallel filling factor for pure H, He gas
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 5
init file "hheonly.ini"
c
c commands controlling geometry =========
filling factor -1
stop efrac 0.05 %stop when 5 percent ionized
sphere
c
c other commands for details =========
iterate
constant temper 7500
no level2
no induced processes (OK)
c
c commands controlling output =========
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "hhe_outppff.ovr" last
punch dr "hhe_outppff.dr" last
c
c helium lines
c hhe_outppff.in
c class limit
c ========================================
c
This is a plane-parallel constant temperature cloud with only
hydrogen and helium. The gas has a filling factor of 0.1.
Induced processes are turned off and a large H atom is used
so that the hydrogen recombination
spectrum will be close to Case B. The calculation stops beyond the
hydrogen ionization front, because of the stop efrac command (it needs
this since this is a constant temperature calculation, so the usual
lower-temperature stopping criterion does not apply). The asserts
confirm that energy is conserved and that the hydrogen spectrum is correct.
hhe_outsp.in
spherical conservation and hydrogenic emission for hydrogen and helium
title spherical conservation and hydrogenic emission for hydrogen and helium
c
c commands controlling continuum =========
ionization -2
blackbody 40000
c
c commands for density & abundances =========
hden 5
init file "hheonly.ini"
c
c commands controlling geometry =========
radius 13
stop efrac 0.05 %stop when 5 percent ionized
sphere static
c
c other commands for details =========
iterate
diffuse outward only
constant temper 7500
no level2
no induced processes (OK)
c
c commands controlling output =========
normalize to "Ca B" 4861 = 1 % normalize to exact result
punch overview "hhe_outsp.ovr" last
punch dr "hhe_outsp.dr" last
c
c helium lines
c hheonlyoutpp.in
c class limit
c ========================================
c
This tests a spherical cloud with only hydrogen and helium.
Diffuse fields are transferred with the outward only approximation.
The asserts check that the ionizing radiation is conserved.
hii_blister.in
conditions similar to Orion nebula hii_blister
title conditions similar to Orion nebula hii_blister
c
c commands controlling continuum =========
c a simple blackbody
blackbody 40000
c the flux of H-ionizing photons
phi(h) 13.0
c add the CMB after hot star shape & luminosity specified
CMB
c
c commands for density & abundances =========
hden 4
init file="ism.ini"
abundances hii region no grains
abundances he1 c1 n1 o1 ne1 mg=.0001 si.01 s1 cl1 ar1
continue fe.001
c
c commands controlling geometry =========
sphere
c
c other commands for details =========
c this is to desaturate Lya to prevent excited state photoionization
c from being important - the other codes did not include this important
c physical process
turbulence 10
c must iterate since optically thick
iterate
c
c commands controlling output =========
normalize to "Ca B" 4861
print line faint .01
print line sum
He 1 5876
totl 2326
C 2 1335
TOTL 1909
N 2 6584
N 2 6548
N 3 57.21m
TOTL 7325
TOTL 3727
O 3 51.80m
O 3 5007
O 3 4959
Ne 2 12.81m
Ne 3 15.55m
Ne 3 3869
Ne 3 3968
S 3 18.67m
S 3 9532
S 3 9069
S 4 10.51m
end of lines
punch overview last "hii_blister.ovr"
punch dr last "hii_blister.dr"
punch results last "hii_blister.rlt"
c
c
c hii_blister.in
c class hii
c ========================================
This is one of the test cases from the Lexington Meeting suite of
nebulae (Ferland et al. 1995). It is a grain-free hii_blister HII region,
similar to inner regions of the Orion Nebula, except for the
absence of grains.
The set of lines entered with the print line sum command lists
the most powerful coolants in this model. This is one of the
tabulated quantities in the Lexington Meeting, and is a
fundamental test of energy conservation in the code.
The ratio of the sum of these lines to Hb is equivalent to the
Stoy ratio, used for determining stellar temperatures.
The "dielec kludge 0" command is to turn off my estimates of the
DR rates for those elements that had none. This was only to allow
comparison with other calculations that did not make similar estimates.
For an actual calculation I would not include this command, since the
guesses are better than nothing.
the turbulence is to stop the balmer lines from becoming optically
thick since few other codes include an actual H atom, but use
case b instead. The Orion HII region does have an observed turbulence
of about 8 km/s.
This calculation stops near the H+ - H0 ionization front, where
the temperature falls below the default lowest temperature of 4000 K.
This model would have continued into the PDR had a lower temperature
been specified with the STOP LOWEST TEMP command.
hii_coolstar.in
cool HII region model from Lexington Meeting
title cool HII region model from Lexington Meeting
c
c commands controlling continuum =========
black body, T = 20000 K
q(h) 49
c add the CMB after hot star shape & luminosity specified
CMB
c
c commands for density & abundances =========
hden = 2
init file="ism.ini"
abund He-1 C-3.6576 N-4.39794 O-3.481146 ne-4.30103 mg-8
continue si-8 s-5.04576 cl-7 ar-8 fe-8
c
c commands controlling geometry =========
radius = 18.477121
sphere
c this is to go deep to pick up all H recombs
stop efract -2
c note that Case B predictions will not be available since this
c temperature overruns the HS table
stop temperature 100
c
c other commands for details =========
c must iterate since fine structure lines are optically thick
iterate
c assume there is neutral gas beyond what we compute
double
c
c commands controlling output =========
plot continuum range .1
print line faint .01
print line sum :this is the sum of lines in table 2 of the Lexington meeting
N 2 6584
N 2 6548
TOTL 3727
Ne 2 12.81m
S 2 6720
S 3 18.67m
S 3 33.47m
S 3 9532
S 3 9069
end of lines
punch overview last "hii_coolstar.ovr"
punch dr last "hii_coolstar.dr"
punch results last "hii_coolstar.rlt"
punch continuum last units microns "hii_coolstar.con"
punch lines, cumulative, "hii_coolstar.cum"
totl 4861
o 3 5007
totl 3727
o 1 6300
end of lines
c
c Hbeta 4.93E36, L(total)4.30xHbeta
c hii_coolstar.in
c class hii
c ========================================
c
This is one of the test cases from the Lexington Meeting suite of nebulae.
It is a grain-free HII region ionized by a very cool star. Hydrogen is
ionized but not helium so this tests the transport of the H Lyman continuum.
The set of lines is entered with the print line sum command to test
energy conservation.
hii_icf.in
HII region with negative He/H ICF
title HII region with negative He/H ICF
c
c commands controlling continuum =========
c this is the result of this command
c table star mihalas 46000
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_mihalas_46000.dat"
ionization parameter -2.
c add the CMB after hot star shape & luminosity specified
CMB
c
c commands for density & abundances =========
hden 3
init file "ism.ini"
abundances ism no grains
grains no qheat single
metals and grains 0.320
c
c commands controlling geometry =========
stop efrac -3
stop temperature 15
sphere
c
c other commands for details =========
failures 1
c
c commands controlling output =========
punch overview "hii_icf.ovr" last
punch dr "icr.dr" last
c
c hii_icf.in
c class hii
c ========================================
This is an example of an H II region irradiated
by a hard stellar continuum - one of the Mihalas NLTE stars.
The hard continuum produces a negative He/H ionization
correction factor, as discussed in Ballantyne, Ferland & Martin (2000).
>>refer HeI icf by Ballantyne,D.R., Ferland, G.J., & Martin, P.G., 2000, ApJ 536, 773-777
hii_paris.in
"New" Paris meeting HII region
title "New" Paris meeting HII region
c "standard" HII region model of the Pequignot Meudon Conferance
c
c commands controlling continuum =========
black body, T = 40000 K radius = 12.113943
c add the CMB after hot star shape & luminosity specified
CMB
c
c commands for density & abundances =========
hden = 2
init file="ism.ini"
abund he -1 C-3.6576 N-4.39794 O-3.481146 ne-4.30103 mg-8
continue si-8 s-5.04576 cl=-7 ar-8 fe-8
c
c commands controlling geometry =========
radius = 18.477121
sphere
c next two to make sure we pick up all possible H recombination
stop temperature 600
stop efrac -2
c
c other commands for details =========
c must iterate since fine structure lines are opticall thick
iterate
c
c commands controlling output =========
normalize to "Ca B" 4861
print column densities
print line optical depths
plot continuum range .1
print line faint 0.01
print line sum
He 1 5876
totl 2326
TOTL 1909
N 2 121.7m
N 2 6584
N 2 6548
N 3 57.21m
TOTL 3727
O 3 5007
O 3 4959
O 3 51.80m
O 3 88.33m
Ne 2 12.81m
Ne 3 15.55m
Ne 3 3869
Ne 3 3968
S 2 6720
S 3 18.67m
S 3 33.47m
S 3 9532
S 3 9069
S 4 10.51m
end of lines
punch overview last "hii_paris.ovr"
punch hydrogen 21 cm last "hii_paris.21cm"
punch results last "hii_paris.rlt"
punch dr last "hii_paris.dr"
punch lines structure "hii_paris.str" last
TOTL 4861
end lines
c
c hii_paris.in
c class hii
c ========================================
c
This is one of the "standard" models computed at the Paris and Lexington
meetings on photoionization and shock calculations. a bable in hazy
compares the predictions of the current version of CLOUDY with predictions
of a few of the other codes. It is necessary to iterate since some
fine structure lines are optically thick. The set of lines entered
with the print line sum command is used to obtain the total luminosity
in detected lines, a measure of the Stoy temperature.
Checks:
- Hb close to case B, Q(H) 4861, intensities.
- Enter answers in Table Error! Reference source not found..
igm_lalpha.in
Ly alpha forest cloud
title Ly alpha forest cloud
c
c commands controlling continuum =========
c ionized by cosmic background
c plus quasar/starburst continuum
cmb z=2
table hm05 z = 2
c
c commands for density & abundances =========
hden -2
metals -1.5 ;reduce all heavy elements by 1.5 dex
c
c commands controlling geometry =========
stop neutral column density 15
double ; mimic two-siced photoionization
c
c other commands for details =========
iterate to convergence; must iterate since optically thin
c
c commands controlling output =========
print line faint -1
print last
punch overview last "igm_lalpha.ovr"
punch dr last "igm_lalpha.dr"
punch results last "igm_lalpha.rlt"
c
c on 05 aug 29 changed from old background command to
c cmb 2 plus hm05 z=2 background - all results changed substantially
c igm_lalpha.in
c class igm
c ========================================
c
This example demonstrates the behavior of the code in the low-density limit.
The ionizing source is the cosmic background at a redshift of z = 2.
The cloud is assumed to be optically thin to ionizing radiation,
in keeping with the Gunn-Peterson test, so the double command is included,
and an iteration is performed to converge the optical depth scale.
Continuum fluorescent excitation of lines is important because the cloud
is otically thin.
igm_primal.in
cloud with primordial abundances exposed to background at Z=10
title cloud with primordial abundances exposed to background at Z=10
c
c commands controlling continuum =========
background, redshift= 10
c background due to agn from Haardt & Madau 1996 ApJ, 461, 20
c this sets both continuum shape and intensity, actually for z=2
table HM96 old
c
c other commands for details =========
c put in a velocity field to stop Lya from causing stability problems
turbulence 20 km/s
c do a second iteration to establish that it is optically thick downstream
iterate
c this sim should run cleanly - do not accept problems
failures 2
c
c commands for density & abundances =========
c this log of the hydrogen density
hden 5
c stored set of primordial abundances
abundances primordial
c
c commands controlling geometry =========
c we want to shielded face to not be exposed to empty space, but rather
c be just very deep in a dense cloud - this prevents lines from freely
c radiating from the shielded face
double optical depths
c set a lower temperature stopping criterion - default is 4000K
stop temperature 30
c
c commands controlling output =========
c don't want to print lots of faint hydrogen lines
print line faint -3
c some output files with info we want to save
punch temperature history "igm_primal_tem.his"
punch ionizing continuum "igm_primal.ion"
punch hydrogen conditions "igm_primal.mol"
punch heating "igm_primal.het"
punch dr "igm_primal.dr"
punch cooling "igm_primal.col"
punch overview "igm_primal.ovr"
punch results last "igm_primal.rlt"
c
c igm_primal.in
c class hii pdr igm
c ========================================
This is a high redshift cloud irradiated by the cosmic background and AGN light.
H- absorption of the z=10 cosmic background is a MAJOR heating agent
This model is very sensitive to treatment of Lya destruction, since resulting
J-bar is major destruction process for H-, which is the H2 formation route.
Lya ots oscillations could develop and are damped by not reevaluating H-
photo rate after first n times, as in hmole.c
igm_z3.in
redshift 1000 recombination epoch
title redshift 1000 recombination epoch
c
c commands controlling continuum =========
c a true blackbody radiation field, in strick thermo equilibrium
blackbody 3000 STE
c
c commands for density & abundances =========
hden 0
abundances primordial
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
atom h-like levels large
case b
iterate
c
c commands controlling output =========
normalize lines to "Ca B" 4861
print ages
print departure coef
c
c igm_z3.in
c class limit igm
c ========================================
c
This is a model of the universe near the recombination epoch, at
a redshift of a thousand. The gas is exposed to a true blackbody at
3000 K, and the abundances are primordial.
ism.in
interstellar cloud irradiated by ism background
title interstellar cloud irradiated by ism background
c
c commands controlling continuum =========
c this uses the Black continuum with no H-ionizing radiation,
c and with background cosmic rays
table ism
extinguish by a column of 22
cosmic rays, background
c
c commands for density & abundances =========
hden 0
init file="ism.ini"
abundances ism
c
c commands controlling geometry =========
c set sphere since matter in all directions
sphere
stop temperature linear 10
stop thickness 0.1 linear parsecs
c
c other commands for details =========
c this prevents continuum pumping from exciting H0
case b
c this is to converge optical depths
iterate
c
c commands controlling output =========
c this model has no heavy element molecules, because H2 does not
c really fully form, since Lyman bands are not self-shielded
normalize to 157.6m "C 2"
print line pump
punch overview "ism.ovr" last
punch dr "ism.dr"
punch heating "ism.het"
punch coolign "ism.col"
c
c ism.in
c class ism
c ========================================
TODO - look at temperature struture - it has jitter at about the
level of convergence, up and down. Temp jitter caused by eden jitter.
This model is nearly isothermal, jitter measures noise in solver,
and is great chance to pin this down.
This is a test of the behavior of the code in the extreme of
photoionization by a relatively hard continuum, at low densities.
The continuum is the galactic background, attenuated by a column density
of 10^22 cm-2. Ionization by galactic background cosmic rays is included.
Case b appears since this region is deep in the ISM, and the Lyman lines
are quite thick. This example checks whether the ionization balance,
thermal balance, and electron density sum, are performed correctly in
this limit.
Checks:
- Numerical stability of solution
- Thickness exact
ism_cosmicray.in
background cosmic ray ionization by suprathermal electrons only
title background cosmic ray ionization by suprathermal electrons only
c
c commands controlling continuum =========
c background CR H0 ionization rate should be 2.5e-17 s^-1,
c Williams et al. ApJ, 503, 689
cosmic rays, background
c this will be VERY faint black 50000
ionization parameter -25
black 50000
c
c commands for density & abundances =========
hden 5
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
iterate
constant temperature 1000 K
c turn off molecules to test only ions
no molecules
c must turn off charge transfer and its heating, which would dominate,
c since we want to test cr heating and ionization
no charge transfer
no ctheat
c must turn of photoionization by recombining species,
c since we want to test cr heating and ionization
no photoionization
c any time photoionization is turned off, it is best to also
c turn off induced processes, which include pumping in the lyman lines
no induced processes
c
c commands controlling output =========
print ages
c
c ism_cosmicray.in
c class limit
c ========================================
c
This test conditions of cosmic ray ionization. Molecules and charge transfer
are disabled so that analytical estimates can be made.
ism_grid.in
interstellar cloud irradiated by ism background
title interstellar cloud irradiated by ism background
c
c commands controlling continuum =========
c this uses the Black continuum with no H-ionizing radiation,
c and with background cosmic rays
table ism
extinguish by a column of 22
cosmic rays, background
c
c commands for density & abundances =========
hden 0 vary
grid from -3 to 0 in 0.5 dex steps
init file="ism.ini"
c this turns on ism abundances and grains
abundances ism
c
c commands controlling geometry =========
c set sphere since matter in all directions
sphere
stop temperature linear 10
stop thickness 0.1 linear parsecs
stop zone 1
c
c other commands for details =========
c this prevents continuum pumping from exciting H0
case b
c
c commands controlling output =========
c this model has no heavy element molecules, because H2 does not
c really fully form, since Lyman bands are not self-shielded
normalize to 157.6m "C 2"
punch grid "ism_grid.grd" last
punch averages "ism_grid.avr" last
temperature hydrogen 1
end of averages
punch overview "ism_grid.ovr" last
punch dr "ism_grid.dr"
punch heating "ism_grid.het"
punch coolign "ism_grid.col"
c
c
c ism_grid.in
c class ism
c ========================================
this shows an S-curve calculation - make plot showing density as X-axis
and gas pressure (nT) as y-axis
ism_hot_brems.in
generate continuum due to hot ism in high Z,z starburst
title generate continuum due to hot ism in high Z,z starburst
c
c TODO the guess of the thickness of the first zone is badly too small,
c because this model is collisionally ionized, and it used Stromgren
c length - better to use collisional balance and dr - as result of this
c the model takes far too many zones
c
c commands controlling continuum =========
c this sets up coronal equilibrium for temperature and continuum
c cloud is predominantly collisionally ionized
coronal equilibrium, t = 1.e6 K
c
c commands for density & abundances =========
c log of hydrogen density (cm-3), by default a constant density model
hden 1
c chemical composition for a well-evolved star cluster with high Z
c the starburst command generates scale factors that will be used to
c multiply the old solar composition, used in ver 84. this is maintain
c record of physical changes
abundances old solar 84
abundances starburst Z=10 Zsun
c
c commands controlling geometry =========
c the log of the total hydrogen column density, cm-2
stop column density 21
c
c other commands for details =========
c continuum is plotted in Hzay, we me must iterate to predict this
iterate
c
c commands controlling output =========
c normalize to O 7 Lya
normalize to "O 7" 21.60
c this shortens the printout somewhat
print line faint 1
c output options
punch continuum last "ism_hot_brems.con" no title, units keV
punch continuum last "ism_hot_brems.con" no title, units keV
punch overview last "ism_hot_brems.ovr"
c this is only to document the problem described above
punch dr last "ism_hot_brems.dr"
c
c ism_hot_brems.in
c class ism
c ========================================
c
This model generates a large column constant density cloud
similar to the hot phase of the interstellar medium.
The continuum is punched to generate one of the figures in Part 2 of Hazy.
There is a strange feature between 1.7e-3A and 2.2e-3A that is the N emission
(head starting at 1.7e-3A with O absorption at 2.2e-3A. This model is
strongly enriched in heavies so many metal edges, esp O, are optically thick.
ism_jura.in
check rate H2 forms on grain surfaces
title check rate H2 forms on grain surfaces
c
c commands controlling continuum =========
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
blackbody, t = 75 K
intensity 2.7 (total)
c this is hot star continuum
black 30000
intensity 1.90 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
cosmic rays, background
c
c commands for density & abundances =========
hden 0
grains ism, abundance log 0.16 no qheat
init file="ism.ini"
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl=-7 ar-8 fe -6.60
c
c commands controlling geometry =========
sphere
c stop when gas is fully neutral
stop efrac -10
c stop when gas is cold
stop temperature 10 linear
stop zone 1
c
c other commands for details =========
turbulence 2.7 km/sec
c set the gas temperature
constant temper 100
case b
c
c commands controlling output =========
normalize to "C 2" 157.6m
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print ages
punch hydrogen 21 cm last "ism_jura.21cm"
c
c ism_jura.in
c class limit
c ========================================
c
This model started out life as the Tielens & Hollenbach 1985 pdr. The
density was set to unity and the incident radiation field adjusted so
that the two default grains have temperatures near 100K. The model
asserts that the H2 formation rate on grain surfaces is close to the
//>>refer H2 grain physics Jura, M., 1975, ApJ, 197, 575
rate.
ism_opacity.in
generate standard ISM opacity curve
title generate standard ISM opacity curve
c
c commands controlling continuum =========
table agn
ionization parameter -3
c
c commands for density & abundances =========
grains no qheat
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
iterate
c
c commands controlling output =========
print last iteration
punch total opacity last "ism_opacity.opc" units kev
c
c ism_opacity.in
c class limit
c ========================================
c
This example creates the file ism_opacity.opc which tabulates
the total opacity of the gas as a function of energy. These
plots are used in ISM studies to understand the transmission
characteristics along a line of sight. The opacity depends on
the dust to gas ratio, the gas phase abundances, and the level
of ionization, all of which can be changed by altering parameters
given above.
The model is of a 1 cm think parcel of gas which is optically thin
in the Lyman continuum and Lyman lines. As a result the hydrogen
emission line spectrum is close to case C. The model iterates
so that the predicted ionization and emission know about this.
ism_set_cr_rate.in
background cosmic ray ionization by suprathermal electrons only
title background cosmic ray ionization by suprathermal electrons only
c
c commands controlling continuum =========
cosmic ray rate -16
c this will be VERY faint black 50000
ionization parameter -25
black 50000
c
c commands for density & abundances =========
hden 5
element hydrogen ionization 1 0.00000001
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
constant temperature 1000 K
c turn off molecules to test only ions
no molecules
c must turn off charge transfer and its heating, which would dominate,
c since we want to test cr heating and ionization
no charge transfer
no ctheat
c must turn of photoionization by recombining species,
c since we want to test cr heating and ionization
no photoionization
c any time photoionization is turned off, it is best to also
c turn off induced processes, which include pumping in the lyman lines
no induced processes
c
c commands controlling output =========
print ages
c
c ism_set_cr_rate.in
c class limit
c ========================================
c
This test conditions of cosmic ray ionization. Molecules and charge transfer
are disabled so that analytical estimates can be made.
limit_compton_hi_t.in
compton exchange in high temper limit
title compton exchange in high temper limit
c
c commands controlling continuum =========
c as hot as STE bb can be on IEEE 32-bit cpu
c equilibrium temperature should also be 2.51e7
c since gas radiated by true blackbody will equilibriate
c at its temperature
blackbody 7.4 STE
c
c commands for density & abundances =========
hden 10
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
set temperature convergence 0.0001
iterate
c
c commands controlling output =========
c
c limit_compton_hi_t.in
c class limit
c ========================================
c
This is the highest Compton temperature that can be computed in
LTE on an IEEE 32-bit processor. This tests the code in the
high-temperature Compton limit. Temperatures as high as 10^10 K
can be computed successfully on CPUs with longer word lengths,
such as a Cray or the new 64 bit processors.
Checks:
- The equilibrium temperature should be exactly 107.4 K (2.51239107 K).
limit_compton_lo_t.in
compton exchange near low temperature limit
title compton exchange near low temperature limit
c
c commands controlling continuum =========
black linear 3 lte
brems 5
ionizat -5
c
c commands for density & abundances =========
hden -10
init file "hheonly.ini"
eden 0
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
set temperature convergence 0.0001
iterate 3
c
c commands controlling output =========
c
asset itrzn < 5
c limit_compton_lo_t.in
c class limit
c ========================================
c
This tests the code in the low temperature Compton limit. The gas is
illuminated by a 3 K blackbody in thermodynamic equilibrium.
The equilibrium temperature should be exactly 3 K.
It is necessary to add an extra component of free electrons to
test the code in this limit with the eden command.
limit_compton_mid_t.in
Compton limit, test continuum partition
title Compton limit, test continuum partition
c test of thermal equil in limit_compton_mid_t limit; temp should EXACTLY equal 2E5K
c check continuum partition;
c energy range, photon densities, luminosities, follow
c 0.25-1. Q=26.6470 L=15.8190 c 1-1.807 Q=26.8867 L=16.3766
c 1.807-4 Q=27.3042 L=17.0948 c 4-20 Q=27.2041 L=17.3179
c 20 -- Q=22.9038 L=22.9038 c total lumin 17.5597
c nufnu(912A) = 1.8029E+16
c
c commands controlling continuum =========
black body t = 2.e5 K lte
c
c commands for density & abundances =========
hden = 6
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
stop lyman continuum optical depth = -6
set dr 0
c
c other commands for details =========
set temperature convergence .0001 % decrease error on heat-cool match
c
c commands controlling output =========
print line faint .1
print departure coef
c
c limit_compton_mid_t.in
c class limit
c ========================================
c
This tests the behavior of the code in the Compton limit.
The incident continuum is a blackbody in strict thermodynamic equilibrium.
Strict thermodynamic equilibrium is expected for all constituents of the gas.
The input stream also lists the expected photon fluxes for the incident
continuum; this tests the normalization of the continuum,
and its distribution. Grains are included to confirm their behavior
in the LTE limit. The set dr command sets the zone thickness to 1 cm.
Checks:
- Luminosity, photon flux, over various energy intervals, 4*J at 912143.
- Electron temperature exactly 2*105 K.
- Grain temperature forced to 2*105 K by radiative processes.
limit_conserve.in
test that energy is limit_conserved
title test that energy is limit_conserved
* trace temperature convergence 0 2
c
c commands controlling continuum =========
c primary continuum is a hot blackbody
black body, temp=5.5
intensity total -4
c add cosmic background
backbround
c
c commands controlling geometry =========
c don't stop the calculation until all energy used up
c in order to check energy conservation
stop temperature 3 linear
c cloud will go very deep with temperature near 3K without this
stop thickness 20
c
c commands for density & abundances =========
no molecules
metals -0.5
hden 1.0
constant pressure
c
c other commands for details =========
iterate
c want a clean calculation
failures 3
c this allows for more molecular gas beyond what we model
double optical depths
c must increase number of levels since so optically thick
atom co levels 30
c
c commands controlling output =========
punch pressure "limit_conserve.pre"
punch pressure history "limit_conserve.his"
punch cooling "limit_conserve.col"
punch heating "limit_conserve.het"
punch dr "limit_conserve.dr"
punch overview "limit_conserve.ovr1"
punch overview last "limit_conserve.ovr"
punch results "limit_conserve.rlt"
punch element carbon "limit_conserve.car"
punch element neon "limit_conserve.ne"
punch element magnesium "limit_conserve.mag"
punch element silicon "limit_conserve.sil"
punch molecules "limit_conserve.mol"
c
asset line "Q(H)" 4861 0.776
c limit_conserve.in
c class hii pdr
c ========================================
This checks that energy is limit_conserved. The code always checks
that it did not radiate more energy than was absorbed. This
calculation extends well past the photo-dissociation zone into
fully molecular gas, so that all of the incident radiation is absorbed.
Grains, CMB, & CRs are not present so that only the incident radiation field
powers the gas.
Small changes can affect this model to surprising extents because
of the presence of a major thermal front at the H0 - H+ transition region.
limit_eden.in
Martin Gaskell's funny model
title Martin Gaskell's funny model
c used to test that electron density convergence is ok
c hydrogen line spectrum strongly pumped by continuum
c all elements > H have HIGH abundances
c ionization of elements Fe, Mg, Si strongly affected
c by charge transfer
c
c commands controlling continuum =========
black 4000
lumin 27.2
c
c commands for density & abundances =========
hden 5.138
abundances all 1000
c
c commands controlling geometry =========
radius 15
stop zone 1
set dr 0
c
c other commands for details =========
constant temper 5500
c
c commands controlling output =========
print line faint -3
c
c limit_eden.in
c class limit
c ========================================
c
This is mainly a test of the ability of the code to converge a model
with a very strange electron density. The electrons are mainly contributed
by heavy elements, and the gas is only slightly ionized.
Ionization of elements Fe, Mg, Si strongly affected by charge transfer
with other heavy elements.
Checks:
- Electron density is correct.
- Hydrogen line spectrum strongly pumped by continuum.
limit_hi_ion.in
test limit of very high ionization parameter
title test limit of very high ionization parameter
c commands controlling continuum =========
table agn
ionization parameter 15
c commands controlling geometry =========
hden 0
stop zone 1
c other commands for details =========
iterate
c
c commands controlling output =========
punch heating "limit_hi_ion.het"
punch cooling "limit_hi_ion.col"
c
c
c func_hi_ion.in
c class functionality
c ========================================
c
This tests a limit of very high ionization
limit_laser_1.in
test of H ionization in optically thin limit
title test of H ionization in optically thin limit
c
c H cross section is 2.09E-18 cm^2, rec coef is 4.18E-14
c answer is neutral fraction 2.00E-4
c also checks that only 3 iterations needed
c
c commands controlling continuum =========
laser 1.5 Ryd
phi(h) 10
c
c commands for density & abundances =========
hden 1
init file="hheonly.ini"
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
constant temperature = 4
iterate to convergence
c
c commands controlling output =========
print last
c
c limit_laser_1.in
c class limit
c ========================================
c
This checks the calculation of the hydrogen photoionization equilibrium.
The continuum is a laser peaked at 1.5 Ryd, where the hydrogen
photoionization cross section is 2.09*10-18 cm-2.
Checks:
- The hydrogen neutral fraction is nearly 2.00*10-4 (not exact
since laser has finite width).
- Hb emissivity close to high density case A. The predicted
TOTL 4861 intensity should be nearly 2.2 times the expected
case B intensity.
limit_laser_2.in
test of H and HeI ionization in optically thin limit
title test of H and HeI ionization in optically thin limit
c H cross section is 0.927E-18 cm^2, rec coef is 4.18E-13
c answer is Ho/H+ = 4.51e-4
c HeI cross section is 6.54E-18 cm^2, rec coef is 4.32E-13
c answer is Heo/He+ = 6.61e-5
c
c commands controlling continuum =========
laser 2.0 Ryd
intensity -0.3604
c
c commands for density & abundances =========
hden 1
init file="hheonly.ini"
abundances all -10
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature = 4
iterate to convergence
c
c commands controlling output =========
print last
c
c limit_laser_2.in
c class limit
c ========================================
c
This checks the calculation of the hydrogen and helium photoionization
equilibrium. The continuum is a laser peaked at 2.0 Ryd, and so can
only ionize hydrogen and atomic helium.
Checks:
- The hydrogen neutral fraction is nearly Ho/H+=4.51*10-4
(not exact since laser has finite width).
- Hb emissivity close to high density case A. The predicted TOTL 4861
intensity should be nearly 2.2 times the expected case B intensity.
- Helium ionization should be Heo/He+ = 6.61*10-4.
limit_laser_200.in
test of ionization in optically thin Auger-dominated limit
title test of ionization in optically thin Auger-dominated limit
c
c commands controlling continuum =========
laser 200 Ryd
phi(h) 10
c
c commands for density & abundances =========
hden 1
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate to convergence
constant temperature = 5
c
c commands controlling output =========
print last
c
c limit_laser_200.in
c class limit
c ========================================
c
This checks the calculation of ionization
equilibrium. The continuum is a laser peaked at 200 Ryd.
It asserts ionization of C, O, and Fe. their ionization is dominated
by the Auger effect.
Checks:
Auger OK
limit_laser_200_low.in
test of ionization in optically thin Auger-dominated limit
title test of ionization in optically thin Auger-dominated limit
c
c commands controlling continuum =========
laser 200 Ryd
phi(h) 5
c
c commands for density & abundances =========
hden 1
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate to convergence
constant temperature = 4
c
c commands controlling output =========
print last
c
c
c
c
c
c
c limit_laser_200_low.in
c class limit
c ========================================
c
This checks the calculation of ionization
equilibrium. The continuum is a laser peaked at 200 Ryd.
It asserts ionization of C, O, and Fe. their ionization is dominated
by the Auger effect.
Checks:
Auger OK
limit_laser_3.in
test of H, HeI, and HeII ionization in optically thin limit
title test of H, HeI, and HeII ionization in optically thin limit
c H cross section is 1.0E-18 cm^2, rec coef is 4.18E-13
c answer is n(Ho)/n(H+)=4.18e-3
c HeI cross section is 1.51E-18 cm^2, rec coef is 4.32e-13
c answer is n(Heo)/n(He+)=2.86e-4, so Heo/He = 4.83e-7
c HeII cross section is 1.30E-18 cm^2, rec coef is 2.20e-12
c answer is n(He+)/n(He2+)=1.69e-3
c
c commands controlling continuum =========
laser 4.3 Ryd
phi(h) 10
c
c commands for density & abundances =========
hden 1
abundances all -10
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate to convergence
constant temperature = 4
c
c commands controlling output =========
print last
c
c limit_laser_3.in
c class limit
c ========================================
c
This checks the calculation of the hydrogen and helium photoionization
equilibrium. The continuum is a laser peaked at 4.3 Ryd, where it
can fully ionize both hydrogen and helium.
Checks:
- The hydrogen neutral fraction is nearly 4.18*10-4
(not exact since laser has finite width).
- Helium ion: The ratio He+/He++ should be 1.69*10-3 and
the ratio Heo/He+ should be 2.86*10-4.
- Hb emissivity should be close to high-density case A. The predicted
TOTL 4861 intensity should be nearly 2.2 times the expected
case B intensity.
limit_lowd0.in
test low density limit, this and lowdm6 should get same results
title test low density limit, this and lowdm6 should get same results
c
c commands controlling continuum =========
table agn
ionization parameter -2
c
c commands for density & abundances =========
hden 0
abundances old solar 84
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate
c
c commands controlling output =========
print line sort wavelength range 100 1m
punch cooling last "limit_lowd0.col"
punch continuum last "limit_lowd0.con"
c
init file="lowd.dat"
c
c limit_lowd0.in
c class limit
c ========================================
c
this test case is paired with lowdm6.in
both tests read in the same set of asserts, those contained in
the file lowd.dat, and they should get exactly the same answer
this is also the test of the print lines intensity command
limit_lowden.in
test optically thin model that extends to very low densities
title test optically thin model that extends to very low densities
c
c commands controlling continuum =========
table agn
ionization parameter 0
c
c commands for density & abundances =========
hden -2 to the -2 power
init file="ism.ini"
c
c commands controlling geometry =========
sphere
c this will take density down to 1e-8, below stated limit
radius 10 13
c
c other commands for details =========
c
c commands controlling output =========
punch overview "limit_lowden.ovr" last
c
c limit_lowden.in
c class limit
c ========================================
c
This model is optically thin, with density falling off as inverse square law,
so ionization and temperature should be nearly constant.
if outer radius increased by 2 dex problem with level3 will appear,
several li seq lines (OVI, NeVIII) will fluctuate when density about
1e-9
limit_lowdm6.in
test low density limit, this and limit_lowdm6 should get same results
title test low density limit, this and limit_lowdm6 should get same results
c
c commands controlling continuum =========
table agn
ionization parameter -2
c
c commands for density & abundances =========
hden -6
abundances old solar 84
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate
c
c commands controlling output =========
print line sort wavelength range 100 1m
punch cooling last "limit_lowdm6.col"
punch continuum last "limit_lowdm6.con"
c
init file="lowd.dat"
c
c limit_lowdm6.in
c class limit
c ========================================
c
this test case is paired with lowd0.in
both tests read in the same set of asserts, those contained in
the file lowd.dat, and they should get exactly the same answer
this also tests the print line sort range command
limit_lowion_low.in
test conditions of very low ionization matrix/simple solver
title test conditions of very low ionization matrix/simple solver
c
c commands controlling continuum =========
blackbody 50000
ionization parameter -30
c
c commands for density & abundances =========
hden 2
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
atom h-like matrix lowt
atom he-like matrix lowt
no molecules
constant temperature 100
c
c commands controlling output =========
c
c limit_lowion_low.in
c class limit
c ========================================
c
The lowion_pops.in and limit_lowion_low.in models form a pair
that have identical boundary conditions but use the two
different hydrogenic level populations solvers. The results
should agree. lowion_pops.in uses the full solution with
the associate matrix inversion. This can fail under conditions
of extreme low ionization due to numerical instabilities and
roundoff. The solver used in limit_lowion_low.in is much simpler
and will work for any conditions.
The model is almost totally molecular.
limit_lowion_pops.in
test conditions of very low ionization matrix/simple solver
title test conditions of very low ionization matrix/simple solver
c
c commands controlling continuum =========
blackbody 50000
ionization parameter -30
c
c commands for density & abundances =========
hden 2
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
atom h-like matrix pops
atom he-like matrix pops
no molecules
constant temperature 100
c
c commands controlling output =========
c
c
c limit_lowion_pops.in
c class limit
c ========================================
c
The limit_lowion_pops.in and lowion_low.in models form a pair
that have identical boundary conditions but use the two
different hydrogenic level populations solvers. The results
should agree. limit_lowion_pops.in uses the full solution with
the associate matrix inversion. This can fail under conditions
of extreme low ionization due to numerical instabilities and
roundoff. The solver used in lowion_low.in is much simpler
and will work for any conditions.
The model is almost totally molecular.
limit_lte_h_t50_cion.in
test collisional ionization only, no excitation, should be in LTE
title test collisional ionization only, no excitation, should be in LTE
c
c commands controlling continuum =========
black body, T = 50000 K
ionization parameter -8
c
c commands for density & abundances =========
hden = 20
set eden 20
init file="hheonly.ini"
c
c commands controlling geometry =========
stop lyman optical depth -6
set dr -10
stop zone 1
c
c other commands for details =========
constant temperature = 50000 K
c want pure collisional model
no induced processes (OK)
atom h-like collisional exciation off
atom h-like collisions l-mixing off
atom h-like levels large
atom he-like collisional exciation off
atom he-like collisions l-mixing off
atom he-like collapsed levels 20
iterate
c
c commands controlling output =========
print h-like departure coefficients
print he-like departure coefficients
print last iteration
print populations h-like
print departure coefficients
c
c limit_lte_h_t50_cion.in
c class limit
c ========================================
c
This is the limiting case pure hydrogen collisional ionization,
There are no excitation or l-mixing collisions, so this tests whether
collisional ionization - three body recombination
works in detailed balance.
limit_lte_h_t50_coll.in
test collisional excitation only, very high density to force H to LTE
title test collisional excitation only, very high density to force H to LTE
c
c commands controlling continuum =========
black body, T = 50000 K
ionization parameter -8
c
c commands for density & abundances =========
hden = 20
init file="hheonly.ini"
c
c commands controlling geometry =========
stop lyman optical depth -6
stop zone 1
set dr 0
set eden 20
c
c other commands for details =========
no induced processes (OK)
atom h-like collisions l-mixing off
atom h-like levels large
atom he-like collapsed levels 20
atom he-like collisions l-mixing off
constant temperature = 50000 K
iterate
c
c commands controlling output =========
print h-like departure coefficients
print h-like populations
print he-like departure coefficients
print populations h-like
c
c limit_lte_h_t50_coll.in
c class limit
c ========================================
c
This checks that the model H atom goes to LTE at high densities.
chng 06 aug 24, had not included collisional ionization, and so he-like
departure coefficients were very large, around 202. comments said there
were problems. turned on collisional ionization, no problems noted
chng 06 jul 22
with RP changes in high-n n-changing collisions the rates are now much
smaller - needed to change density to be far higher and several
quantities changed. at lower density (1e18 cm-3) the populations are very
unphysical and runaway maser now occurs. this is only a homework problem
and intended to only test n-changing collisions. with higher density
this test is done.
limit_lte_he1_coll.in
test hei atom at high densities
title test hei atom at high densities
c
c commands controlling continuum =========
laser 2
ionization -3
c
c commands for density & abundances =========
hden 18
set eden 18
element helium abundance -1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
set dr -10
c
c other commands for details =========
atom he-like levels 10
atom he-like collapsed levels 20
constant temper 20000 K
iterate
c
c commands controlling output =========
print he-like departure helium
normalise to "He 1" 4471
print line faint -2
c
c limit_lte_he1_coll.in
c class limit
c ========================================
c
test whether he-like ion populations go to lte in high density limit.
The level populations should be in LTE, and the departure coefficients
should be unity.
limit_lte_he1_ste.in
test He I atom LTE at high densities
title test He I atom LTE at high densities
c
c commands controlling continuum =========
blackbody 50000 STE
c
c commands for density & abundances =========
hden 18
* set eden 18
c a high he abundance so that it dominates
element helium abundance 1
init file "hheonly.ini"
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
atom he-like levels 7
atom he-like collapsed levels 30
iterate
c
c commands controlling output =========
print he-like departure helium
print h-like departure hydrogen
normalise to "He 1" 4471
print line faint -2
punch heating "limit_lte_he1_ste.het"
punch cooling "limit_lte_he1_ste.col"
c
c limit_lte_he1_ste.in
c class limit
c ========================================
c
test whether a gas dominated by He goes to LTE in high-density limit.
The level populations should be in LTE, the departure coefficients
should be unity, and the temperature equal to the BB temp.
limit_lte_hhe_coll_t50.in
high electron density approach to LTE
title high electron density approach to LTE
c test from Ferland and Rees 88, collisions drive H to LTE
c collisions should drive all departure coef to unity
c
c commands controlling continuum =========
black body t = 50000
ionization parameter -5
c
c commands for density & abundances =========
hden = 19
init file="hheonly.ini"
set eden 20
c
c commands controlling geometry =========
set dr -10
stop zone 1
stop lyman optical depth -6
c
c other commands for details =========
constant temperature = 50000
atom h-like levels large
iterate
c
c commands controlling output =========
print last
print departure coefficients
print populations h-like
punch dr "limit_lte_hhe_coll_t50.dr" last
c
c limit_lte_hhe_coll_t50.in
c class limit
c ========================================
c
This model is a test of the behavior of hydrogen and helium
in the high density, collision dominated, limit.
The temperature is preset, the hydrogen density is set
to a very high value, and the ionization parameter is very low.
The resulting model is collision dominated, so this case
checks that the collision physics occurs in detailed balance.
The predicted departure coefficients should all equal unity.
The set dr command sets the zone thickness to 1 cm.
Checks:
- Hydrogen departure coefficients exactly unity.
- Helium departure coefficients near unity.
(Density not high enough to bring helium departure coefficients
exactly to unity.)
- H-, H2, H2+ H3+, and HeH+ departure coefficients exactly unity.
`
limit_lte_hhe_induc.in
half H, He gas with induced BF processes dominate, go to LTE
title half H, He gas with induced BF processes dominate, go to LTE
c
c commands controlling continuum =========
black body, t = 50000 LTE
c
c commands for density & abundances =========
hden 4
init file="hheonly.ini"
c
c commands controlling geometry =========
stop lyman continuum optical depth -6
set dr -10
stop zone 1
c
c other commands for details =========
iterate
atom h-like levels large
set temperature convergence 0.002
c
c commands controlling output =========
print populations h-like
print heating
print departure coef
print line faint -1
c
c limit_lte_hhe_induc.in
c class limit
c ========================================
c
This is a H, He-only gas that is optically thin in the Lyman continuum. It is
irradiated by a blackbody in strict thermodynamic equilibrium. The tests confirm
that the gas temperature equilibriates close to the black body temperature.
limit_lte_hhe_ste.in
thermal equil black body LTE limit from Ferland and Rees 1988
title thermal equil black body LTE limit from Ferland and Rees 1988
c this tests whether thermal processes go to LTE
c
c commands controlling continuum =========
black body, t = 50000 LTE
c
c commands for density & abundances =========
hden 10
init file="hheonly.ini"
c
c commands controlling geometry =========
stop lyman optical depth -6
set dr -10
stop zone 1
c
c other commands for details =========
atom h-like levels large
iterate
c
c commands controlling output =========
print heating
print populations h-like
print departure coef h-like
print departure coef he-like
print last
c
c limit_lte_hhe_ste.in
c class limit
c ========================================
c
This is the ultimate test of the behavior of the code in the
strict thermodynamic equilibrium limit. The temperature is not
held constant, so the resulting equilibrium temperature determines
whether cooling processes are treated properly in the detailed balance limit.
The equilibrium temperature should be exactly 5*104 K, and all
departure coefficients should equal unity. A small amount of
grains are included to check that the grain thermal balance is
handled properly in this limit.
Checks:
- Electron temperature exactly 5*104 K.
- Departure coefficients unity.
limit_lte_hminus.in
hminus test of LTE
title hminus test of LTE
c
c commands controlling continuum =========
blackbody 5000 STE
c
c commands for density & abundances =========
hden 10
c
c commands controlling geometry =========
stop zone 1
set dr 0
c
c other commands for details =========
iterate
constant temperature 5000
c
c commands controlling output =========
print populations h-like hydrogen
print departure coef
c
c limit_lte_hminusin
c
c class limit
c ========================================
c
This checks that the negative hydrogen ion goes
to thermodynamic equilibrium when irradiated by
a blackbody in thermodynamic equlibrium. It was
originally presented in
>>refer H- test Ferland, G. J., & Persson, S. E. 1989, ApJ, 347, 656
limit_lte_metal.in
LTE (actually strict thermodynamic equilibrium) with metals
title LTE (actually strict thermodynamic equilibrium) with metals
c this tests whether thermal processes go to LTE
c
c commands controlling continuum =========
black body, t = 20000 LTE
c
c commands for density & abundances =========
hden 10
abundances starburst 5
c
c commands controlling geometry =========
stop zone 1
stop lyman optical depth -6
set dr -10
c
c other commands for details =========
atom h-like levels large
iterate
c
c commands controlling output =========
print departure coef
print populations h-like
print last
punch heating "limit_lte_metal.het" last
punch cooling "limit_lte_metal.col" last
c
c limit_lte_metal.in
c class limit
c ========================================
c
This checks that the code goes to strict thermodynamic equilibrium
for the case of a metal rich gas exposed to a true black body.
The many heavy element lines should dominate cooling, so this is
a test that the multilevel atoms go to LTE in the radiation-dominated limit.
Checks:
- Temperature should equilibrate at 20000 K.
- Departure coefficients should equal unity.
limit_recoil_ion.in
test compton recoil ionization of hydrogen
title test compton recoil ionization of hydrogen
c
c commands controlling continuum =========
laser 2000 ryd
ionization parameter -5
c
c commands for density & abundances =========
hden 0
init "honly.ini"
c
c commands controlling geometry =========
set dr 14
stop zone 10
c
c other commands for details =========
iterate
no secondary ionization
constant temperature 1000
c
c commands controlling output =========
c
c
c limit_recoil_ion.in
c class limit
c ========================================
c
H ionization is totally due to recoil ionization in this model.
The assert checks the final hydrogen ionization.
limit_strom.in
check pure hydrogen Stromgren sphere
title check pure hydrogen Stromgren sphere
c the answer is R(Stromgren) = 4.16E17 cm
c
c commands controlling continuum =========
blackbody 50000 K
q(h) 49
c
c commands for density & abundances =========
hden 4
init file "honly.ini"
abundances all -10
c
c commands controlling geometry =========
sphere static
radius 16
stop efrac -0.5
c
c other commands for details =========
no level2
turbulence 20 km/s
constant temper 7500
iterate
c
c commands controlling output =========
print last iteration
punch overview last "limit_strom.ovr"
punch results last "limit_strom.rlt"
punch dr last "limit_strom.dr"
c
c limit_strom.in
c class limit
c ========================================
c
This case checks that the code computes the geometry and emissivity
correctly for a pure hydrogen spherical shell. The low temperature
is chosen to avoid collisional ionization. The model stops at the
Ho-H+ ionization front. The turbulence is to prevent the Balmer lines
from becoming optically thick.
Checks
- Outer radius should be 4.16391017 cm.
- Predicted Hb, case B Hb, and Q(H) Hb, all agree.
limit_supra.in
test very high levels of seciondary ionization, like SN envelope
title test very high levels of seciondary ionization, like SN envelope
c
c commands controlling continuum =========
c this continuum will be VERY faint
table agn
ionization parameter -25
c secondary ioinzations will dominate
set csupra -5
c
c commands for density & abundances =========
hden 5
abundances old solar 84
c
c commands controlling geometry =========
stop temperature 20
stop zone 10
set dr 0
c
c other commands for details =========
constant temperature 1000 K
iterate convergence
c
c commands controlling output =========
print ages
punch heating "limit_supra.het"
punch ionizing continuum "limit_supra.ion" last
punch molecules "limit_supra.mol"
c
c limit_supra.in
c class limit
c ========================================
c
This model computes the ionization within cool gas that is totally ionized
by suprathermal secondary electrons.
limit_vbhum.in
test against Van Blerkom and Hummer, fig 4
title test against Van Blerkom and Hummer, fig 4
c
c test hydrogen ground state rec effic against vb+h exact results
c this is their case e) - "zero condition"
c their answer for H0/Htot at the illuminated edge is approx 5.8E-4,
c and a Stromgren radius of approximately 7.7E16 cm
c
c commands controlling continuum =========
blackbody 50000 K
phi(h) 12.30103
c commands for density & abundances =========
hden 4
init "honly.ini"
abundances all -10
c
c commands controlling geometry =========
stop efrac -0.5
c
c other commands for details =========
diffuse ots
no level2 lines
constant temper 4
iterate
c
c commands controlling output =========
punch overview last "limit_vbhum.ovr"
punch dr last "limit_vbhum.dr"
c
c
c limit_vbhum.in
c
c class limit
c ========================================
c
This is a test of the treatment of the diffuse fields, their transfer,
and their effects on the ionization structure of a nebula. The comparison
is made against the exact calculation published by Van Blerkom and
Hummer (1967). The geometry is open, that is, similar to that assumed
in most BLR calculations.
>>refer H ionization Van Blerkom, D., & Hummer, D. G. 1967, MNRAS, 137, 353
The diffuse ots command is entered in order to reproduce the
Van Blerkom and Hummer results. The default assumption, outward only,
does not agree as well. I changed the default from OTS to outward only
to be in better agreement with predictions by Harrington and Rubin
at the Lexington meeting. They have not checked whether their codes
are in agreement with the Van Blerkom and Hummer paper.
Checks:
- Neutral fraction at illuminated face 5.8*10-4.
- Location of ionization front at 7.8*1016 cm.
- 34TOTL 486134 and 34CA B 486134 agree; both slightly lower than 34Q(H) 486134.
- Answers with OTS agree with 1967 results.
limit_veryveryfast.in
very fast model for running with debuggers
title very fast model for running with debuggers
c
c commands controlling continuum =========
table agn
ionization parameter -2
c
no file opacity
set dr 0
c
c commands for density & abundances =========
hden 9
init file "hheonly.ini"
element oxygen on
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
no level 2
constant temperature 4
atom h-like levels very small
atom he-like levels very small
c
c commands controlling output =========
print last iteration
punch dr last "limit_veryveryfast.dr"
c
c
c limit_veryveryfast
c class limit
c ========================================
c
This is meant to be a very fast calculation to
use when running extensive debug-enabled runtimes.
lines.in
create output file with list of lines
title create output file with list of lines
atom h2 levels large
atom feii
test
c
c commands controlling continuum =========
c commands for density & abundances =========
c commands controlling geometry =========
c other commands for details =========
c
c commands controlling output =========
c this creates the file "lines.lab that contains a list
c of all emission labels
punch line labels "lines.lab"
c this creates the file "lines.lis" that contains a list
c of the emission lines that are transported.
punch line data "lines.lis"
c
c
c lines.in
c class function
c ========================================
c
This runs the standard "test" case, and then creates
the line data and labels files. Test by itself includes
many asserts, so no further asserts are needed here.
The file lines.lab is a useful list of all lines predicted by the code.
Cut and paste this into other places when you need to find a particular
emission line.
The lines lines.dat gives atomic data for all the lines, and their
critical density at 10000 K.
nlr_lex00.in
Hagai's nlr_lex00 model for Lexington Meeting
title Hagai's nlr_lex00 model for Lexington Meeting
c
c commands controlling continuum =========
phi(h) 12.47712
c shape of continuum is interpolated table
interpolate (0 -10) (0.08 -10) (0.1 1) (3676 -4.935) (4700 -11)
continue (4750 -20) (7.4e6 -30)
c
c commands controlling geometry =========
stop column 22
c
c commands for density & abundances =========
hden 4
init file="c84.ini"
abundances -1 c-3.52 n-4 o-3.097 ne-4 na-9
continue mg-4.523 al-8 si-4.523 s-4.824 a-9 ca-8 fe-8 ni-9
c
c other commands for details =========
iterate
c
c commands controlling output =========
* normalize to "totl" 4861 of 0.923
print line faint .01
print diffuse continuum
c add sum of these lines to printout for Stoy ratio
print line sum
TOTL 1216
He 1 5876
TOTL 1909
TOTL 1549
N 2 6584
N 2 6548
N 3 57.21m
totl 1750
totl 1486
o 1 6300
o 1 63.17m
TOTL 3727
totl 1665
O 3 51.80m
O 3 5007
O 3 4959
totl 1402
Ne 3 15.55m
Ne 3 3869
Ne 3 3968
Ne 4 2424
Ne 5 3426
Ne 5 3346
totl 2798
si 2 34.81m
s 2 6720
S 3 18.67m
S 3 9532
S 3 9069
S 4 10.51m
end of lines
c
punch overview last "nlr_lex00.ovr"
punch dr last "nlr_lex00.dr"
punch results last "nlr_lex00.rlt"
c
table lines "LineList_NLR.dat"
c nlr_lex00.in
c class nlr_lex00
c ========================================
This is one of the test cases from the Lexington (1993) Meeting suite of nebulae.
It is a grain-free NLR model.
nlr_liner.in
nlr_liner model
title nlr_liner model
c a constant pressure (gas+radiation) model of a nlr_liner cloud
c in the spirit of Ferland and Netzer 1983
c
c commands controlling continuum =========
table agn
ionization parameter -3
c
c commands for density & abundances =========
hden 6
init file "ism.ini"
abundances old solar 84
metals 0.3
c
c commands controlling geometry =========
stop thickness 16.3
stop column density 23
constant pressure
c
c other commands for details =========
iterate to convergence ;must iterate to get radiation pressure correctly
c
c commands controlling output =========
punch overview last "nlr_liner.ovr"
punch element nitrogen last "nlr_liner.nit"
punch dr last "nlr_liner.dr"
punch results last "nlr_liner.rlt"
c
c nlr_liner.in
c class nlr
c ========================================
This is a model somewhat like the Liner parameters proposed by
Ferland and Netzer (1983). A second iteration is performed to
allow the calculation of the line radiation pressure.
>>refer nlr_liner model Ferland, G. J., & Netzer, H. 1983, ApJ, 264, 105
nlr_liner_grains.in
liner model with grains
title liner model with grains
c
c a constant pressure (gas+radiation) model of a liner cloud
c in the spirit of Ferland and Netzer 1983
c
c commands controlling continuum =========
ionization parameter -3
table agn
c
c commands for density & abundances =========
hden 6
c this turns on both grains and depleted abundances
abundances ism no grains
grains ism no qheat single
init file "ism.ini"
c
c commands controlling geometry =========
stop column density 23
c
c other commands for details =========
constant pressure
iterate to convergence ;must iterate to get radiation pressure correctly
c
c commands controlling output =========
punch overview last "nlr_liner_grains.ovr"
punch dr last "nlr_liner_grains.dr"
punch continuum last "nlr_liner_grains.con" units microns
punch results last "nlr_liner_grains.rlt"
c
c nlr_liner_grains.in
c class nlr
c ========================================
c
This is a model somewhat like the Liner parameters proposed by
Ferland and Netzer (1983). A second iteration is performed to
allow the calculation of the line radiation pressure.
It includes grains and is so more realistic.
>>refer liner model Ferland, G. J., & Netzer, H. 1983, ApJ, 264, 105
nlr_paris.in
paris meeting NLR model
title paris meeting NLR model
c
c commands controlling continuum =========
interpolate (0 -10) (0.073 -10) (0.074 1) (7352 -4)
continue (7400 -15) (7.353e6 -20)
ionization parameter -2
c
c commands controlling geometry =========
stop lyman continuum optical depth 4
c
c commands for density & abundances =========
hden 3
init file="c84.ini"
abund -1 C-3.5229 N-4. O-3.22185 ne-3.82391 na=-8 mg-4.5229 al=-8
continue si-4.5229 s-4.82391 ar-8 ca=-8 fe-7 ni-8
c
c other commands for details =========
c this large group of lines are turned off so that the model runs
c 30 percent faster - lines are mainly fuv and x-ray, which we will
c not track in this model
no level2
iterate
c
c commands controlling output =========
punch overview last "nlr_paris.ovr"
punch dr last "nlr_paris.dr"
punch results last "nlr_paris.rlt"
c
c nlr_paris.in
c class nlr
c ========================================
This is the NLR model presented in the Meudon meeting on model nebulae.
The init file is entered to make the code behave more like version 84.
Checks:
- init file works
nova_dqher.in
(roughly) Ferland et al. DQ Her model
title (roughly) Ferland et al. DQ Her model
c model of room temperature ionized cloud arund old nova DQ Her
c roughly that of
c tests behaviour of code at very low temperatures
c
c commands controlling continuum =========
c flat continuum, followed by nu^-2 power law
interpolate (0 0) (0.3 0) (8.e6 -14.8)
luminosity total 34
c
c commands for density & abundances =========
hden = 2.
abundances nova
grains no qheat single
c
c commands controlling geometry =========
filling factor 0.667
covering factor 0.667
radius 16.5682 , thickness = 16.14613
c need this since gas is so cool
stop temp = 100
c
c other commands for details =========
# >>chng 06 aug 02, add cosmic rays, should always have been included
cosmic rays background
iterate
turbulence 300
age 70 years %it exploded in the 1930's
c
c commands controlling output =========
normalize to 4861 label="CaBo" intensity 1.169
print line all
punch overview last "nova_dqher.ovr"
punch results last "nova_dqher.rlt"
punch dr last "nova_dqher.dr"
punch transmitted continuum "nova_dqher.trn" last
c
c nova_dqher.in
c class nova
c ========================================
c
This tests the code39s behavior in the limit posed by the
metal rich low density nebula surrounding DQ Her (Ferland et al. 1984).
Checks:
- Thickness exact
- Thermal stability
nova_photos.in
dense nova_photos shell
title dense nova_photos shell
c
c commands controlling continuum =========
c this is the result of this command
c table star kurucz 35000; use old atlas atmosphere
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_kurucz_35000.dat"
absolute bolometric magnitude -8.1 ; set luminosity of central star
c
c commands for density & abundances =========
hden 10 ;log of total hydrogen density (cm-3) at illuminated face
abundances old solar 84 ; use old solar abundances for record keeping
c
c commands controlling geometry =========
sphere ; model is a full sphere covering the central star
radius 12 to 14; log of inner and outer radii (cm)
c
c other commands for details =========
iterate to convergence
c
c commands controlling output =========
print line faint -1
punch overview "nova_photos.ovr" last
punch results "nova_photos.rlt" last
punch dr "nova_photos.dr" last
punch heat "nova_photos.het" last
punch cool "nova_photos.col" last
c
c nova_photos.in
c class nova_photos
c ========================================
c
this model is intensely affected by continuum pumping of atoms. The
hydrogen ionizaiton is by lyman line pumping, followed by photoionization
from excited states.
o8_caseb.in
o8_caseb O VIII case B
title o8_caseb O VIII case B
c
c commands controlling continuum =========
table agn
ionization parameter 4
c
c commands for density & abundances =========
hden = 8
init file="ism.ini"
element oxygen abundance 1
set eden 9
c
c commands controlling geometry =========
c set following so the n^2 dr is unity
set dr -18
stop zone 1
c
c other commands for details =========
atom h-like oxygen levels large
case b hummer no photoionization
no induced processes (OK)
constant temperature = 1.e5
c
c commands controlling output =========
punch results "o8_caseb.rlt"
c
c caseb_c6.in
c class limit
c ========================================
c
This test case compares the predictions of the multi-level H0like O VIII atom
with the Storey and Hummer (1995) results. The set dr command sets the
zone thickness to 1 cm. The case b command sets Lyman line optical depths
to very large values.
Checks:
- Neutral fractions
- H* emissivity
- Relative line intensities
optimal.in
test optimizer limits on stellar atmosphere grids
title test optimizer limits on stellar atmosphere grids
c
c commands controlling continuum =========
TABLE STAR COSTAR AGE= 6.302320 LOG , MASS= 36.000000
ionization par -3
c
c commands for density & abundances =========
HDEN=5.002525
init file "hheonly.ini"
element oxygen on
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
optimize tolerance 0.03
CONStant TEMP 3.999435
optimize phymir sequential
optimize iteration 500
c the following spectrum was derived by running this sim at hden 5, t=4,
c and log(age)=6.3
c >>chng 06 dec 01, change line intensities due to new Badnell DR rates
optimize lines
he 2 303.8 28.5845
o ii 3729 0.0252
o ii 3726 0.0713
o ii 7323 0.0928
o ii 7332 0.0748
o 3 5007 2.4914
totl 4363 0.0235
end of lines
c
c commands controlling output =========
print line faint -2
c
c
c optimize_phymir.in
c class optimizer
c ========================================
c
optimize_phymir.in
test optimizers, spectrum computed with hden 5, temp 4
title test optimizers, spectrum computed with hden 5, temp 4
c subplex, phymir, powell
c
c commands controlling continuum =========
black 40000 K
ionization par -3
c
c commands for density & abundances =========
hden 4 vary
init file "hheonly.ini"
element oxygen on
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
optimize tolerance 0.03
constant temper 4.3 vary
optimize phymir
optimize iteration 100
c the following spectrum was derived by running this sim at hden 5 and t=4
c >>chng 06 apr 03, backdate to old [OII] As and rederive spectrum
optimize lines
o ii 3729 0.0918
o ii 3726 0.2565
o ii 7323 0.3152
o ii 7332 0.2317
o 3 5007 1.537
totl 4363 0.0148
end of lines
c
c commands controlling output =========
print line faint -2
c
c
c optimize_phymir.in
c class optimizer
c ========================================
c
This checks whether the optimizer can recover a known solution. The line
spectrum was calculated at T = 1e4 K and and n_H=1e5 cm^-3,
and resulted in the given electron
density. There are four models that are copies of this file,
amoeba.in, phymir.in, powell.in, and subplex.in
optimize_subplex.in
test optimizers, spectrum computed with hden 5, temp 4
title test optimizers, spectrum computed with hden 5, temp 4
c subplex, phymir, amoeba, powell
c
c commands controlling continuum =========
black 40000 K
ionization par -3
c
c commands for density & abundances =========
hden 4 vary
init file "hheonly.ini"
element oxygen on
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
constant temper 4.3 vary
optimize iteration 100
optimize tolerance 0.01
c the following spectrum was derived by running this sim at hden 5 and t=4
c >>chng 06 apr 03, backdate to old [OII] As and rederive spectrum
optimize lines
o ii 3729 0.0918
o ii 3726 0.2565
o ii 7323 0.3152
o ii 7332 0.2317
o 3 5007 1.537
totl 4363 0.0148
end of lines
c
c commands controlling output =========
print line faint -2
c
c
c optimize_subplex.in
c class optimizer
c ========================================
c
This checks whether the optimizer can recover a known solution. The line
spectrum was calculated at T = 1e4 K and and n_H=1e5 cm^-3,
and resulted in the given electron
density. There are four models that are copies of this file,
amoeba.in, phymir.in, powell.in, and subplex.in
orion_hii_dist_grn.in
conditions similar to Orion nebula blister
title conditions similar to Orion nebula blister
c orion_hii_dist_grn and orion_hii_single_grn form a pair
c which are identical expect for resolved grains
c
c commands controlling continuum =========
c this is the result of this command
c table star kurucz 39600 K
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_kurucz_39600.dat"
phi(h) 13.0
c
c
c this is to include velocity field in radiative transfer, so not in pressure
turbulence 8 km/sec no pressure
c
c commands for density & abundances =========
hden 4
abundances hii region no grains
c want to leave quantum heating on for a few cases,
c this is one, and we check the continuum
grains orion distributed
c
c commands controlling geometry =========
sphere
constant pressure
c this is the account for continued matter beyond i front
double optical depths
c
c other commands for details =========
iterate 3
c
c commands controlling output =========
c print lines as surface brightness
print lines surface brightness arcsec
print diffuse continuum
c want to print line intensities are surface brightness, per arcsec^2
print lines surface brightness arcsec
print diffuse continuum
print line faint -1.5
punch overview last "orion_hii_dist_grn.ovr"
punch dr last "orion_hii_dist_grn.dr"
punch continuum last "orion_hii_dist_grn.con" units microns
punch overview last "orion_hii_dist_grn.ovr"
c
c orion_hii_dist_grn.in
c class hii
c ========================================
This is a model similar in spirit to the blister geometry HII region
model computed by Baldwin et al. (1991). Orion grains are turned on
with the abundances command. The constant pressure command does a
hydrostatic equilibrium structure. The predicted emission line spectrum is
affected by the reddening of the internal grains. The resulting t2
analysis produces artificial results as a result.
Checks:
- Pressure convergence
- Helium neutral fraction
orion_hii_open.in
conditions similar to Orion nebula blister
title conditions similar to Orion nebula blister
c orion hii region with open geometry
c
c commands controlling continuum =========
c this is the result of this command
c table star kurucz 39600 K
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_kurucz_39600.dat"
phi(h) 13.0
c
c
c this is to include velocity field in radiative transfer
turbulence 8 km/sec no pressure
c
c commands for density & abundances =========
hden 4
abundances hii region no grains
grains orion
c
c commands controlling geometry =========
c sphere is NOT included in this one
constant pressure
c this is the account for continued matter beyond i front
double optical depths
c
c other commands for details =========
iterate 2
c
c commands controlling output =========
print line faint -5
c print lines as surface brightness
print lines surface brightness arcsec
print line optical depths
print diffuse continuum
c want to print line intensities are surface brightness, per arcsec^2
print lines surface brightness arcsec
print diffuse continuum
punch overview last "orion_hii_open.ovr"
punch dr last "orion_hii_open.dr"
punch continuum last "orion_hii_open.con" units microns
punch transmitted continuum last "orion_hii_open.trn"
punch overview last "orion_hii_open.ovr"
punch fine opacity last "orion_hii_open.opc" range 0.7 0.9 ryd
c
c orion_hii_open.in
c class hii
c ========================================
This is a model similar in spirit to the blister geometry HII region
model computed by Baldwin et al. (1991). Orion grains are turned on.
The constant pressure command does a
hydrostatic equilibrium structure. The predicted emission line spectrum is
affected by the reddening of the internal grains. The resulting t2
analysis produces artificial results as a result.
this has an open geometry
Checks:
- Pressure convergence
- Helium neutral fraction
orion_hii_pdr.in
constant pressure orion into pdr
title constant pressure orion into pdr
c
c commands controlling continuum =========
c the incident continuum is two parts
c star and flux of photons striking it
c this is the photosphere of the OVI star, its temperature and Q(H)
c this is the result of this command
c table star kurucz 39600 K
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_kurucz_39600.dat"
Q(H) 49
c this adds the observed hot brems, its temperature and the flux of
c photons striking the cloud
brems 6
phi(h) 10
c add cosmic rays, which are important in the pdr
cosmic rays, background
c
c commands controlling geometry =========
c this sets the separation from the star and the face of the cloud
radius 17.4507
c this turns off the stop temperature opption
c the sim will not stop due to temperature
stop temperature off
c this is typical of a gmc, larson 1981
stop total hydrogen column density 22
c this sets the thickness of the HII region & PDR
stop thickness 0.5 linear parsec
c assume constant pressure
c constant gas pressure
constant pressure
c this will result in a milli gauss B-field in molecular region
magnetic field -5 gauss
c the orion environment full covers the star, so turn on sphere
sphere
c we have a spherical geometry but want to simulate observing
c through a spectrometer's slit. use the aperture
c command for this
aperture beam
c
c other commands for details =========
c mimic existance of unmodeled molecular gas
double
c iterate since lines optically thick
iterate
c the observed microturbulence, partially a flow, so not included in pressure
turbulence 8 km/sec no pressure
c set the line width so lines appear on the punch continuum
set punchLwidth 10 km/s
c there should be no failures at all
failures 3
c
c this model goes quite deep, needs lots of zones
set nend 2200
c
c commands for density & abundances =========
c this is the log of the initial H density, cm-3
hden 4
c this will speed up the calculation a bit
init file="ism.ini"
c this uses HII region abundances, but no grains
abundances hii region no grains
c this turns on orion grains, but leaves quantum heating off to save time
grains orion no qheat
c turn on PAHs, with an abundance that depends on H0 fraction,
c as suggested by long-slit observations of Orion bar,
c with an abundance 3x larger than default built into the code
grains pah function 3 no qheat
c
c commands controlling output =========
normalize to "O 1" 63.17m
c print lots of faint CO lines
print line faint -4
punch overview last "orion_hii_pdr.ovr"
punch heating "orion_hii_pdr.het"
punch cooling "orion_hii_pdr.col"
punch dr last "orion_hii_pdr.dr"
punch results last "orion_hii_pdr.rlt"
punch continuum last "orion_hii_pdr.con" units microns
punch hydrogen 21 cm last "orion_hii_pdr.21cm"
punch hydrogen lya last "orion_hii_pdr.lya"
punch element silicon last "orion_hii_pdr.sil"
punch grain extinction last "orion_hii_pdr.grnext"
punch grain charge last "orion_hii_pdr.grnchr"
punch grain potential last "orion_hii_pdr.grnpot"
punch grain temperature last "orion_hii_pdr.grntem"
punch grain charge last "orion_hii_pdr.grnchr"
c
c
c orion_hii_pdr.in
c class hii pdr
c ========================================
This model is as close as I can get to the Tielens & Hollenbach (1985) pdr,
except that it starts at the face of the hii region.
orion_hii_pdr_fast.in
constant gas pressure orion into pdr
title constant gas pressure orion into pdr
* trace molecules
c
c commands controlling continuum =========
c the incident continuum is two parts
c star and flux of photons striking it
c this is the result of this command
c table star kurucz 39600 K
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_kurucz_39600.dat"
Q(H) 49
c add cosmic rays, which are important at depth
cosmic rays, background
c plus hot brems
brems 6
phi(h) 10
c
c commands for density & abundances =========
hden 4
init file="fast.ini"
abundances hii region no grains
grains orion no qheat single
c
c commands controlling geometry =========
radius 17.4507
sphere
c let it go into the molecular cloud
stop temperature 10 linear
c this is roughly the thickness of OMC1
stop thickness 1 linear parsec
c mimic existance of unmodeled molecular gas with double command
double
c
c other commands for details =========
failures 3
iterate
c make constant gas pressure since fast turned off line radiation pressure
constant gas pressure
c the observed microturbulence
turbulence 8 km/sec
c we have a spherical geometry but want to simulate observing
c through a spectrometer's slit. use the aperture
c command for this
aperture beam
c
c commands controlling output =========
normalize to "CaBo" 4861
print line faint -4
punch overview last "orion_hii_pdr_fast.ovr"
punch heating "orion_hii_pdr_fast.het"
punch cooling "orion_hii_pdr_fast.col"
punch dr last "orion_hii_pdr_fast.dr"
punch results last "orion_hii_pdr_fast.rlt"
punch continuum last "orion_hii_pdr_fast.con" units microns
punch ionizing continuum last "orion_hii_pdr_fast.ion"
punch hydrogen 21 cm last "orion_hii_pdr_fast.21cm"
punch hydrogen lya last "orion_hii_pdr_fast.lya"
punch element silicon last "orion_hii_pdr_fast.sil"
punch element sulphur last "orion_hii_pdr_fast.sul"
c
c
c orion_hii_pdr_fast.in
c class hii pdr
c ========================================
c
This model is as close as I can get to the Tielens & Hollenbach (1985) pdr,
except that it starts at the face of the hii region.
orion_hii_pdr_pp.in
the Orion HII Region / PDR / Molecular cloud with an open geometry
title the Orion HII Region / PDR / Molecular cloud with an open geometry
c
c commands controlling continuum =========
c the incident continuum has two parts
c kurucz continuum with flux of photons striking cloud
c this is the the OVI star, its temperature and phi(H)
c this is the result of this command
c table star kurucz 39600 K
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_kurucz_39600.dat"
phi(H) 13
c
c this adds the observed hot brems
c its temperature (as log of T) and the flux of
c photons striking the cloud
brems 6
phi(h) 10
c
c cosmic rays are important for pdr chemistry
cosmic rays, background
c
c commands controlling geometry =========
c this turns off the stop temperature option
c so the sim will not stop due to temperature
stop temperature off
c this sets the thickness of the HII region & PDR
stop thickness 0.5 linear parsec
c this is typical of a gmc, larson 1981
stop total hydrogen column density 22
c this will result in a milli gauss B-field in molecular region
magnetic field -5 gauss
c assume constant pressure
constant pressure
set nend 2000
c
c other commands for details =========
failures 3
c mimic existance of unmodeled molecular gas
double
c iterate since lines optically thick
iterate
c set microturbulence in equipartition with B field
turbulence equipartition
c set the line width so lines appear on the punch continuum
set punchLwidth 10 km/s
c
c commands for density & abundances =========
c this is the log of the initial H density, cm-3
hden 4
c this will speed up the calculation a bit
init file="ism.ini"
c this uses HII region abundances, but no grains
abundances hii region no grains
c this uses orion grains
grains orion
c turn on PAHs, with an abundance that depends on H0 fraction,
c as suggested by long-slit observations of Orion bar,
c with an abundance 3x larger than default built into the code
grains pah function 3
c
c commands controlling output =========
c print lots of faint CO lines
print line faint -6
c normalize to Ha
normalize to "H 1" 6563
punch overview last "orion_hii_pdr_pp.ovr"
punch temperature last "orion_hii_pdr_pp.tem"
punch overview "orion_hii_pdr_pp.ovr1"
punch molecules last "orion_hii_pdr_pp.mol"
punch molecules "orion_hii_pdr_pp.mol1"
punch heating "orion_hii_pdr_pp.het"
punch cooling "orion_hii_pdr_pp.col"
punch dr last "orion_hii_pdr_pp.dr"
punch results last "orion_hii_pdr_pp.rlt"
punch continuum last "orion_hii_pdr_pp.con" units microns
punch hydrogen 21 cm last "orion_hii_pdr_pp.21cm"
punch hydrogen lya last "orion_hii_pdr_pp.lya"
punch element silicon last "orion_hii_pdr_pp.sil"
punch pressure last "orion_hii_pdr_pp.pre"
punch pressure "orion_hii_pdr_pp.pre1"
punch grain extinction last "orion_hii_pdr_pp.grnext"
punch grain charge last "orion_hii_pdr_pp.grnchr"
punch grain potential last "orion_hii_pdr_pp.grnpot"
punch grain temperature last "orion_hii_pdr_pp.grntem"
punch grain temperature "orion_hii_pdr_pp.grntem1"
punch grain charge last "orion_hii_pdr_pp.grnchr"
c
c
c orion_hii_pdr_pp.in
c class hii pdr
c ========================================
orion_hii_single_grn.in
conditions similar to Orion nebula blister
title conditions similar to Orion nebula blister
c
c orion_hii_dist_grn and orion_hii_single_grn form a pair
c which are identical expect for resolved grains
c
c commands controlling continuum =========
c this is the result of this command
c table star kurucz 39600 K
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files
table read "star_kurucz_39600.dat"
phi(h) 13.0
c
c this is to include velocity field in radiative transfer
c but don't include in pressure
turbulence 8 km/sec no pressure
c
c commands for density & abundances =========
hden 4
abundances hii region no grains
c fastest treatment of grain physics, single size bin
c for each species, no quantum heating
grains orion no qheat single
c
c commands controlling geometry =========
constant pressure
sphere
c this is the account for continued matter beyond i front
double optical depths
c
c other commands for details =========
cosmic rays, background
iterate 3
c
c commands controlling output =========
c print lines as surface brightness
print lines surface brightness arcsec
print diffuse continuum
c want to print line intensities are surface brightness, per arcsec^2
print lines surface brightness arcsec
print diffuse continuum
print line faint -1.5
punch overview last "orion_hii_single_grn.ovr"
punch hydrogen 21 cm last "orion_hii_single_grn.21cm"
punch dr last "orion_hii_single_grn.dr"
punch continuum last "orion_hii_single_grn.con" units microns
c
c orion_hii_single_grn.in
c class hii
c ========================================
c
This is a model similar in spirit to the blister geometry HII region
model computed by Baldwin et al. (1991). Orion grains are turned on
with the abundances command. The constant pressure command does a
hydrostatic equilibrium structure. The predicted emission line spectrum is
affected by the reddening of the internal grains. The resulting t2
analysis produces artificial results as a result.
Checks:
- Pressure convergence
- Helium neutral fraction
pdr_co_fully.in
test case where H2 and CO go into fully molecular limit
title test case where H2 and CO go into fully molecular limit
c this is a pair with pdr_co_fully_noneq - that tests non equilibrium chem
c
c commands controlling continuum =========
blackbody 5000
luminosity total solar linear 2
brems 6
luminosity total solar log -2.7
c
c commands for density & abundances =========
c set the density and composition, ism with grains
hden 10
abundances ism no qheat
c
c other commands for details =========
c this is an important ionization source
cosmic ray background
c code will stop at 4000 K by default, lower this temp so we
c can do fully molecular gas
stop temperature 10K linear
c
c commands controlling geometry =========
c the radius and thickness
radius 15.8
stop thickness 6
turbulence 3 km / s
c
c commands controlling output =========
punch overview "pdr_co_fully.ovr"
punch molecules "pdr_co_fully.mol"
punch dr "pdr_co_fully.dr"
punch heating "pdr_co_fully.het"
punch cooling "pdr_co_fully.col"
punch grain charge "pdr_co_fully.grnchr"
punch grain drift "pdr_co_fully.grndft"
punch grain temperature "pdr_co_fully.grntem"
c
c
c gas has more he++ than He+ due to fast charge transfer with co
c neutralizing he+
c pdr_co_fully.in
c class pdr
c ========================================
test code in fully molecular limit
this is a pair with pdr_co_fully_noneq - that tests non equilibrium chem
pdr_co_fully_noneq.in
test case where H2 and CO go into fully molecular limit
title test case where H2 and CO go into fully molecular limit
c this is a pair with pdr_co_fully.in, which does not include non-equil chem
c
c commands controlling continuum =========
blackbody 5000
luminosity total solar linear 2
brems 6
luminosity total solar log -2.7
c
c commands for density & abundances =========
c set the density and composition, ism with grains
hden 10
abundances ism no qheat
c
c other commands for details =========
c this is an important ionization source
cosmic ray background
c code will stop at 4000 K by default, lower this temp so we
c can do fully molecular gas
stop temperature 10K linear
c
c commands controlling geometry =========
c the radius and thickness
radius 15.8
stop thickness 6
turbulence 3 km / s
c this turns on non-equilibrium chemistry as suggested by Federman et al. 1996
set chemistry non equilibrium
c
c commands controlling output =========
punch overview "pdr_co_fully_noneq.ovr"
punch molecules "pdr_co_fully_noneq.mol"
punch dr "pdr_co_fully_noneq.dr"
punch heating "pdr_co_fully_noneq.het"
punch cooling "pdr_co_fully_noneq.col"
punch grain charge "pdr_co_fully_noneq.grnchr"
punch grain drift "pdr_co_fully_noneq.grndft"
punch grain temperature "pdr_co_fully_noneq.grntem"
c
c
c gas has more he++ than He+ due to fast charge transfer with co
c neutralizing he+
c pdr_co_fully_noneq.in
c class pdr
c ========================================
test code in fully molecular limit with Federman non-equilibrium chem
this is a pair with pdr_co_fully.in, which does not include non-equil chem
pdr_dense_persei.in
dense phase model of zeta persei cloud, Table 3
title dense phase model of zeta persei cloud, Table 3
c
c density and abundances ==========
c density of dense phase from Table 2 of Le Petit paper ========
hden 4.3
c
c abundandces from Table 1 of Le Petit paper ==========
element carbon abundance 0.000132 linear
element helium abundance 0.10 linear
element oxygen abundance 0.00032 linear
element nitrogen abundance 0.000075 linear
element sulphur abundance 0.0000186 linear
element silicon abundance 0.000029 linear
c set abundance of all other elements to zero ====================
element copper off
element magnesium off
element manganese off
element sodium off
element chlorine off
element vanadium off
element potassium off
element phosphorous off
element calcium off
element iron off
element zinc off
element neon off
element argon off
element fluorine off
element aluminum off
element boron off
element lithium off
element beryllium off
element scandium off
element nickel off
element titanium off
element chromium off
element cobalt off
c do not use Federman rates for this model ==============
set federman chemistry off
c use standard ism grain size distribution ============
grains ism
c
c command controlling the continuum, for this model is Draine 1978 field ====
table draine 0.5 linear
c fix the temperature to 20 K =============
constant temperature 20
c stop at a radius of 4.3e-4 parsecs ============
stop thickness 0.00043 parsecs linear
c Le Petit model does not consider ices, so turn this off ============
no grain molecules
c turn on cosmic rays =========
cosmic rays background
c Set cosmic ray ionization rate to Table 2 of Le Petit paper ========
set csupra -15.6
c Allow calculation to go extend into cold environment ============
stop temperature linear 3
c
c commands controlling output ============
punch overview "pdr_dense_persei.ovr"
punch dr "pdr_dense_persei.dr"
punch molecules "pdr_dense_persei.mol"
punch heating "pdr_dense_persei.het"
punch cooling "pdr_dense_persei.col"
c
c pdr_dense_persei.in
c class pdr
c ====================================
This is the dense phase model presented by LePetit, Roueff, and Herbst
in order to reproduce C2 and C3 column densities observed along the line
of sight to zera persei. This is our attempt at reproducing their calculation.
This is the dense molecular phase, not the phase that produces H3+
// >>refer model pdr Le Petit, F., Roueff, E., & Herbst, E. 2004,
// >>refercon A&A, 417, 993
pdr_HTT91.in
Hollenbach et al. 1991 low-density PDR
title Hollenbach et al. 1991 low-density PDR
c >>refer test model Hollenbach, D., Takahashi, T., &
c >>refcon Tielens, A.G.G.M., 1991, ApJ,377, 192-209
c
c commands controlling continuum =========
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
c this is hot star continuum
black 30000
intensity 0.2 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
cosmic rays, background
c
c commands for density & abundances =========
hden 3
init file="ism.ini"
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl-7 ar-8 fe -6.60
grains orion, abundance log 0.16 no qheat
c
c commands controlling geometry =========
sphere
c stop when gas is fully neutral
stop efrac -10
c this is to stop at an intended place, since results would be
c very dependent on details if we stop on temperature
stop thickness 18.954
c stop when gas is cold
stop temperature 10 linear
c add this to mimic unmodelled neutral gas
double optical depths
c
c other commands for details =========
iterate 2
case b
turbulence 1.5 km/sec
c
c commands controlling output =========
normalize to 157.6m "C 2"
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print ages
punch overview last "pdr_HTT91.ovr"
punch dr "pdr_HTT91.dr "
punch molecules last "pdr_HTT91.mol"
punch results last "pdr_HTT91.rlt"
punch continuum units microns last "pdr_HTT91.con"
punch heating last "pdr_HTT91.het"
punch cooling last "pdr_HTT91.col"
c
c HTT91.in
c class pdr
c ========================================
c
This is the Hollenbach et al 1991 Low-density PDR
The case b command appears because the Lyman lines are
vastly optically thick in this environment.
If the command is removed the Lyman lines will be optically thin,
which will result in fast fluorescent excitation of the
(nearly totally neutral) hydrogen atom.
there is very little CO in this model since it is not thick enough
for the UV pumping lines to become optically thick
pdr_leiden_f1.in
model 1 as defined in e-mail
title model 1 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
c hydrogen density
hden 3.
grains ism 1.16 no qheat
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch overview "pdr_leiden_f1.ovr"
punch continuum "pdr_leiden_f1.con"
punch leiden lines "pdr_leiden_f1.lin"
punch leiden "pdr_leiden_f1.lei"
punch dr "pdr_leiden_f1.dr"
punch molecules "pdr_leiden_f1.mol"
punch grain physical conditions "pdr_leiden_f1.grn"
punch chemistry rates CO "pdr_leiden_f1.rat"
c
c
c pdr_leiden_f1.in
c class pdr
c ========================================
c
This sim has some interesting properties. The grain temp is so low that
O freezes onto grains as H2O. This is so efficient that CO does not become
well formed - the CO/C ratio only reaches 0.3.
pdr_leiden_f2.in
model 2 as defined in e-mail
title model 2 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use
c half the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c hydrogen density
hden 3.
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch overview "pdr_leiden_f2.ovr"
punch continuum "pdr_leiden_f2.con"
punch molecules "pdr_leiden_f2.mol"
punch leiden lines "pdr_leiden_f2.lin"
punch leiden "pdr_leiden_f2.lei"
punch dr "pdr_leiden_f2.dr"
punch grain physical conditions "pdr_leiden_f2.grn"
c
c
c pdr_leiden_f2.in
c class pdr
c ========================================
c
pdr_leiden_f3.in
* print arrays only chlorine
set Badnell dielectronic recombination on
model 3 as defined in e-mail
title model 3 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch overview "pdr_leiden_f3.ovr"
punch continuum "pdr_leiden_f3.con"
punch leiden lines "pdr_leiden_f3.lin"
punch leiden "pdr_leiden_f3.lei"
punch dr "pdr_leiden_f3.dr"
punch grain physical conditions "pdr_leiden_f3.grn"
c
c
c
c pdr_leiden_f3.in
c class pdr
c ========================================
c
pdr_leiden_f4.in
model 4 as defined in e-mail
title model 4 as defined in e-mail
c
c commands controlling continuum =========
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c
c commands controlling geometry =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c hydrogen density
hden 5.5
c
c commands controlling output =========
punch overview "pdr_leiden_f4.ovr"
punch continuum "pdr_leiden_f4.con"
punch leiden lines "pdr_leiden_f4.lin"
punch leiden "pdr_leiden_f4.lei"
punch dr "pdr_leiden_f4.dr"
punch grain physical conditions "pdr_leiden_f4.grn"
c
c
c
c pdr_leiden_f4.in
c class pdr
c ========================================
c
pdr_leiden_hack_f1.in
model 1 as defined in e-mail
title model 1 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
c hydrogen density
hden 3.
grains ism 1.16 no qheat
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_f1.lin"
punch leiden "pdr_leiden_hack_f1.lei"
punch dr "pdr_leiden_hack_f1.dr"
punch grain physical conditions "pdr_leiden_hack_f1.grn"
punch overview "pdr_leiden_hack_f1.ovr"
c
c
c
c
c pdr_leiden_f1_hack.in
c class pdr
c ========================================
c
pdr_leiden_hack_f2.in
model 2 as defined in e-mail
title model 2 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use
c half the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c hydrogen density
hden 3.
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_f2.lin"
punch leiden "pdr_leiden_hack_f2.lei"
punch dr "pdr_leiden_hack_f2.dr"
punch grain physical conditions "pdr_leiden_hack_f2.grn"
punch overview "pdr_leiden_hack_f2.ovr"
c
c
c pdr_leiden_hack_f2.in
c class pdr
c ========================================
c
pdr_leiden_hack_f3.in
model 3 as defined in e-mail
title model 3 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_f3.lin"
punch leiden overview "pdr_leiden_hack_f3.ovr"
punch leiden "pdr_leiden_hack_f3.lei"
punch dr "pdr_leiden_hack_f3.dr"
punch grain physical conditions "pdr_leiden_hack_f3.grn"
c
c
c pdr_leiden_hack_f3.in
c class pdr
c ========================================
c
pdr_leiden_hack_f4.in
model 4 as defined in e-mail
title model 4 as defined in e-mail
c
c commands controlling continuum =========
c
c commands for density & abundances =========
grains ism 1.16 no qheat
c
c commands controlling geometry =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c This command defines the grain temperature to be a constant 20 Kelvin
constant grain temperature 20
c This sets the temperature to a constant 50 Kelvin
constant temperature 50 linear
c hydrogen density
hden 5.5
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_f4.lin"
punch leiden "pdr_leiden_hack_f4.lei"
punch dr "pdr_leiden_hack_f4.dr"
punch grain physical conditions "pdr_leiden_hack_f4.grn"
punch overview "pdr_leiden_hack_f4.ovr"
c
c
c pdr_leiden_hack_f4.in
c class pdr
c ========================================
c
pdr_leiden_hack_v1.in
model 5 as defined in e-mail
title model 5 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c is half the requested value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
c hydrogen density
hden 3.
grains ism 1.16 no qheat
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
* element nitrogen on
* element nitrogen scale -1
* element silicon on
* element silicon scale -1
* element sulphur on
* element sulphur scale -1
* element chlorine on
* element chlorine scale -1
c Define the abundance of PAHs so we can reproduce the same photoelectric
c heating curve as BT94
c We have also scaled the grain abundance to achieve the desired AV/NH ratio
* grains PAH no qheat 8
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_v1.lin"
punch leiden "pdr_leiden_hack_v1.lei"
punch dr "pdr_leiden_hack_v1.dr"
punch grain temperature "pdr_leiden_hack_v1.grn"
punch overview "pdr_leiden_hack_v1.ovr"
c
c
c pdr_leiden_hack_v1.in
c class pdr
c ========================================
c
pdr_leiden_hack_v2.in
model 6 as defined in e-mail
title model 6 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c commands for density & abundances =========
c Define the abundance of PAHs so we can reproduce the same photoelectric
c heating curve as BT94
c We have also scaled the grain abundance to achieve the desired AV/NH ratio
* grains PAH no qheat 40
grains ism 1.16 no qheat
c hydrogen density
hden 3.
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_v2.lin"
punch leiden "pdr_leiden_hack_v2.lei"
punch dr "pdr_leiden_hack_v2.dr"
punch grain temperature "pdr_leiden_hack_v2.grn"
punch overview "pdr_leiden_hack_v2.ovr"
punch heating "pdr_leiden_hack_v2.het"
punch cooling "pdr_leiden_hack_v2.col"
punch h2 rates "pdr_leiden_hack_v2.h2rat"
c
c
c pdr_leiden_hack_v2.in
c class pdr
c ========================================
c
pdr_leiden_hack_v3.in
model 7 as defined in e-mail
title model 7 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
c Define the abundance of PAHs so we can reproduce the same photoelectric
c heating curve as BT94
c We have also scaled the grain abundance to achieve the desired AV/NH ratio
* grains PAH no qheat 5.7
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
c this sim has more than one thermal solution - force into lowest one,
c the one found by the rest of the group - without this we will find a
c soln at 1400 K
force temperature to 60K
failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_v3.lin"
punch leiden "pdr_leiden_hack_v3.lei"
punch dr "pdr_leiden_hack_v3.dr"
punch grain temperature "pdr_leiden_hack_v3.grn"
punch overview "pdr_leiden_hack_v3.ovr"
punch heating "pdr_leiden_hack_v3.het"
punch cooling "pdr_leiden_hack_v3.col"
c
c
c pdr_leiden_hack_v3.in
c class pdr
c ========================================
c
pdr_leiden_hack_v4.in
model 8 as defined in e-mail
title model 8 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c commands for density & abundances =========
c Define the abundance of PAHs so we can reproduce the same photoelectric
c heating curve as BT94
c We have also scaled the grain abundance to achieve the desired AV/NH ratio
* grains PAH no qheat 40
* grains ism 0.936 no qheat
grains ism 1.16 no qheat
set PAH Bakes
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
* failures 3
c use leiden initialization file
init file="pdr_leiden_hack.ini"
c
c commands controlling output =========
punch leiden lines "pdr_leiden_hack_v4.lin"
punch leiden "pdr_leiden_hack_v4.lei"
punch dr "pdr_leiden_hack_v4.dr"
punch grain temperature "pdr_leiden_hack_v4.grn"
punch overview "pdr_leiden_hack_v4.ovr"
punch heating "pdr_leiden_hack_v4.het"
punch cooling "pdr_leiden_hack_v4.col"
punch transmitted continuum "pdr_leiden_hack_v4.con"
c
c
c pdr_leiden_hack_v4.in
c class pdr
c ========================================
c
pdr_leiden_v1.in
model 5 as defined in e-mail
title model 5 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c is half the requested value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
c hydrogen density
hden 3.
grains ism 1.16 no qheat
c add PAHs
grains PAH no qheat 3 function
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c
c commands controlling output =========
punch overview "pdr_leiden_v1.ovr"
punch leiden lines "pdr_leiden_v1.lin"
punch leiden "pdr_leiden_v1.lei"
punch dr "pdr_leiden_v1.dr"
punch grain temperature "pdr_leiden_v1.grn"
c
c
c
c pdr_leiden_v1.in
c class pdr
c ========================================
c
pdr_leiden_v2.in
model 6 as defined in e-mail
title model 6 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c commands for density & abundances =========
c add PAHs and grains
grains PAH no qheat 3 function
grains ism 1.16 no qheat
c hydrogen density
hden 3.
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c
c commands controlling output =========
punch overview "pdr_leiden_v2.ovr"
punch leiden lines "pdr_leiden_v2.lin"
punch leiden "pdr_leiden_v2.lei"
punch dr "pdr_leiden_v2.dr"
punch grain temperature "pdr_leiden_v2.grn"
c
c
c pdr_leiden_v2.in
c class pdr
c ========================================
c
pdr_leiden_v3.in
model 7 as defined in e-mail
title model 7 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 5
extinguish 24
c
c commands for density & abundances =========
c add PAHs and grains
grains PAH no qheat 3 function
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c
c commands controlling output =========
punch overview "pdr_leiden_v3.ovr"
punch leiden lines "pdr_leiden_v3.lin"
punch leiden "pdr_leiden_v3.lei"
punch dr "pdr_leiden_v3.dr"
punch grain temperature "pdr_leiden_v3.grn"
punch heating "pdr_leiden_v3.het"
punch cooling "pdr_leiden_v3.col"
c
c
c pdr_leiden_v3.in
c class pdr
c ========================================
c
pdr_leiden_v4.in
model 8 as defined in e-mail
title model 8 as defined in e-mail
c
c commands controlling continuum =========
c Use the Draine 1978 field, for a semi-infinite slab we have to use half
c the requested value, so the actual value
c they want for the model is actually twice the value below
table draine 50000
extinguish 24
c
c commands for density & abundances =========
grains PAH no qheat 3 function
grains ism 1.16 no qheat
c hydrogen density
hden 5.5
c
c commands controlling geometry =========
c
c other commands for details =========
failures 3
c use leiden initialization file
init file="pdr_leiden.ini"
c
c commands controlling output =========
punch overview "pdr_leiden_v4.ovr"
punch leiden lines "pdr_leiden_v4.lin"
punch leiden "pdr_leiden_v4.lei"
punch dr "pdr_leiden_v4.dr"
punch grain temperature "pdr_leiden_v4.grn"
punch transmitted continuum "pdr_leiden_v4.con"
c
c
c pdr_leiden_v4.in
c class pdr
c ========================================
c
pdr_orion_veil.in
model like Orion's veil
title model like Orion's veil
c
c commands controlling continuum =========
q(h) 50.000000
table star rubin
c AGN & CMB background
background z=0
cosmic rays background
c as a test this was tried - atomic temperature raised fm 100 to 140
c but nothing really wrong
* cosmic rays equipartition
c the observed magnetic field, roughly 100 micro Gauss
magnetic field -4
c
c commands for density & abundances =========
hden 3.1
init "ism.ini"
abundances orion no grains
grains orion no qheat
c
c commands controlling geometry =========
stop neutral hydrogen column density 21.6
radius 19.1
stop temperature linear 5
c this will not be hit, but put it in to insure that we stop if
c gas becomes totally molecular
stop av 10
c
c other commands for details =========
failures 3
* uncomment to include large H2 molecule - results will CHANGE!
* atom h2
turbulence 0.5 km/sec
sphere
c c this simulates Lyman lines in the stellar continuum
c atom h-like lyman pumping off
c
c commands controlling output =========
print line optical depths
punch PDR "pdr_orion_veil.pdr"
punch feii column densities "pdr_orion_veil.fe2col"
punch hydrogen 21 cm "pdr_orion_veil.21cm"
punch molecules "pdr_orion_veil.mol"
punch overveiw "pdr_orion_veil.ovr"
punch column densities "pdr_orion_veil.col"
punch heating "pdr_orion_veil.het"
punch cooling "pdr_orion_veil.col"
punch dr "pdr_orion_veil.dr"
punch H2 rates "pdr_orion_veil.h2rat"
punch pressure "pdr_orion_veil.pre"
punch wind "pdr_orion_veil.wnd"
print line optical depths
c
c pdr_orion_veil.in
c class pdr
c ========================================
c
This is Nick Abel's model of Orion's veil, the layer of gas in front
of the Orion Nebula.
pdr_th85ism.in
Tielens and Hollenbach pdr model with ism grains, Table 2, paper b
title Tielens and Hollenbach pdr model with ism grains, Table 2, paper b
c
c commands controlling continuum =========
c cosmic radiation background
background z=0
c cosmic ray background ionization and heating
cosmic rays, background
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
blackbody, t = 75 K
intensity 2.7 (total)
c this is the attenuated hot stellar continuum
black 30000
intensity 2.204 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c
c commands controlling geometry =========
c simulate effects of gas we do not model
double
c their turbulence
turbulence 2.7 km/sec
c stopping criteria
c this should be the one actually used
stop AV 33.2 extended
c stop when gas is fully neutral
stop efrac -10
c stop when gas is cold
stop temperature 10 linear
c stop at thickness of 16.5 so that mole limit does not
c stop this calculation - that would make results very
c detail dependent
stop thickness 19
c
c commands for density & abundances =========
hden 5.362
c this turns off some elements, and processes we don't need
c so needs to come before the elements
init file="ism.ini"
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl-7 ar-8 fe -6.60
c use ism abundances but turn of quantum heating
grains ism, abundance log 0.16 no qheat
c
c other commands for details =========
c must iterate since many FIR MIR lines are optically thick
iterate
c this is a closed geometry, in Orion, veil covers 2pi sr
sphere
c make Lyman lines very optically thick, which stops pumping by
c Balmer continuum
case b
c this should run cleanly - turn down number of allow conv fails
failures 3
c
c commands controlling output =========
normalize to "C 2" 157.6m
c want the spectrum to be relative to this [C II] line
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print ages
punch overview last "pdr_th85ism.ovr"
punch pdr last "pdr_th85ism.pdr"
punch molecules last "pdr_th85ism.mol"
punch dr last "pdr_th85ism.dr"
punch hydrogen physical conditions last "pdr_th85ism.hyd"
punch results last "pdr_th85ism.rlt"
punch continuum units microns last "pdr_th85ism.con"
punch heating last "pdr_th85ism.het"
punch cooling last "pdr_th85ism.col"
punch hydrogen 21 cm last "pdr_th85ism.21cm"
punch element magnesium last "pdr_th85ism.mag"
punch grain extinction last "pdr_th85ism.grnext"
punch grain potential last "pdr_th85ism.grnpot"
punch grain temperature last "pdr_th85ism.grntem"
punch grain charge last "pdr_th85ism.grnchr"
c
table lines "LineList_PDR.dat"
c pdr_th85ism.in
c class pdr
c ========================================
c
This is the Tielens and Hollenbach (1985a, b) standard model of
the Orion photodissociation region (PDR). Specifically,
this is my attempt at their Table 2 of paper 2, ApJ 291, p749.
The case b command appears because the Lyman lines are
vastly optically thick in this environment.
If the command is removed the Lyman lines will be optically thin,
which will result in fast fluorescent excitation of the
(nearly totally neutral) hydrogen atom.
// >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985a, ApJ, 291, 722
// >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985b, ApJ, 291, 746
pdr_th85ism_cgto.in
Tielens and Hollenbach pdr with ism grains and C > O
title Tielens and Hollenbach pdr with ism grains and C > O
c
c commands controlling continuum =========
c cosmic radiation background
background z=0
c cosmic ray background ionization and heating
cosmic rays, background
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
blackbody, t = 75 K
intensity 2.7 (total)
c this is the attenuated hot stellar continuum
black 30000
intensity 2.204 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c
c commands controlling geometry =========
c simulate effects of gas we do not model
double
c their turbulence
turbulence 2.7 km/sec
c stopping criteria
c this should be the one actually used
stop AV 33.2 extended
c stop when gas is fully neutral
stop efrac -10
c stop when gas is cold
stop temperature 10 linear
c stop at thickness of 16.5 so that mole limit does not
c stop this calculation - that would make results very
c detail dependent
stop thickness 19
c
c commands for density & abundances =========
hden 5.362
c this turns off some elements, and processes we don't need
c so needs to come before the elements
init file="ism.ini"
c this is the TH85 mix except C is raised by 0.2 dex,
c and O lowered byt his amount, so that C/O > 1
abundances he -1.01 c -3.4 n-8 o-3.50 ne-8 mg-5.89
continue si -6.10 s -5.10 cl-7 ar-8 fe -6.60
c use ism abundances but turn of quantum heating
grains ism, abundance log 0.16 no qheat
c
c other commands for details =========
c must iterate since many FIR MIR lines are optically thick
iterate
c this is a closed geometry, in Orion, veil covers 2pi sr
sphere
c make Lyman lines very optically thick, which stops pumping by
c Balmer continuum
case b
c this should run cleanly - turn down number of allow conv fails
failures 3
c
c commands controlling output =========
normalize to "C 2" 157.6m
c want the spectrum to be relative to this [C II] line
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print ages
punch overview last "pdr_th85ism_cgto.ovr"
punch pdr last "pdr_th85ism_cgto.pdr"
punch molecules last "pdr_th85ism_cgto.mol"
punch dr last "pdr_th85ism_cgto.dr"
punch hydrogen physical conditions last "pdr_th85ism_cgto.hyd"
punch results last "pdr_th85ism_cgto.rlt"
punch continuum units microns last "pdr_th85ism_cgto.con"
punch heating last "pdr_th85ism_cgto.het"
punch cooling last "pdr_th85ism_cgto.col"
punch hydrogen 21 cm last "pdr_th85ism_cgto.21cm"
punch element magnesium last "pdr_th85ism_cgto.mag"
punch grain extinction last "pdr_th85ism_cgto.grnext"
punch grain potential last "pdr_th85ism_cgto.grnpot"
punch grain temperature last "pdr_th85ism_cgto.grntem"
punch grain charge last "pdr_th85ism_cgto.grnchr"
c
c pdr_th85ism_cgto.in
c class pdr
c ========================================
c
This is the Tielens and Hollenbach (1985a, b) standard model of
the Orion photodissociation region (PDR) but with the C abundance
increased by 0.2 dex, and O lowered by this amount, so that C/O>1.
Specifically, this is their Table 2 of paper 2, ApJ 291, p749.
The case b command appears because the Lyman lines are
vastly optically thick in this environment.
If the command is removed the Lyman lines will be optically thin,
which will result in fast fluorescent excitation of the
(nearly totally neutral) hydrogen atom.
// >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985a, ApJ, 291, 722
// >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985b, ApJ, 291, 746
pdr_th85orion.in
Tielens and Hollenbach pdr model with orion grains, Table 2, paper b
title Tielens and Hollenbach pdr model with orion grains, Table 2, paper b
c
c commands controlling continuum =========
c cosmic background
background
c galactic cosmic ray background
cosmic rays, background
c first continuum is FIR hot grain continuum produced in
c unmodeled HII Region
blackbody, t = 75 K
intensity 2.7 (total)
c this is the attenuated hot stellar continuum
black 30000
intensity 2.204 range 0.4412 to 1 Ryd
c this will remove all ionizing radiation
extinguish 24 0
c
c commands for density & abundances =========
hden 5.362
init file="ism.ini"
grains orion, abundance log 0.16 no qheat
abundances he -1.01 c -3.52 n-8 o-3.30 ne-8 mg-5.89
continue si -6.10 s -5.10 cl-7 ar-8 fe -6.60
c
c commands controlling geometry =========
c simulate effects of gas we do not model
double
sphere
c stopping criteria
c this should be the one actually used
c >>chn 04 mar 13, from 10 to 36.3, bug in AV had stopped at this depth,
c so change it to this to keep model the same
stop AV 36.3 extended
c stop when gas is fully neutral
stop efrac -10
c stop when gas is cold
stop temperature 10 linear
c stop at thickness so that would make results very
c detail dependent
stop thickness 19
c
c other commands for details =========
turbulence 2.7 km/sec
iterate
failures 3
c
c this is done to not allow pumping and subsequent photoionization
c of H by Balmer continuum in keeping with std pdr assumptions
case b
c
c commands controlling output =========
normalize to "C 2" 157.6m
c uv lines are strongly pumped by stellar continuum, break out contribution
print line pump
print line optical depths
print ages
c
punch overview "pdr_th85orion.ovr" no hash
punch pdr last "pdr_th85orion.pdr"
c this only outputs major molecules
punch molecules last "pdr_th85orion.mol"
punch dr "pdr_th85orion.dr"
c this is to get all of H molecules
punch hydrogen physical conditions last "pdr_th85orion.hyd"
punch results last "pdr_th85orion.rlt"
punch continuum units microns last "pdr_th85orion.con"
punch ionizing continuum last "pdr_th85orion.ion"
punch heating last "pdr_th85orion.het"
punch cooling last "pdr_th85orion.col"
punch hydrogen 21 cm last "pdr_th85orion.21cm"
c these are for the grains
punch grain extinction last "pdr_th85orion.grnext"
punch grain potential last "pdr_th85orion.grnpot"
punch grain temperature last "pdr_th85orion.grntem"
punch grain charge last "pdr_th85orion.grnchr"
c
c
c pdr_th85orion.in
c class pdr
c ========================================
This is the Tielens and Hollenbach (1985a, b) standard model of
the Orion photodissociation region (PDR). Specifically,
this is my attempt at their Table 2 of paper 2, ApJ 291, p749.
The case b command appears because the Lyman lines are
vastly optically thick in this environment.
If the command is removed the Lyman lines will be optically thin,
which will result in fast fluorescent excitation of the
(nearly totally neutral) hydrogen atom.
// >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985a, ApJ, 291, 722
// >>refer model pdr Tielens, A. G. G. M., & Hollenbach, D. 1985b, ApJ, 291, 746
pn_fluc.in
Paris meeting Planetary nebula with density fluctuations
title Paris meeting Planetary nebula with density fluctuations
c
c commands for density & abundances =========
fluctuations 16.5 4 3
c this sets up the code like version c84, with fewer chemical elements
init file="c84.ini"
c this sets the abundances of the elements that we will use
abund -1 C-3.523 N-4. O-3.222 ne-3.824 na=-10 mg-4.523 al=-10
continue si-4.523 s-4.824 ar-10 ca=-10 fe-10 ni=-10
c
c commands controlling continuum =========
black body, T = 150000 K, radius = 10
black 5
luminosity total 38
c
c commands controlling geometry =========
sphere
radius 17
c
c other commands for details =========
c
c commands controlling output =========
punch overview last "pn_fluc.ovr"
punch element nitrogen last "pn_fluc.nit"
punch dr last "pn_fluc.dr"
punch results last "pn_fluc.rlt"
c
c pn_fluc.in
c class pn function
c ========================================
The boundary conditions are similar to those for the Paris meeting PN,
a homogeneous grain-free PN. This model has density fluctuations,
with values chosen so that the mean ne np V are the same
for the nebulae with and without fluctuations.
Checks:
- Check that the fluctuations command works.
- How do results compare with homogeneous Paris pn?
pn_ots.in
Paris meeting Planetary nebula with ots
title Paris meeting Planetary nebula with ots
c recompute "standard" PN model of the Pequignot Meudon Conferance
c
c commands controlling continuum =========
black body, T = 150000 K, radius = 10
c
c commands for density & abundances =========
hden = 3.4771213
init file="c84.ini"
abund -1 C-3.523 N-4. O-3.222 ne-3.824 na=-10 mg-4.523 al=-10
continue si-4.523 s-4.824 ar-10 ca=-10 fe-10 ni=-10
c
c commands controlling geometry =========
radius = 17
sphere
c
c other commands for details =========
diffuse ots
no level2
iterate
c
c commands controlling output =========
punch overview last "pn_ots.ovr"
punch element nitrogen last "pn_ots.nit"
punch results last "pn_ots.rlt"
punch dr last "pn_ots.dr"
c
c pn_ots.in
c class pn
c ========================================
c
This is the on-the-spot version of the Paris Planetary Nebula.
Checks:
- Q(H) total 4861 luminosity is close to expected value
- Line spectrum similar to that predicted by default conditions.
pn_paris.in
pn_paris.in Meudon Planetary nebula
title pn_paris.in Meudon Planetary nebula
c standard" PN model of the Pequignot Meudon Conferance
c
c commands controlling continuum =========
black body, T = 150000 K, radius = 10
c
c commands for density & abundances =========
hden = 3.4771213
init file="ism.ini"
abund he -1 C-3.523 N-4. O-3.222 ne-3.824 mg-4.523
continue si-4.523 s-4.824 cl-7 ar-10 fe-10
c
c commands controlling geometry =========
radius = 17
sphere
c
c other commands for details =========
age 10000 years
c
c commands controlling output =========
normalize to "Ca b" 4861
print ages
print column densities
punch overview "pn_paris.ovr"
punch transmitted continuum "pn_paris.trn"
punch continuum units microns "pn_paris.con"
punch element carbon "pn_paris.car"
punch physical conditions "pn_paris.phy"
punch element nitrogen "pn_paris.nit"
punch dr "pn_paris.dr"
punch temperature "pn_paris.tem"
punch lines emissivity "pn_paris.ems"
h 1 6563
s 2 6720
end of lines
punch averages "pn_paris.avr" last
temperature hydrogen 1
end of averagers
c
c pn_paris.in
c class pn
c ========================================
c
This is one of the "standard" models computed at the Paris meeting on
photoionization and shock calculations. A table in Hazy compares the
predictions of the current version of CLOUDY with predictions of a
few of the other codes.
The model is meant to represent a planetary nebula ionized by a very
hot central object. As a result there is a great deal of He++ and
the associated line and continuum emission is very important.
An age of 10000 years was set. This tests logic in the age command.
pn_paris_fast.in
pn_paris_fast.in Meudon Planetary nebula
title pn_paris_fast.in Meudon Planetary nebula
c standard" PN model of the Pequignot Meudon Conferance
c
c commands controlling continuum =========
black body, T = 150000 K, radius = 10
c
c commands for density & abundances =========
hden = 3.4771213
init file="fast.ini"
abund he -1 C-3.523 N-4. O-3.222 ne-3.824 mg-4.523
continue si-4.523 s-4.824 ar-10
c
c commands controlling geometry =========
sphere
radius = 17
c
c other commands for details =========
c
c commands controlling output =========
normalize to "Ca b" 4861
punch overview "pn_paris_fast.ovr"
punch physical conditions "pn_paris_fast.phy"
punch element nitrogen "pn_paris_fast.nit"
punch dr "pn_paris_fast.dr"
c
c pn_paris_fast.in
c class pn
c ========================================
c
This is one of the "standard" models computed at the Paris meeting on
photoionization and shock calculations. A table in Hazy compares the
predictions of the current version of CLOUDY with predictions of a
few of the other codes.
The model is meant to represent a planetary nebula ionized by a very
hot central object. As a result there is a great deal of He++ and
the associated line and continuum emission is very important.
pn_sqrden.in
test with density falling as R^-2, and filling factor
title test with density falling as R^-2, and filling factor
c this is an optically thin model, which will have a broad range
c of density but the same ionization parameter across the model.
c the physical condtions should not vary by much across the structure
c
c commands controlling continuum =========
c continuum is hot blackbody at ~Ledd for solar mass
blackbody, t=5 luminosity=38
c
c commands for density & abundances =========
c the density will fall off as the inverse square of the radius, hence the name
hden 5 -2
init file "ism.ini"
c use default PN abundances but turn off quantum heating since
c we are not going to save the predicted near IR continuum
abundances planetary no qheat
c
c commands controlling geometry =========
c small filling factor so that region will be optically thin
filling factor -2
c set inner and outer radius
radius 16 18
sphere
c
c other commands for details =========
iterate
c
c commands controlling output =========
print line faint 1
punch overview "pn_sqrden.ovr" last
punch dr "pn_sqrden.dr" last
c
c pn_sqrden.in
c class pn
c ========================================
Checks:
- Zone thickness budgeting handled OK.
- physical conditions nearly constant across computed structure
stars_atlas.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a hot star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the new Kurucz model option.
c I choose the T_eff and log(g) values arbitrarily. It implies a star
c with luminosity 32300 L_sun, radius 6.4745 R_sun, and mass 3.745 M_sun.
table star atlas 30400.0 4.2
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_atlas.ovr"
punch results last "stars_atlas.rlt"
punch transmitted continuum last "stars_atlas.trn" no header
punch continuum last "stars_atlas.con" units microns
c
c stars_atlas.in
c class stars
c ========================================
This is a test that the code can correctly access the
large block of Atlas model atmosphere continua described by Kurucz (1991).
Kevin Volk provided it as part of his original coding of these
stellar atmosphere files.
Constant temperature is assumed since this tests shape of stellar continuum
rather than thermal physics.
Checks:
table star atlas command works.
stars_atlas_3d.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a hot star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the atlas 3dim model option.
c I choose the T_eff and log(g) values arbitrarily.
table star atlas 3-dim 38400.0 4.82 -1.42
luminosity 4.509202522 solar
table star atlas odfnew 3-dim 38400.0 4.82 0.42
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_atlas_3d.ovr"
punch results last "stars_atlas_3d.rlt"
punch transmitted continuum last "stars_atlas_3d.trn" no header
punch continuum last "stars_atlas_3d.con" units microns
c
c stars_atlas_3d.in
c class stars
c ========================================
This is a test that the code can correctly access the
large 3-dim block of Atlas model atmosphere continua.
Checks:
table star atlas 3dim command works.
stars_atlas_all.in
Test access to various Atlas grids
title Test access to various Atlas grids
c
c commands controlling continuum =========
c Include stars with a wide range of metallicities to test access to these grids.
table star atlas Z+1.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z+0.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z+0.3 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z+0.2 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z+0.1 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z+0.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-0.1 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-0.2 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-0.3 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-0.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-1.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-1.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-2.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-2.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-3.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-3.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-4.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-4.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas Z-5.0 38400.0 4.82
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -2.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_atlas_all.ovr"
punch results last "stars_atlas_all.rlt"
punch transmitted continuum last "stars_atlas_all.trn" no header
punch continuum last "stars_atlas_all.con" units microns
c
c stars_atlas_all.in
c class stars
c ========================================
This is a test that the code can correctly access the
2D grids of various metallicity of Atlas model atmosphere continua.
Checks:
access to the various atlas grids
stars_atlas_odfnew_all.in
Test access to various Atlas ODFNEW grids
title Test access to various Atlas ODFNEW grids
c
c commands controlling continuum =========
c Include stars with a wide range of metallicities to test access to these grids.
table star atlas odfnew Z+0.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas odfnew Z+0.2 38400.0 4.82
luminosity 4.509202522 solar
table star atlas odfnew Z+0.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas odfnew Z-0.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas odfnew Z-1.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas odfnew Z-1.5 38400.0 4.82
luminosity 4.509202522 solar
table star atlas odfnew Z-2.0 38400.0 4.82
luminosity 4.509202522 solar
table star atlas odfnew Z-2.5 38400.0 4.82
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -2.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_atlas_odfnew_all.ovr"
punch results last "stars_atlas_odfnew_all.rlt"
punch transmitted continuum last "stars_atlas_odfnew_all.trn" no header
punch continuum last "stars_atlas_odfnew_all.con" units microns
c
c stars_atlas_odfnew_all.in
c class stars
c ========================================
This is a test that the code can correctly access the
Atlas ODFNEW model atmosphere continua.
Checks:
access to the various atlas odfnew grids
stars_bstar2006.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a B star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the new Tlusty model option.
c I choose the T_eff and log(g) values arbitrarily.
table star tlusty bstar 20400.0 4.2
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
c
c in this model the Lyman lines are pumped significantly by the stellar continuum;
c disable this process to avoid critical dependence on the precise shape of the continuum.
no induced processes
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_bstar2006.ovr"
punch results last "stars_bstar2006.rlt"
punch transmitted continuum last "stars_bstar2006.trn" no header
c
c stars_bstar2006.in
c class stars
c ========================================
This is a test that the code can correctly access the
large block of Tlusty B-star model atmosphere continua.
Checks:
table star tlusty bstar command works.
stars_bstar2006_3d.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a B star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the new Tlusty model option.
c I choose the T_eff, log(g) and log(Z) values arbitrarily.
table star tlusty bstar 3-dim 20400.0 4.2 -0.573
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
c
c in this model the Lyman lines are pumped significantly by the stellar continuum;
c disable this process to avoid critical dependence on the precise shape of the continuum.
no induced processes
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_bstar2006_3d.ovr"
punch results last "stars_bstar2006_3d.rlt"
punch transmitted continuum last "stars_bstar2006_3d.trn" no header
c
c stars_bstar2006_3d.in
c class stars
c ========================================
This is a test that the code can correctly access the
large 3-dim block of Tlusty B-star model atmosphere continua.
Checks:
table star tlusty bstar 3-dim command works.
stars_bstar2006_all.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of 6 B stars ionizing their surroundings.
c Assume H density of 1e+6 per cc. This run tests all the various
c metallicities contained in the BSTAR 2006 grids. I choose the T_eff
c and log(g) values arbitrarily.
table star tlusty bstar Z+0.3 20000.0 4.5
luminosity 4.509202522 solar
table star tlusty bstar Z+0.0 21400.0 4.5
luminosity 4.509202522 solar
table star tlusty bstar Z-0.3 15000.0 4.2
luminosity 4.509202522 solar
table star tlusty bstar Z-0.7 20400.0 4.2
luminosity 4.509202522 solar
table star tlusty bstar Z-1.0 27500.0 3.0
luminosity 4.509202522 solar
table star tlusty bstar Z-INF 18000.0 4.75
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -2.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
c
c in this model the Lyman lines are pumped significantly by the stellar continuum;
c disable this process to avoid critical dependence on the precise shape of the continuum.
no induced processes
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_bstar2006_all.ovr"
punch results last "stars_bstar2006_all.rlt"
punch transmitted continuum last "stars_bstar2006_all.trn" no header
c
c stars_bstar2006_all.in
c class stars
c ========================================
This is a test that the code can correctly access the
large block of Tlusty B-star model atmosphere continua.
Checks:
table star tlusty bstar metallicity option works.
stars_costar1.in
test costar continuum, interpolation mode 1
title test costar continuum, interpolation mode 1
c
c commands controlling continuum =========
table star costar, 40000 K
ionization parameter -2
table star costar, 31000 K, index 4
ionization parameter -2
c
c commands for density & abundances =========
hden 2
grains orion single no qheat
abundances ism no grains
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
* >>chng 06 jul 01, add this to stars set since want to test shape
* of continuum and resulting ion dist, not thermal equilibrium of nebula
constant temperature 6000K
iterate
c
c commands controlling output =========
print last iteration
punch continuum "stars_costar1.con" units microns last
c
c stars_costar1.in
c class stars
c ========================================
c
This is a test of whether the code can read the CoStar stellar atmospheres.
It checks the resulting ionization to confirm that the continuum
shape was ok.
Constant temperature is assumed since this tests shape of stellar continuum
rather than thermal physics.
stars_costar2.in
test costar continuum, interpolation mode 2
title test costar continuum, interpolation mode 2
c
c commands controlling continuum =========
table star costar 35575.4 4.2367
ionization parameter -2
table star costar 34468.5 3.86765
ionization parameter -2
table star costar 41741.2 4.02
ionization parameter -2
table star costar 45000 4.1593
ionization parameter -2
table star costar 27500 3.3
ionization parameter -2
c
c commands for density & abundances =========
hden 2
grains orion single no qheat
abundances ism no grains
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
iterate
c
c commands controlling output =========
print last iteration
punch continuum "stars_costar2.con" units microns last
c
c stars_costar2.in
c class stars
c ========================================
c
This is a test of whether the code can read the CoStar stellar atmospheres.
It checks the resulting ionization to confirm that the continuum
shape was ok.
stars_costar3.in
test costar continuum, interpolation mode 3
title test costar continuum, interpolation mode 3
c
c commands controlling continuum =========
table star costar zams 40.0 Msol 2.85e6
ionization parameter -2
table star costar zams 40.0 Msol 4e6
ionization parameter -2
table star costar zams 50.9 Msol 2.85e6
ionization parameter -2
table star costar zams 42.9 Msol 1.6e6
ionization parameter -2
table star costar zams 70.5 Msol 2.1e6
ionization parameter -2
c
c commands for density & abundances =========
hden 2
grains orion single no qheat
abundances ism no grains
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
iterate
c
c commands controlling output =========
print last iteration
punch continuum "stars_costar3.con" units microns last
c
c stars_costar3.in
c class stars
c ========================================
c
This is a test of whether the code can read the CoStar stellar atmospheres.
It checks the resulting ionization to confirm that the continuum
shape was ok.
stars_costar4.in
test costar continuum, interpolation mode 4
title test costar continuum, interpolation mode 4
c
c commands controlling continuum =========
table star costar age 2.85e6 40.0
ionization parameter -2
table star costar age 4e6 40.0
ionization parameter -2
table star costar age 2.85e6 50.9
ionization parameter -2
table star costar age 1.6e6 42.9
ionization parameter -2
table star costar age 2.1e6 70.5
ionization parameter -2
c
c commands for density & abundances =========
hden 2
grains orion single no qheat
abundances ism no grains
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
iterate
c
c commands controlling output =========
print last iteration
punch continuum "stars_costar4.con" units microns last
c
c stars_costar4.in
c class stars
c ========================================
c
This is a test of whether the code can read the CoStar stellar atmospheres.
It checks the resulting ionization to confirm that the continuum
shape was ok.
stars_costarhalo.in
test costar halo abundances continuum
title test costar halo abundances continuum
c
c commands controlling continuum =========
table star costar, halo abundances, 40000 K
ionization parameter -2
c
c commands for density & abundances =========
hden 2
abundances ism
c
c commands controlling geometry =========
set dr 0
stop zone 1
c
c other commands for details =========
* >>chng 06 jul 01, add this to stars set since want to test shape
* of continuum and resulting ion dist, not thermal equilibrium of nebula
constant temperature 6900K
iterate
c
c commands controlling output =========
print last iteration
punch continuum last "stars_costarhalo.con" units microns
c
c
c stars_costarhalo.in
c class stars
c ========================================
c
This test checks that the code can read the CoStar
stellar atmospheres. It checks the resulting
ionization to make sure that the continuum shape is ok.
Constant temperature is assumed since this tests shape of stellar continuum
rather than thermal physics.
stars_kurucz79.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a hot star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the Kurucz model option.
c I choose the T_eff and log(g) values arbitrarily. It implies a star
c with luminosity 32300 L_sun, radius 6.4745 R_sun, and mass 3.745 M_sun.
table star kurucz 30400.0
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_kurucz79.ovr"
punch results last "stars_kurucz79.rlt"
punch transmitted continuum last "stars_kurucz79.trn" no header
punch continuum last "stars_kurucz79.con" units microns
c
c stars_kurucz79.in
c class stars
c ========================================
This is a test that the code can correctly access the
Kurucz model atmosphere continua described by Kurucz (1979).
Kevin Volk provided it as part of his original coding of these
stellar atmosphere files.
Constant temperature is assumed since this tests shape of stellar continuum
rather than thermal physics.
Checks:
table star kurucz command works.
stars_mihalas.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a hot star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the new Kurucz model option.
c I choose the T_eff and log(g) values arbitrarily. It implies a star
c with luminosity 32300 L_sun, radius 6.4745 R_sun, and mass 3.745 M_sun.
table star mihalas 30400.0
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_mihalas.ovr"
punch results last "stars_mihalas.rlt"
punch transmitted continuum last "stars_mihalas.trn" no header
punch continuum last "stars_mihalas.con" units microns
c
c stars_mihalas.in
c class stars
c ========================================
This is a test that the code can correctly access the
small Mihalas grid of NLTE model atmospeheres
Constant temperature is assumed since this tests shape of continuum
not thermal physics.
Checks:
table star mihalas command works.
stars_optimize1.in
test optimizer limits on stellar atmosphere grids
title test optimizer limits on stellar atmosphere grids
c
c commands controlling continuum =========
table star tlusty ostar 3-dim 31200 3.6 -1.12 vary
ionization par -3
c
c commands for density & abundances =========
hden 4 vary
init file "hheonly.ini"
element oxygen on
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
optimize tolerance 0.03
constant temper 4.3 vary
optimize phymir sequential
optimize iteration 500
c the following spectrum was derived by running this sim at hden 5, t=4,
c and T_star=51200
c >>chng 06 dec 01, change line intensities due to new Badnell DR rates
optimize lines
he 2 303.8 102.3658
o ii 3729 0.0350
o ii 3726 0.0990
o ii 7323 0.1286
o ii 7332 0.1037
o 3 5007 6.0385
totl 4363 0.0570
end of lines
c
c commands controlling output =========
print line faint -2
c
c
c optimize_phymir.in
c class optimizer
c ========================================
c
This checks whether the optimizer sets corrects limits for Teff. The line
spectrum was calculated at T = 1e4 K, n_H=1e5 cm^-3, T_star = 51200 K
and resulted in the given electron density.
stars_optimize2.in
test optimizer limits on stellar atmosphere grids
title test optimizer limits on stellar atmosphere grids
c
c commands controlling continuum =========
table star costar 36200 vary
ionization par -3
c
c commands for density & abundances =========
hden 4 vary
init file "hheonly.ini"
element oxygen on
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
optimize tolerance 0.03
constant temper 4.3 vary
optimize phymir sequential
optimize iteration 500
c the following spectrum was derived by running this sim at hden 5, t=4,
c and T_star=51200
c >>chng 06 dec 01, change line intensities due to new Badnell DR rates
optimize lines
he 2 303.8 132.5159
o ii 3729 0.0249
o ii 3726 0.0705
o ii 7323 0.0918
o ii 7332 0.0740
o 3 5007 5.3344
totl 4363 0.0504
end of lines
c
c commands controlling output =========
print line faint -2
c
c
c optimize_phymir.in
c class optimizer
c ========================================
c
This checks whether the optimizer sets corrects limits for Teff. The line
spectrum was calculated at T = 1e4 K, n_H=1e5 cm^-3, T_star = 51200 K
and resulted in the given electron density.
stars_optimize3.in
test optimizer limits on stellar atmosphere grids
title test optimizer limits on stellar atmosphere grids
c
c commands controlling continuum =========
table star costar age 5.3 log 36 msol vary
ionization par -3
c
c commands for density & abundances =========
hden 4 vary
init file "hheonly.ini"
element oxygen on
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
optimize tolerance 0.03
constant temper 4.3 vary
optimize phymir sequential
optimize iteration 500
c the following spectrum was derived by running this sim at hden 5, t=4,
c and log(age)=6.3
c >>chng 06 dec 01, change line intensities due to new Badnell DR rates
optimize lines
he 2 303.8 28.5845
o ii 3729 0.0252
o ii 3726 0.0713
o ii 7323 0.0928
o ii 7332 0.0748
o 3 5007 2.4914
totl 4363 0.0235
end of lines
c
c commands controlling output =========
print line faint -2
c
c
c optimize_phymir.in
c class optimizer
c ========================================
c
This checks whether the optimizer sets corrects limits for log(age). The line
spectrum was calculated at T = 1e4 K, n_H=1e5 cm^-3, log(age/yr) = 6.3
and resulted in the given electron density.
stars_ostar2002.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a hot star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the new Tlusty model option.
c I choose the T_eff and log(g) values arbitrarily. It implies a star
c with luminosity 32300 L_sun, radius 6.4745 R_sun, and mass 3.745 M_sun.
table star tlusty ostar 30400.0 4.2
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
c
c in this model the Lyman lines are pumped significantly by the stellar continuum;
c disable this process to avoid critical dependence on the precise shape of the continuum.
no induced processes
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_ostar2002.ovr"
punch results last "stars_ostar2002.rlt"
punch transmitted continuum last "stars_ostar2002.trn" no header
c
c stars_ostar2002.in
c class stars
c ========================================
This is a test that the code can correctly access the
large block of Tlusty O-star model atmosphere continua.
Checks:
table star tlusty ostar command works.
stars_ostar2002_3d.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a hot star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the new Tlusty model option.
c I choose the T_eff, log(g) and log(Z) values arbitrarily.
table star tlusty ostar 3-dim 30400.0 4.2 -1.573
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -3.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
c
c in this model the Lyman lines are pumped significantly by the stellar continuum;
c disable this process to avoid critical dependence on the precise shape of the continuum.
no induced processes
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_ostar2002_3d.ovr"
punch results last "stars_ostar2002_3d.rlt"
punch transmitted continuum last "stars_ostar2002_3d.trn" no header
c
c stars_ostar2002_3d.in
c class stars
c ========================================
This is a test that the code can correctly access the
large block of Tlusty O-star model atmosphere continua.
Checks:
table star tlusty ostar 3-dim command works.
stars_ostar2002_all.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of 10 hot stars ionizing their surroundings.
c Assume H density of 1e+6 per cc. This run tests all the various
c metallicities contained in the OSTAR 2002 grids. I choose the T_eff
c and log(g) values such that the selection algorithm gets a good workout.
table star tlusty ostar Z+0.3 30000.0 4.5
luminosity 4.509202522 solar
table star tlusty ostar Z+0.0 31400.0 4.5
luminosity 4.509202522 solar
table star tlusty ostar Z-0.3 45000.0 4.2
luminosity 4.509202522 solar
table star tlusty ostar Z-0.7 30400.0 4.2
luminosity 4.509202522 solar
table star tlusty ostar Z-1.0 27500.0 3.0
luminosity 4.509202522 solar
table star tlusty ostar Z-1.5 55000.0 4.75
luminosity 4.509202522 solar
table star tlusty ostar Z-1.7 27500.0 4.75
luminosity 4.509202522 solar
table star tlusty ostar Z-2.0 55000.0 3.00
luminosity 4.509202522 solar
table star tlusty ostar Z-3.0 44400.0 3.2
luminosity 4.509202522 solar
table star tlusty ostar Z-INF 53000.0 4.75
luminosity 4.509202522 solar
table star tlusty ostar Z+0.0 48400.0 3.87
luminosity 4.509202522 solar
table star tlusty ostar Z-0.7 48400.0 3.55
luminosity 4.509202522 solar
table star tlusty ostar Z+0.3 43400.0 3.65
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -2.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
c
c in this model the Lyman lines are pumped significantly by the stellar continuum;
c disable this process to avoid critical dependence on the precise shape of the continuum.
no induced processes
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_ostar2002_all.ovr"
punch results last "stars_ostar2002_all.rlt"
punch transmitted continuum last "stars_ostar2002_all.trn" no header
c
c stars_ostar2002_all.in
c class stars
c ========================================
This is a test that the code can correctly access the
large block of Tlusty O-star model atmosphere continua.
Checks:
table star tlusty ostar metallicity option works, also test selection algorithm.
stars_rauch_3d.in
Test Rauch 3-dim option
title Test Rauch 3-dim option
c
c commands controlling continuum =========
table star rauch 3-dim T = 150000 K, g = 6.5 , log(Z) = -0.6
luminosity 4.5 solar
table star rauch old 3-dim T = 233000 K, g = 7.5 , log(Z) = -0.32
luminosity 4.5 solar
c
c commands controlling geometry =========
sphere
radius 18
stop zone 1
set dr 0
c
c commands for density & abundances =========
abundances planetary
hden 3.0
c
c other commands for details =========
iterate
no level2
c
c commands controlling output =========
print last iteration
plot continuum range .1
punch overview last "stars_rauch_3d.ovr"
punch results last "stars_rauch_3d.rlt"
punch continuum last "stars_rauch_3d.con" units microns
c
c stars_rauch_3d.in
c class stars
c ========================================
This model checks that the code is able to access the 3-dim versions
of the Rauch grids of stellar atmosphere models.
Constant temperature is set since this test is intended to test
shape of ionizing continuum, not thermal physics.
stars_rauch_h+he.in
hot PN model
title hot PN model
c
c commands controlling continuum =========
c this tests the various interpolation modes of the 3d grid
table star rauch h+he T = 150000 K, g = 6.5 , f(He) = 0.2
luminosity 4.5 solar
table star rauch h+he T = 155000 K, g = 6.5 , f(He) = 0.2
luminosity 4.5 solar
table star rauch h+he T = 150000 K, g = 6.25 , f(He) = 0.2
luminosity 4.5 solar
table star rauch h+he T = 150000 K, g = 6.5 , f(He) = 0.25
luminosity 4.5 solar
table star rauch h+he T = 155000 K, g = 6.25 , f(He) = 0.2
luminosity 4.5 solar
table star rauch h+he T = 155000 K, g = 6.5 , f(He) = 0.25
luminosity 4.5 solar
table star rauch h+he T = 150000 K, g = 6.25 , f(He) = 0.25
luminosity 4.5 solar
table star rauch h+he T = 155000 K, g = 6.25 , f(He) = 0.25
luminosity 4.5 solar
c
c commands controlling geometry =========
sphere
radius 18
stop zone 1
set dr 0
c
c commands for density & abundances =========
abundances planetary
hden 3.0
c
c other commands for details =========
iterate
no level2
c
c commands controlling output =========
print last iteration
plot continuum range .1
punch overview last "stars_rauch_h+he.ovr"
punch results last "stars_rauch_h+he.rlt"
punch continuum last "stars_rauch_h+he.con" units microns
c
c stars_rauch_h+he.in
c class stars
c ========================================
This is a model of a very hot planetary nebula, and checks that the code
is able to access the H+He version of the Rauch grid of stellar
atmosphere models.
stars_rauch_h-ca.in
very hot PN model
title very hot PN model
c NB
c this uses the H-Ca rauch atmospheres
c
c commands controlling continuum =========
table star old rauch T = 250000 K, g = 7.5
luminosity 4.5 solar
c
c commands controlling geometry =========
sphere
radius 18
stop zone 1
set dr 0
c
c other commands for details =========
no level2
c
c commands for density & abundances =========
hden 3.0
abundances planetary
c
c commands controlling output =========
print last iteration
iterate
plot continuum range .1
punch overview last "stars_rauch_h-ca.ovr"
punch results last "stars_rauch_h-ca.rlt"
punch continuum last "stars_rauch_h-ca.con" units microns
c
c stars_rauch_h-ca.in
c class stars
c ========================================
This is a model of a very hot planetary nebula, and checks that the code
is able to access the H-Ca Rauch grid of stellar atmosphere models.
stars_rauch_h-ni.in
hot PN model
title hot PN model
c
c commands controlling continuum =========
table star rauch T = 150000 K, g = 6.5
luminosity 4.5 solar
c
c commands controlling geometry =========
sphere
radius 18
stop zone 1
set dr 0
c
c commands for density & abundances =========
abundances planetary
hden 3.0
c
c other commands for details =========
iterate
no level2
c
c commands controlling output =========
print last iteration
plot continuum range .1
punch overview last "stars_rauch_h-ni.ovr"
punch results last "stars_rauch_h-ni.rlt"
punch continuum last "stars_rauch_h-ni.con" units microns
c
c stars_rauch_h-ni.in
c class stars
c ========================================
This is a model of a very hot planetary nebula, and checks that the code
is able to access the H-Ni version of the Rauch grid of stellar
atmosphere models.
Constant temperature is set since this test is intended to test
shape of ionizing continuum, not thermal physics.
stars_rauch_helium.in
hot PN model
title hot PN model
c
c commands controlling continuum =========
table star rauch helium T = 240000 K, g = 8.5
luminosity 4.5 solar
c
c commands controlling geometry =========
sphere
radius 18
stop zone 1
set dr 0
c
c commands for density & abundances =========
abundances planetary
hden 3.0
c
c other commands for details =========
iterate
no level2
c
c commands controlling output =========
print last iteration
plot continuum range .1
punch overview last "stars_rauch_helium.ovr"
punch results last "stars_rauch_helium.rlt"
c
c stars_rauch_helium.in
c class stars
c ========================================
This is a model of a very hot planetary nebula, and checks that the code
is able to access the pure helium version of the Rauch grid of stellar
atmosphere models.
stars_rauch_hydr.in
hot PN model
title hot PN model
c
c commands controlling continuum =========
table star rauch hydr T = 240000 K, g = 8.5
luminosity 4.5 solar
c
c commands controlling geometry =========
sphere
radius 18
stop zone 1
set dr 0
c
c commands for density & abundances =========
abundances planetary
hden 3.0
c
c other commands for details =========
iterate
no level2
c
c commands controlling output =========
print last iteration
plot continuum range .1
punch overview last "stars_rauch_hydr.ovr"
punch results last "stars_rauch_hydr.rlt"
c
c stars_rauch_hydr.in
c class stars
c ========================================
This is a model of a very hot planetary nebula, and checks that the code
is able to access the pure hydrogen version of the Rauch grid of stellar
atmosphere models.
stars_rauch_pg1159.in
hot [WR] type PN model
title hot [WR] type PN model
c
c commands controlling continuum =========
table star rauch pg1159 T=150000K g=6.5
luminosity 4.5 solar
c
c commands controlling geometry =========
sphere
radius 18
stop zone 1
set dr 0
c
c commands for density & abundances =========
abundances planetary
hden 3.0
c
c other commands for details =========
iterate
no level2
c
c commands controlling output =========
print last iteration
plot continuum range .1
punch overview last "stars_rauch_pg1159.ovr"
punch results last "stars_rauch_pg1159.rlt"
c
c stars_rauch.in
c class stars
c ========================================
This is a model of a very hot planetary nebula, and checks that the code
is able to access the PG1159 version of the Rauch grid of stellar
atmosphere models.
stars_starburst99.in
demonstrate use of Starburst 99 spectrum
title demonstrate use of Starburst 99 spectrum
c
c commands controlling continuum =========
c read the table starburst and interploate to an age of 1e8 years
table star log age=8 "starburst99.mod"
ionization -2
c add background continuum - starburs 99 does not extend
c into fir so code would complain about zero continuum
c adding the cosmic background stops this complaint
background
c
c commands for density & abundances =========
hden 5
c
c commands controlling geometry =========
stop zone 1
c
c other commands for details =========
c iterate since needed for accurate punch continuum
iterate
constant temper 4
c
c commands controlling output =========
punch continuum "stars_starburst99.con" last
punch continuum last "stars_starburst99.con" units microns
c
c
c stars_starburst99.in
c class stars
c ========================================
c
this is a demonstration of the use of a Starburst 99 spectrum. It was
provided to me by anand Srianand, and we used it in Srianand et al. 2003.
A constant temperature is set since this is to test shape of continuum
not thermal physics. It also tests access to user-supplied grids as the
mechanism is the same.
stars_werner.in
test run with Werner stellar atmosphere
title test run with Werner stellar atmosphere
c
c commands controlling continuum =========
table star werner 190000 g=7.5
luminosity total 38
c
c commands for density & abundances =========
hden 4
abundances planetary no grains
grains agm no qheat single
c
c commands controlling geometry =========
radius 17
sphere
set dr 0
stop zone 1
c
c other commands for details =========
* >>chng 06 jul 01, add this to stars set since want to test shape
* of continuum and resulting ion dist, not thermal equilibrium of nebula
constant temperature 4.373
iterate 3
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_werner.ovr"
punch results last "stars_werner.rlt"
punch transmitted continuumlast "stars_werner.trn"
c
c stars_werner.in
c class stars
c ========================================
c
This checks that the code can access Kevin Volk39s Werner atmospheres.
Checks:
- table star Werner atmosphere works.
constant temperature is set to that this tests shape of stellar continuum
stars_wmbasic.in
Model of a Compact HII Region
title Model of a Compact HII Region
c
c commands controlling continuum =========
c Approximate model of a hot star ionizing its surroundings.
c Assume H density of 1e+6 per cc. This run tests the WMBASIC model option.
c I choose the T_eff, log(g), and log(Z) values arbitrarily.
table star wmbasic 30000.0 4.0 -0.3
luminosity 4.509202522 solar
table star wmbasic 32000.0 4.0 0.0
luminosity 4.509202522 solar
table star wmbasic 40000.0 3.6 -0.3
luminosity 4.509202522 solar
table star wmbasic 35000.0 4.0 -0.15
luminosity 4.509202522 solar
table star wmbasic 57000.0 4.0 0.0
luminosity 4.509202522 solar
table star wmbasic 37000.0 3.4 -0.3
luminosity 4.509202522 solar
table star wmbasic 32000.0 4.0 -0.1
luminosity 4.509202522 solar
table star wmbasic 36000.0 3.5 -0.12
luminosity 4.509202522 solar
c
c commands for density & abundances =========
abundances old solar 84
hden 6.0
grains ism
c
c commands controlling geometry =========
radius -2.0 parsec
sphere
stop zone 1
set dr 0
c
c other commands for details =========
c
c in this model the Lyman lines are pumped significantly
c by the stellar continuum; disable this process to avoid critical
c dependence on the precise shape of the continuum.
no induced processes
constant temperature 4
iterate
c
c commands controlling output =========
print last iteration
c
punch overview last "stars_wmbasic.ovr"
punch results last "stars_wmbasic.rlt"
punch transmitted continuum last "stars_wmbasic.trn" no header
c
c stars_wmbasic.in
c class stars
c ========================================
This is a test that the code can correctly access the
WMBASIC O-star model atmosphere continua.
Checks:
table star wmbasic command works.