MAIN
# This section is quite free form: parameters may appear in any order,
# lines begining with #, % or blank space are for comments, otherwise
# from column 1 to 10 an identifier of the parameter (case sensitive)
# from 11 to 20 its value (always real with explicit decimal point)
# from 21 on comments.
#2345678901234567890123456789
#         
#istarfo=0  SFR from che_evo code
#istarfo=1  analytical SFR (parameters in section ANALY).
#istarfo=2  SFR from a chemical evolution code. When required the 
# residual gas fraction can be imposed through gastmp and igastmp in
# section MAIN
#istarfo=3  burst superimposed to the SFR provided by a chemical 
# evolution code. Parameters in section BURST
#istarfo=4  SFR from semianalytical hierarchical clustering models
#rsfmin, dtcurmax, dtcurmin are specific of istarfo 2 or 3.
#their value should not be changed
tgal      13.00      age of the model
# NOTE: when tgal < 0 it is ignored and read instead in age.par.
# this is a trick to be used in programs running GRASIL for the same
# model at different ages. See for instance program AGETAB.
# tau_eff = sqrt(tau_a*tau_e) is used in the diffuse component,
# if 0.0 tau_a is used
gastmp    1.5e-4    residual gas fraction is forced to gastmp when
# istarfo=2 and igastmp=1
igastmp   0.0       see gastmp
mmolfraz  0.0       fraction of molecular over total gas
# NOTE: if mmolfraz<0 no gas, useful for bare EPS
mcloud    1.0e6     Mass (Msun) of molecular clouds (MCs)
rcloud    16.0      radius (pc) of molecular clouds
etastart  0.1       Gyr see below
# fraction f of SSP energy inside MCs depends on its age t according to
# f=1 is t<etastart
# f=1-vout(t-etastart) if etastart<t<etastart+1/vout
# f=0 otherwise
flutflag  0.0       if 0.0 no computation of T fluctuations of small
# grains (speeds up), if 1.0 fluctuations are computed
dsug      0.4545    
idz       1.0            
# if idz=0.0 dsug is simply the dust/gas mass ratio, if idz=1.0
# (not meaningful for istarfo=1)
# then dsug becomes the proportionality factor between Z and dust/gas.
# in this case dsug=0.45 gives the standard DL dsug=1/110 for Zsun
tgra      400.0     maximum T of graphite grains in MCs
tsil      400.0     maximum T of silicate grains in MCs
igeo      1.0       sets the geometry of stars and gas
#igeo=0 homogeneous
#igeo=1 King profiles: core radii given by rcstar and rcdiff
#igeo=2 Exponentials of distances from polar axis and from equator.
#scale lengths given by (rdstar,rddiff) and  (zdstar,zddiff) 
rcstar    0.5       [kpc] igeo=1
rcdiff    6.0       [kpc] igeo=1  (not used for istarfo=4)
rdstar    3.0       [kpc] igeo=2  (not used for istarfo=4)
rddiff    6.0       [kpc] igeo=2  (or rddiff/rdstar for istarfo=4)
zdstar    1.0       [kpc] igeo=2  (or zdstar/rdstar for istarfo=4)
zddiff    2.0       [kpc] igeo=2  (or zddiff/rdstar for istarfo=4)
#
#
SPECTRA
Directions in degrees from the polar axis in which the SEDs are computed:
1  number of direction
0.0  degrees
45.0 degrees
90.0 degrees
#
#
SSP
Names of SSPs to be used. List ended by a comment line
First column put 1 if ssp is to be used, 0 otherwise. 
1    'Z0001.sal'
1    'Z0004.sal'
1    'Z004.sal'
1    'Z008.sal'
1    'Z02.sal'
1    'Z05.sal'
1    'Z1.sal'
#
che_evo must be run with an IMF identical to that used to generate the SSP
files (sal=Salpeter, ken=Kennicut, sca=Scalo, msc=Miller & Scalo). The only
thing that can be (slightly) different is the lower mass limit (in the SSP it
is always 0.15), as long as the value used in che_evo affects only the
normalization (i.e. M/L) but not the shape of the SSPs.
#
MIXDIF
# Parameters setting the dust mixture in the diffuse component
# This section is quite free form: parameters may appear in any order,
# lines begining with #, % or blank space are for comments, otherwise
# from column 1 to 10 an identifier of the parameter (case sensitive)
# from 11 to 20 its value (always real with explicit decimal point)
# from 21 on comments.
#2345678901234567890123456789
#parameters for graphite follow  
amin_gr   0.0008e-4 minimum radius
amax_gr   0.25e-4   maximum radius
ab_gr     0.0050e-4 break radius between the two power laws
beta1_gr  -3.5      power law index from amax to ab
beta2_gr  -4.05     pawer law index from ab to amin
corgra    1.00      correction factor with respect to abundance in DL model
# can be used to vary the ratio between graphite and silicate
# parameters for silicate follow  
amin_si   0.0050e-4 minimum radius
amax_si   0.25e-4   maximum radius
ab_si     0.0003e-4 break radius between the two power laws
beta1_si  -3.5      power law index from amax to ab
beta2_si  -3.95     pawer law index from ab to amin
corsil    1.00      correction factor with respect to abundance in DL model
# can be used to vary the ratio between graphite and silicate
# parameters for PAHs follow
xx        3.03e-25  normalization power law PAHs
xh        0.2       hydro coverage PAHs
#
MIXMOL
# Parameters setting the dust mixture in the molecular clouds
# This section is quite free form: parameters may appear in any order,
# lines begining with #, % or blank space are for comments, otherwise
# from column 1 to 10 an identifier of the parameter (case sensitive)
# from 11 to 20 its value (always real with explicit decimal point)
# from 21 on comments.
#2345678901234567890123456789
#parameters for graphite follow  
amin_gr   0.0008e-4 minimum radius
amax_gr   0.25e-4   maximum radius
ab_gr     0.0050e-4 break radius between the two power laws
beta1_gr  -3.5      power law index from amax to ab
beta2_gr  -4.05     pawer law index from ab to amin
corgra    1.00      correction factor with respect to abundance in DL model
# can be used to vary the ratio between graphite and silicate
# parameters for silicate follow  
amin_si   0.0050e-4 minimum radius
amax_si   0.25e-4   maximum radius
ab_si     0.0003e-4 break radius between the two power laws
beta1_si  -3.5      power law index from amax to ab
beta2_si  -3.95     pawer law index from ab to amin
corsil    1.00      correction factor with respect to abundance in DL model
# can be used to vary the ratio between graphite and silicate
# parameters for PAHs follow
xx        3.03e-26  normalization power law PAHs
xh        0.2       hydro coverage PAHs
#
# is better to leave some comment line at the end