title test generating line ratios in a grid run c c commands controlling continuum ===== blackbody 4e4 K ionization parameter -2 c c commands for density & abundances ===== c these are to speed up the calculation, only do H, O, and Ne init "honly.ini" element oxygen on element neon on element sulphur on element oxygen ionization 1 1 1 0.01 element neon ionization 1 1 1 0.01 element sulphur ionization 1 1 1 0.01 c vary the hydrogen density hden 4 vary grid 2 6.1 1 c c other commands for details ======= c these are constant temperature models, vary T constant temperature 4 vary grid 4000 17000 3000 linear stop zone 1 c c commands controlling output ===== save overview "func_grid_line_ratios.ovr" save performance "func_grid_line_ratios.per" save line list "func_grid_line_ratios.pun" "func_grid_line_ratios.dat" ratio no hash save grid "func_grid_line_ratios.grd" c c commands giving (lack of) assert ========== monitor nothing 0 c c func_grid_line_line_ratios.in c class function c =================================================== This uses the grid command to compute line ratios for a wide range of density and temperature. The ionization is set to a uniform value and only a few elements are included. this makes the calculation faster and prevents recombination [O III] 4363 from becoming important (there is no O+3). These are the line ratios mentioned as limits in the Johnstone et al. Spitzer cooling flow filament paper (2007).