# 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