param.dat

	###############################
# Input parameters for HELIOS #
###############################

### GENERAL ###                                                                 ### INFO ###

name =                                  0                                      (name of the output subdirectory and file prefix.)
output directory =                      ./output/                              (main path to output)
precision =                             double                                 (double or single. but actually only double makes sense.)
realtime plotting =                     no                                     (yes, no, interval value)
———

### GRID ###

isothermal layers =                     no                                     (set 'no' for iterative runs and 'yes' for post-processing only)
number of layers =                      105                                     (10 layers per magnitude of pressure as rule of thumb)
TOA pressure [10^-6 bar] =              1e-1                                   (typically 1e0 or 1e-1 for common applications)
BOA pressure [10^-6 bar] =              1e9                                    (bottom pressure. or surface pressure. for gas planets usually 1e3 or 1e2.)
planet type =                           gas                                    (rocky, gas)
———

### ITERATION ###

post-processing only =                  no                                     (yes: pure post-processing of given TP profile, no: iterative run)
path to temperature file =              ./output/0/0_tp.dat
temperature file format & P unit =      helios                                 (helios, TP (bar), PT (bar))
adaptive interval =                     20                                     (numerical iteration speed. fast: '6', slow & stable: '20')
TP profile smoothing =                  no                                     (use 'no' by default. smoothes sharp kinks in TP profile at the cost of energy balance.)
———

### RADIATION ###

direct irradiation beam =               no                                    ('yes': star is a point-source. 'no': isotropic irradiation assumed.)
scattering =                            yes
improved two-stream correction =        no                                     ('yes', 'no')
asymmetry factor g_0 =                  0                                      (values: [-1, 1]. 0: isotropic scattering. 1: pure forward, -1: pure backward.)
path to opacity file =                  ./input/opac_sample.h5 
diffusivity factor =                    2                                     (usually set to 2)
second Eddington coefficient =          0.5                                    (modulates the scattering of the direct irradiation beam. per default set to 0.5)
f factor =                              0.667                                   (0.667: dayside only, no redistribution, 0.25: global redistribution, 1: substellar point)
stellar zenith angle [deg] =            60                                     (values: 0 - 89, with 0 being vertical irradiation. exclusive with setting the f factor.)
geometric zenith angle correction =     no                                    (recommended: 'no' if zenith angle < 70. 'yes' if angle >= 70)
internal temperature [K] =              100
surface/BOA temperature [K] =           0                                      (value only relevant for post-processing. prioritizing value from TP-file, if exists.)
surface/BOA albedo =                    0.0                                    (set to '0' for no surface. maximum value < 1.0)
energy budget correction =              yes                                     (set to 'yes' for iterative run and 'no' for post-processing only.)
———

### CONVECTIVE ADJUSTMENT ###

convective adjustment =                 yes
kappa value =                           0.2857                                 (file, number)
entropy/kappa file path =               insert_path_here
damping parameter =                     auto                                     (default choice: 'auto'. or any integer for debugging.)
———

### ASTRONOMICAL PARAMETERS ###

stellar spectral model =                blackbody                   (blackbody, HDF5 data set in below file)
path to stellar spectrum file =         ./input/stellar_sample.h5

——pre-tabulated——
planet =                                HD189733b                              (manual, name in the file)
path to planet data file =              ./input/planet_data.dat
——manual input——
surface gravity [cm s^-2] =             3.05
orbital distance [AU] =                 0.0115
radius planet [R_Jup] =                 0.1
radius star [R_Sun] =                   0.2
temperature star [K] =                  5000			
———

############## DANGER ZONE ############## PROCEED UPON OWN RISK ##############

### EXPERIMENTAL ###

debugging feedback =                    no
I2S transition point =                  0.1                                       (transition in w0. only above this value I2S is used.)
clouds =                                no
path to cloud file =                    ./data/clouds/cloudy.h5
total cloud opacity [cm^2 g^-1] =       1e-1
cloud center pressure [10^-6 bar] =     1e0
cloud width (std.dev.) [dex] =          1
number of run-in timesteps =            0
artificial shortw. opacity =            0
varying timestep =                      yes                                    
timestep [s] =                          1e3                                    
use f approximation formula =           no
———

### VULCAN COUPLING ###

VULCAN coupling =                       no
path to individual opacities =          ./input/r3000/
species file =                          ./input/r3000/species.dat
mixing ratios file =                    ../helios_mix.txt
———