# This example 2 of an input file for STOKES represents a system of an
# electron-dominated emitting and scattering disk and a pair of polar
# scattering cones. The double-cone is highly ionized closer in and
# dominated by dust farther out. It represents a radial outflow.
# Continuum emission from the disk and an H-alpha line are modeled. A
# broad line component comes from a cylindrical region surrounding the
# disk and a narrow line component is produced in the outer part of
# the double-cone.
# set data file name, model is run on only 1 CPU
OutputFile Disk+Cones 1
# number of photons sampled
PhotonNum 50000000
# wavelength range considered
LambdaMin 5500
LambdaMax 7500
# viewing directions:
#
# 5 directions in theta
# 1 direction in phi (axis-symmetric system)
# symmetry with respect to the xy-plane
# 50 channels of spectral resolution
#
ThetaViewAng 5
PhiViewAng 1
PlaneSym yes
SpectRes 50
# continuum emission region: electron disk
#
# inner radius = R-b = 0 pc
# outer radius = R+b = 8e-4 pc
# disk height = 2*a = 6.5e-7 pc
# spectral index = alpha = 1
#
# R a b alpha rad. velocity
ContSource 0.0004 3.25e-7 0.0004 1 0
# broad line emission region: cylindrical
#
# inner radius = R-b = 0.002 pc
# outer radius = R+b = 0.012 pc
# cylinder height = 2*a = 0.01 pc
# line centroid = lambda0 = 6563 Angstroem
# line width = Gamma = 400 Angstroem
# rel. strength = 30%
#
# R a b lambda0 Gamma strength rad. velocity
BLSource 0.007 0.005 0.005 6563 400 30 0
# narrow line emission region: double-conical
#
# inner radius = r1 = 25 pc
# outer radius = r2 = 75 pc
# half-op. angle = theta = 30 deg
# line centroid = lambda0 = 6563 Angstroem
# line width = Gamma = 40 Angstroem
# rel. strength = 20%
# rad. velocity = 1000 km/s
#
# r1 r2 theta lambda0 Gamma strength rad. velocity
NLSource 25 75 30 6563 40 20 1000
# scattering region 1: electron disk (same as the continuum emission region)
#
# inner radius = R-b = 0 pc
# outer radius = R+b = 8e-4 pc # The optical depth from the
# disk height = 2*a = 6.5e-7 pc # equatorial plane to the
# elec. density = 4.5e12/ccm # disk surface is 3.
#
# R a b dust/elec. density rad. velocity
Cylinder 0.0004 3.25e-7 0.0004 electrons 4.5e12 0
# scattering region 2: optically thin electron cones
#
# inner radius = r1 = 1 pc
# outer radius = r2 = 20 pc # optically thin double-cone
# half-op. angle = theta = 30 deg # with an optical depth of 0.5
# elec. density = 12800/ccm # along the polar direction
# rad. velocity = 6000 km/s
#
# r1 r2 theta dust/elec. density rad. velocity
Double-cone 1 20 30 electrons 12800 6000
# scattering region 3: optically thick dusty cones (same as narrow line region)
#
# inner radius = r1 = 25 pc
# outer radius = r2 = 75 pc # optically thick dusty cones
# half-op. angle = theta = 30 deg
# dust density = 0.001/ccm
# rad. velocity = 1000 km/s
#
# r1 r2 theta dust/elec. density rad. velocity
Double-cone 25 75 30 dust 0.001 1000
# dust composition:
#
# dust with 37.5% graphite and 62.5% astronomical silicate
# grain radii ranging from 0.005 to 0.250 micrometer
# grain size index of -3.5
# computed at 200 wavelength and 200 grain radii
# saved to dust model "MilkyWay"
#
DustComp 12.5 25 62.5 # The optical symmetry of graphite
GrainRadMin 0.005 # grains requires twice as much
GrainRadMax 0.250 # ParaGraphite as OrthoGraphite
GrainSizeInd -3.5 # The percentages of all components
GrainLambdaNum 200 # must sum up to 100.
GrainRadNum 200
DustModel MilkyWay
end