# 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
# no imaging
#
ThetaViewAng 5
PhiViewAng 1
PlaneSym yes
SpectRes 50
ResX 1
ResY 1
# continuum emission region: electron disk
#
# inner radius = R-b = 0 pc (no funnel),
# outer radius = R+b = 8e-4 pc
# disk height = 2*a = 6.5e-7 pc
# no vertical offset
# spectral index = alpha = 2
# no velocity field
#
# R a b c alpha spherical velocity
ContSource 0.0004 3.25e-7 0.0004 0. 2 0 0. 0. 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 = 40 Angstroem
# rel. strength = 30%
# rad. velocity = -500 km/s
# azim. velocity = 3000 km/s
#
# R a b lambda0 Gamma strength cylindrical velocity
BLSource 0.007 0.005 0.005 6563 400 30 1 -500. 3000. 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
# azim. velocity = 100 km/s
#
# r1 r2 theta lambda0 Gamma strength spherical velocity
NLSource 25 75 30 6563 40 20 0 1000. 100. 0.
# 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
# no vertical offset
# elec. density = 4.5e12/ccm # disk surface is 3.
# no velocity field
#
# R a b c dust/elec. density spherical velocity
Cylinder 0.0004 3.25e-7 0.0004 0. electrons 4.5e12 0 0. 0. 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
# azim. velocity = 1000 km/s
#
# r1 r2 theta dust/elec. density spherical velocity
Double-cone 1 20 30 electrons 12800 0 6000. 1000. 0.
# 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
# azim. velocity = 100 km/s
#
# r1 r2 theta dust/elec. density spherical velocity
Double-cone 25 75 30 dust 0.001 0 1000. 100. 0.
# 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