hden 4 abundances ISM radius 18.5 cm blackbody 4.36e4 K Q(H) 49.34 stop zone 1 save total opacity "sopacity.opc" units microns