# This example 3 of an input file for STOKES represents a system of a
# point source emitting a narrow [0III] 5007 emission line. The source
# is shielded by a rapidly expanding electron-dominated shell. A few
# spherical dust clouds are located further out. The clouds are
# distributed asymmetrically, hence the polarization is evaluated at
# several azimuthal angles.


# set data file name, model is run on only 1 CPU
OutputFile Shell+Cloud 1


# number of photons sampled
PhotonNum 30000000


# wavelength range considered
LambdaMin 4500
LambdaMax 7000


# viewing directions:
#
#    6 directions in theta
#    12 directions in phi
#    no symmetry with respect to the xy-plane
#    50 channels of spectral resolution
#
ThetaViewAng 6
PhiViewAng 12
PlaneSym no
SpectRes 50


# narrow line emission region: point-like
#
#    inner radius    = r1 = 0 pc
#    outer radius    = r2 = 0 pc
#    cylinder height = theta = 0 deg
#    line centroid   = lambda0 = 5007 Angstroem
#    line width      = Gamma = 25 Angstroem
#    rel. strength   = 100%
#
#           r1   r2   theta   lambda0    Gamma    strength   rad. velocity
NLSource    0    0      0      5007       25        100            0  


# scattering region 1: optically thin electron shell
#
#    inner radius   = r1 = 1 pc            
#    outer radius   = r2 = 3 pc            # expanding spherical shell
#    half-op. angle = theta = 90 deg       # with an optical depth of 1
#    elec. density  = 244000/ccm
#    rad. velocity  = 5000 km/s
#
#             r1    r2    theta    dust/elec.   density   rad. velocity
Double-cone   1     3      90      electrons    244000        5000


# scattering region 2: optically thick dust cloud
# 
#    position       = (x0,y0,z0) = (10 pc, 10pc, 10pc)   # optically thick
#    radius         = r = 3 pc                           # dust sphere    
#    dust density   = 0.01/ccm
#
#             x0   y0   z0   r    dust/elec.   density   rad. velocity
Sphere        10   10   10   3       dust       0.01           0


# scattering region 3: optically thick dust cloud
# 
#    position       = (x0,y0,z0) = (-10 pc, 8pc, 3pc)    # optically thick
#    radius         = r = 3 pc                           # dust sphere    
#    dust density   = 0.01/ccm
#
#             x0   y0   z0   r    dust/elec.   density   rad. velocity
Sphere        -10  8    3    3       dust       0.01           0


# scattering region 4: optically thick dust cloud
# 
#    position       = (x0,y0,z0) = (5 pc, -10pc, 0pc)    # optically thick
#    radius         = r = 3 pc                           # dust sphere    
#    dust density   = 0.01/ccm
#
#             x0   y0   z0   r    dust/elec.   density   rad. velocity
Sphere        5    -10  0    3       dust       0.01           0


# scattering region 5: optically thick dust cloud
# 
#    position       = (x0,y0,z0) = (1 pc, 0 pc, -8 pc)   # optically thick
#    radius         = r = 3 pc                           # dust sphere    
#    dust density   = 0.01/ccm
#
#             x0   y0   z0   r    dust/elec.   density   rad. velocity
Sphere        1    0    -8   3       dust       0.01           0


# Galactic dust model read from a file (see example 2)
DustModel MilkyWay

end