[pyrat]
# Pyrat Bay run mode: [tli atmosphere spectrum opacity retrieval radeq]
runmode = atmosphere

# Output file name and verbosity
logfile = WASP39b_hybrid_atmosphere.log
verb = 2

# Pressure profile
ptop = 1.0e-8 bar
pbottom = 100.0 bar
nlayers = 81

# Temperature profile [isothermal guillot madhu]
tmodel = isothermal
tpars = 900.0

# Chemistry model and composition [free equilibrium]
chemistry = equilibrium
species =
    H  He  C  N  O  F  Na  Mg  Si  S  K  Ti  V Fe
    H2  H2O  CH4  CO  CO2  HCN  NH3  N2  OH  C2H2  C2H4 OCS CS CS2
    S2  SH  H2S  SO2  SO  TiO  VO  TiO2  VO2   SiO  SiH  SiS  SiH4
    e-  H-  H+  H2+  He+  Na-  Na+  Mg+ K-  K+  Fe+  Ti+  V+ SiH+ Si- Si+

# Scale independently the elemental abundances of carbon, oxygen, and
# all other metals (dex units, relative to solar)
vmr_vars =
    [M/H]     1.0
    [C/H]     0.8
    [O/H]     1.1
    log_SO2  -5.0