[pyrat]

# Pyrat Bay run mode [tli atmosphere opacity spectrum radeq retrieval]
runmode = radeq
# Number of radiative-equilibrium iterations
nsamples = 250

# Output file names
logfile = radeq_WASP69b_1x_solar_metallicity.log

# Verbosity level [1--3]
verb = 2

# Observing geometry, select between: [transit eclipse emission_two_stream]
rt_path = emission_two_stream

# Wavelength sampling
wl_low = 0.15 um
wl_high = 33.0 um
wl_thinning = 2

# System parameters (Bonomo et al. 2017)
rstar = 0.813 rsun
mstar = 0.826 msun
starspec = inputs/phoenix_WASP69b_4700K_logg4.5.dat

rplanet = 1.057 rjup
mplanet = 0.250 mjup
smaxis = 0.04527 au
refpressure = 1e-2 bar

# Atmospheric model
ptop = 1.0e-9 bar
pbottom = 100 bar
nlayers = 81

# Initial thermal profile for radiative-equilibrium runs
tmodel = isothermal
tpars = 900.0

# Chemistry and composition
chemistry = equilibrium
species =
    H  He  C  N  O  Na  K  Mg  Si  P  S  Cl Ti  V Fe
    H2  H2O  CH4  CO  CO2  HCN  NH3  N2  OH  C2H2  C2H4
    S2  SH  H2S  SO2  SO  TiO  VO  TiO2  VO2   SiO  SiH  SiS  SiH4
    OCS PH PH3 PN PO KOH CS CS2 NaH  NaCl  HCl  KCl NaOH
    e-  H-  H+  H2+  He+  Na-  Na+  Mg+ K-  K+  Fe+  Ti+  V+  SiH+ Si- Si+
# Elemental composition
vmr_vars =
    [M/H] 0.0
    C/O 0.59

# Hydrostatic equilibrium
radmodel = hydro_m

# Stellar irradiation scaling factor, beta = (1-A)*f
beta_irr = 0.66

# Internal heat (K)
tint = 100.0

# Line-by-line sampled cross sections
sampled_cross_sec =
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_H2O_exomol_pokazatel.npz
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_CO_hitemp_2019.npz
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_CO2_ames_ai3000k.npz
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_CH4_exomol_mm.npz
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_H2S_exomol_ayt2.npz
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_HCN_exomol_harris_larner.npz
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_NH3_exomol_coyute.npz
    inputs/cross_section_0.15-33.0um_0200-5000K_R025K_SiO_exomol_siouvenir.npz

# Alkali models, select from: [sodium_vdw potassium_vdw]
alkali =
    sodium_vdw
    potassium_vdw

# Continuum cross sections
continuum_cross_sec =
    {ROOT}/pyratbay/data/CIA/CIA_Borysow_H2H2_0060-7000K_0.6-500um.dat
    {ROOT}/pyratbay/data/CIA/CIA_Borysow_H2He_0050-7000K_0.5-031um.dat

# Rayleigh models, select from: [rayleigh_H rayleigh_He rayleigh_H2]
rayleigh =
    rayleigh_H2
    rayleigh_He
    rayleigh_H

# H- continous absorption
h_ion = h_ion_john1988