[pyrat]

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

# Output spectrum file:
logfile = ./transmission_spectrum_tutorial.log

# Radiative-transer geometry, select from
# [transit eclipse eclipse_two_stream emission emission_two_stream f_lambda]
rt_path = transit

# Atmospheric model:
ptop = 1e-7 bar
pbottom = 100 bar
nlayers = 81

tmodel = guillot
tpars = -5.4  -0.68  0.0  0.0  1240.0  100.0

chemistry = free
species = H2   He    H2O  Na   K
uniform_vmr = 0.85 0.145 1e-4 1e-6 1e-7

# Spectrum boundaries, sampling rate, and oversampling factor:
wl_low  = 0.4 um
wl_high = 8.0 um
resolution = 25_000

# System parameters:
rstar   = 1.0 rsun
rplanet = 1.35 rjup
mplanet = 0.70 mjup
refpressure = 0.1 bar

# Line-by-line opacities:
tlifile = ./HITRAN_H2O_0.4-8.0um.tli

# 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_1000-7000K_0.5-400um.dat

# Radius-profile model, select from: [hydro_m hydro_g]
radmodel = hydro_m

# Cloud models: [lecavelier deck]
clouds = lecavelier 0.0 -4.0

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

# Number of parallel processors:
ncpu = 7

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