Input variables
damocles-code edited this page Aug 3, 2020
·
3 revisions
Input variables are described in the input files themselves and descriptions should be self-explanatory. The variables are summarised below:
All general variables relating to the models are defined in the input.in
file.
========================GENERAL OPTIONS================================================
false ! create a unique folder to store all data? (false will overwrite current output files)
false ! using observed data?
'input/line.in' ! data file containing observed line
'input/exclude.in' ! data file containing details of regions to exclude from chi calculation
false ! doublet?
true ! include multiple scattering events?
false ! use line of sight?
false ! use multiple lines of sight? (complete coverage)
false ! include e- scattering?
'input/e_scat.in' ! e- scattering file
1100 ! day no (used to calculate R_out)
hg ! dust scattering (enter isotropic or hg (Henyey-Greenstein))
========================GEOMETRY OPTIONS===============================================
false ! decouple dust and gas? (If not decoupled, gas distribution will follow the dust distribution)
shell ! dust geometry? Enter 'shell' or 'torus' or 'arbitrary'
'input/dust.in' ! dust geometry file (required - specifies parameters to create grid and specify dust mass)
'input/species.in' ! name of file detailing dust species and grain sizes
'input/dust_grid.in' ! dust grid file neme for use in arbitary geometries
shell ! gas geometry? Enter 'shell' or 'torus' or 'arbitrary'
'input/gas.in' ! gas geometry file (required - - specifies parameters to create grid or is the grid file)
=========================GRID==========================================================
40 ! number of cells in each direction
20 ! number of divisions in phi and theta to generate lines of sight covering whole ejecta
=======================PHOTON PACKETS==================================================
100000 ! number of photons
1000 ! number of frequency bins for output energy distribution
8 ! number of threads to run on
=======================VELOCITY LAW====================================================
true ! use a velocity law that is independent of radius?
2200 ! if true, max vel
300 ! if true, min vel
0.0 ! if true, power-law vel distribution
All dust variables except those controlling species and grain radius distributions.
================DUST OPTIONS=========================================
0.0e-5 ! total dust mass (Msun)
================DUST CLUMPING OPTIONS================================
0.0 ! mass fraction in clumps (set to 0 for no clumping) [remainder will follow smooth power-law as specified]
0.2 ! clump volume filling factor (volume in clumps)
2.0 ! clump number density distribution
================DUST GEOMETRY========================================
1900 ! max dust velocity in km/s
0. ! Rout in 1e15cm. Set to 0 to couple to maximum velocity (r=vt).
0.0001 ! Rin/Rout ratio
1.0 ! value of l where velocity profile is v~r^l
1.2 ! value of where dust density profile is rho~r^-q
2
NOTE THAT RELATIVE SPECIES WEIGHTS ARE FRACTIONAL CROSS_SECTIONAL AREAS - OTHER WEIGHTINGS CALCULATED WITHIN THE CODE
id filename weight amin amax power no of divisions
1 'dustData/amC-zb1.nk' 0.8 0.01 0.5 1.0 100
2 'dustData/sil-dlee.nk' 0.2 0.01 1.0 1.0 100
========================LINE===========================================================
500.7 ! wavelength (nm) of line or first doublet component
495.9 ! wavelength (nm) of second doublet component (enter 0.0 if n/a)
4.12e2 ! total observed luminosity (e40 ergs-1) [not required unless modelling absolute fluxes]
9.8e3 ! total flux of line to scale model to (in e15)
2.98 ! flux ratio of line 1 to line 2 (0 if n/a) [must be ratio of first line to second line]
========================SHELL GEOMETRY=================================================
0 ! fraction of emission emanated from clumps (currently only 0 or 1)
0 ! gas clump filling factor
0 ! gas clump power distribution
=======================================================================================
2000. ! max gas velocity in km/s
0. ! Rout in 1e15cm. Set to 0 to couple to velocity (r=vt).
0.05 ! Rin/Rout ratio for gas shell
1.0 ! value of l where gas velocity profile is v~r^l
1.2 ! value of where gas density profile is rho~r^-q
2.0 ! value of b where emissivity~r^-qb (applicable only if dust and gas coupled)
=====================ELECTRON SCATTERING PARAMETERS==============================
0.2 ! total Halpha luminosity in e40 ergs-1 (for electron densities)
10000 ! electron scattering temperature (set to 0 if no es)