=76.2, respectively
X = A; # start infiltration rate; here, A was taken from the soil table, since parameters SU_PAR01 ... SU_PAR10 are mapped to internal variables A to J, see description below. Also possible: X = 65.2; but then no variation for different soil types would be possible
Y = B*K^C*(L*100)^D; # end infiltration rate
C1 = (100*L)^F; # F = SU_PAR06 C1 will be stored in a new internal variable
C2 = K^G; # G = SU_PAR07 C2 will be stored in a new internal variable
C3 = (O+0.001)^H; # H = SU_PAR08 C3 will be stored in a new internal variable
Z = (O<=0) + (O>0)*(E*C1*C2*C3); # E = SU_PAR05, O = time since last soil tillage, see below
# attention: only starting with V= and following following expressions will be called after internal update of Q. V must be set only after this internal update, but any other expression may be placed herunder for preparation of the V-call. However, they will be called before the internal update of Q, so they should not touch any of the variables needed for EKIN update
V = ((X-Y)*exp(-Z*Q)+Y)*R/60; # potential infiltration, will be limited internally by available precipitation.
}
# Short description of the expression parser and the expression list syntax for method 3
# - Expressions can be defined following algebraic rules:
# - Each line contains a single expression which must be closed with a semi colon.
# - Each assignment (e.g. A = 15) results in creating or updating a value in the internal variable list.
# - A number of values is already defined by WaSiM (as interface from the calling module), and WaSiM expects some other values to be defined after all expressions were called
# The expression parser is based on the source code of the expression parser used in SpeQ Mathematics (http://www.speqmath.com/tutorials/expression_parser_cpp/index.html),
# written by Jos de Jong, 2007. It was adopted to the usage in WaSiM by simplifying the error handling (exceptions are to be handled by WaSiM),
# extracting the variable list as an external class (to be handled by WaSiM) and some other minor technical changes
# Operators (ascending precedence per line, no precedence within a line):
# & | << >> (AND, OR, BITSHIFTLEFT, BITSHIFTRIGHT)
# = <> < > <= >= (EQUAL, UNEQUAL, SMALLER, LARGER, SMALLEREQ, LARGEREQ)
# + - (PLUS, MINUS)
# * / % || (MULTIPLY, DIVIDE, MODULUS, XOR)
# ^ (POW)
# ! (FACTORIAL)
# Functions (must be used with brackets):
# Abs(arg), Exp(arg), Sign(arg), Sqrt(arg), Log(arg), Log10(arg)
# Sin(arg), Cos(arg), Tan(arg), ASin(arg), ACos(arg), ATan(arg)
# Factorial(arg)
# Variables:
# Pi, Euler (not only e, e is a predefined variable used by WaSiM to deliver a value to the expression parser interface)
# you can define your own variables, even with with more than one significant character length, e.g. Inf0 or Help etc.
# there is no distinction between upper and lower case in function names and variables.
# Other:
# Scientific notation supported
#
# ====> what values WaSiM defines forinput (can be used in any expression)
# A to J: values as used in soiltable with names SU_PAR01 to SU_PAR10
# K: grain size distribution Dg, internally calculated after
# double FClay = log004+log2;
# double FSilt = 0.3326 * (log2+log6_3) + 0.3348 *(log6_3+log20) + 0.1704 * (log20+log36) + 0.1622 * (log36+log63) ;
# double FSand = 0.1336 * (log63+log100) + 0.2005 *(log100+log200) + 0.3318 *(log200+log630) + 0.3341 *(log630+log2000);
# double FStones1 = (log2000+log6300);
# double FStones2 = (log6300+log20000);
# double FStones3 = (log20000+log63000);
# double FStones4 = (log63000+log200000);
# double dg = (FClay*dFractionClay + FSilt*dFractionSilt + FSand*dFractionSand + FStones1*dFractionStones1 + FStones2*dFractionStones2 + FStones3*dFractionStones3 + FStones4*dFractionStones4)/2.0;
# with fractions of each grain size class taken from the soil table
# L: fraction of sand
# M: fraction of clay
# N: fraction of silt
# O: t_cult, time since last soil cultivation (in days)
# P: rain intensity in mm/h, taken from precipitation input
# Q: e_kin: accumulated cinetic energy: for all expressions resulting in W, X, Y or Z: result value of the last time time step; for V: value of the actual time step
# R: internal time step in minutes
# ====> What WaSiM expects for output: (ranging from Z downwards, will be used by WaSiM when going ahead)
# Z: silting up disposition SDISP
# Y: end infiltration rate i_inf
# X: start infiltration rate i0
# W: actual cinetic energy
# V: potential infiltration rate inf_pot, depending on energy, siting up disposition, inf_start and inf_infinite
# order of expressions evanulated by WaSiM:
# expressions returning W, X, Y and Z are independently of each other.
# expression V must be called as last call in any case, since WaSiM will update EKIN internally using the energy-result (in W) and V depends on all the other results W to Z
# other expressions for storing intermediate results may be defined at any position in the expression list before the results will be used in another expression
[SurfaceRoutingModel]
0 # 0=ignore this module, 1 = run the module
$time # duration of a time step in minutes
2 # method: 1=MultipleFlowPaths for diverging areas, 2=single flowpaths (nearest direction as given by aspect)
$outpath//qdsr//$grid//.//$code//$year $hour_mean # direct discharge from surface routing module
$outpath//qisr//$grid//.//$code//$year $hour_mean # interflow from surface routing module
$outpath//qbsr//$grid//.//$code//$year $hour_mean # baseflow from surface routing module
$outpath//qgsr//$grid//.//$code//$year $hour_mean # total discharge from surface routing module
$outpath//$surfspeed_grid # grid with actual flow velocity of surface flow in m/s
$Writegrid # writegrid for this grid
$outpath//$surfflux_grid # grid with actual flow amounts of surface flow in m^3/s
$Writegrid # writegrid for this grid
0.001 # maximum wake lenght iteration difference (if Delta_A_nl < this value, iteration for a_NL stops)
40 # maximum number of iterations for a_NL
0.0001 # maximum flow velocity iteration difference (if Delta v is less than this value, iteration stops)
40 # maximum number of iterations for v
30 # shortest sub-time step in seconds
3600 #longest allowed sub time step (even if flow travel times are longer, the time step is subdivided into sub timesteps of this lenght) be careful: tracers are mixed much faster when multiple sub time steps are applied
0.02 # minimum water depth for regarding roughenss of crops in m (shallower sheet flow: only roughness of bare soil will be regarded)
2.0 # ConcentrationFactor takes into account the micro scale concentration of flow pathes, flow will take place on a fraction of the cell only, so the amount flowing per meter width will be multiplied by this factor (1..n)
$readgrids # readgrid code 0 do not read, 1 = read grids
$outpath//sfstsr//$grid//.//$code//$year $hour_mean # statistics for surface storage in mm per sub catchment
[lake_model]
0 # 0=ignore this module, 1 = run the module
2 # method for recalculating DHM,
# 1 = do not change the DHM, it refects already the ground surface of the lakes,
# 2 = use mean_pond_grid to calculate dhm corrections
# max_pond_grid will be used for mapping the cells pond content to a lake during model runs - so the lake level may well rise above the normal surface
0.1 # Albedo_OpenWater (will be used only, when the pond is filled with water when calculating potential evaporation -> otherwise, the normal landuse for this cell is referenced for this parameter)
# 10 # rsc for water (usage as above)
0.4 # z0 for water (usage as above)
# 10.0 # LAI_OpenWater (usage as above)
# 1.0 # VCF_OpenWater (usage as above)
$readgrids # readgrid code 0 do not read, 1 = read grids -->
# if 0, the initial valte for the POND-grid as Volume of Lakes and Reservoirs is set by V0 from the routing description,
# if readgrids=1, no initialization in done (POND-Grid is read in) but the Vakt-Value is set by the various grids
[unsatzon_model]
1 # 0=ignore this module, 1 = run the module
$time # duration of a time step in minutes
3 # method, 1=simple method (will not work anymore from version 7.x), 2 = FDM-Method 3 = FDM-Method with dynamic time step down to 1 secound
2 # controlling interaction with surface water: 0 = no interaction, 1 = exfitration possible 2 = infiltration and exfiltration possible
0 # controlling surface storage in ponds: 0 = no ponds, 1 = using ponds for surface storage (pond depth as standard grid needed -> height of dams oround fields)
0 # controlling artificial drainage: 0 = no artificial drainage 1 = using drainage (drainage depth and horizontal pipe distances as standard grids needed!)
0 # controlling clay layer: 0 = no clay layer, 1 = assuming a clay layer in a depth, specified within a clay-grid (declared as a standard grid)
5e-8 # permeability of the clay layer (is used for the clay layer only)
4 # parameter for the initialization of the gw_level (range between 1..levels (standard: 4))
$outpath//qdra//$grid//.//$code//$year $hour_mean # results drainage discharge in mm per zone
$outpath//gwst//$grid//.//$code//$year $hour_mean # results groundwater depth
$outpath//gwn_//$grid//.//$code//$year $hour_mean # results mean groundwater recharge per zone
$outpath//sb05//$grid//.//$code//$year $hour_mean # results rel. soil moisture within the root zone per zone
$outpath//sb1_//$grid//.//$code//$year $hour_mean # results rel. soil moisture within the unsat. zone (0m..GW table) per zone
$outpath//wurz//$grid//.//$code//$year $hour_mean # results statistic of the root depth per zone
$outpath//infx//$grid//.//$code//$year $hour_mean # results statistic of the infiltration excess
$outpath//pond//$grid//.//$code//$year $hour_mean # results statistic of the ponding water storage content
$outpath//qdir//$grid//.//$code//$year $hour_mean # results statistic of the direct discharge
$outpath//qifl//$grid//.//$code//$year $hour_mean # results statistic of the interflow
$outpath//qbas//$grid//.//$code//$year $hour_mean # results statistic of the baseflow
$outpath//qges//$grid//.//$code//$year $hour_mean # results statistic of the total discharge
$outpath//gwin//$grid//.//$code//$year $hour_mean # statistic of the infiltration from surface water into groundwater (from rivers and lakes)
$outpath//gwex//$grid//.//$code//$year $hour_mean # statistic of the exfiltration from groundwater into surface water (into rivers and lakes)
$outpath//macr//$grid//.//$code//$year $hour_mean # statistic of infiltration into macropores
$outpath//qinf//$grid//.//$code//$year $hour_mean # statistic of total infiltration into the first soil layer
$outpath//$SB_1_grid # grid with actual soil water content for the root zone
$Writegrid # writegrid for this grid
$outpath//$SB_2_grid # grid with actual soil water content for the entire unsaturated zone
$Writegrid # writegrid for this grid
$outpath//$ROOTgrid # grid with root depth
$Writegrid # writegrid for this grid
$outpath//$Thetastack # stack, actual soil water content for all soil levels
$Writestack # Writecode for this stack
$outpath//$hydraulic_heads_stack # stack, contaiing hydraulic heads
$Writestack # Writecode for this stack
$outpath//$geodetic_altitude_stack # stack, containig geodaetic altitudes of the soil levels (lower boudaries)
$Writestack # Writecode for this stack
$outpath//$flowstack # stack, containing the outflows from the soil levels
$Writestack # Writecode for this stack
$outpath//$GWdepthgrid # grid with groudwaterdepth
$Writegrid # writegrid for this grid
$outpath//$GWthetagrid # grid with theta in GWLEVEL
$Writegrid # writegrid for this grid
$outpath//$GWNgrid # grid with groundwater recharge
$Writegrid # writegrid for this grid
$outpath//$GWLEVELgrid # grid with level index of groundwater surface (Index der Schicht)
$Writegrid # writegrid for this grid
$outpath//$QDRAINgrid # grid with the drainage flows
$Writegrid # writegrid for this grid
$outpath//$SATTgrid # grid with code 1=saturation at interval start, 0 no sat.
$Writegrid # writegrid for this grid
$outpath//$INFEXgrid # grid with infiltration excess in mm (surface discharge)
$Writegrid # writegrid for this grid
$outpath//$QDgrid # grid with direct discharge
$Writegrid # writegrid for this grid
$outpath//$QIgrid # grid with Interflow
$Writegrid # writegrid for this grid
$outpath//$QBgrid # grid with baseflow
$Writegrid # write code for baseflow
$outpath//$GWINgrid # grid with infiltration from rivers into the soil (groundwater)
$Writegrid # writegrid for re-infiltration
$outpath//$GWEXgrid # grid with exfiltration (baseflow) from groundwater (is only generated, if groundwater module is active, else baseflow is in QBgrid)
$Writegrid # writegrid for exfiltration
$outpath//$act_pond_grid # grid with content of ponding storge
$Writegrid #
$outpath//$UPRISEgrid # grid with amount of capillary uprise (mm)
$Writegrid #
$outpath//$PERCOLgrid # grid with amount of percolation (mm)
$Writegrid #
$outpath//$MACROINFgrid # grid with amount of infiltration into macropores (mm)
$Writegrid #
$outpath//$irrig_grid # grid with irrigation amount (will be written when irrigation is used, only)
$Writegrid # writegrid for this grid (however: will be written when irrigation is used, only)
36 50 # coordinates of control plot, all theta and qu-values are written to files (qu.dat, theta.dat in the directory, from which the model is started)
$outpath//qbot//$grid//.//$code//$year # name of a file containing the flows between the layers of the control point
$outpath//thet//$grid//.//$code//$year # name of a file containing the soil moisture as theta values of the layers of the control point
$outpath//hhyd//$grid//.//$code//$year # name of a file containing the hydraulic head of the layers of the control point
$outpath//otherdata//$grid//.//$code//$year # name of a file containing some other water balance data of the control point (non layer data)
$outpath//etrd//$grid//.//$code//$year # name of a file containing the hydraulic head of the layers of the control point
$outpath//intd//$grid//.//$code//$year # name of a file containing the hydraulic head of the layers of the control point
1 2 3 4 5 6 7 8 9 10 11 12 13 # range for subbasin codes
3 3 3 2 0.5 0.5 0.5 1 2 6 0.5 2 3 # kelsqd
15 11 11 9 1 1.1 7 6 6 9 1 15 15 # kelsqi
18 9 12 35 80 100 100 25 40 21 60 8 18 # drainage density
0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 # k in qb = Q0 * exp(-k/z) with z = depth to groundwater
0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 # Q0 in the above formula
0.05 0.05 0.05 0.05 0.05 0.2 0.2 0.05 0.05 0.05 0.05 0.15 0.05 # fraction of surface runoff on snow melt
# values for RHB optimized, all others not yet
#3 5 5 3 3 3 3 3 3 6 3 5 3 # kelsqd
#15 25 20 15 8 8 10 10 8 9 8 20 15 # kelsqi
#18 18 18 18 16 25 18 20 13 21 18 18 18 # drainage density
##0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 # k in qb = Q0 * exp(-k/z) with z = depth to groundwater
#0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 # Q0 in the above formula
#0.08 0.05 0.05 0.08 0.08 0.10 0.12 0.12 0.12 0.05 0.12 0.15 0.08 # fraction of surface runoff on snow melt
$readgrids # meanings are extended now! read the follwing comments
$outpath//storage_richards.ftz # if readgrids = 1, then this file contains the contents of the flow travel time zones for interflow and surface flow and for the tracers
300 # minimum dynamic time step in secounds. the smaller this number, the longer the model runs but the results will be more accurate due to a maintained Courant condition
$outpath//step//$grid//.//$code//$year $hour_mean # results statistic of the number of substeps
$outpath//$SUBSTEPSgrid # grid with number of substeps --> a good idea is to use writecode 5x (e.g. 53) to get the average number of substeps per cell for the model run
$Writegrid # for substeps, the areal distribution is of interest for the annual average value. This is code 6 as first digit in 2-digit codes. Or use 5 for the entire model run
[ExternalCoupling]
0 # 0 = no coupling, 1=coupling
$exchngpath//wasim.inf # name of the semaphore file to inform wasim that all grids written by the groundwater model are available now
50 # wait interval in ms between scanning the directory for the new semaphore file (wasim for windows will use a second thread to minimize CPU time, whereas wasim as
# console application will use 100% CPU time while waiting for the output file of the groundwater model. The wait time is then used to minimize disk access
# A follwing version will use a DLL with a memory pointer to the required grid and a flag, which is used by both programs to couple the models.
# But this is music for the future yet...
H # Coupling mode: I=each interval, H=each hour, D=each Day, M=each month, Y=each year
60 # time interval in minutes, the external model uses. This is important to convert changes in groundwater level into fluxes as used by WaSiM
1 # number of following grid names which must be available once the semaphore file was written. Each following row (1..n) will contain a symbolic name
# (grid names from modconst.h) Thus, any grid may be read in, even for other sub models like the boundary conditions as gw_boundary_fix_h_1 for the first aquifer or a changed landuse grid a.s.o.
$exchngpath//gwtable38.grd GWTableExtern 1 0 # the first value is the file name, the second the internal grid name (see English-section in modconst.h ), the third parameter is the fillMissings-parameter (0=no fill, 1=fill with nearest neighbors value), the last ine is the rename(1)/delete(0) parameter
#$exchngpath//bh.grd gw_boundary_fix_h_1 0 0 # the first value is the file name, the second the internal grid name (see English-section in modconst.h ), the third parameter is the fillMissings-parameter (0=no fill, 1=fill with nearest neighbors value), the last ine is the rename(1)/delete(0) parameter
3 # number of grids (each matching one of the following rows) which should be written when the next synchronisation is due
$exchngpath//gwn.grd groundwater_recharge Mean # hier als Beispiel das Grid mit der Grundwasserneubildung
$exchngpath//gwstand.grd groundwater_distance Last # hier als Beispiel das Grid mit der Grundwasserneubildung
$exchngpath//balance.grd Balance Sum # hier als Beispiel das Grid mit der Bilanz aller Wasserinhaltsänderungen durch die Kopplung (sollte 0 sein)
2 # number of subbasin correlated statistics (mean values) which should be written as table (in ASCII-Format) (this is actually limited to directflow and interflow)
$exchngpath//qdir.table direct_discharge Sum # direct flow per subbasin/zone in mm
$exchngpath//difl.table Interflow Sum # interflow per subbasin/zone in mm
$exchngpath//geofim.inf # name of the semaphore file wasim will write after all of the output above was written
geofim # content of the semaphore file
[irrigation]
0 # 0=ignore this module, 1 = run the module
$time # duration of a time step in minutes
$outpath//irgw//$grid//.//$code//$year $hour_mean # statistic of the irrigation water from groundwater
$outpath//irsw//$grid//.//$code//$year $hour_mean # statistic of the irrigation water from surface water
[groundwater_flow]
1 # 0=ignore the module, 1 = run the module
$time # duration of a time step in minutes; doen't change the value unless you have strong reasons to do so!!
1 # solving method: 1=Gauss-Seidel-iteration (using alpha for control wether it is explicite, partly or fully implicite), 2=PCG (not yet implemented
100 # if iterative solving method (1): max.numberof iterations
0.0005 # if iterative solving method (1): max. changes between two iterations
0.0 # Alpha for estimation of central differences 0.5 = Crank-Nicholson Method, 0 = fully explicite, 1 = fully implicite
1.0 # factor for relaxing the iteration if using iterativemethod (successive over[/under] relaxation)
$readgrids # 1=read grids for heads from disk, 0=do not read but initialize with gw-level from unsaturated zone
1 # number of layers
36 50 # coordinates of a control point for all fluxes and for each layer : q0..q4, leakage up and down
$outpath//glog//$grid//.//$code//$year # name of a file containing the flows between of the control point
1 # use Pond Grid -> this enables the model to use the hydraulic head of a pond in addition to the groundwater itself 0=use traditional method without pond (default), 1=use ponds
$outpath//$head1grid # (new) grid for hydraulic heads for layer 1
$Writegrid # writecode for hydraulic heads for layer 1
$outpath//$flowx1grid # (new) grid for fluxes in x direction for layer 1
$Writegrid # writecode for flux-x-grid in layer 1
$outpath//$flowy1grid # (new) grid for fluxes in y direction for layer 1
$Writegrid # writecode for flux-y-grid in layer 1
$outpath//$GWbalance1grid # (new) grid for balance (difference of storage change vs. balance of fluxes -> should be 0 or the amount of in-/outflows by boundary conditions
13 # writecode for balance control grid in layer 1 (should be at least one sum grid per year --> Code = 20 or 23 (if old grids must be read in)
$outpath//$head2grid # (new) grid for hydraulic heads for layer 2
$Writegrid # writecode for hydraulic heads for layer 2
$outpath//$flowx2grid # (new) grid for fluxes in x direction for layer 2
$Writegrid # writecode for flux-x-grid in layer 2
$outpath//$flowy2grid # (new) grid for fluxes in y direction for layer 2
$Writegrid # writecode for flux-y-grid in layer 2
$outpath//$GWbalance2grid # (new) grid for balance (difference of storage change vs. balance of fluxes -> should be 0 or the amount of in-/outflows by boundary conditions
13 # writecode for balance control grid in layer 2 (should be at least one sum grid per year --> Code = 20 or 23 (if old grids must be read in)
$outpath//$head3grid # (new) grid for hydraulic heads for layer 3
$Writegrid # writecode for hydraulic heads for layer 3
$outpath//$flowx3grid # (new) grid for fluxes in x direction for layer 3
$Writegrid # writecode for flux-x-grid in layer 3
$outpath//$flowy3grid # (new) grid for fluxes in y direction for layer 3
$Writegrid # writecode for flux-y-grid in layer 3
$outpath//$GWbalance3grid # (new) grid for balance (difference of storage change vs. balance of fluxes -> should be 0 or the amount of in-/outflows by boundary conditions
13 # writecode for balance control grid in layer 3 (should be at least one sum grid per year --> Code = 20 or 23 (if old grids must be read in)
# this paragraph is not needed for WaSiM-uzr but for the WaSiM-version with the variable saturated area approach (after Topmodel)
[soil_model]
1 # 0=ignore this module, 1 = run the module
$time # duration of a time step in minutes
1 # method, 1 = without slow baseflow, 2 = with slow baseflow (not recommended)
$outpath//$sat_def_grid # (new) saturation deficite-grid (in mm)
$Writegrid # writegrid for this grid
$outpath//$SUZgrid # (new) storage grid for unsat. zone
$Writegrid # writegrid for this grid
$outpath//$SIFgrid # (new) storage grid for interflow storage
$Writegrid # writegrid for this grid
$outpath//$SBiagrid # (new) grid for soil moisture in the inaktive soil storage
$Writegrid # Writegrid for inaktive soil moisture
$outpath//$fcia_grid # (new) grid for plant available field capacity in the inaktiven soil storage
$Writegrid # writegrid for this grid
$outpath//$SSPgrid # (new) grid for the relative fraction of the soil storages, which is in contact with ground water
$Writegrid # writegrid for this grid
$outpath//$QDgrid # (new) grid for surface runoff
$Writegrid # writegrid for this grid
$outpath//$QIgrid # (new) grid for Interflow
$Writegrid # writegrid for this grid
$outpath//$Peakgrid # (new) grid for Peakflow (maximum peakflow for the entire model time)
$outpath//qdir//$grid//.//$code//$year $hour_mean # statistic of the surfeca discharge
$outpath//qifl//$grid//.//$code//$year $hour_mean # statistic of the Interflows
$outpath//qbas//$grid//.//$code//$year $hour_mean # statistic of the base flow
$outpath//qbav//$grid//.//$code//$year $hour_mean # statistic of the slow base flow
$outpath//qges//$grid//.//$code//$year $hour_mean # statistic of the total discharge
$outpath//sb__//$grid//.//$code//$year $hour_mean # soil storage in mm per zone
$outpath//suz_//$grid//.//$code//$year $hour_mean # drainage storage in mm per zone
$outpath//sifl//$grid//.//$code//$year $hour_mean # interflow storage in mm per zone
$outpath//sd__//$grid//.//$code//$year $hour_mean # saturation deficite per zone in mm
1 2 3 4 5 6 7 8 9 10 11 12 13 # Codes der Teilgebiete im Zonengrid
0.025 0.045 0.042 0.022 0.032 0.052 0.027 0.040 0.031 0.015 0.030 0.013 0.025 # Rezessionsparameterter m fuer Saettigungsflaechenmodell in Metern
5.0 5.0 5.0 7.0 3.0 3.0 3.0 3.0 8.0 40.0 3.0 20.0 5.0 # Korrekturfaktor fuer Transmissivitaeten
2.0 2.0 2.0 2.0 2.0 2.0 1.0 2.0 2.0 8.0 2.0 8.0 2.0 # Korrekturfaktor fuer K-Wert (vertikale Versickerung), Modell erwartet k in m/s
6.0 7.0 7.0 6.0 4.0 4.0 6.0 3.0 6.0 6.0 6.0 8.0 6.0 # Speicherrueckgangskonstante Direktabflus ELS in h
25.0 25.0 25.0 20.0 15.0 15.0 18.0 15.0 15.0 0.0 15.0 0.0 25.0 # Saettigungsdefizit, bei dessen Unterschreitung lokaler Interflow gebildet wird
80.0 80.0 80.0 60.0 70.0 40.0 50.0 40.0 60.0 1.0 60.0 1.0 80.0 # Speicherrueckgangskonstante Interflow ELS in h
3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 3600.0 # Rueckgangskonstante verzoegerter Basisabfluss in h
0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 # maximale Tiefenversickerungsrate bei Saettigung in mm/h
0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 # Anfangswert QBB
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # Anfangsfuellung des SUZ-Speichers in n*nFK
0.75 0.75 0.75 0.6 0.45 0.45 0.45 0.4 0.45 0.45 0.3 0.7 0.75 # Anfangssaettigungsdefizit in n*nFK, beeinflusst den ersten Basisabfluss
3.0 3.0 4.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 5.0 3.0 3.0 # Anspringpunkt fuer Makroporenabfluss (in mm/h!, bezogen auf Stundenniederschlag!), alles darueber geht direkt in den Drainspeicher!
0.9 0.9 0.9 0.9 0.9 0.9 1.0 0.9 0.9 0.9 0.9 0.9 0.9 # Reduktionsfaktor fuer Auffuellung von Verdunstungsverlusten aus dem Grundwasser und aus dem Interflowspeicher
0.1 0.1 0.1 0.4 0.4 0.5 0.4 0.5 0.4 0.4 0.35 0.35 0.1 # Anteil an der effektiven Schneeschmelze, der bei geschlossener Schneedecke direkt abfliesst und nicht in den Boden gelangen kann
$readgrids # 1=read grids from disk, else generate internal
$outpath//storage_topmodel.ftz # if readgrids = 1, then this file contains the contents of the flow travel time zones for interflow and surface flow and for the tracers
#
# now the routing descriptions follow
# TG (Teilgebiet) is the code for subbasin, OL (Oberlieger) is the code for a upper subbasin (tributary) , flowing into the actual subbasin TG
# AE and AErel (in TG (...)) are the area of the subbasin (TG) and the relative fraction of the modelled
# compared to the real subbasin area (normally 1.0, but if inflows have to be considered or only parts of a subbasin are modelled, AErel can be smaller)
#
# example:
# TG 9 (AE= 484.0, AErel=0.8125) # TG 11 is external inflow ZL 1:
# from OL 10 (kh=0.4, kv=0.4, Bh=3.5, Bv=20.0, Th=0.5, Mh=25.0, Mv=15.0, I=0.0128, L=10408.3, AE=3.75 )
# and ZL 1 (modus = extern $outpath//abstract11.dat 4 5 , kh=0.4, kv=0.4, Bh=7.0, Bv=50.0, Th=1.0, Mh=25.0, Mv=15.0, I=0.0066, L=10838.5, AE=90.75 )
# and OL 12 (kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=36339.6, AE=81.5 )
# TG 9 (AE= 484.0, AErel=0.6435)
# from OL 10 (kh=0.4, kv=0.4, Bh=3.5, Bv=20.0, Th=0.5, Mh=25.0, Mv=15.0, I=0.0128, L=10408.3, AE=3.75 ) # TG 9 with external inflow from OL 11 and 12
# and ZL 2 (modus = extern $outpath//abstract11.dat 4 5, kh=0.4, kv=0.4, Bh=7.0, Bv=50.0, Th=1.0, Mh=25.0, Mv=15.0, I=0.0066, L=10838.5, AE=90.75 )
# and ZL 1 (modus = extern $outpath//abstract12.dat 4 5, kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=36339.6, AE=81.5 )
#
# kh, kv are storage coefficients for mean channel and flood plain flow; measure in h
# Mh, Mv are Manning roughness for mean channel (Mh) and flood plains (Mv); measure in m^(1/3)/s
# Bh, Bv are width of mean channel resp. floodplains; measure in m
# Th is the water depth in the mean channel, measure in m
# I is the slope angle in flow direction in m/m, L is flow length in m, AE is the area of the sub-basin in km^2
#
# for abstractions: modus = intern or extern (if extern, a file name must follow),
# then 3 numeric parameters: min. threshould, fraction of abstracted water, max. capacity of abstraction channel,
# then the measure (standard is mm/h, other poss.: m3/s or m^3/s; all other measures are interpreted as mm/h)
# example:
# and AL 1 (modus = intern 0.0 1.0 1 m^3/s ) -> there must be an corresponding internal inflow ZL 1 !
# and AL 2 (modus = extern $outpath//abstrac2.dat 0 1 15 m^3/s )
#
# for inflows: like normal routings except that the first parameters are modus = intern / extern,
# if modus = extern, a filename and two integer parameters indicating the start row/start column for the given file must follow
# example:
# from ZL 1 (modus = intern , kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=10, AE=81.5 )
# or: from ZL 1 (modus = extern $outpath//abstrac1.dat 4 5, kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=36339.6, AE=81.5 )
#
# for reservoirs: a filename for the output of reservoir content/tracer concentrations, an initial filling and initial tracer concentrations must be given
# additionally for each reservoir a paragraph containing the abstraction rule has to be found within the control file
# example:
# and SP 1 ( file = $outpath//spv___01.//$code//$year , V0 = 1E8, C0 = 1.0 0.1 3.0 0.0 0.0 0.0 0.0 0.0 0.0 ) # C0: initial concentration for all of max. 9 tracers in the reservoir water
# ..
# [abstraction_rule_reservoir_1]
# 2 # number of following reservoir-content - abstraction points
# 0 20.0 # first value: reservoir content (storage), second value: abstraction(outflow) at this filling level
# 1e8 20.0
#
#
$set $ManningH = 30.0
[routing_model]
1 # 0=ignore this module, 1 = run the module, 2=run the module with observed inflows into the routing channels (from discharge files)
$time # duration of a time step in minutes
5 1200 90 24 # minimum/maximum specific discharge (l/s/km^2), number of log. fractions of the range, splitting of the timeintervall (24= 1 hour-intervalls are splitted into 24 Intervalls each of 2.5 min. duration)
# 5 1200 90 24 means: the flow velocities are calculated befor the model start for a range of 5 to 1200 l/s/km^2 with a logarithmic splitting into 90 flow range intervals. each timestep is also splitted into 24 time sub-intervals
$outpath//qgko//$grid//.//$code//$year $routing_code # name of the statistic file with routed discharges
$inpath//spende.//$year # name of the file with observed discharges (l/s/km^2)
13 # number of following column descripotr (which column in the spec. disch. file corresponds to which subbasin
1 1 # first number: subbasin, second: column index
2 2
3 3
4 4
5 5
6 6
7 7
8 8
9 9
10 10
11 11
12 12
13 0 # this is not required, if there is no column in the file for a subbasin, then this is indicated by an absent entry for this subbasin or by setting column index to 0
48 # timeoffset (for r-square calculation. intervals up to this parameter are not evaluated in r-square calculation. e.g. 12: first 12 intervals are neglected )
TG 2 (AE=210.250, AErel=1.0)
from OL 3 (kh=0.1, kv=0.4, Bh= 9.1, Bv= 36.4, Th= 0.91, Mh= $ManningH , Mv=10.0, I=0.0080, L=10363.9, AE=75.750)
TG 5 (AE=269.750, AErel=1.0)
from OL 7 (kh=0.1, kv=0.4, Bh= 8.2, Bv= 32.7, Th= 0.93, Mh= $ManningH , Mv=10.0, I=0.0129, L=13778.1, AE=76.500)
and OL 6 (kh=0.1, kv=0.4, Bh= 6.4, Bv= 25.8, Th= 0.86, Mh= $ManningH , Mv=10.0, I=0.0178, L=10156.8, AE=66.250)
TG 9 (AE=491.750, AErel=1.0)
from OL 10 (kh=0.1, kv=0.4, Bh= 4.7, Bv= 18.9, Th= 0.55, Mh= $ManningH , Mv=10.0, I=0.0138, L=8828.4, AE=2.500)
and OL 11 (kh=0.1, kv=0.4, Bh=10.3, Bv= 41.3, Th= 1.03, Mh= $ManningH , Mv=10.0, I=0.0085, L=8535.5, AE=110.250)
and OL 12 (kh=0.1, kv=0.4, Bh= 9.1, Bv= 36.3, Th= 0.91, Mh= $ManningH , Mv=10.0, I=0.0099, L=31263.3, AE=77.250)
TG 4 (AE=1086.000, AErel=1.0)
from OL 5 (kh=0.1, kv=0.4, Bh=15.7, Bv= 62.7, Th= 1.57, Mh= $ManningH , Mv=10.0, I=0.0055, L=22071.0, AE=269.750)
and OL 9 (kh=0.1, kv=0.4, Bh=21.2, Bv= 84.7, Th= 2.12, Mh= $ManningH , Mv=10.0, I=0.0036, L=22520.7, AE=491.750)
and OL 8 (kh=0.1, kv=0.4, Bh= 4.9, Bv= 19.5, Th= 0.49, Mh= $ManningH , Mv=10.0, I=0.0097, L=22142.0, AE=15.750)
TG 13 (AE=1592.000, AErel=1.0)
from OL 2 (kh=0.1, kv=0.4, Bh=19.6, Bv= 78.4, Th= 1.96, Mh= $ManningH , Mv=10.0, I=0.0010, L=1000.0, AE=210.250)
and OL 4 (kh=0.1, kv=0.4, Bh=30.2, Bv=120.7, Th= 3.02, Mh= $ManningH , Mv=10.0, I=0.0027, L=30263.3, AE=1086.000)
TG 1 (AE=1697.000, AErel=1.0)
from OL 13 (kh=0.1, kv=0.4, Bh=41.9, Bv=167.4, Th= 4.19, Mh= $ManningH , Mv=10.0, I=0.0010, L=18535.5, AE=1592.000)
# abstration rules are defined this way:
# first row: number of following columns, followed by the julian days for which rules will be established
# the Julian day describes the LAST day, the rule is valid for, so the year doesn't have to begin with 1
# but may begin with 31 instead to indicate, that rule one is valid for the entire January.
# Also, the last JD doesn't have to be 366 - when no other rule follows the actual rule, the last rule
# is valid until the end of the year
# other rows: discharge (m^3/s), followed by the abstraction valid for this discharge (m^3/s)
[abstraction_rule_abstraction_4]
12 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD)
7 0 0 0 0 0 0 0 0 0 0 0 0
8 1 0 0 0 0 0 0 0 0 0 0 1
9 2 1 0 0 0 0 0 0 0 0 1 2
10 3 2 1 0 0 0 0 0 0 1 2 3
11 4 3 2 1 0 0 0 0 1 2 3 4
14 7 6 5 4 3 3 3 3 4 5 6 7
15 7 7 6 5 4 4 4 4 5 6 7 7
16 7 7 7 6 5 5 5 5 6 7 7 7
17 7 7 7 7 6 6 6 6 7 7 7 7
18 7 7 7 7 7 7 7 7 7 7 7 7
27 7 7 7 7 7 7 7 7 7 7 7 7
27 8 7 8 7 8 7 8 7 8 7 8 7
TargetCap = 8 8 8 8 8 8 8 8 8 8 8 8
[abstraction_rule_abstraction_1]
12 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD)
7 0 0 0 0 0 0 0 0 0 0 0 0
8 1 0 0 0 0 0 0 0 0 0 0 1
9 2 1 0 0 0 0 0 0 0 0 1 2
10 3 2 1 0 0 0 0 0 0 1 2 3
11 4 3 2 1 0 0 0 0 1 2 3 4
14 7 6 5 4 3 3 3 3 4 5 6 7
15 7 7 6 5 4 4 4 4 5 6 7 7
16 7 7 7 6 5 5 5 5 6 7 7 7
17 7 7 7 7 6 6 6 6 7 7 7 7
18 7 7 7 7 7 7 7 7 7 7 7 7
27 7 7 7 7 7 7 7 7 7 7 7 7
27 8 7 8 7 8 7 8 7 8 7 8 7
TargetCap = 160 160 160 160 160 160 160 160 160 160 160 160
[abstraction_rule_abstraction_2]
12 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD)
14 0 0 0 0 0 0 0 0 0 0 0 0
16 2 0 0 0 0 0 0 0 0 0 0 2
18 4 2 0 0 0 0 0 0 0 0 2 4
20 6 4 2 0 0 0 0 0 0 2 4 6
22 8 6 4 2 0 0 0 0 2 4 6 8
28 14 12 10 8 6 6 6 6 8 10 12 14
30 14 14 12 10 8 8 8 8 10 12 14 14
32 14 14 14 12 10 10 10 10 12 14 14 14
34 14 14 14 14 12 12 12 12 14 14 14 14
36 14 14 14 14 14 14 14 14 14 14 14 14
54 14 14 14 14 14 14 14 14 14 14 14 14
54 16 16 16 16 16 16 16 16 16 16 16 16
TargetCap = 16 16 16 16 16 16 16 16 16 16 16 16
[abstraction_rule_reservoir_1] #AE ca 1640 km^2 = ca. 41 m3^s bei 800 mm/Jahr Abfluss
6 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD)
0 0 0 0 0 0 0 0 0 0 0 0 0 # = leer
2.7625e08 0 0 0 0 0 0 0 0 0 0 0 0 # = leer -> kein Abfluss bis 17 m Wasserstand
2.925e08 2 2 2 2 2 2 2 2 2 2 2 2 # = 18m -> sehr wenig Abfluss bis 18m Wasserstand
3.0875e08 10 8 5 3 10 8 5 3 10 8 5 3 # = 19m
3.25e08 40 30 20 10 40 30 20 10 40 30 20 10 # = 20m -> mittlerer Abfluss 41m^3/s = ca. 800mm/Jahr bei 1600km^2 EZG
3.43e08 200 150 100 70 200 150 100 70 200 150 100 70 # = 21m -> incl. Flachwasserzone im Norden (7 Zellen mit max 1m Tiefe bis 21 m Seewasserstand)
[abstraction_rule_reservoir_2] #AE ca. 715 km^2, = ca 18 m^3/s bei 800 mm/Jahr Abfluss
6
0 0 # = leer
3.675e07 0 # = 7 m
4.200e07 2.6 # = 8 m
4.725e07 20 # = 9 m
5.250e07 50 # = 10 m #-> mittlerer Abfluss 18 m^3/S
5.8125e07 120 # = 11 m #incl. 3 Zellen Flachwasserzone im Nordwesten (0.5m)
[abstraction_rule_reservoir_3] #AE ca. 556 km^2, = ca 14 m^3/s bei 800 mm/Jahr Abfluss
6
0 0
1.000e07 0 # = 4 m
1.250e07 2.4 # = 5 m
1.500e07 12 # = 6 m
1.750e07 40 #F = 7 m Normalwasserstand
2.0125e07 140 # = 8 m incl. 1 Zelle 0.5 m Flachwasserzone im Nordosten
[abstraction_rule_reservoir_4] #AE ca. 111 km^2, = ca 2.81 m^3/s bei 800 mm/Jahr Abfluss
6
0 0
8.750e05 0 # = 3.5 m
1.000e06 0.1 # = 4 m
1.125e06 2.8 # = 4.5 m
1.250e06 8 # = 5m Normalwasserstand
1.375e06 40 # = 5.5 m
[abstraction_rule_reservoir_5] # AE ca. 110 km^2, = ca 2.8 m^3/s bei 800 mm/Jahr Abfluss
6
0 0
2.100e07 0 # = 12 m
2.275e07 0.1 # = 13 m
2.450e07 3.8 # = 14 m
2.625e07 12 # = 15 m Normalwasserstand
2.800e07 50 # = 16 m
#[abstraction_rule_reservoir_1]
#9 # number of points, reading now as x-y-values
#1e4 1.0
#1e5 1.0
#1.0001e5 10.0
#1e6 10.0
#1.0001e6 50.0
#1e7 50.0
#1.0001e7 100.0
#1e8 100.0
#1.0001e8 500.0
# the following section defines combinations of single landuse types to combinations of them.
# e.g. a landuse type deciduous forest may contain of oaks, bushes, and herbs, so each of those three components
# must be parameterised in the traditional landuse table. Example: oaks = code 1, bushes = code 2, herbs = code 3
# here, the combination of oaks, bushes and herbs will be parameterised like: 1 deciduous_forest { layers = 1, 2, 3;}
# The VCF (vegetation covered fraction) of each landuse will define the amount of water and radiation (except diffuse
# radiation which will go through the canopy layer) reaching the next layer. The uppermost layer must be listed first,
# the next layer follows then a.s.o.
# All multilayer-landuses must have an equal number of layers. Missing layers can be filled up from the end of the
# list using landuse code 9999, e.g. grassland would be defined in a 3-layer configuration by "2 grass {layers = 4, 9999, 9999;}
# When the multilayer_landuse table is used, the codes of the LANDUSE-Grid are referring no longer to the landuse_table
# anymore but to the multilayer_table following. The codes in the old landuse table are reffering to the entries in the
# multilayer_landuse table
[multilayer_landuse]
10 # count of multilayer landuses
1 water { Landuse_Layers = 1 ; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
2 settlements { Landuse_Layers = 2; k_extinct = 0.3; LAI_scale = 20;}
3 pine_forest { Landuse_Layers = 3; k_extinct = 0.3; LAI_scale = 20;}
4 decidous_forest { Landuse_Layers = 4; k_extinct = 0.3; LAI_scale = 20;}
5 mixed_forest { Landuse_Layers = 5; k_extinct = 0.3; LAI_scale = 20;}
6 agriculture { Landuse_Layers = 6; k_extinct = 0.3; LAI_scale = 20;}
7 grass_variable { Landuse_Layers = 7; k_extinct = 0.3; LAI_scale = 20;}
8 bushes { Landuse_Layers = 8; k_extinct = 0.3; LAI_scale = 20;}
15 rock { Landuse_Layers = 15; k_extinct = 0.3; LAI_scale = 20;}
19 horticulture { Landuse_Layers = 19; k_extinct = 0.3; LAI_scale = 20;}
1 water { Landuse_Layers = 1, 1, 1; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
2 settlements { Landuse_Layers = 2, 2, 2; k_extinct = 0.3; LAI_scale = 20;}
3 pine_forest { Landuse_Layers = 3, 3, 3; k_extinct = 0.3; LAI_scale = 20;}
4 decidous_forest { Landuse_Layers = 4, 4, 4; k_extinct = 0.3; LAI_scale = 20;}
5 mixed_forest { Landuse_Layers = 5, 5, 5; k_extinct = 0.3; LAI_scale = 20;}
6 agriculture { Landuse_Layers = 6, 6, 6; k_extinct = 0.3; LAI_scale = 20;}
7 grass_variable { Landuse_Layers = 7, 7, 7; k_extinct = 0.3; LAI_scale = 20;}
8 bushes { Landuse_Layers = 8, 8, 8; k_extinct = 0.3; LAI_scale = 20;}
15 rock { Landuse_Layers = 15,15,15; k_extinct = 0.3; LAI_scale = 20;}
19 horticulture { Landuse_Layers = 19,19,19; k_extinct = 0.3; LAI_scale = 20;}
1 water { Landuse_Layers = 1, -9999, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
2 settlements { Landuse_Layers = 2, -9999, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
3 pine_forest { Landuse_Layers = 3, 8, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
4 decidous_forest { Landuse_Layers = 4, 8, 7; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
5 mixed_forest { Landuse_Layers = 5, 8, 7; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
6 agriculture { Landuse_Layers = 6, -9999, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
7 grass_variable { Landuse_Layers = 7, -9999, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
8 bushes { Landuse_Layers = 8, 7, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
15 rock { Landuse_Layers = 15, -9999, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
19 horticulture { Landuse_Layers = 19, -9999, -9999; # references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
k_extinct = 0.3; # extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
LAI_scale = 20; # Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
}
# declaring some common variables for vegetation period dependent grid-writing
# default (if not used in land use table at all) is JDVegReset = 1 and JDVegWrite = 365
$set $JDVegReset = 3
$set $JDVegWrite = 365
[landuse_table]
10 # number of following land use codes, per row one use
#Co name of the albe- surface resistances rsc as monthly values julian day for LAI (eff. veget. height) Veg.covering root depth [m] Param. root theta-value for beginning
#de Landuse type do 1 2 3 4 5 6 7 8 9 10 11 12 the param.-sets 1 2 3 4 z01 2 3 4 1 2 3 4 1 2 3 4 distribution. etp-reduction
#-- -------------- ----- ----------------------------------------------- ---------------- ------------- --------------- --------------- ------------------ ------------- ------------
1 water { method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter. The count of days doesn't matter.
Albedo = 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05; # Albedo (snow free)
rsc = 20 20 20 20 20 20 20 20 20 20 20 20; # leaf surface resistance in s/m
rs_interception = 20 20 20 20 20 20 20 20 20 20 20 20; # INTERCEPTION surface resistance in s/m
LAI = 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0; # Leaf Area Index (1/1)
Z0 = 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01; # Roughness length in m
VCF = 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
2 settlements { method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.2; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter. The count of days doesn't matter.
Albedo = 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 ; # Albedo (snow free)
rsc = 100 100 100 100 100 100 100 100 100 100 100 100; # leaf surface resistance in s/m
rs_interception = 100 100 100 100 100 100 100 100 100 100 100 100; # INTERCEPTION surface resistance in s/m
LAI = 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0; # Leaf Area Index (1/1)
Z0 = 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0; # Roughness length in m
VCF = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
3 pine_forest { method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.4; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
Albedo = 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12; # Albedo (snow free)
rsc = 80 80 75 65 55 55 55 55 55 75 80 80; # leaf surface resistance in s/m
# rs_interception = 21 21 21 21 21 21 21 21 21 21 21 21 ; # INTERCEPTION surface resistance in s/m
rs_interception = 80 80 75 65 55 55 55 55 55 75 80 80; # INTERCEPTION surface resistance in s/m
rs_evaporation = 80 80 75 65 55 55 55 55 55 75 80 80; # SOIL surface resistance in s/m (for evaporation only)
LAI = 8.0 8.0 8.0 9.0 12.0 12.0 12.0 12.0 11.0 10.0 8.0 8.0; # Leaf Area Index (1/1)
Z0 = 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0; # Roughness length in m
VCF = 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2 1.2; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
#4 decidous_forest { method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
# RootDistr = 1.0; # parameter for root density distribution
# TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
# LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
# HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
# IntercepCap = 0.3; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
# JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
# Albedo = 0.17 0.17 0.17 0.17 0.17 0.17 0.17 0.17 0.17 0.17 0.17 0.17; # Albedo (snow free)
# rsc = 100 100 95 75 65 65 65 65 65 85 100 100; # leaf surface resistance in s/m
## rs_interception = 22 22 22 22 22 22 22 22 22 22 22 22; # INTERCEPTION surface resistance in s/m
# rs_interception = 100 100 95 75 65 65 65 65 65 85 100 100; # INTERCEPTION surface resistance in s/m
# rs_evaporation = 220 220 220 220 220 220 220 220 220 220 220 220; # SOIL surface resistance in s/m (for evaporation only)
# LAI = 0.5 0.5 0.5 3 8 8 8 8 8 3 0.5 0.5; # Leaf Area Index (1/1)
# Z0 = 0.3 0.3 0.3 3.0 8.00 10.0 10.0 10.0 10.0 3.0 0.5 0.3; # Roughness length in m
# VCF = 0.7 0.7 0.7 0.8 0.95 0.95 0.95 0.95 0.9 0.8 0.7 0.7; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
# RootDepth = 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4 1.4; # Root depth in m
# AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
# }
# for DynamicPhenology_1, one julian day in spring must be named "-1" which means the TStart-day. The next julian day will be interpreted as Delta to the Tstart-day
# the "Tstart"-day will be the day, the F* calculation starts from (i.e. end of Dormation period), the Delta1 will be the end of the leaves unfolding-time.
# all other julian days before and after this two days will still be handled in the usual way (as in Variable Day Count)
# if a julian day "moves" into the time-span of the two special days, it will be ignored.
# Tstart and the next julian day will not be altitude dependent, the other days will.
# There will be no interpolation of Parameters between the Julian day before Tstart and Tstart (because we don't know Tstart in advance)
#
# to use the DynamicPhenology_1, F*, t1_dorm and TBf must be given as Parameters
# Also, a number of grids must be hold in memory (and they will be modelled after each time step), like:
# - actual sum of forcing units (per layer) (also called F in the equations)
# - TStart as julian day (per layer), if Sum(F)>F* --> Phase was started already (otherwise: -9999)
# Ablauf: in InterpolateTVarParams muss festgestellt werden, ob der aktuelle Tag < Tstart ist.
# Wenn ja: Alte Werte weiterverwenden (keine Interpolation)
# Wenn nein: prüfen, ob TStart-Grid schon einen Wert hat.
# TStart schon bekannt (und auch Delta): mit diesen Werten weiterrechnen, falls Julian Day innerhalb TStart bis TStart+Delta liegt,
# sonst ebenfalls normal weitermachen mit den festgelegten Werten
# TStart nicht bekannt: Forcing units berechnen und im F-Grid aufsummieren (nach Formel 12b)
# -> weiter die alten Werte verwenden von vor Tstart
# Für die Berechnung der Forcing units wird die Temperatur benötigt sowie die fixen Parameter F* und T_Bf
#
# Für die interne Abarbeitung:
# - beim EINlesen der Tabelle wird für den Julian-Day, an welchem "-1" steht, die Nummer der SPalte wird als Tstart gespeichert
# (neues Member-Feld der LU-Tabelle. Der Zugriff auf den Julian Day erfolgt dann zwischen Tstart und Delta indirekt über das Tstart-Array
# - der nächste Wert wird als Delta ebenfalls in neuem Member gespeichert -> nachdem Tstart feststeht, wird er intern durch Tstart+Delta ersetzt und
# wird dann normal wie alle anderen Julian Days behandelt
# - beim Beginn eines neuen Jahres wird (bzw. bei Erreichen des Julian days JDReset_TStart) das Tstart-Grid -1 gesetzt, FORC wird auf 0 gesetzt.
#
# internal used: TStart as julian day for t_start
# Tstart_col: welcher der x julian days ist eigentlich tstart (gewöhnlich der zweite, muss aber nicht sein)
#
4 decidous_forest { method = DynamicPhenology_1; # valid methods: "VariableDayCount" with variable number of fix points, DynamicPhenology_1 which estimates the begin of the diurnal cycle in spring dependent on daily temperature sums only (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.3; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
StressFactorDynPhen = 1.5; # optional: specifying the maximum scaling factor for Forcing-Values dependent on soil moisture. Range 0..+infinity, use values between 0.25 and 1 to reduce growth for dry soils and values between 1 and 3 to enforce growth under drying soil conditions
F* = 175; # "Temperatursum" which must be exceeded for starting the phenological cycle (unfolding leaves)
DP1_t1_dorm = 60; # starting day (julian day number), forcing units will be summed up after this day of year
DP1_T_Bf = 0; # threshold temperatur for a positive forcing unit after Model 12b (thermal time model)
JDReset_TStart = 1; # Julian Day when TStart is reset to -1 and Forcing untis are reset to 0 for a new vegetation period
maxStartJDforDP1 = 150; # latest start day for the model run to use DynamicPhenology_1. If start date is after this date, then TStart is set to maxStartJDforDP1 minus the delta of the next column (e.g. 150 - 18 = 132), so we assume that this start date meets a fully developed vegetation. If start day is even after DP2_t0_dorm, then the next year will use DP1 only
StartVegetationPeriodForBalance = 2 ; # the sampling point in the following JD-Table when the vegetation period starts, default = 0 (start of model run)
StopVegetationPeriodForBalance = 6 ; # the sampling point in the following JD-Table when the vegetation period ends, default = n+1 (end of model run)
JDVegetationResetForBalance = $JDVegReset ; # Julian day, when vegetetaion start and vegetation stop grids are re-initialized to -1 (northern hemisphere: usually day 1)
JDVegetationWriteForBalance = $JDVegWrite ; # Julian day, when vegetetaion period dependent grids should be written (usually just before JDVegetationResetForBalance, e.g. 365). Attention: this Value should be identical for all land uses, since grids cannot be written for specific land uses only
JulDays = 1 -1 +17 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
Albedo = 0.17 0.17 0.17 0.17 0.17 0.17 0.17; # Albedo (snow free)
rsc = 100 100 65 65 85 100 100; # leaf surface resistance in s/m
# rs_interception = 22 22 22 22 22 22 22; # INTERCEPTION surface resistance in s/m
rs_interception = 100 100 65 65 85 100 100; # INTERCEPTION surface resistance in s/m
rs_evaporation = 100 100 65 65 85 100 100; # SOIL surface resistance in s/m (for evaporation only)
LAI = 0.5 0.5 8 8 3 0.5 0.5; # Leaf Area Index (1/1)
Z0 = 0.3 0.3 8.00 10.0 3.0 0.5 0.3; # Roughness length in m
VCF = 0.7 0.7 0.95 0.9 0.8 0.7 0.7; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 1.4 1.4 1.4 1.4 1.4 1.4 1.4; # Root depth in m
AltDep = 0.0 0.0 0.0 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben). Parameter beziehen sich auf 400m.ü.NN
}
#5 mixed_forest { method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points, DynamicPhenology_1 which estimates the begin of the diurnal cycle in spring dependent on daily temperature sums only (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
# RootDistr = 1.0; # parameter for root density distribution
# TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
# LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
# HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
# IntercepCap = 0.35; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
# JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
# Albedo = 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15 0.15; # Albedo (snow free)
# rsc = 90 90 85 70 60 60 60 60 60 80 90 90; # leaf surface resistance in s/m
## rs_interception = 23 23 23 23 23 23 23 23 23 23 23 23; # INTERCEPTION surface resistance in s/m
# rs_interception = 90 90 85 70 60 60 60 60 60 80 90 90; # INTERCEPTION surface resistance in s/m
# rs_evaporation = 230 230 230 230 230 230 230 230 230 230 230 230; # SOIL surface resistance in s/m (for evaporation only)
# LAI = 2 2 2 4 8 10 10 10 8 5 3 3; # Leaf Area Index (1/1)
# Z0 = 3.0 3.0 3.0 5.0 8.0 10.0 10.0 10.0 9.0 5.0 3.0 3.0; # Roughness length in m
# VCF = 0.8 0.8 0.8 0.9 0.92 0.92 0.92 0.92 0.9 0.8 0.8 0.8; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
# RootDepth = 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3 1.3; # Root depth in m
# AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
# }
#
# valid methods:
# VariableDayCount with variable number of fix points, the following methods are extension to VriableDayCounts: some key days will be estimated dynamically
# DynamicPhenology_1 which estimates the begin of the diurnal cycle in spring dependent on daily temperature sums (Thermal Time Model 12b),
# DynamicPhenology_2 uses the Sequential Model 2 (24b).
# DynamicPhenology_3 uses the thermal time model for multiple/subsequent sample phases
# other methods will follow --> old method: if the table is structured like the old ones, they are still valid
5 mixed_forest { method = DynamicPhenology_3;
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.35; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
StressFactorDynPhen = 1.5; # optional: specifying the maximum scaling factor for Forcing-Values dependent on soil moisture. Range 0..+infinity, use values between 0.25 and 1 to reduce growth for dry soils and values between 1 and 3 to enforce growth under drying soil conditions
F* = 175.2; # used for DynamicPhenology_1 and _2!: "Temperatursum" which must be exceeded for starting the phenological cycle (unfolding leaves) (if the model period starts between t0_dorm and t1_dorm, then F* will not calculated by sequential model (24b) but by thermal time model (12b)
DP1_t1_dorm = 60; # used for DynamicPhenology_1 and _2!: starting day (julian day number), forcing units will be summed up after this day of year until F* is reached
DP1_T_Bf = 0; # used for DynamicPhenology_1 and _2!: threshold temperatur for a positive forcing unit after Model 12b (thermal time model)
DP2_t0_dorm = 244; # used for DynamicPhenology_2 only: starting day (julian day number), chilling units will be summed up after this day of year until t1_dorm_DP2 is reached
DP2_t1_dorm = 110; # used for DynamicPhenology_2 only: starting day (julian day number), forcing units will be summed up after this day of year
DP2_T_Bf = 0; # used for DynamicPhenology_2 only: threshold temperatur for a positive forcing unit after Model 24b (sequential model 2)
DP2_T_Bc = 11.1; # used for DynamicPhenology_2 only: threshold temperatur for a chilling unit after Model 24b (sequential model 2)
DP2_Par_a = 303.2; # used for DynamicPhenology_2 only: Parameter a in F*=a*exp(bC*) after Model 24b (sequential model 2)
DP2_Par_b = -0.019; # used for DynamicPhenology_2 only: Parameter b in F*=a*exp(bC*) after Model 24b (sequential model 2)
DP2_Offset_1 = -3.4; # used for DynamicPhenology_2 only: value for z1 in R_c(T_i)=(T_i-z1)/(T_Bc-z1) when z1 < T_i < T_Bc
DP2_Offset_2 = 10.4; # used for DynamicPhenology_2 only: value for z2 in R_c(T_i)=(T_i-z2)/(T_Bc-z2) when T_Bc < T_i < z2
JDReset_TStart = 1; # used for DynamicPhenology_1 and _2!: Julian Day when TStart is reset to -1 and Forcing untis are reset to 0 for a new vegetation period
maxStartJDforDP1 = 150; # latest start day for the model run to use DynamicPhenology_1. If start date is after this date, then TStart is set to maxStartJDforDP1 minus the delta of the next column (e.g. 150 - 18 = 132), so we assume that this start date meets a fully developed vegetation. If start day is even after DP2_t0_dorm, then the next year will use DP1 only
StartVegetationPeriodForBalance = 2 ; # the sampling point in the following JD-Table when the vegetation period starts
StopVegetationPeriodForBalance = 6 ; # the sampling point in the following JD-Table when the vegetation period ends
JDVegetationResetForBalance = $JDVegReset ; # Julian day, when vegetetaion start and vegetation stop grids are re-initialized to -1 (northern hemisphere: usually day 1)
JDVegetationWriteForBalance = $JDVegWrite ; # Julian day, when vegetetaion period dependent grids should be written (usually just before JDVegetationResetForBalance, e.g. 365). Attention: this Value should be identical for all land uses, since grids cannot be written for specific land uses only
# dynphen_2 JulDays = 1 -1 +10 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
(max) JulDays = 1 120 150 258 288 319 366 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! For DynamicPhenology_3 these days mark the latest allowed day (when ForcingThreshold was not stepped over, the corresponding julian day will be taken automatically
ForcingThreshold = -1 100 455 2300 -1 -1 -1 ; # Forcing units as Rf=T-DP1_T_Bf, summed up starting from DP1_t1_dorm (not using the functions for model 12b or 24b, pure thermal time model after model 11 or 12a!)
Albedo = 0.15 0.15 0.15 0.15 0.15 0.15 0.15; # Albedo (snow free)
rsc = 90 60 60 60 80 90 90; # leaf surface resistance in s/m
rs_interception = 23 23 23 23 23 23 23; # INTERCEPTION surface resistance in s/m
rs_interception = 120 120 120 80 80 80 120; # INTERCEPTION surface resistance in s/m
rs_evaporation = 90 60 60 60 80 90 90; # SOIL surface resistance in s/m (for evaporation only)
LAI = 3 3 10 8 5 3 3; # Leaf Area Index (1/1)
Z0 = 3.0 8.0 10.0 9.0 5.0 3.0 3.0; # Roughness length in m
VCF = 0.8 0.92 0.92 0.9 0.8 0.8 0.8; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 1.3 1.3 1.3 1.3 1.3 1.3 1.3; # Root depth in m
AltDep = 0.0 0.0 0.0 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
6 agriculture {method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.3; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
StartVegetationPeriodForBalance = 3 ; # the sampling point in the following JD-Table when the vegetation period starts, default = 0 (start of model run)
StopVegetationPeriodForBalance = 10 ; # the sampling point in the following JD-Table when the vegetation period ends, default = n+1 (end of model run)
JDVegetationResetForBalance = $JDVegReset ; # Julian day, when vegetetaion start and vegetation stop grids are re-initialized to -1 (northern hemisphere: usually day 1)
JDVegetationWriteForBalance = $JDVegWrite ; # Julian day, when vegetetaion period dependent grids should be written (usually just before JDVegetationResetForBalance, e.g. 365). Attention: this Value should be identical for all land uses, since grids cannot be written for specific land uses only
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
Albedo = 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25; # Albedo (snow free)
rsc = 80 80 75 75 65 65 65 65 65 75 90 90; # leaf surface resistance in s/m
# rs_interception = 24 24 24 24 24 24 24 24 24 24 24 24 ; # INTERCEPTION surface resistance in s/m
rs_interception = 80 80 75 75 65 65 65 65 65 75 90 90; # INTERCEPTION surface resistance in s/m
rs_evaporation = 240 240 240 240 240 240 240 240 240 240 240 240; # SOIL surface resistance in s/m (for evaporation only)
LAI = 1 1 2 3 4 5 5 4 3 2 1 1; # Leaf Area Index (1/1)
Z0 = 0.05 0.05 0.1 0.25 0.5 0.5 0.5 0.4 0.3 0.15 0.05 0.05; # Roughness length in m
VCF = 0.3 0.3 0.3 0.5 0.8 0.8 0.8 0.7 0.6 0.3 0.3 0.3 ; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 0.05 0.05 0.1 0.2 0.3 0.4 0.4 0.4 0.3 0.2 0.05 0.05; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
7 grass_variable {method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.2; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
Albedo = 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25; # Albedo (snow free)
rsc = 90 90 75 65 50 55 55 55 60 70 90 90; # leaf surface resistance in s/m
# rs_interception = 25 25 25 25 25 25 25 25 25 25 25 25; # INTERCEPTION surface resistance in s/m
rs_interception = 90 90 75 65 50 55 55 55 60 70 90 90; # INTERCEPTION surface resistance in s/m
rs_evaporation = 250 250 250 250 250 250 250 250 250 250 250 250; # SOIL surface resistance in s/m (for evaporation only)
LAI = 2 2 2 2 3 4 4 4 4 2 2 2; # Leaf Area Index (1/1)
Z0 = 0.15 0.15 0.15 0.15 0.3 0.4 0.36 0.33 0.3 0.15 0.15 0.15; # Roughness length in m
VCF = 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
8 bushes {method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 0.5; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.35; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
Albedo = 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20; # Albedo (snow free)
rsc = 80 80 70 70 50 50 50 55 55 70 70 80; # leaf surface resistance in s/m
# rs_interception = 26 26 26 26 26 26 26 26 26 26 26 26; # INTERCEPTION surface resistance in s/m
rs_interception = 80 80 70 70 50 50 50 55 55 70 70 80; # INTERCEPTION surface resistance in s/m
rs_evaporation = 260 260 260 260 260 260 260 260 260 260 260 260; # SOIL surface resistance in s/m (for evaporation only)
LAI = 3.0 3.0 3.0 4.0 5.0 5.0 5.0 5.0 5.0 4.0 3.0 3.0; # Leaf Area Index (1/1)
Z0 = 1.5 1.5 1.5 2.0 2.5 2.5 2.5 2.5 2.0 2.0 1.5 1.5; # Roughness length in m
VCF = 0.9 0.9 0.9 0.9 0.95 0.95 0.95 0.95 0.95 0.9 0.9 0.9; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
15 rock {method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 0.5; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.1; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
Albedo = 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12; # Albedo (snow free)
rsc = 250 250 250 250 250 250 250 250 250 250 250 250; # leaf surface resistance in s/m
# rs_interception = 27 27 27 27 27 27 27 27 27 27 27 27; # INTERCEPTION surface resistance in s/m
rs_interception = 250 250 250 250 250 250 250 250 250 250 250 250; # INTERCEPTION surface resistance in s/m
rs_evaporation = 400 400 400 400 400 400 400 400 400 400 400 400; # SOIL surface resistance in s/m (for evaporation only)
LAI = 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0; # Leaf Area Index (1/1)
Z0 = 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05; # Roughness length in m
VCF = 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
19 horticulture {method = VariableDayCount; # valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
SoilTillage = 90 240; # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
RootDistr = 1.0; # parameter for root density distribution
TReduWet = 0.95; # relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
LimitReduWet = 0.5; # minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
HReduDry = 3.45; # hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
IntercepCap = 0.35; # optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
Albedo = 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25; # Albedo (snow free)
rsc = 100 100 90 70 60 60 60 60 60 80 100 100; # leaf surface resistance in s/m
# rs_interception = 28 28 28 28 28 28 28 28 28 28 28 28; # INTERCEPTION surface resistance in s/m
rs_interception = 100 100 90 70 60 60 60 60 60 80 100 100; # INTERCEPTION surface resistance in s/m
rs_evaporation = 280 280 280 280 280 280 280 280 280 280 280 280; # SOIL surface resistance in s/m (for evaporation only)
LAI = 0.5 0.5 0.5 2 4 5 5 5 4 3 0.5 0.5; # Leaf Area Index (1/1)
Z0 = 0.4 0.40 0.40 1.00 2.5 3.0 3.0 3.0 3.0 1.0 0.4 0.4; # Roughness length in m
VCF = 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75; # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
RootDepth = 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8; # Root depth in m
AltDep = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
}
#original landuse table up to WaSiM 6.4
1 water 0.05 20 20 20 20 20 20 20 20 20 20 20 20 110 150 250 280 1 1 1 1 .01 .01 .01 .01 .1 .1 .1 .1 0.01 0.01 0.01 0.01 1.0 3.45
2 settlements 0.1 100 100 100 100 100 100 100 100 100 100 100 100 110 150 250 280 1 1 1 1 10 10 10 10 .5 .5 .5 .5 0.2 0.2 0.2 0.2 1.0 3.45
3 pine_forest 0.12 80 80 75 65 55 55 55 55 55 75 80 80 110 150 250 280 8 12 12 8 10 10 10 10 .9 .9 .9 .9 1.2 1.2 1.2 1.2 1.0 3.45
4 decidous_forest 0.17 100 100 95 75 65 65 65 65 65 85 100 100 110 150 250 280 .5 8 8 .5 .3 10 10 .3 .7 .95 .95 .7 1.4 1.4 1.4 1.4 1.0 3.45
5 mixed_forest 0.15 90 90 85 70 60 60 60 60 60 80 90 90 110 150 250 280 2 10 10 2 3 10 10 3 .8 .92 .92 .8 1.3 1.3 1.3 1.3 1.0 3.45
6 agriculture 0.25 80 80 75 75 65 65 65 65 65 75 90 90 110 150 250 255 1 5 3 1 .05 .5 .2 .05 .3 .8 .7 .3 0.05 0.4 0.3 0.05 1.0 3.45
7 grass 0.25 90 90 75 65 50 55 55 55 60 70 90 90 110 150 250 280 2 4 4 2 .15 .4 .3 .15 .95 .95 .95 .95 0.4 0.4 0.4 0.4 1.0 3.45
8 bushes 0.20 80 80 70 70 50 50 50 55 55 70 80 80 110 150 250 280 3 5 5 3 1.5 2.5 2.5 1.5 .9 .95 .95 .9 0.5 0.5 0.5 0.5 1.0 3.45
15 rock 0.12 250 250 250 250 250 250 250 250 250 250 250 250 110 150 250 280 1 1 1 1 .05 .05 .05 .05 .8 .8 .8 .8 0.1 0.1 0.1 0.1 1.0 0.5
19 horticulture 0.25 100 100 90 70 60 60 60 60 60 80 100 100 110 150 250 280 .5 5 5 .5 .4 3 3 .4 .75 .75 .75 .75 0.8 0.8 0.8 0.8 1.0 3.45
#[soil_table]
#14 # number of following entries
#Code name FC(Vol.%) mSB(Vol.%) ksat(m/s) suction. parameter #k theta-values rel.k-values for the given thetas #h theta-values for the following h-values h-values for the left Thetas [m] #+thickns. maxratio for each layer (downwards) one entry for k_sat [m/s] has to follow
# [mm] 1=par 2=table of layers ko_rel/ku_rel
#---- ---------------- ---------- --------- -------- -------- ------------- ----- --------------------------------------------------------- ------------------------------------------------------------- ----- --------------------------------------------------------------------------------- ------------------------------------------------- --------- ------------- -----------------------------------------------------------------------------------------------------------------------------------------------------------------------
1 Sand_(S) 6.21 38.5 8.25E-5 385 2 10 .43 .4286 .3537 .2143 .1071 .0650 .0513 .0470 .0455 .0450 1 .9221 .2543 .0212 4.8E-04 7.1E-6 9.9E-8 1.3E-9 1.1E-11 0.0 14 0.43 .4286 .3537 .2143 .1071 .0650 .0513 .0470 .0455 .0452 .0450 .0450 .0450 .0450 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 8.16E-5 7.86E-5 7.38E-5 6.85E-5 6.36E-5 5.90E-5 5.48E-5 5.09E-5 4.72E-5 4.39E-5 4.07E-5 3.78E-5 3.51E-5 3.26E-5 3.03E-5 2.81E-5 2.61E-5 2.42E-5 2.25E-5 2.09E-5 1.94E-5
2 loamy_sand_(LS) 10.91 37.3 4.05E-5 373 2 10 .43 .4282 .3740 .2736 .1661 .1043 .0767 .0651 .0600 .0570 1 .8651 .2679 .0422 .0022 7.3E-5 2.1E-6 5.6E-8 1.0E-09 0.0 14 0.43 .4282 .3740 .2736 .1661 .1043 .0767 .0651 .0600 .0583 .0573 .0571 .0570 .0570 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 4.01E-5 3.86E-5 3.62E-5 3.37E-5 3.12E-5 2.90E-5 2.69E-5 2.50E-5 2.32E-5 2.15E-5 2.00E-5 1.86E-5 1.72E-5 1.60E-5 1.49E-5 1.38E-5 1.28E-5 1.19E-5 1.10E-5 1.02E-5 9.52E-6
3 sandy_loam_(SL) 12.28 34.5 1.23E-5 345 2 10 .41 .4088 .3871 .3431 .2659 .1878 .1339 .1026 .0841 .0657 1 .8097 .3595 .1269 .0207 1.7E-3 1.1E-4 6.0E-6 2.4E-07 0.0 14 0.41 .4088 .3871 .3431 .2659 .1878 .1339 .1026 .0841 .0753 .0690 .0668 .0660 .0657 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 1.21E-5 1.17E-5 1.10E-5 1.02E-5 9.47E-6 8.79E-6 8.16E-6 7.57E-6 7.03E-6 6.53E-6 6.06E-6 5.63E-6 5.22E-6 4.85E-6 4.50E-6 4.18E-6 3.88E-6 3.60E-6 3.34E-6 3.10E-6 2.88E-6
4 silty_loam_(SIL) 22.58 38.3 1.25E-6 383 2 10 .45 .4496 .4458 .4392 .4239 .3936 .3469 .2928 .2373 .1039 1 .6382 .3764 .2441 .1251 4.5E-2 1.1E-2 1.9E-3 2.2E-04 0.0 14 0.45 .4496 .4458 .4392 .4239 .3936 .3469 .2928 .2373 .1966 .1513 .1249 .1106 .1039 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 1.24E-6 1.19E-6 1.12E-6 1.04E-6 9.64E-7 8.95E-7 8.31E-7 7.71E-7 7.16E-7 6.65E-7 6.17E-7 5.73E-7 5.32E-7 4.94E-7 4.58E-7 4.26E-7 3.95E-7 3.67E-7 3.41E-7 3.16E-7 2.93E-7
5 loam_(L) 12.9 35.2 2.89E-6 352 2 10 .43 .4293 .4217 .4074 .3754 .3223 .2608 .2071 .1633 .0884 1 .7131 .3875 .2154 .0811 1.8E-2 2.7E-3 3.0E-4 2.4E-05 0.0 14 0.43 .4293 .4217 .4074 .3754 .3223 .2608 .2071 .1633 .1361 .1101 .0972 .0910 .0884 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 2.86E-6 2.75E-6 2.58E-6 2.40E-6 2.23E-6 2.07E-6 1.92E-6 1.78E-6 1.65E-6 1.54E-6 1.43E-6 1.32E-6 1.23E-6 1.14E-6 1.06E-6 9.83E-7 9.13E-7 8.48E-7 7.87E-7 7.31E-7 6.78E-7
6 sandy_clay_(SC) 19.43 28.0 3.33E-7 280 2 10 .38 .3794 .3758 .3706 .3606 .3437 .3206 .2943 .2656 .1704 1 .3191 .1423 .0797 .0357 1.2E-2 3.1E-3 6.6E-4 1.0E-04 0.0 14 0.38 .3794 .3758 .3706 .3606 .3437 .3206 .2943 .2656 .2422 .2117 .1906 .1772 .1704 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 3.30E-7 3.18E-7 2.98E-7 2.77E-7 2.57E-7 2.39E-7 2.21E-7 2.06E-7 1.91E-7 1.77E-7 1.65E-7 1.53E-7 1.42E-7 1.32E-7 1.22E-7 1.13E-7 1.05E-7 9.78E-8 9.08E-8 8.43E-8 7.83E-8
7 silty_clay_(SIC) 27.65 29.0 5.56E-8 290 2 10 .36 .3599 .3596 .3591 .3581 .3562 .3526 .3465 .3363 .2665 1 .1439 .0803 .0569 .0369 2.1E-2 1.0E-2 4.2E-3 1.2E-03 0.0 14 0.36 .3599 .3596 .3591 .3581 .3562 .3526 .3465 .3363 .3245 .3042 .2865 .2737 .2665 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 5.49E-8 5.29E-8 4.97E-8 4.61E-8 4.28E-8 3.98E-8 3.69E-8 3.43E-8 3.18E-8 2.95E-8 2.74E-8 2.55E-8 2.36E-8 2.19E-8 2.04E-8 1.89E-8 1.76E-8 1.63E-8 1.51E-8 1.40E-8 1.30E-8
8 clay_(C) 29.12 31.2 5.56E-7 312 2 10 .38 .3799 .3792 .3784 .3767 .3735 .3679 .3592 .3461 .2707 1 .1244 .0641 .0429 .0258 1.3E-2 5.8E-3 2.0E-3 5.1E-04 0.0 14 0.38 .3799 .3792 .3784 .3767 .3735 .3679 .3592 .3461 .3324 .3101 .2915 .2782 .2707 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 5.49E-7 5.29E-7 4.97E-7 4.61E-7 4.28E-7 3.98E-7 3.69E-7 3.43E-7 3.18E-7 2.95E-7 2.74E-7 2.55E-7 2.36E-7 2.19E-7 2.04E-7 1.89E-7 1.76E-7 1.63E-7 1.51E-7 1.40E-7 1.30E-7
9 moor_(M) 47.31 75.0 8.25E-5 750 2 10 .80 .7995 .7965 .7920 .7821 .7617 .7248 .6694 .5937 .2868 1 .4269 .2497 .1733 .1049 5.2E-2 2.0E-2 5.8E-3 1.1E-03 0.0 14 0.80 .7995 .7965 .7920 .7821 .7617 .7248 .6694 .5937 .5231 .4248 .3545 .3099 .2868 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 8.16E-5 7.86E-5 7.38E-5 6.85E-5 6.36E-5 5.90E-5 5.48E-5 5.09E-5 4.72E-5 4.39E-5 4.07E-5 3.78E-5 3.51E-5 3.26E-5 3.03E-5 2.81E-5 2.61E-5 2.42E-5 2.25E-5 2.09E-5 1.94E-5
10 settlements_rock_(R) 14.00 15.0 5.56E-8 150 2 10 .20 .1999 .1996 .1992 .1984 .1969 .1942 .1900 .1837 .1474 1 .1244 .0641 .0429 .0258 1.3E-2 5.8E-3 2.0E-3 5.1E-04 0.0 14 0.20 .1999 .1996 .1992 .1984 .1969 .1942 .1900 .1837 .1771 .1664 .1575 .1510 .1474 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 5.49E-8 5.29E-8 4.97E-8 4.61E-8 4.28E-8 3.98E-8 3.69E-8 3.43E-8 3.18E-8 2.95E-8 2.74E-8 2.55E-8 2.36E-8 2.19E-8 2.04E-8 1.89E-8 1.76E-8 1.63E-8 1.51E-8 1.40E-8 1.30E-8
11 clay_loam_(CL) 21.24 31.5 7.22E-7 315 2 10 .41 .4096 .4067 .4021 .3920 .3729 .3434 .3074 .2675 .1496 1 .5005 .2721 .1720 .0882 3.4E-2 9.3E-3 1.9E-3 2.7E-04 0.0 14 0.41 .4096 .4067 .4021 .3920 .3729 .3434 .3074 .2675 .2357 .1968 .1717 .1569 .1496 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 7.14E-7 6.88E-7 6.46E-7 6.00E-7 5.57E-7 5.17E-7 4.80E-7 4.46E-7 4.14E-7 3.84E-7 3.56E-7 3.31E-7 3.07E-7 2.85E-7 2.65E-7 2.46E-7 2.28E-7 2.12E-7 1.97E-7 1.83E-7 1.70E-7
12 silt_(SI) 28.17 42.6 6.94E-7 426 2 10 .46 .4596 .4565 .4511 .4386 .4129 .3702 .3157 .2549 .0901 1 .6139 .3711 .2503 .1384 5.7E-2 1.6E-2 3.3E-3 4.3E-04 0.0 14 0.46 .4596 .4565 .4511 .4386 .4129 .3702 .3157 .2549 .2075 .1519 .1181 .0991 .0901 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 6.87E-7 6.62E-7 6.21E-7 5.77E-7 5.35E-7 4.97E-7 4.61E-7 4.28E-7 3.98E-7 3.69E-7 3.43E-7 3.18E-7 2.95E-7 2.74E-7 2.55E-7 2.36E-7 2.19E-7 2.04E-7 1.89E-7 1.76E-7 1.63E-7
13 silty_clay_loam_(SICL) 28.16 34.1 1.94E-7 341 2 10 .43 .4298 .4284 .4264 .4218 .4126 .3958 .3706 .3362 .1967 1 .4269 .2497 .1733 .1049 5.2E-2 2.0E-2 5.8E-3 1.1E-03 0.0 14 0.43 .4298 .4284 .4264 .4218 .4126 .3958 .3706 .3362 .3041 .2594 .2275 .2072 .1967 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 1.92E-7 1.85E-7 1.74E-7 1.61E-7 1.50E-7 1.39E-7 1.29E-7 1.20E-7 1.11E-7 1.03E-7 9.60E-8 8.91E-8 8.27E-8 7.68E-8 7.13E-8 6.62E-8 6.15E-8 5.71E-8 5.30E-8 4.92E-8 4.57E-8
14 sandy_clay_loam_(SCL) 13.35 29.0 3.64E-6 290 2 10 .39 .3886 .3760 .3566 .3221 .2772 .2335 .1976 .1680 .1112 1 .5525 .2157 .0922 .0256 4.6E-3 6.4E-4 7.6E-5 6.7E-06 0.0 14 0.39 .3886 .3760 .3566 .3221 .2772 .2335 .1976 .1680 .1489 .1294 .1189 .1136 .1112 0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150 21 .3333 200 3.60E-6 3.47E-6 3.25E-6 3.02E-6 2.81E-6 2.60E-6 2.42E-6 2.24E-6 2.08E-6 1.93E-6 1.80E-6 1.67E-6 1.55E-6 1.44E-6 1.33E-6 1.24E-6 1.15E-6 1.07E-6 9.91E-7 9.20E-7 8.54E-7
[soil_table]
14 # number of following entries
#Code name FC(Vol.%) mSB(Vol.%) ksat(m/s) suction. parameter Theta_sat Theta_res alpha n layer thick maxratio k-recession
# [mm] 1=par 2=tab 1/1 1/1 1/m [m] ko_rel/ku_rel per m ku/ko
#---- ---------------- ---------- --------- -------- ------- ----------- --------- --------- ----- ----- ----- ----- ------------- -----------
# 1 Sand_(S) 6.21 38.5 8.25E-5 385 1 .43 .045 14.5 2.68 31 .3333 90 .4
# 2 loamy_sand_(LS) 10.91 37.3 4.05E-5 373 1 .43 .057 7.00 1.70 31 .3333 90 .4
## 2 loamy_sand_(LS) 10.91 37.3 4.05E-5 373 1 .43 .057 12.40 2.28 31 .3333 90 .4
# 3 sandy_loam_(SL) 12.28 34.5 1.23E-5 345 1 .41 .065 7.50 1.89 31 .3333 90 .4
# 4 silty_loam_(SIL) 22.58 38.3 1.25E-6 383 1 .45 .067 2.00 1.41 31 .3333 90 .4
# 5 loam_(L) 12.9 35.2 2.89E-6 352 1 .43 .078 3.60 1.56 31 .3333 90 .4
# 6 sandy_clay_(SC) 19.43 28.0 3.33E-7 280 1 .38 .100 2.70 1.23 31 .3333 90 .4
# 7 silty_clay_(SIC) 27.65 29.0 5.56E-8 290 1 .36 .070 0.50 1.09 31 .3333 90 .4
# 8 clay_(C) 29.12 31.2 5.56E-7 312 1 .38 .068 0.80 1.09 31 .3333 90 .4
# 9 moor_(M) 47.31 75.0 8.25E-5 750 1 .80 .200 4.00 1.23 31 .3333 90 .4
# 10 settlements_rock_(R) 14.00 15.0 1.00E-9 50 1 .20 .040 8.00 1.80 31 .3333 90 .
# 11 clay_loam_(CL) 21.24 31.5 7.22E-7 315 1 .41 .095 1.90 1.31 31 .3333 90 .4
# 12 silt_(SI) 28.17 42.6 6.94E-7 426 1 .46 .034 1.60 1.37 31 .3333 90 .4
# 13 silty_clay_loam_(SICL) 28.16 34.1 1.94E-7 341 1 .43 .089 1.00 1.23 31 .3333 90 .4
# 14 sandy_clay_loam_(SCL) 13.35 29.0 3.64E-6 290 1 .39 .010 5.90 1.48 31 .3333 90 .4
#
#co- name of the
#de soil profile
#-- ---------------
1 sand_(S) { method = MultipleHorizons;
FCap = 6.21; mSB = 38.5; ksat_topmodel = 8.25E-5; suction = 385; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.75 0.1 0.05 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
PMacroThresh = 1000; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 0 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 1.0 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 0.0 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = Sand1m Sand2m ; # short descriptions
ksat = 8.25e-4 6.25e-4 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.43 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.045 0.045 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 14.5 14.5 ; # van Genuchten Parameter Alpha
Par_n = 2.68 2.68 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
2 loamy_sandsand_(S) {method = MultipleHorizons;
FCap = 10.91; mSB = 37.3; ksat_topmodel = 4.05E-5; suction = 373; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.55 0.15 0.2 0.05 0.03 0.02 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
PMacroThresh = 150 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 5 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.9 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 1.5 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = LS1m LS2m ; # short descriptions
ksat = 4.05E-5 3.05E-5 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.43 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.057 0.057 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 7.00 7.00 ; # van Genuchten Parameter Alpha
Par_n = 1.70 1.70 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
3 sandy_loam_(SL) {method = MultipleHorizons;
FCap = 12.28; mSB = 34.5; ksat_topmodel = 1.23E-5; suction = 345; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.4 0.4 0.2 0.0 0.0 0.0 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
PMacroThresh = 200 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 3 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.5 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 1.0 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = SL1m SL2m ; # short descriptions
ksat = 1.23e-5 1.03e-5 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.41 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.065 0.065 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 7.50 7.50 ; # van Genuchten Parameter Alpha
Par_n = 1.89 1.89 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
4 silty_loam_(SIL) {method = MultipleHorizons;
FCap = 22.58; mSB = 38.3; ksat_topmodel = 1.25E-6; suction = 383; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.3 0.4 0.2 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 100 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 4 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 1.0 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 2.0 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = SIL10m SIL10m ; # short descriptions
ksat = 1.25e-6 0.95e-6 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.45 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.067 0.067 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 2.0 2.0 ; # van Genuchten Parameter Alpha
Par_n = 1.41 1.41 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
maxratio = 100; # maximum proportion of effektive k-values between two layers is limited by this factor, using the higher value as reference
}
5 loam_(L) {method = MultipleHorizons;
FCap = 12.9; mSB = 35.2; ksat_topmodel = 2.89E-6; suction = 352; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.2 0.35 0.35 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 100 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 4 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 1.0 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 0.8 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = L10m L10m ; # short descriptions
ksat = 2.89e-6 2.29e-6 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.43 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.078 0.078 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 3.6 3.6 ; # van Genuchten Parameter Alpha
Par_n = 1.56 1.56 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
6 sandy_clay_(SC) {method = MultipleHorizons;
FCap = 19.43; mSB = 28.0; ksat_topmodel = 3.33E-7; suction = 280; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.3 0.3 0.3 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 100 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 3 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.4 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 0.5 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = SC10m SC10m ; # short descriptions
ksat = 3.33e-7 2.63e-7 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.38 0.30 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.1 0.1 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 2.7 2.7 ; # van Genuchten Parameter Alpha
Par_n = 1.23 1.23 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
7 silty_clay_(SIC) {method = MultipleHorizons;
FCap = 27.65; mSB = 29.0; ksat_topmodel = 5.56E-8; suction = 290; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.1 0.35 0.45 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 80 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 2 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.4 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 0.6 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = SIC10m SIC10m ; # short descriptions
ksat = 5.56e-8 4.56e-8 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.36 0.30 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.07 0.07 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 0.5 0.5 ; # van Genuchten Parameter Alpha
Par_n = 1.09 1.09 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
8 clay_(C) {method = MultipleHorizons;
FCap = 29.12; mSB = 31.2; ksat_topmodel = 5.56E-7; suction = 312; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.05 0.1 0.75 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 80 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 3 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.5 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 0.7 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = C10m C10m ; # short descriptions
ksat = 5.56e-8 4.56e-8 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.38 0.30 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.068 0.068 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 0.8 0.8 ; # van Genuchten Parameter Alpha
Par_n = 1.09 1.09 ; # van Genuchten Parameter n
Par_tau = 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
9 Moor_(M) {method = MultipleHorizons;
FCap = 47.31; mSB = 75.0; ksat_topmodel = 8.25E-5; suction = 750; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.7 0.1 0.1 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 38 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 12 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.8 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 1.6 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = M10m M10m ; # short descriptions
ksat = 8.e-4 6.e-4 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.8 0.7 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.2 0.2 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 4.0 4.0 ; # van Genuchten Parameter Alpha
Par_n = 1.2 1.2 ; # van Genuchten Parameter n
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
10 Settlement_Rock_(R) {method = MultipleHorizons;
FCap = 14.00; mSB = 15.0; ksat_topmodel = 1E-9; suction = 50; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.1 0.1 0.7 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 100 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 1 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 1.0 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 2.0 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = R10m R10m ; # short descriptions
ksat = 1.e-3 0.9e-3 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.2 0.18 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.04 0.04 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 8.0 8.0 ; # van Genuchten Parameter Alpha
Par_n = 1.8 1.8 ; # van Genuchten Parameter n
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
11 clay_loam_(CL) {method = MultipleHorizons;
FCap = 21.24; mSB = 31.5; ksat_topmodel = 7.22E-7; suction = 315; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.1 0.5 0.3 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 120 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 3 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.5 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 1.2 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = CL10m CL10m ; # short descriptions
ksat = 7.22e-7 5.22e-7 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.41 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.095 0.095 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 1.9 1.9 ; # van Genuchten Parameter Alpha
Par_n = 1.31 1.31 ; # van Genuchten Parameter n
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
12 silt_(SI) {method = MultipleHorizons;
FCap = 28.17; mSB = 42.6; ksat_topmodel = 6.94E-7; suction = 426; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.1 0.7 0.1 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 150 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 4 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 1.0 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 1.5 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = SI10m SI10m ; # short descriptions
ksat = 6.94e-7 5.94e-7 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.46 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.034 0.034 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 1.6 1.6 ; # van Genuchten Parameter Alpha
Par_n = 1.37 1.37 ; # van Genuchten Parameter n
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
13 silty_clay_(SICL) {method = MultipleHorizons;
FCap = 28.16; mSB = 34.1; ksat_topmodel = 1.94E-7; suction = 341; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.2 0.25 0.45 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 150 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 4 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 1.0 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 1.5 ; # maximum depth of the macropores
horizon = 1 2 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = SICL10m SICL10m ; # short descriptions
ksat = 1.94e-7 1.44e-7 ; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 ; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.43 0.40 ; # saturated water content (fillable porosity in 1/1)
theta_res = 0.089 0.089 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 1.00 1.00 ; # van Genuchten Parameter Alpha
Par_n = 1.23 1.23 ; # van Genuchten Parameter n
thickness = 0.16666 0.3333 ; # thickness of each single numerical layer in this horizon in m
layers = 2 59 ; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
14 profile_1 {method = MultipleHorizons;
FCap = 13.35; mSB = 29.0; ksat_topmodel = 3.64E-6; suction = 290; # optional parameters which are needed for Topmodel only
GrainSizeDist = 0.4 0.3 0.2 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4 must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
PMacroThresh = 100 ; # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
MacroCapacity = 5 ; # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
CapacityRedu = 0.8 ; # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter
MacroDepth = 1.0 ; # maximum depth of the macropores
horizon = 1 2 3; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
Name = Sand04m Clay01m Loam10m; # short descriptions
ksat = 4.0e-5 3.3e-7 3.0e-6; # saturated hydraulic conductivity in m/s
k_recession = 0.8 0.9 0.9; # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
theta_sat = 0.43 0.38 0.40; # saturated water content (fillable porosity in 1/1)
theta_res = 0.057 0.10 0.078; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
alpha = 7.00 2.70 3.60; # van Genuchten Parameter Alpha
Par_n = 1.70 1.23 1.56; # van Genuchten Parameter n
thickness = 0.10 0.05 0.355; # thickness of each single numerical layer in this horizon in m
layers = 4 2 55; # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
}
# allowed keywords for substance transport (without ""-chars):
# "radioactive" resp. "non_radioactive"
# "evaporating" resp. "non_evaporating"
# "half_time" with its unit "d"
# "min_conc" and "max_conc"
# measures: "mg/l", "g/l", "kg/kg", "Kg/Kg"; all other units will be interpreted as kg/kg (relative concentration)
[substance_transport]
0 # number of tracers to be considered (max. 9)
#
# name radioact. or not half time in days evapor. or not minim. concentr. max.conc. with unit initial initial output code writecode output path output extension output extension
#3chars if no: -9999 mg/l g/l kg/kg conc. in soil conc. in gr.w statfiles for grids with closing "\" for stat-files for grid files
#------ ---------------- ------------------- --------------- ----------------- ---------------------- ------------- ------------- --------------------- --------------------- ---------------- ---------------------- ----------------
3H radioactive half_time = 4500 d evaporating min_conc = 0 max_conc = 3500 kg/kg soilini = 1.0 gwini = 2.0 statcode = $hour_mean gridcode = $Writegrid path = $outpath statext = $code//$year gridext = $suffix
18O non_radioactive half_time = -9999 d evaporating min_conc = -9999 max_conc = -9999 kg/kg soilini = 1.0 gwini = 1.0 statcode = $hour_mean gridcode = $Writegrid path = $outpath statext = $code//$year gridext = $suffix
NACL non_radioactive half_time = -9999 d non_evaporating min_conc = 0 max_conc = 0.35 kg/kg soilini = 0.01 gwini = 0.01 statcode = $hour_mean gridcode = $Writegrid path = $outpath statext = $code//$year gridext = $suffix
# irrigation descriptions
# method 1: count MM1 DD1 amount1 MM2 DD2 amount2 MM3 DD3 amount3 MM4 DD4 amount4 MM5 DD5 amount5 MM6 DD6 amount6 MM7 DD7 amount7 MM8 DD8 amount8 MM9 DD9 amount9 MM10 DD10 amount10
# method 2a: "starting from MM DD with XX mm to MM DD with YY mm every ZZ days" here, the start end end date are explicitly given
# method 2b: "starting from MM DD with XX mm YY times every ZZ days" Here, the number of irrigation events is given explicitly */
# method 3: by demand: without additional parameters
# method 4: when ETR