configuration_weekly.py

import pathlib

# number of timesteps to be run 
#numberOfTimeSteps=10000 * 52   # long run
numberOfTimeSteps=15000   # test run

# option to fix both the regolith and the vegetation, not typically used
# in normal simulations
fixedStates=False

# option to call the methods that change the geomorphology
# this is typically on
changeGeomorphology=True

# number of realizations generated
# output is written to subdirectories 1, 2, 3, ..
nrOfSamples = 1

# calculation of early warning signals
intervalForStatsCalculated=52
#intervalForStatsCalculated=1

# interval for map reports
#reportInterval=50
reportInterval=1

# for calculation of early warning signals
# does not work anymore
variances=False

# option to define some parameters as realizations from random
# functions, typically False
createRealizations=False

theDurationOfRainstorm=2.0

# option to define the variables to report to disc
# these are defined in the modules, it is typically either 'full' or 'filtering'
setOfVariablesToReport = 'none'
#setOfVariablesToReport = 'some'
#setOfVariablesToReport = 'full'

# option to define if stats are calculated and reported for zones
# across the map
calculateStatsForZones = False

# option to define if stats are calculated for the central
# part of the map (almost the whole map), i.e. one value
calculateStatsAverageOverMap = False

# calculate and store ad hoc report of timeseries as tss file
reportAdHocTimeseries = False


# calculate upstream totals (with accuflux) in subsurfacewateronelayer module
# and interceptionuptomaxstore module
# may be needed for some reports
# and possibly for budget checks (if one needs these)
# normal use: set to False
calculateUpstreamTotals = False

# when True, a particle filtering run is done
# first time users should have this False
filtering = False

# introduce random jumps in biomass and regolith thickness
# for normal runs this should be False
jumpsInRegolithAndBiomass = False

#
# Reporting for the model components
#
if setOfVariablesToReport == 'full':
  interception_report_rasters = ["Vo", "Vi", "Vgf", "Vms"] # reports of totals (Vot) only make sense if calculateUpstreamTotals is True
  infiltration_report_rasters = ["Ii", "Is", "Iks"]
  runoff_report_rasters = ["Rq", "Rqs"]
  subsurface_report_rasters = ["Gs", "Go"]  # reports of totals (Gxt, Got) only make sense if calculateUpstreamTotals is True
  surfacestore_report_rasters = ["Ss", "Sc"]
  rainfalleventsfromgammadistribution_report_rasters = ["Pf"]
  exchange_report_rasters = ["Xrc"]
  soilwashMMF_report_rasters = ["Wde", "Wdm", "Wfl"]
  regolith_report_rasters = ["Ast"]
  bedrockweathering_report_rasters = ["Cwe"]
  evapotranspirationsimple_report_rasters = ["Ep", "Ea"]
  biomassmodifiedmay_report_rasters = ["Xs"]
  baselevel_report_rasters = ["Ll"]
  creep_report_rasters = ["Ds"]
elif setOfVariablesToReport == 'some':
  interception_report_rasters = [] # reports of totals (Vot) only make sense if calculateUpstreamTotals is True
  infiltration_report_rasters = []
  runoff_report_rasters = []
  subsurface_report_rasters = []  # reports of totals (Gxt, Got) only make sense if calculateUpstreamTotals is True
  surfacestore_report_rasters = []
  rainfalleventsfromgammadistribution_report_rasters = []
  soilwashMMF_report_rasters = []
  exchange_report_rasters = []
  soilwashMMF_report_rasters = []
  regolith_report_rasters = ["Ast", "Ade"]
  bedrockweathering_report_rasters = []
  evapotranspirationsimple_report_rasters = []
  biomassmodifiedmay_report_rasters = ["Xs"]
  baselevel_report_rasters = []
  creep_report_rasters = []
elif setOfVariablesToReport == 'none':
  interception_report_rasters = [] # reports of totals (Vot) only make sense if calculateUpstreamTotals is True
  infiltration_report_rasters = []
  runoff_report_rasters = []
  subsurface_report_rasters = []  # reports of totals (Gxt, Got) only make sense if calculateUpstreamTotals is True
  surfacestore_report_rasters = []
  rainfalleventsfromgammadistribution_report_rasters = []
  soilwashMMF_report_rasters = []
  exchange_report_rasters = []
  regolith_report_rasters = []
  bedrockweathering_report_rasters = []
  evapotranspirationsimple_report_rasters = []
  biomassmodifiedmay_report_rasters = []
  baselevel_report_rasters = []
  creep_report_rasters = []