set dr 0
stop zone 1
set save prefix "M16"
radius 15
#
# just do a slab 1 cm thick for speed
set dr 0
stop zone 1
# <== this is a comment
#
# the cloud is ionized by a nearby star cluster, NGC 6611.  The brightest
# star is  the O3 - O5 V star W205 which is about 2kpc distant
# McLeod+15 http://adsabs.harvard.edu/abs/2015MNRAS.450.1057M
#
# 50 000 K blackbody, roughly an O3 - O5 V star
blackbody, T=4.87e4 K # the AGN3 Table 2.3 entry for O4 V
Q(H)  49.70
#
# RADIUS gives the separation between the star and the cloud.
# units are log cm.  The projected separation between star 
# and nebula is about 2 pm according to McLeod.  This is
# slightly more than 2 pc (log 2 pc 18.78 cm)
#
# this is the log of the hydrogen density
hden 3
#
# this uses a standard set of H II region 
# abundances, including grains
abundances HII region
#
# not important in the H II region, 
# but will be critical when we
# extend it to the PDR
cosmic ray background
#
# cosmic microwave background at z=0
CMB
#
# we must iterate at least one time
# to establish line optical depths
iterate
# we only want the output for the last iteration
print last iteration
#
# 
# this will save the temperature and ionization of the cloud
save overview ".ovr" last
#
# save emission from our unit cell
save emitted continuum ".con" units microns last
#