# -*- coding: utf-8 -*-
"""
Created on Fri Sep 16 15:19:15 2022

@author: s7970488
"""

#
#
#Custom filename tags to distinguish from other runs
project = "Test"
runVers = "PP"

#Default plot

makePlotStruc = False
makePlotSED = False
makePlotSpec = False
makePlotLDC = True
makePlotFT = False
makePlotTLA = False
makePlotTrans = True
makePlotPPress = False
plotSpec = 'H'

#Spectrum synthesis mode
# - uses model in Restart.py with minimal structure calculation
specSynMode = False
whichLineList = "NIST"
ifSphere = False

#Model atmosphere

##Canonical grid-point Sun
#teff = 5800.0  #,    K
#logg = 4.5 #,      cgs
#log10ZScale = 0.0     # [A/H]
#massStar = 1.0 #,      solar masses

##Neilson & Lester 2013 Fig. 1 parameters:
#teff = 5000.0  #,    K
#logg = 2.0 #,      cgs
#log10ZScale = 0.0     # [A/H]
#massStar = 10.0 #,      solar masses    

#Neilson & Lester 2012 Fig. 6 parameters:
teff = 5000.0  #,    K
logg = 2.0 #,      cgs
log10ZScale = 0.0     # [A/H]
massStar = 5.0 #,      solar masses    

##Neilson et al. 2017 Fig. 3 parameters:
#    # Similar to Figs. 1 and 2 parameters
#teff = 5700.0  #,    K
#logg = 4.5 #,      cgs
#log10ZScale = 0.0     # [A/H]
#massStar = 1.0 #,      solar masses    

##Fieldus, Lester, Rogers, 1990
#teff = 10000.0  #,    K
#logg = 1.5 #,      cgs
#log10ZScale = 0.0     # [A/H]
#massStar = 2.0 #,      solar masses

##Red giant - Arcturus-like
#teff = 4200.0  #,    K
#logg = 4.5 #,      cgs
#log10ZScale = 0.0     # [A/H]
#massStar = 0.5 #,      solar masses

xiT = 1.0  #,       km/s
logHeFe = 0.0  #,   [He/Fe]
logCO = 0.0  #,   [C/O]
logAlphaFe = 0.0   #,   [alpha-elements/Fe]


#Spectrum synthesis
lambdaStart = 588.5  #,       nm    
lambdaStop = 590.0  #,     nm
##Halpha
#lambdaStart = 655.5  #,       nm    
#lambdaStop = 657.5  #,     nm
#lambdaStart = 700.0  #,       nm    
#lambdaStop = 701.0  #,     nm

fileStem = project + "-"\
 + str(round(teff, 7)) + "-" + str(round(logg, 3)) + "-" + str(round(log10ZScale, 3))\
 + "-" + str(round(lambdaStart, 5)) + "-" + str(round(lambdaStop, 5))\
 + "-" + runVers  

lineThresh = -3.0  #,    min log(KapLine/kapCnt) for inclusion at all - areally, being used as "lineVoigt" for now
voigtThresh = -3.0  #,     min log(KapLine/kapCnt) for treatment as Voigt - currently not used - all lines get Voigt
logGammaCol = 0.0
logKapFudge = 0.0
macroV = 0.0  #,     km/s
rotV = 0.0  #,   km/s
rotI = 90.0 #,    degrees
RV = 0.0 #,   km/s
vacAir = "vacuum"    
sampling = "fine"
ifSweep = True #sweep wavelength grid
ifVoigt = True

#Performance vs realism
nOuterIter = 12   #,     no of outer Pgas(HSE) - EOS - kappa iterations
nInnerIter = 12   #,    no of inner (ion fraction) - Pe iterations
ifMols = 1   #,     where to include TiO JOLA bands in synthesis 

#Gaussian filter for limb darkening curve, fourier transform
diskLambda = 500.0  #,      nm
diskSigma = 0.01  #,     nm

#Two-level atom and spectral line
userLam0 = 589.592  #,   nm 
userA12 = 6.24  #,    A_12 logarithmic abundance = log_10(N/H_H) = 12
userLogF = -0.495   #,  log(f) oscillaotr strength // saturated line
userStage = 0  #,   ionization stage of user species (0 (I) - 3 (IV)
userChiI1 = 5.139   #,  ground state chi_I, eV
userChiI2 = 47.29   #,  1st ionized state chi_I, eV
userChiI3 = 71.62   #,  2nd ionized state chi_I, eV
userChiI4 = 98.94   #,  3rd ionized state chi_I, eV
userChiL = 0.0  #,   lower atomic E-level, eV
userGw1 = 2   #,  ground state state. weight or partition fn (stage I) - unitless
userGw2 = 1   #,  ground state state. weight or partition fn (stage II) - unitless
userGw3 = 1   #,   ground state state. weight or partition fn (stage III) - unitless
userGw4 = 1  #,  ground state state. weight or partition fn (stage IV) - unitless
userGwL = 2  #,   lower E-level state. weight - unitless
userMass = 22.9   #,  amu
userLogGammaCol = 1.0   #,  log_10 Lorentzian broadening enhancement factor

#Planetary transit parameters for transit light curve modelling
# ** Also depends on rotI chosen above!!
# Oribital period is not a free parameter - it is set by the 
# size of the orbit and the planet's mass by basic form of 
#Kepler's 3rd law
rJupiter = 11.21 # Earth radii - handy reference
rSunAU = 0.00465047  #radius of the Sun in AU - handy reference

ifTransit = False
#Data source:  Wikipedia for now...
rOrbit = 1.0  #AU
rPlanet = 1.0 #Earth radii
#mPlanet = 1.0 #Earth masses #not needed (yet?)
