# This example 1 of an input file for STOKES represents a system of an
# electron-dominated emitting and scattering disk inside the funnel of
# a dusty torus. The disk emits a continuum in the optical/UV.


# set data file name, model is run on only 1 CPU
OutputFile Disk+Torus 1


# number of photons sampled
PhotonNum 50000000


# wavelength range considered
LambdaMin 2000
LambdaMax 8000


# viewing directions:
#
#    1O directions in theta
#    1 direction in phi (axis-symmetric system)
#    symmetry with respect to the xy-plane
#    30 channels of spectral resolution
#
ThetaViewAng 10
PhiViewAng 1
PlaneSym yes
SpectRes 30


# 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
#    spectral index = alpha = 2
#
#             R        a       b    alpha   rad. velocity
ContSource 0.0004  3.25e-7  0.0004    2           0


# scattering region 1: electron disk (same as the 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 thick dusty torus
#
#    inner radius = R-b = 0.25 pc
#    outer radius = R+b = 100 pc        # optically thick torus with
#    torus height = 2*a = 17.3206 pc    # an opening angle of 30 deg
#    dust density = 0.001/ccm     
#
#           R        a        b    dust/elec.   density   rad. velocity
Torus    50.125   8.6603   49.875     dust       0.001          0


# dust composition:
#
#    dust with 15% graphite and 85% astronomical silicate
#    grain radii ranging from 0.005 to 0.200 micrometer
#    grain size index of -2.05
#    computed at 200 wavelengths and 200 grain radii
#    saved to dust model "AGN"
#
DustComp 5 10 85                        # The optical symmetry of graphite
GrainRadMin 0.005                       # grains requires twice as much 
GrainRadMax 0.200                       # ParaGraphite as OrthoGraphite.   
GrainSizeInd -2.05                      # The percentages of all components
GrainLambdaNum 200                      # must sum up to 100.
GrainRadNum 200
DustModel AGN

end