Configuration-Section Options

[eqparam], [rk4input] and [initials]

The options of these sections are always the same and are always required! After the Configuration.importer function processes the .ini-file, these options are distributed to the functions of the model algorithm. For more information about the parameters, for [eqparam] see ModelEquation, for [rk4input] see variable_importer and for [initials] see variable_importer:

[eqparam]
C_ao=70*4.2e6

[rk4input]
number_of_integration=365*10
stepsize_of_integration=60*60*24
spatial_resolution=0
both_hemispheres=True
latitudinal_circle=True
latitudinal_belt=False

eq_condition=False
eq_condition_length=100
eq_condition_amplitude=1e-3

data_readout=1
number_of_externals=0

[initials]
time=0
zmt=273+15
gmt=273+15
initial_temperature_cosine=False
initial_temperature_amplitude=30
initial_temperature_noise=True
initial_temperature_noise_amplitude=5

[func]

The options for a [func] section are all physical functions defined in Functions which describe an Energy flux. There are four classes which contain such functions, Functions.flux_down, Functions.flux_up, Functions.transfer and Functions.forcing.

The functions and examples of their required parameters are listed here after their class.

Important

There may only be one option each of Functions.flux_down, Functions.flux_up and Functions.transfer be used while multiple functions of Functions.forcing might be used.

flux_down Options

There is only one option, Insolation:

[func0]
func=flux_down.insolation
q=1
m=1
dq=0

albedo=albedo.dynamic_bud
albedoread=True
albedoparam=[273.15-5,273.15-15,0.32,0.5,0.75]

noise=False
noiseamp=342*0.03
noisedelay=1
seed=True
seedmanipulation=0

solarinput=True
convfactor=1
timeunit='annualmean'
orbital=False
orbitalyear=0

flux_up Options

Option 1, Budyko clear sky:

[func1]
func=flux_up.budyko_noclouds
activation=True
a=210
b=2.1

Option 2, Budyko cloudy sky:

[func1]
func=flux_up.budyko_clouds
activation=True
a=230.31
b=2.2274
a1=3.0*15.91
b1=0.1*15.91
fc=0.5

Option 3, Stefan-Boltzmann radiation:

[func1]
func=flux_up.planck
activation=True
grey=0.612
sigma=const.sigma

Option 4, Sellers:

[func1]
func=flux_up.sellers
activation=True
m=0.5
sig=const.sigma
gamma=1.9*10**(-15)
k=1

transfer Options

Option 1, Budyko transfer:

[func2]
func=transfer.budyko
beta=3.18
read=True
activated=True

Option 2, Sellers transfer:

[func2]
func=transfer.sellers
readout=True
activated=True
k_wv=10**5
k_h=10**6
k_o=10**2
g=9.81
a=2/100
eps=0.622
p=1000
e0=1700
L=const.Lhvap/1000
rd=const.Rd/1000
dy=1.11*10**6
dp=800
cp=const.cp
dz=2000
l_cover=0.5
radius=const.a
cp_w=4182
dens_w=998
factor_wv=1
factor_air=1
factor_oc=1
factor_kwv=1
factor_kair=1

forcing Options

Important

If you use multiple forcing.random and forcing.predefiend you have to increase the value of the option forcingnumber by 1, this will create an additional space in the output-array and an internal counter of the forcings.

Option 1, Random forcing:

[func3]
func=forcing.random
forcingnumber=0
start=1958
stop=2018
steps=1/365
timeunit='year'
strength=10
frequency='rare'
behaviour='exponential'
lifetime=365
seed=None
sign='negative'

Option 2, Imported predefined forcing:

[func3]
func=forcing.predefined
forcingnumber=0
datapath="../Config/Data/"
name="EVA_800_2000.txt"
delimiter=","
header=1
footer=0
col_time=0
col_forcing=1
timeunit='year'
bp=False
time_start=0
k_output=1
m_output=0
k_input=1
m_input=0

Option 3, Imported 1D predefined forcing:

[func3]
func=forcing.predefined1
forcingnumber=0
datapath="../Config/Data/"
name="Forcingdata1D.csv"
delimiter=","
header=1
footer=0
col_time=0
colrange_forcing=[1,19]
timeunit='year'
bp=False
time_start=0
k_output=1
m_output=0
k_input=1
m_input=0

Option 4, Imported CO2 forcing after Myhre:

[func3]
func=forcing.co2_myhre
A=5.35
C_0=280
CO2_base=280
datapath="../Config/Data/"
name="CO2data.csv"
delimiter=","
header=0
footer=0
col_time=3
col_forcing=8
timeunit='year'
bp=False
time_start=0

Option 5, Imported orbital parameter data:

[func3]
func=forcing.orbital
datapath="../Config/Data/"
name="Orbitaldata.csv"
delimiter=","
header=0
footer=0
col_time=0
col_ecc=1
col_per=2
col_obl=3
timeunit='year'
bp=False
time_start=0
initial={'ecc': 0.017236, 'long_peri': 281.37, 'obliquity': 23.446}
perishift=180

Option 6, Imported total solar irradiance:

[func3]
func=forcing.solar
datapath="../Config/Data/"
name="tsi_SBF_11yr.txt"
delimiter=""
header=4
footer=0
col_time=0
col_forcing=1
timeunit='year'
bp=False
time_start=0
k_output=1
m_output=0
k_input=1
m_input=0

Option 7, Imported AOD forcing:

[func3]
func=forcing.aod
datapath="../Config/Data/"
name="AODdata.csv"
delimiter=","
header=0
footer=0
col_time=0
col_forcing=1
timeunit='year'
bp=False
time_start=0
k_output=1
m_output=0
k_input=1
m_input=0