[pyrat]

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

# Output files
logfile = transmission_spectrum_tutorial.log

# Radiative-transer observing geometry, select from: [transit emission]
rt_path = transit

# Wavelength sampling:
wl_low  = 0.3 um
wl_high = 5.0 um

# System parameters:
rstar = 1.27 rsun
tstar = 5800.0
smaxis = 0.045 au
mplanet = 0.6 mjup
rplanet = 1.0 rjup
# Reference pressure level at rplanet:
refpressure = 0.1 bar

# Atmospheric model:
ptop = 1.0e-8 bar
pbottom = 100.0 bar
nlayers = 81

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

# Chemistry model [free equilibrium]
chemistry = free
species =       H2    He   Na  H2O  CH4   CO  CO2  NH3  HCN   N2
uniform_vmr = 0.85 0.149 5e-6 2e-3 1e-4 5e-4 1e-6 1e-5 1e-9 2e-4

# Radius-profile model [hydro_m hydro_g]
radmodel = hydro_m


# Line-sampled cross sections
sampled_cross_sec =
    cross_section_0.15-33.0um_0200-5000K_R025K_H2O_exomol_pokazatel.npz

# Continuum cross-section files:
continuum_cross_sec =
    {ROOT}/pyratbay/data/CIA/CIA_Borysow_H2H2_0060-7000K_0.6-500um.dat
    {ROOT}/pyratbay/data/CIA/CIA_Borysow_H2He_0050-3000K_0.3-030um.dat

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

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

# Cloud models, select from: [lecavelier deck ccsgray]
clouds =
    deck       -0.5
    lecavelier  0.0 -4.0

# Verbosity level (<0:errors, 0:warnings, 1:headlines, 2:details, 3:debug):
verb = 2

# If defined, plot x-axis in log scale and set ticks at logxticks locations:
logxticks = 0.3 0.5 0.7 1.0 2.0 3.0 5.0