Options¶
Effectively this are (almost 1 2) all capabilities of the SCOPE model.
This is an input structure that controls the workflow.
The values have binary (or tertiary) logic thus equal to 0 or 1 (or 2).
Influence on the output files is highlighted in the corresponding section Output files
Note
Not all combinations can bring to the desired result
Initialized¶
SCOPE.m
: read from setoptions.csv
Rules of input reading¶
simulation
¶
Defines rules of input reading
Switch in SCOPE.m
(multiple)
- 0
individual run(s): specify one value for fixed input parameters, and an equal number (> 1) of values for all parameters that vary between the runs.
- 1
time series (uses text files with meteo input as time series from “./input/dataset X” with files similar to ./input/dataset for_verification specified in the
filenames.csv
- 2
Lookup-Table: specify a number of values in the row of input parameters. All possible combinations of inputs will be used.
Let us illustrate what the difference is in details.
It is possible to specify several values in a row on input_data.csv
. Suppose we have an the following combination of input parameters. Notice, we provide two values for Cab and Cca parameters.
If individual run(s) (options.simulation == 0
) was chosen the given combination will end up in two simulations:
Cab=80, Cca=20
Cab=40, Cca=10
If Lookup-Table (options.simulation == 2
) was chosen the given combination will end up in four simulations:
Cab=80, Cca=20
Cab=80, Cca=10
Cab=40, Cca=20
Cab=40, Cca=10
Variations in input¶
soilspectrum
¶
Calculate soil reflectance or use from a file in ./input/soil_spectra
Switch in SCOPE.m
0
use soil spectrum from the file with soil.spectrumdefault file issoilnew.txt
, can be changed on thefilenames
sheetsoil_file
cellvariable name isrsfile
- 1
simulate soil spectrum with the BSM model (
BSM()
) parameters are fixed in code
soil_heat_method
¶
Method of ground heat flux (G) calculation. In soil_heat_method 0 and 1 soil thermal inertia (GAM) is calculated from inputs.
Switch in SCOPE.m
, select_input()
, ebal()
0
standard calculation of thermal inertia from soil characteristic
- 1
- empirically calibrated formula from soil moisture content
Soil_Inertia1()
inselect_input()
- 2
- as constant fraction (0.35) of soil net radiation
calc_rss_rbs
¶
soil resistance for evaporation from the pore space (rss) and soil boundary layer resistance (rbs)
Switch in select_input()
0
use resistance rss and rbs as provided in inputdata soil
- 1
calculate rss from soil moisture content and correct rbs for LAI
calc_rssrbs()
mSCOPE
¶
Switch in SCOPE.m
0
traditional single layer SCOPE
- 1
multilayer mSCOPE
Variations in output¶
lite
¶
Switch in RTMo()
0
Normal SCOPE execution with [13 x 36 x nlayers] sunlit leaves
1
Lite SCOPE execution with [nlayers x 1] sunlit leaves, sunlit leaf inclinations are not accounted for
calc_planck
¶
Calculate spectrum of thermal radiation with spectral emissivity instead of broadband
Warning
only effective with calc_ebal == 1
Switch in SCOPE.m
, calc_brdf()
0
RTMt_sb()
- broadband brightness temperature is calculated in accordance to Stefan-Boltzman’s equation.
- 1
- Calculation is done per each wavelength thus takes more time than Stefan-Boltzman.
calc_directional
¶
Calculate BRDF and directional temperature for many angles specified in the file: directional.
Warning
only effective with
calc_ebal == 1
Be patient, this takes some time
Switch in SCOPE.m
, calc_brdf()
0
- 1
- struct directional is loaded from the file directional
calc_brdf()
is launched inSCOPE.m
calc_xanthophyllabs
¶
Calculate dynamic xanthopyll absorption (zeaxanthin) for simulating PRI (photochemical reflectance index)
Warning
only effective with
calc_ebal == 1
Switch in SCOPE.m
0
- 1
RTMz()
is launched inSCOPE.m
andcalc_brdf()
(ifcalc_directional
)
calc_vert_profiles
¶
Calculation of vertical profiles (per 60 canopy layers).
Corresponding structure profiles
Switch in SCOPE.m
, RTMo()
and ebal()
0
Profiles are not calculated
- 1
calc_fluor
¶
Calculation of fluorescence
Switch in SCOPE.m
, calc_brdf()
0
No fluorescence
- 1
- total emitted fluorescence is calculated by
SCOPE.m
Fluorescence_model
¶
Fluorescence model
Switch in ebal()
- 0
empirical, with sigmoid for Kn:
biochemical()
(Berry-Van der Tol)- 1
biochemical_MD12()
(von Caemmerer-Magnani)
applTcorr
¶
correct Vcmax and rate constants for temperature
Warning
only effective with Fluorescence_model == 0
i.e. for biochemical()
Switch in ebal()
- 0
- 1
correction in accordance to Q10 rule
MoninObukhov
¶
Switch in ebal()
- 0
do not apply Monin-Obukhov atmospheric stability correction
- 1
apply Monin-Obukhov atmospheric stability correction
For users’ comfort¶
verify
¶
verify the results (compare to saved ‘standard’ output) to test the code for the first time
Switch in SCOPE.m
- 0
- 1
runs
output_verification()
saveCSV
¶
Switch in SCOPE.m
, bin_to_csv()
- 0
leave .bin files in output folder
- 1
convert .bin files to .csv with
bin_to_csv()
, delete .bin files
save_spectral
¶
Save files with full spectrum. May reach huge sizes in long time-series.
Switch in create_output_files_binary()
- 0
do not save
- 1
save