root/trunk/tsuite/auto/agn_S_curve_grid.in

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 will be varied
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 will be varied
c create file with list of grid parameters
punch grid "agn_S_curve_grid.grd" last
c save average H+ temperature
punch average "agn_S_curve_grid.avr" last
temperature hydrogen 1 
end of temperatures 
punch asserts "agn_S_curve_grid.asr"
c save heating and cooling
punch heating "agn_S_curve_grid.het" last
punch cooling "agn_S_curve_grid.col" last
punch iron element "agn_S_curve_grid.fe" last 
punch xspec atable "agn_S_curve_grid.fit" 
c
c commands giving the asserts    =========
// >>chng 07 jan 25, use only Gu uta
// asse hydro 2 temper 4.512 err0.05 grid 4.600 4.720 4.899 5.206 5.438 5.642 5.997 6.219 6.359 6.555 6.779 6.927
// >>chng 08 oct 18, from 5.519 to 5.543, include Badnell UTA data
asse hydro 2 temper 4.512 err0.05 grid 4.600 4.746 4.964 5.543 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 computes a series of models that check the temperature
through the S curve in the Fields et al. three-phase model of
ISM stability.