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
no level 2 lines // slight speed up - here for consistency with old ini file
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
save overview "agn_S_curve_grid.ovr"
save performance "agn_S_curve_grid.per"
save grid "agn_S_curve_grid.grd" last
c save average H+ temperature
save average "agn_S_curve_grid.avr" last
temperature hydrogen 1 
end of temperatures 
save monitors "agn_S_curve_grid.asr"
c save heating and cooling
save heating "agn_S_curve_grid.het" last
save cooling "agn_S_curve_grid.col" last
save iron element "agn_S_curve_grid.fe" last 
save xspec atable "agn_S_curve_grid.fit" 
save neon element "agn_S_curve_grid.ne" last separate
c
c commands giving the monitors    =========
// >>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
// >>chng 11 jul 13, several changes due to collisional suppression of DR
// >>chng 12 feb 27, from 5.152 to 5.127; upgrade to Chianti 7
monitor hydro 2 temper grid "agn_S_curve_grid.dat" err 0.05
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.