=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 evaluated 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 $once_per_interval # direct discharge from surface routing module $outpath//qisr//$grid//.//$code//$year $once_per_interval # interflow from surface routing module $outpath//qbsr//$grid//.//$code//$year $once_per_interval # baseflow from surface routing module $outpath//qgsr//$grid//.//$code//$year $once_per_interval # 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 $once_per_interval # 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 max_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) 0.4 # z0 for water (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 # kd --> recession constant for single linear reservoir for direct runoff $set $kd1 = 5 # ki --> recession constant for single linear reservoir for interflow $set $ki1 = 12 # ki # dr --> drainage density (interflow generation parameter) $set $dr1 = 30 # dr # sdf --> Snow melt: Direct Runoff fraction $set $sdf1 = 0.2 # sdf [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 1 # controlling interaction with surface water: 0 = no interaction, 1 = exfiltration 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) 12 # parameter for the initialization of the gw_level (range between 1..levels (standard: 4)) $outpath//qdra//$grid//.//$code//$year $once_per_interval # results drainage discharge in mm per zone $outpath//gwst//$grid//.//$code//$year $once_per_interval # results groundwater depth $outpath//gwn_//$grid//.//$code//$year $once_per_interval # results mean groundwater recharge per zone $outpath//sb05//$grid//.//$code//$year $once_per_interval # results rel. soil moisture within the root zone per zone $outpath//sb1_//$grid//.//$code//$year $once_per_interval # results rel. soil moisture within the unsat. zone (0m..GW table) per zone $outpath//wurz//$grid//.//$code//$year $once_per_interval # results statistic of the root depth per zone $outpath//infx//$grid//.//$code//$year $once_per_interval # results statistic of the infiltration excess $outpath//pond//$grid//.//$code//$year $once_per_interval # results statistic of the ponding water storage content $outpath//qdir//$grid//.//$code//$year $once_per_interval # results statistic of the direct discharge $outpath//qifl//$grid//.//$code//$year $once_per_interval # results statistic of the interflow $outpath//qbas//$grid//.//$code//$year $once_per_interval # results statistic of the baseflow $outpath//qges//$grid//.//$code//$year $once_per_interval # results statistic of the total discharge $outpath//gwin//$grid//.//$code//$year $once_per_interval # statistic of the infiltration from surface water into groundwater (from rivers and lakes) $outpath//gwex//$grid//.//$code//$year $once_per_interval # statistic of the exfiltration from groundwater into surface water (into rivers and lakes) $outpath//macr//$grid//.//$code//$year $once_per_interval # statistic of infiltration into macropores $outpath//qinf//$grid//.//$code//$year $once_per_interval # statistic of total infiltration into the first soil layer $outpath//$SB_1_grid # grid with actual soil water content for the root zone $Writegrid # Writecode for this grid $outpath//$SB_2_grid # grid with actual soil water content for the entire unsaturated zone $Writegrid # Writecode for this grid $outpath//$ROOTgrid # grid with root depth $Writegrid # Writecode for this grid $outpath//$Thetastack # stack, actual soil water content for all soil levels $Writegrid # 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 # Writecode for this grid $outpath//$GWthetagrid # grid with theta in GWLEVEL $Writegrid # Writecode for this grid $outpath//$GWNgrid # grid with groundwater recharge $Writegrid # Writecode for this grid $outpath//$GWLEVELgrid # grid with level index of groundwater surface (Index der Schicht) $Writegrid # Writecode for this grid $outpath//$QDRAINgrid # grid with the drainage flows $Writegrid # Writecode for this grid $outpath//$SATTgrid # grid with code 1=saturation at interval start, 0 no sat. $Writegrid # Writecode for this grid $outpath//$INFEXgrid # grid with infiltration excess in mm (surface discharge) $Writegrid # Writecode for this grid $outpath//$QDgrid # grid with direct discharge 1//$Writegrid # Writecode for this grid $outpath//$QIgrid # grid with Interflow 1//$Writegrid # Writecode for this grid $outpath//$QBgrid # grid with baseflow 1//$Writegrid # Writecode for this grid $outpath//$GWINgrid # grid with infiltration from rivers into the soil (groundwater) 1//$Writegrid # Writecode for this grid $outpath//$GWEXgrid # grid with exfiltration (baseflow) from groundwater (is only generated, if groundwater module is active, else baseflow is in QBgrid) 1//$Writegrid # Writecode for this grid $outpath//$act_pond_grid # grid with content of ponding storge $Writegrid # Writecode for this grid $outpath//$UPRISEgrid # grid with amount of capillary uprise (mm) 1//$Writegrid # Writecode for this grid $outpath//$PERCOLgrid # grid with amount of percolation (mm) 1//$Writegrid # writegrid for this grid $outpath//$MACROINFgrid # grid with amount of infiltration into macropores (mm) 1//$Writegrid # Writecode for this grid $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) 3 3 # 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 withdrawal of soil water for each layer for the control point (due to transpiration) $outpath//intd//$grid//.//$code//$year # name of a file containing the interflow for the soil layers of the control point 10 # codes of the subbasins (in the subbasin grid) $kd1 # recession parameters QD (h) $ki1 # recession parameters QI (h) $dr1 # flow density (for Interflow, channels per km) 0.4 # recession parameters k for Base discharge (in QB = Q0*exp(-k/z)) with z = depth to groundwater 0.1 # correction of transmissivities Q0 for Baseflow in QB = Q0 * exp(-k/z) $sdf1 # fraction of snow melt, which is direct flow (no infiltration) $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 100 # 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 $once_per_interval # 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 5//$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 # the following section for heat transfer can be used with WaSiM version 9.0 ff [heat_transfer] 1 # 0 = do not model heat transfer, 1 = heat transfer is modelled 11 # vertical 1D heat transfer in the unsaturated zone (0=no, 1=yes, only heat diffusion, 2 = yes, heat diffusion and advection (by infiltrating water, not yet available)), 11=heat diffusion 1D with implicit solution method (recommended) 0 # vertical heat transfer in snow cover (not yet available) 0 # 2D lateral heat transfer by advection (coupled to water transport) in groundwater (not yet available) #parameters # the lower boundary condition for temperature may either be defined by a grid with the internal name _T_Lower_Boundary_Condition_ or created by using the annual temperature and the lapse rate as defined in the next two lines -10.0 # used when no grid "_T_Lower_Boundary_Condition_" was read in only: mean annual air temperature reduced to sea level to be used as lower boundary condition (e.g. 5°C) --> used for definition of the lower boundary condition at lower soil boundary, if no grid with lower boundary condition was read in -0.007 # used when no grid "_T_Lower_Boundary_Condition_" was read in only: temperature gradient (e.g. -0.007 K/m) for defining the lower boundary condition (used if no grid with lower boundary condition was found) # default soil "constants": can be changed in the soil table (using DryHeatCapacity, DryDensity and DryThermalConduct as parameter names) 800 # default heat capacity of dry soil in J/(Kg*K), default 800 --> value may be given in detail for each soil type in the soil table 1500 # default density of dry soil in Kg/m^3 , default 1500 --> value may be given in detail for each soil type in the soil table 0.58 # default thermal conductivity for dry soil in J/(m*s*K) or W/(m*K): default: 0.58 --> value may be given in detail for each soil type in the soil table 1e-12 # reduced k_sat (minimum hydraulic conductivity for fully frozrn soils) # thermodynamic constants of water and ice (not for calibration! these are constants giving only marginal room for variations) 0.5562 # thermal conductivity of liquid water 2.33 # thermal conductivity of ice (0°C...-20°C) 4187 # heat capacity of water in J/(Kg*K) 1940 # heat capacity of ice at -20°C in J/(Kg*K) 2090 # heat capacity of ice at 0°C in J/(Kg*K) 334000 # latent heat of freezing in J/Kg 1000 # density of water in Kg/m^3 # other parameters (not for calibrating, but there is no clear lioterature value) 1.22 # scaling factor (solution of the clapeyron equation, literature gives values of 1.8 up to 123, bhut this may be measure dependent. Theoretical value is dH/T_m = 1.22 J/(Kg*K)) 300 # minimum sub time step allowed for heat transfer model (if the required time step would be shorter, numeric errors like extrem temperature fluctuations are possible). Recommendation: for soil layers of 5cm: 3...180, 1cm layers: 3...30 1200 # maximum sub time step allowed for heat transfer model (to avoid instabilites induced by the nonlinearity of the processes) Recommendation: for soil layers of 5cm: 180; 1cm layers: 30 1.0 # n-factor for freezing (factor applied to the air temperature to get the temperature at the soil surface as upper boundary condition when temperatures are negative 1.0 # n-factor for thawing (factor applied to the air temperature to get the temperature at the soil surface as upper boundary condition when temperatures are positive # this value ranges from 0.01 to 0.99 with 0.01 defining beginning freezing (1% ice) and 0.99 defining complete freezing (99% ice, only smallest pores may contain water # output grids and statistics $outpath//ts_loc//$grid//.//$code//$year # results soil temperature for control point $outpath//ts_avg//$grid//.//$code//$year $once_per_interval # results soil temperature thaw depth or active layer thickness as average value for subbasins $outpath//$Temperaturestack # stack, actual soil water content for all soil levels $Writegrid # Writecode for this stack $outpath//$ThawDepthGridTMod # grid containing the active layer thickness relative to the soil surface (deepest thawing front in the soil profile) $Writegrid # Writecode for this stack $readgrids # --> could be used from soil model? [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//gwtable.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 2 # 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 $once_per_interval # statistic of the irrigation water from groundwater $outpath//irsw//$grid//.//$code//$year $once_per_interval # 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 1000 # if iterative solving method (1): max.numberof iterations 0.0001 # 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.20 # 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 3 3 # 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 $Writegrid # 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) # this paragraph is not needed for WaSiM-uzr but for the WaSiM-version with the variable saturated area approach (after Topmodel) [soil_model] 0 # 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 $once_per_interval # statistic of the surfeca discharge $outpath//qifl//$grid//.//$code//$year $once_per_interval # statistic of the Interflows $outpath//qbas//$grid//.//$code//$year $once_per_interval # statistic of the base flow $outpath//qbav//$grid//.//$code//$year $once_per_interval # statistic of the slow base flow $outpath//qges//$grid//.//$code//$year $once_per_interval # statistic of the total discharge $outpath//sb__//$grid//.//$code//$year $once_per_interval # soil storage in mm per zone $outpath//suz_//$grid//.//$code//$year $once_per_interval # drainage storage in mm per zone $outpath//sifl//$grid//.//$code//$year $once_per_interval # interflow storage in mm per zone $outpath//sd__//$grid//.//$code//$year $once_per_interval # saturation deficite per zone in mm 10 # Codes der Teilgebiete im Zonengrid 0.015 # Rezessionsparameterter m fuer Saettigungsflaechenmodell in Metern 40.0 # Korrekturfaktor fuer Transmissivitaeten 8.0 # Korrekturfaktor fuer K-Wert (vertikale Versickerung), Modell erwartet k in m/s 6.0 # Speicherrueckgangskonstante Direktabflus ELS in h 0.0 # Saettigungsdefizit, bei dessen Unterschreitung lokaler Interflow gebildet wird 1.0 # Speicherrueckgangskonstante Interflow ELS in h 3600 # Rueckgangskonstante verzoegerter Basisabfluss in h 0.03 # maximale Tiefenversickerungsrate bei Saettigung in mm/h 0.01 # Anfangswert QBB 0.0 # Anfangsfuellung des SUZ-Speichers in n*nFK 0.45 # Anfangssaettigungsdefizit in n*nFK, beeinflusst den ersten Basisabfluss 3.0 # Anspringpunkt fuer Makroporenabfluss (in mm/h!, bezogen auf Stundenniederschlag!), alles darueber geht direkt in den Drainspeicher! 0.9 # Reduktionsfaktor fuer Auffuellung von Verdunstungsverlusten aus dem Grundwasser und aus dem Interflowspeicher 0.4 # 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 [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 1 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) $outpath//qgko//$grid//.//$code//$year $routing_code # name of the statistic file with routed discharges $inpath//spend_//$year//.dat # name of the file with observed discharges (mm/Timestep or m^3/s) 1 # number of following collumn descriptor 10 1 # if the first code would be a 7, then it would mean, that the modeled discharge of subbasin 1 (or lowest subbasin code) would communicate with the data column 7 in the specific discharge data file (date-columns are not counted!) 720 # 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 16 (AE=93.060, AErel=1.0) from OL 17 (kh=0.1, kv=0.4, Bh= 9.3, Bv= 37.1, Th= 0.93, Mh=25.0, Mv=10.0, I=0.0058, L=15046.7, AE=45.820) TG 8 (AE=333.080, AErel=1.0) from OL 11 (kh=0.1, kv=0.4, Bh=16.1, Bv= 64.3, Th= 1.61, Mh=25.0, Mv=10.0, I=0.0032, L=1772.8, AE=147.230) and OL 10 (kh=0.1, kv=0.4, Bh=13.7, Bv= 55.0, Th= 1.37, Mh=25.0, Mv=10.0, I=0.0079, L=9366.9, AE=153.390) TG 7 (AE=704.950, AErel=1.0) from SUMTRIB 8&9 (kh=0.1, kv=0.4, Bh=20, Bv= 80, Th= 2.00, Mh=25.0, Mv=10.0, I=0.0025, L=23150.1, AE=498.82) TG 6 (AE=940.580, AErel=1.0) from OL 7 (kh=0.1, kv=0.4, Bh=28.3, Bv=113.1, Th= 2.83, Mh=25.0, Mv=10.0, I=0.0036, L=20340.8, AE=704.950) TG 5 (AE=1388.380, AErel=1.0) from SUMTRIB 16&18&19 (kh=0.1, kv=0.4, Bh=12, Bv= 50, Th= 1.2, Mh=25.0, Mv=10.0, I=0.0010, L=20000, AE=135.91) and OL 6 (kh=0.1, kv=0.4, Bh=40.0, Bv=160.0, Th= 4.00, Mh=25.0, Mv=10.0, I=0.0010, L=20549.7, AE=940.580) and SP 1 ( file = $outpath//Lake__01.//$code//$year , V0 = 2.215E09, C0 = 0 0 0 0 0 0 0 0 0 ) TG 4 (AE=1547.030, AErel=1.0) from OL 5 (kh=0.1, kv=0.4, Bh=45.9, Bv=183.5, Th= 4.59, Mh=25.0, Mv=10.0, I=0.0010, L=16460.9, AE=1388.380) TG 20 (AE=1579.340, AErel=1.0) from OL 4 (kh=0.1, kv=0.4, Bh=48.2, Bv=192.8, Th= 4.82, Mh=25.0, Mv=10.0, I=0.0010, L=100.0, AE=1547.030) TG 21 (AE=1579.460, AErel=1.0) from OL 20 (kh=0.1, kv=0.4, Bh=48.6, Bv=194.3, Th= 4.86, Mh=25.0, Mv=10.0, I=0.0010, L=200.0, AE=1579.340) and AL 1 ( modus = intern_with_rule ) TG 13 (AE=180.840, AErel=1.0) from OL 15 (kh=0.1, kv=0.4, Bh=11.8, Bv= 47.2, Th= 1.18, Mh=25.0, Mv=10.0, I=0.0056, L=14456.8, AE=85.870) TG 12 (AE=370.080, AErel=1.0) from SUMTRIB 13&14 (kh=0.05, kv=0.4, Bh=15, Bv= 60, Th= 1.5, Mh=25.0, Mv=10.0, I=0.0037, L=27000, AE=271.730) TG 22 (AE=1955.980, AErel=1.0) from OL 21 (kh=0.1, kv=0.4, Bh=36.5, Bv=146.0, Th= 3.65, Mh=25.0, Mv=10.0, I=0.0046, L=2938.5, AE=1579.460) and OL 12 (kh=0.1, kv=0.4, Bh=20.5, Bv= 81.8, Th= 2.05, Mh=25.0, Mv=10.0, I=0.0055, L=3431.4, AE=370.080) TG 23 (AE=1956.530, AErel=1.0) from OL 22 (kh=0.1, kv=0.4, Bh=52.6, Bv=210.5, Th= 5.26, Mh=25.0, Mv=10.0, I=0.0010, L=282.8, AE=1955.980) and AL 2 ( modus = intern_with_rule ) TG 3 (AE=1960.040, AErel=1.0) from OL 23 (kh=0.1, kv=0.4, Bh=48.8, Bv=195.1, Th= 4.88, Mh=25.0, Mv=10.0, I=0.0015, L=1165.7, AE=1956.530) TG 24 (AE=1976.010, AErel=1.0) from OL 3 (kh=0.1, kv=0.4, Bh=43.8, Bv=175.0, Th= 4.38, Mh=25.0, Mv=10.0, I=0.0027, L=8184.0, AE=1960.040) and ZL 1 ( modus = intern , kh=0.4, kv=0.4, Bh=3.0, Bv=10.0, Th=2.0, Mh=25.0, Mv=15.0, I=0.0066, L=2000.5, AE=1 ) TG 99 (AE=1976.260, AErel=1.0) from OL 24 (kh=0.1, kv=0.4, Bh=43.8, Bv=175.0, Th= 4.38, Mh=25.0, Mv=10.0, I=0.0027, L=8184.0, AE=1976.010) and ZL 2 ( modus = intern , kh=0.4, kv=0.4, Bh=10.0, Bv=20.0, Th=2.0, Mh=25.0, Mv=15.0, I=0.0066, L=5000, AE=1 ) TG 25 (AE=1976.460, AErel=1.0) from OL 99 (kh=0.1, kv=0.4, Bh=32.2, Bv=128.9, Th= 3.44, Mh=25.0, Mv=10.0, I=0.0114, L=200, AE=1976.260) TG 26 (AE=2107.290, AErel=1.0) from OL 25 (kh=0.1, kv=0.4, Bh=44.0, Bv=176.1, Th= 4.40, Mh=25.0, Mv=10.0, I=0.0026, L=14371.0, AE=1976.460) TG 27 (AE=2107.840, AErel=1.0) from OL 26 (kh=0.1, kv=0.4, Bh=54.1, Bv=216.5, Th= 5.41, Mh=25.0, Mv=10.0, I=0.0010, L=241.4, AE=2107.290) TG 2 (AE=2215.900, AErel=1.0) from OL 27 (kh=0.05, kv=0.4, Bh=43.4, Bv=173.5, Th= 4.34, Mh=25.0, Mv=10.0, I=0.0033, L=5835.5, AE=2107.840) TG 1 (AE=2255.60, AErel=1.0) from OL 2 (kh=0.1, kv=0.4, Bh=46.6, Bv=186.2, Th= 4.66, Mh=25.0, Mv=10.0, I=0.0025, L=15588.2, AE=2215.900) # 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_1] 4 20 0 20 7 27 7 27 8 TargetCap = 8 [abstraction_rule_abstraction_2] 12 60 91 121 182 213 244 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) # 28.02. 31.03. 30.04. 30.06. 31.07. 31.08. 31.12. # 7 7.5 10.5 12.5 11 8.5 7 # Restwassermengen in m3/s 7 0 0 0 0 0 0 0 7.5 0.5 0 0 0 0 0 0.5 8.5 1.5 1 0 0 0 0 1.5 10.5 3.5 3 0 0 0 2 3.5 11 4 3.5 0.5 0 0 2.5 4 12.5 5.5 5 2 0 1.5 4 5.5 60 53 52.5 50.5 47.5 49 51.5 53 60.5 53 53 51 48 49.5 52 53 61.5 53 53 52 49 50.5 53 53 63.5 53 53 53 51 52.5 53 53 64 53 53 53 51.5 53 53 53 65.5 53 53 53 53 53 53 53 TargetCap = 60 60 60 60 60 60 60 [abstraction_rule_reservoir_1] 6 0 0 8.750e05 0 1.000e06 0.1 1.125e06 2.8 1.250e06 8 1.375e06 40 # 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] 3 # count of multilayer landuses 2 settlements { Landuse_Layers = 2, -9999, -9999; k_extinct = 0.3; LAI_scale = 20;} 4 mixed_forest { Landuse_Layers = 5, 8, -9999; k_extinct = 0.3; LAI_scale = 20;} 8 grassland { Landuse_Layers = 7, -9999, -9999; k_extinct = 0.3; LAI_scale = 20;} # 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 = 1 $set $JDVegWrite = 365 [landuse_table] 8 # number of following land use codes 16 water {method = VariableDayCount; RootDistr = 1; TReduWet = 1; LimitReduWet = 1; HReduDry = 150; IntercepCap = 0; JulDays = 365; Albedo = 0.1; rsc = 0.1; rs_interception = 0; rs_evaporation = 0; LAI = 0; Z0 = 0.3; VCF = 0; RootDepth = 0; AltDep = 0; } 2 settlements {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.2; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; Albedo = 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 ; rsc = 100 100 100 100 100 100 100 100 100 100 100 100 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 1 1 1 1 1 1 1 1 1 1 1 1 ; Z0 = 1 1 1 1 1 1 1 1 1 1 1 1 ; VCF = 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 ; 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 ; 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; } 3 pine_forest {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.6; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; 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 ; rsc = 80 80 75 65 55 55 55 55 55 75 80 80 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 ; LAI = 6 6 8 8 10 10 10 10 8 8 6 6 ; Z0 = 3 3 3 3 3 3 3 3 3 3 3 3 ; 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 ; 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 ; 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; } 4 decidous_forest {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.6; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; 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 ; rsc = 100 100 95 75 65 65 65 65 65 85 100 100 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 1500 1500 1500 1500 1500 1500 1500 1500 1500 1500 1500 1500 ; LAI = 1 1 4 4 6 7 7 6 5 4 1 1 ; Z0 = 2 2 2 2 2 2 2 2 2 2 2 2 ; 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 ; 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 ; 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; } 5 mixed_forest {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.6; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; 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 ; rsc = 90 90 85 70 60 60 60 60 60 80 90 90 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 1200 1200 1200 1200 1200 1200 1200 1200 1200 1200 1200 1200 ; LAI = 3 3 3 6 8 8 8 8 8 6 3 3 ; Z0 = 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 ; 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 ; 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 ; 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; } 6 agriculture {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.4; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; 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 ; rsc = 80 80 75 75 65 55 55 55 65 75 90 90 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 1 1 2 3 4 5 5 4 3 2 1 1 ; Z0 = 0.03 0.03 0.03 0.04 0.05 0.05 0.05 0.05 0.04 0.03 0.03 0.03 ; VCF = 0.3 0.3 0.3 0.7 0.8 0.95 0.95 0.8 0.7 0.3 0.3 0.3 ; RootDepth = 0.15 0.15 0.2 0.4 0.5 0.5 0.5 0.5 0.4 0.2 0.15 0.15 ; 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; } 7 extensive_grassland {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.4; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; 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 ; rsc = 90 90 80 70 60 55 50 55 60 70 90 90 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 600 600 600 600 600 600 600 600 600 600 600 600 ; LAI = 2 2 2 2 3 3 3 3 3 2 2 2 ; Z0 = 0.03 0.03 0.03 0.04 0.04 0.04 0.04 0.04 0.04 0.03 0.03 0.03 ; VCF = 0.8 0.8 0.8 0.9 0.9 0.9 0.9 0.9 0.8 0.8 0.8 0.8 ; 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 ; 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; } 8 forest_grass {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.4; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; 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 ; rsc = 90 90 80 75 70 65 60 65 70 80 90 90 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 ; LAI = 2 2 2 2 2 2 2 2 2 2 2 2 ; Z0 = 0.03 0.03 0.03 0.04 0.04 0.04 0.04 0.04 0.04 0.03 0.03 0.03 ; VCF = 0.7 0.7 0.7 0.8 0.8 0.8 0.8 0.8 0.7 0.7 0.7 0.7 ; 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 ; 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; } 9 bushes {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.6; JulDays = 15 46 74 105 135 166 196 227 258 288 319 349 ; Albedo = 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 ; rsc = 80 80 70 70 60 50 50 60 60 70 70 80 ; rs_interception = 5 5 5 5 5 5 5 5 5 5 5 5 ; rs_evaporation = 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 ; LAI = 3 3 3 4 5 5 4 4 3 3 3 3 ; Z0 = 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 ; 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 ; 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 ; 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; } # 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 = 0; # 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 = 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 } # 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 = 0; # 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 (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 = 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) } $set $e3 = e-2 $set $e4 = e-3 $set $e5 = e-4 $set $e6 = e-5 $set $e7 = e-6 $set $e8 = e-7 $set $e9 = e-8 $set $e10 = e-9 [soil_table] 2 # number of following entries 7 silty_clay_(SIC) {method = MultipleHorizons; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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 DryHeatCapacity = 810 # dry heat capacity in J/(Kg*K) DryDensity = 1450 # dry density in m^3/m^3 DryThermalConduct = 0.57 # dry thermal conductivity in W/(m*K) (or J/(m*s*K) --> 1J = 1Ws) KMinFrozenSoil = 1e-15 # minimum hydraulic conductivity in m/s when the soil is coimpletely frozen horizon = 1 2 3 ; # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row Name = Peat SIC SIC ; # short descriptions ksat = 1e-4 1e-5 1e-7 ; # saturated hydraulic conductivity in m/s k_recession = 0.9 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.8 0.40 0.30 ; # saturated water content (fillable porosity in 1/1) theta_res = 0.05 0.07 0.07 ; # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation) alpha = 1.5 1.5 1.5 ; # van Genuchten Parameter Alpha (small values for silt of e.g. 0.5 to high values of e.g. 4 to 8 for sand, gravel and peat Par_n = 2 2 2 ; # van Genuchten Parameter n (1.5 for silt to 4.5 for gravel and peat) Par_tau = 0.5 0.5 0.5 ; # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5) thickness = 0.05 0.1 0.1 ; # thickness of each single numerical layer in this horizon in m layers = 4 6 91 ; # 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) } 1 sand_(S) {method = MultipleHorizons; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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.05 0.1 ; # thickness of each single numerical layer in this horizon in m layers = 2 99 ; # 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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 = 2.0 2.0 ; # 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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 = 1.0 1.0 ; # van Genuchten Parameter Alpha Par_n = 1.9 1.9 ; # 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.5 0.1 ; # thickness of each single numerical layer in this horizon in m layers = 2 99 ; # 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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.49 1.49 ; # 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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; # max 11 default parameters SU_PAR00 to SU_PAR10 for silting up module. These parameters can be used when defining expressions for inf_infinite (Y), energy (W), SDISP (Z), inf_start(X), and inf_pot (V) in method 3 SU_PAR01 = 65.0; # for method 3: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) SU_PAR02 = 12.2; # for method 3: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) SU_PAR03 = 0.52; # for method 3: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) SU_PAR04 = -0.64; # for method 3: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) SU_PAR05 = 0.013; # for method 3: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) SU_PAR06 = -1.03; # for method 3: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) SU_PAR07 = 0.7; # for method 3: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) SU_PAR08 = -0.19; # for method 3: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) 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 #------ ---------------- ------------------- --------------- ----------------- ---------------------- ------------- ------------- --------------------- --------------------- ---------------- ---------------------- ---------------- 18O non_radioactive half_time = -9999 d evaporating min_conc = -9999 max_conc = -9999 kg/kg soilini = 1.0 gwini = 1.0 statcode = $once_per_interval 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 = $once_per_interval gridcode = $Writegrid path = $outpath statext = $code//$year gridext = $suffix 3H radioactive half_time = 4493 d evaporating min_conc = 0 max_conc = 3500 kg/kg soilini = 3.0 gwini = 3.0 statcode = $once_per_interval 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 [irrigation_table] 3 # number of following irrigation codes, per row one use # #Code name method from control by # (0=no irr, (1=GW demand: table: # 1=table1, 2=river) psi[m] [mm] [mm] [mm] [mm] [mm] [mm] [mm] [mm] [mm] [mm] # 2=table2) start stop 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 # 3=demand or # 4=ETR/ETP