title cool atomic ISM with Si grain quantum heating
c
c commands controlling continuum =========
table draine
extinguish 20
c
c commands for density & abundances =========
hden 0
c turns on grains with default abundance (log abundance ratio to ism = 0)
c and with quantum heating enabled by default
c (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     =========
cosmic rays background
iterate
c
c commands controlling output    =========
normalize to "C  2" 157.6m 
print line faint -2
print line column
save overview "grains_qheat.ovr"
save performance "grains_qheat.per"
save monitors last "grains_qheat.asr"
save grains temperature last "grains_qheat.qht"
save continuum  "grains_qheat.con" units microns last
save two photon continuum "grains_qheat.2nu" units microns
save qheat last "grains_qheat.qhd"
c
c commands giving the monitors    =========
// >>chng 10 sep 29, from -26.77 to -26.66, use atom_leveln for all temps
// >>chng 12 apr 09, from -26.66 to -26.22, update CR backbround to Indriolo+07, ~8x larger
monitor line intensity "C  2"  157.6m -26.22
//
// >>chng 12 apr 09, from -28.7 to -27.63, update CR backbround to Indriolo+07, ~8x larger
monitor line intensity "o  1"  63.17m -27.63 error 0.1
//
// >>chng 12 apr 09, from 64.4 to 136.3, update CR backbround to Indriolo+07, ~8x larger
monitor temperature hydrogen 2 136.3
// 
monitor line intensity "nFnu" 122.4m -23.976
monitor line intensity "nFnu" 60.83m -24.774
monitor line intensity "nFnu" 41.22m -25.713
monitor line intensity "nFnu" 30.87m -26.596
c
c grains_qheat.in
c class limit 
c ========================================
c 

This sim produces dust emission with a Wein trail that is
dominated by quantum heating emission.