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