# comment lines are indicated by an "#". They are allowed after the entries or at the beginning of the lines # # paragraphs of the control file: # # [output_list] # [output_interval] # [coordinates] (geogr. coordinates) # [region_transition_distance] for multiple regression regions # [elevation_model] (name of the elevation model) # [zonengrid] (name of the zone grid) # [standardgrids] (name of other static grids like slope angle, slope aspect, topogr. factor) # [variable_grids] (names of albedo and soil storage - used by more than one modules) # [model_time] (start end end-dates of model time) # [meteo_data_count] (number of mete data to interpolate) # [meteo_names] (names of meteo data to interpolate - each name is later the headline of a paragraph for interpolation # [precipitation_correction] (paragraph for parameters of the prec.-correction) # [radiation_correction] (paragraph with parameters for radiation correction) # [evapotranspiration] (paragraph with parameters for evapotranspiration) # [snow_model] (paragraph with parameters for the snow model) # [ice_firn] (paragraph with parameters for the glacier model) # [interzeption_model] (paragraph with parameters for the interzeption modell) # [infiltration_model] (paragraph with parameters for the infiltration modell) # [soil_model] (paragraph with parameters for the soilmodel) # [SiltingUpModel] (paragraph with parameters for the silting up model) # [unsatzon_model] (paragraph with parameters for the unsaturated zone model) # [SurfaceRoutingModel] (paragraph with parameters for the surface routing model) # [irrigation] (paragraph with parameters for irrigation model) # [groundwater_flow] (paragraph with parameters for groundwater model) # [ExternalCoupling] (paragraph with parameters for external coupling) # [routing_model] (paragraph with Parametern for discharge routing) # [landuse_table] (paragraph with land use (vegetation) parameters) # [soil_table] (paragraph with soil properties) # [substance_transport] (paragraph with tracer properties) # [abstraction_rule_reservoir_n] rules for routing submodel, each reservoir has its own rule # [multilayer_landuse] new multilayer landuse table definition # # symbol definitions begin with the set command. $set $mainpath = c:\_sim_\web\wasim_9-1-0\ $set $InitialStateDirectory = $mainpath//StateIni\ $set $DefaultOutputDirectory = $mainpath//output\ $set $inpath_grid = $mainpath//input_grid\ $set $inpath_meteo = $mainpath//input_meteo\ $set $inpath_hydro = $mainpath//input_hydro\ $set $inpath_ini = $mainpath//input_ini\ # it is important to set $outpath to an empty string in order to activate $DefaultOutputDirectory $set $outpath = # readgrids : 1 = read storage grids (as SI, SSNOW,SLIQ...) from hard disk, 0=generate and initialize with 0 $set $readgrids = 1 # read grids for dynamic phenology -> usually chilling grid should be read in if availabe because otherwise thermal time method will be applied and not the sequential model $set $DPreadgrids = 1 $set $time = 60.0 $set $year = 10 $set $starthour = 1 $set $startday = 1 $set $startmonth = 1 $set $startyear = 2010 $set $endhour = 24 $set $endday = 31 $set $endmonth = 12 $set $endyear = 2010 $set $grid = r500 $set $stack = s500 $set $suffix = grd $set $code = s # variables for standardgrids # first section: grids, which differ for different subdivisions of the basin $set $zone_grid = $grid//.ezg $set $subcatchments = $grid//.ezg $set $flow_time_grid = $grid//.fzs $set $river_links_grid = $grid//.lnk $set $regio_grid = $grid//.reg #second section: grids, which doesn't depend on subdivision (only pixel-values are of interest) $set $elevation_model = $grid//.dhm $set $RelCellArea_grid = $grid//.rca $set $CellSizeX_grid = $grid//.csx $set $CellSizeY_grid = $grid//.csy $set $slope_grid = $grid//.slp $set $FlowDirection_grid = $grid//.fld $set $aspect_grid = $grid//.exp $set $land_use_grid = $grid//.use $set $ice_firn_grid = $grid//.ice $set $field_capacity_grid = $grid//.nfk $set $ATBgrid = $grid//.atb $set $hydr_cond_grid = $grid//.k $set $soil_types = $grid//.soil $set $sky_view_factor_grid = $grid//.hor $set $drain_depth_grid = $grid//.drn $set $drain_distance_grid = $grid//.dis $set $irrigationcodes = $grid//.irr $set $max_pond_grid = $grid//.maxpond $set $clay_depth_grid = $grid//.cly $set $river_depth_grid = $grid//.dep $set $river_width_grid = $grid//.wit $set $tracer_1 = $grid//.c1 $set $tracer_2 = $grid//.c2 $set $tracer_3 = $grid//.c3 $set $tracer_4 = $grid//.c4 $set $tracer_5 = $grid//.c5 $set $tracer_6 = $grid//.c6 $set $tracer_7 = $grid//.c7 $set $tracer_8 = $grid//.c8 $set $tracer_9 = $grid//.c9 $set $kolmationsgrid = $grid//.kol $set $gw_kx_1_grid = $grid//.kx1 $set $gw_kx_2_grid = $grid//.kx2 $set $gw_kx_3_grid = $grid//.kx3 $set $gw_ky_1_grid = $grid//.ky1 $set $gw_ky_2_grid = $grid//.ky2 $set $gw_ky_3_grid = $grid//.ky3 $set $gw_bound_h_1_grid = $grid//.bh1 $set $gw_bound_h_2_grid = $grid//.bh2 $set $gw_bound_h_3_grid = $grid//.bh3 $set $gw_bound_q_1_grid = $grid//.bq1 $set $gw_bound_q_2_grid = $grid//.bq2 $set $gw_bound_q_3_grid = $grid//.bq3 $set $aquiferthick1 = $grid//.aq1 $set $aquiferthick2 = $grid//.aq2 $set $aquiferthick3 = $grid//.aq3 $set $gw_storage_coeff_1 = $grid//.s01 $set $gw_storage_coeff_2 = $grid//.s02 $set $gw_storage_coeff_3 = $grid//.s03 $set $gw_kolmation_1 = $grid//.gk1 $set $gw_kolmation_2 = $grid//.gk2 $set $gw_kolmation_3 = $grid//.gk3 $set $lake_grid = $grid//.lak $set $taucrit_grid = $grid//.tau $set $ThawCoeffPermaFrost = $grid//.alpha $set $T_lower_boundary_grid = $grid//.tlowbdry $set $debris_on_glaciers = $grid//.debris # grids for surface hydrology modules $set $forcingunitsgrid1 = forc1//$grid//.//$suffix $set $TStartPhenoGrid1 = phen1//$grid//.//$suffix $set $chillingunitsgrid1 = chill1//$grid//.//$suffix $set $FStargrid1 = fstar1//$grid//.//$suffix $set $forcingunitsgrid2 = forc2//$grid//.//$suffix $set $TStartPhenoGrid2 = phen2//$grid//.//$suffix $set $chillingunitsgrid2 = chill2//$grid//.//$suffix $set $FStargrid2 = fstar2//$grid//.//$suffix $set $forcingunitsgrid3 = forc3//$grid//.//$suffix $set $TStartPhenoGrid3 = phen3//$grid//.//$suffix $set $chillingunitsgrid3 = chill3//$grid//.//$suffix $set $FStargrid3 = fstar3//$grid//.//$suffix $set $albedo = albe//$grid//.//$suffix $set $soilstoragegrid = sb__//$grid//.//$suffix $set $throughfall = qi__//$grid//.//$suffix $set $snowcover_outflow = qsno//$grid//.//$suffix $set $melt_from_snowcover = qsme//$grid//.//$suffix $set $days_snow = sday//$grid//.//$suffix $set $snow_age = sage//$grid//.//$suffix $set $snow_rate = snow//$grid//.//$suffix $set $rain_rate = rain//$grid//.//$suffix $set $firn_melt = qfir//$grid//.//$suffix $set $ice_melt = qice//$grid//.//$suffix $set $preci_grid = prec//$grid//.//$suffix $set $preci_grid1 = prec1//$grid//.//$suffix $set $preci_grid2 = prec2//$grid//.//$suffix $set $irrig_grid = irri//$grid//.//$suffix $set $etr2etpgrid = er2ep//$grid//.//$suffix $set $tempegrid = temp//$grid//.//$suffix $set $tempegrid1 = temp1//$grid//.//$suffix $set $tempegrid2 = temp2//$grid//.//$suffix $set $windgrid = wind//$grid//.//$suffix $set $sunshinegrid = ssd_//$grid//.//$suffix $set $radiationgrid = rad_//$grid//.//$suffix $set $humiditygrid = humi//$grid//.//$suffix $set $vaporgrid = vapo//$grid//.//$suffix $set $ETPgrid = etp_//$grid//.//$suffix $set $EIPgrid = eip_//$grid//.//$suffix $set $ETRgrid = etr_//$grid//.//$suffix $set $EVAPgrid = evap//$grid//.//$suffix $set $EVARgrid = evar//$grid//.//$suffix $set $ETRSgrid = etrs//$grid//.//$suffix $set $SSNOgrid = ssno//$grid//.//$suffix $set $SLIQgrid = sliq//$grid//.//$suffix $set $SSTOgrid = ssto//$grid//.//$suffix $set $sat_def_grid = sd__//$grid//.//$suffix $set $SUZgrid = suz_//$grid//.//$suffix $set $SIFgrid = sif_//$grid//.//$suffix $set $EIgrid = ei__//$grid//.//$suffix $set $SIgrid = si__//$grid//.//$suffix $set $ExpoCorrgrid = exco//$grid//.//$suffix $set $Tcorrgrid = tcor//$grid//.//$suffix $set $Shapegrid = shap//$grid//.//$suffix $set $INFEXgrid = infx//$grid//.//$suffix $set $SATTgrid = satt//$grid//.//$suffix $set $Nagrid = na__//$grid//.//$suffix $set $SSPgrid = ssp_//$grid//.//$suffix $set $Peakgrid = peak//$grid//.//$suffix $set $SBiagrid = sbia//$grid//.//$suffix $set $fcia_grid = nfki//$grid//.//$suffix $set $tavg_grid = tavg//$grid//.//$suffix # now variables for unsaturated zone model $set $SB_1_grid = sb05//$grid//.//$suffix $set $SB_2_grid = sb1_//$grid//.//$suffix $set $ROOTgrid = wurz//$grid//.//$suffix $set $QDgrid = qd__//$grid//.//$suffix $set $QIgrid = qifl//$grid//.//$suffix $set $GWdepthgrid = gwst//$grid//.//$suffix $set $GWthetagrid = gwth//$grid//.//$suffix $set $GWNgrid = gwn_//$grid//.//$suffix $set $UPRISEgrid = uprs//$grid//.//$suffix $set $PERCOLgrid = perc//$grid//.//$suffix $set $GWLEVELgrid = gwlv//$grid//.//$suffix $set $QDRAINgrid = qdrn//$grid//.//$suffix $set $QBgrid = qb__//$grid//.//$suffix $set $GWINgrid = gwin//$grid//.//$suffix $set $GWEXgrid = gwex//$grid//.//$suffix $set $act_pond_grid = pond//$grid//.//$suffix $set $MACROINFgrid = macr//$grid//.//$suffix $set $SUBSTEPSgrid = step//$grid//.//$suffix $set $SnowFreeDaysGrid = sfre//$grid//.//$suffix $set $SnowCoverDaysGrid = scov//$grid//.//$suffix $set $ThawDepthGrid = thdp//$grid//.//$suffix $set $ThawDepthGridTMod = thaw//$grid//.//$suffix # variables for groundwater modeling $set $flowx1grid = gwx1//$grid//.//$suffix $set $flowx2grid = gwx2//$grid//.//$suffix $set $flowx3grid = gwx3//$grid//.//$suffix $set $flowy1grid = gwy1//$grid//.//$suffix $set $flowy2grid = gwy2//$grid//.//$suffix $set $flowy3grid = gwy3//$grid//.//$suffix $set $head1grid = gwh1//$grid//.//$suffix $set $head2grid = gwh2//$grid//.//$suffix $set $head3grid = gwh3//$grid//.//$suffix $set $GWbalance1grid = gwbalance1//$grid//.//$suffix $set $GWbalance2grid = gwbalance2//$grid//.//$suffix $set $GWbalance3grid = gwbalance3//$grid//.//$suffix # result grids for surface routing model $set $surfspeed_grid = sfcv//$grid//.//$suffix $set $surfflux_grid = sflx//$grid//.//$suffix # some new stacks and grids for the dynamic glacier model $set $firn_WE_stack = glfirn//$stack//.//$suffix $set $GlacierMassBalance = glmb//grid//.//$suffix $set $OldGlacierMassBalance = glmb_old//grid//.//$suffix $set $glacierizedCells_grid = glc_//$grid//.//$suffix $set $glacier_codes_grid = glid//$grid//.//$suffix # result-stacks for Unsatzonmodel $set $Thetastack = teth//$stack//.//$suffix $set $hydraulic_heads_stack = hhyd//$stack//.//$suffix $set $geodetic_altitude_stack = hgeo//$stack//.//$suffix $set $flowstack = qu__//$stack//.//$suffix $set $concstack = conc//$stack//.//$suffix # result-stacks for temperatures in Unsatzonmodel $set $Temperaturestack = tsoil//stack//.//$suffix # parameters for interpolation of meteorological input data $set $SzenUse = 0 $set $IDWmaxdist = 20000 $set $IDWweight = 2 $set $Anisoslope = 0.0 $set $Anisotropie = 1.0 # explanation of writegrid and outputcode some lines below $set $Writegrid = 3 $set $Writestack = 3 $set $once_per_interval = 2001 $set $avrg_per_24Invs = 2024 $set $sum_per_24Invs = 4024 $set $routing_code = 5001 # Writegrid : max. 4 digits (nnnn) # # only if writegrid >= 1000: 1. digit (1nnn, or 2nnn) # 0 = no vegetation period based grid is written # 1 = sum grid is written for vegetation period (summing up each value as long as this cells vegetation period is active) # 2 = average value grid is written for vegetation period (summing up each value as long as this cells vegetation period is active) # only if writegrid >= 100: 2. digit (n1nn, or n2nn or n3nn or 1nn..3nn -> leading digits may be omitted)) # 0 = no minimum or maximum grid is written # 1 = minimum grid is written (minimum value for each of the grid cells over the entire model period) # 2 = maximum grid is written (maximum value for each of the grid cells over the entire model period) # 3 = both grids are written (minimum and maximum value for each of the grid cells over the entire model period) # only if Writegrid >= 10: 3rd digit: sums or means (1n ... 8n or n1n..n8n or nn1n..nn8n -> leading digits may be omitted)) # 0 = no sum grid will be written # 1 = one sum grid will be written at the end of the model run # 2 = one sum grid per model year # 3 = one sum grid per model month # 4 = one sum grid per day (only, if timestep < 1 day) # 5 = one mean value grid at the end of the model run # 6 = one mean value grid per model year # 7 = one mean value grid per month # 8 = one mean value grid per day # last digit (nnn1 .. nnn5 or nn1..nn5 or n1..n5 or 1..5 -> leading digits may be omitted) (for actual values, not for Sums or means) # 1 = (over)write each timestep into the same grid (for security in case of model crashs) # 2 = write grids each timestep to new files, the name is build from the first 4 letters # of the regular grid name and then from the number of month, day and hour (hoer as file extension). # example: tempm500.grd will become prec0114.07 for 14.January, 7:00. # 3 = only the last grid of the model run will be stored # 4 = the grid from the last hour of each day (24:00) will be stored (for each day the same file will be overwritten) # 5 = like 4, but each day a new grid file is created (like for code 2) # 6 = actual grid at the end of each month # 7 = actual grid at the end of each year # 8 = write immediately after reading the grid from file and filling missing values. This is used for an automated filling of missing values only. Should not be used productive # # outputcode (for statistic files for zones or subcatchments) # # the Codes behind the names of the statistic files have the meaning of: # <1000 : no output # 1 : spatial mean values for the entire basin, averaged in time over intervals (timesteps) # 2 : spatial mean values for all zones (subbasin) and for the entire basin, averaged in time over intervals (timesteps) # 3 : spatial means for the entire basin, added up in time over intervals (timesteps) # 4 : spatial means for all zones (subbasin) and for the entire basin, added up in time over intervals (timesteps) # 5 : spatial means for the entire basin and for those subbasins which are specified in the output-list, averaged in time over intervals # 6 : spatial means for the entire basin and for those subbasins which are specified in the output-list, added up in time over intervals # # example: # 2001 = per timestep for all subcatchments (and for the entire basin) one (spatially averaged) value, # 2004 = each 4 time steps one averaged value over the last 4 time steps for all subcatchments and for the entire basin, # 4024 = Sums of the mean subcatchment/entire basin values of the timesteps over 24 timesteps (e.g. daily rain sums for subcatchments), # 3120 = averaged values (over 120 time steps!) only for the entire basin (spatially averaged) # 5012 = averaged values (over 12 timesteps) as spatial averages for the entire basin and for each of the subbasins specified in the output-list [output_list] 1 # number of subbasins which are scheduled for output (is only of interest, if the code for the statistic files are >5000) 10 [output_interval] 24 # increment of time steps until an output to the screen is done (24 = each day one output, if time steo = 1h) 1 # warning level for interpolation (no station within search radius) 0 # unit of routed discharge (0=mm/timestep, 1=m3/s) 0 # minutes from the hour-entry in the input data files until the end # of the time step is reached 0 if the end of time step is given like "84 01 01 01", # but it should be $time if the begin is given like in "84 01 01 00" WriteAsciiGrids = 1 # 0 if grids should be written in WaSiM native format, 1 if in ESRI ASCII format InitialStateDirectory = $InitialStateDirectory # if using this parameter, all state grids as well as the storage_richards.ftz file will be expected in that directory for reading DefaultOutputDirectory = $DefaultOutputDirectory # this is the default output directory, all output is written to unless the given filename contains an absolute path (starting with either / or ~ for UNIX or a drive letter and :\ for Windows # there are some exceptions, though: for external coupling no default output path is used # relative pathnames may be used as well. # for compatibility reasins with older control files and WaSiM versions, both directories will only be used if the given filename has no absolute path, # so in order to use the new features, all $outpath uses should be reviewed and removed if necessary (or the variable should be set to an empty string) [coordinates] 47.4 # geogr. latitude (center of the basin -> for radiation calculations) 9.2 # geogr. longitude (center of the basin) 15.0 # meridian according to the official time (middle europe: 15)(east: 0 ... +180 degree, west: 0 ... -180 (or 360 ... 180) 1 # time shift of Meteo-data-time with respect to the true local time (mean sun time) # e.g.: if meteo-data are stored in UTC-time and the time meridian is 15 east (central europe), # than the local time is 1 hour later than the time in the meteo-data-file, so 1 hour has to be added to the time from this file # this is important for calculation of sunshine duration and radiation [region_transition_distance] 10000 # in m [soil_surface_groundwater_substeps]. 1 # number of sub time steps for the module group surface routing, unsaturated zone model and groundwater model (and accumulation of real evapotranspiration) # Values to start with are 1 (default), 2 (half of the common time step), 3 etc. # Please be carefull to set too high values here since the model performance will go down dramatically, since unsatzonmodel and surface routing are called each time! [elevation_model] $inpath_grid//$elevation_model # grid with the digital elevation data [zone_grid] $inpath_grid//$zone_grid # grid with Zone codes $set $lai_grid = lai_//$grid//.//$suffix $set $z0_grid = z0_//$grid//.//$suffix $set $root_grid = root_//$grid//.//$suffix $set $rse_grid = rse_//$grid//.//$suffix $set $rsi_grid = rsi_//$grid//.//$suffix $set $rsc_grid = rsc_//$grid//.//$suffix $set $albedo_grid = albedo_//$grid//.//$suffix $set $vcf_grid = vcf_//$grid//.//$suffix $set $lai_stat = lai_//$grid//.//$code//$year $set $z0_stat = z0_//$grid//.//$code//$year $set $root_stat = root_//$grid//.//$code//$year $set $rse_stat = rse_//$grid//.//$code//$year $set $rsi_stat = rsi_//$grid//.//$code//$year $set $rsc_stat = rsc_//$grid//.//$code//$year $set $albedo_stat = albedo_//$grid//.//$code//$year $set $vcf_stat = vcf_//$grid//.//$code//$year # there is a simple possibility starting with WaSiM 8.10.03 to do the nearest neighbor filling permanently: simply set the writecode for the standardgrid to 8 and the grid # will be writen to the default output directory with it's original name but an additional suffix "filled". Once thsi grid is written, it can be converted to binary optionally and #used as input grid (without fillcode = 1 then). [standard_grids] 20 # number of standard grids # path # identification # fillcode 0=no, 1=yes (fill missing values with values of nearest neighbor) # $inpath_grid//$RelCellArea_grid RelCellArea 1 # grid with land use data # $inpath_grid//$CellSizeX_grid CellSizeX 1 # grid with cellsize in x-direction for each cell (in meter) # $inpath_grid//$CellSizeY_grid CellSizeY 1 # grid with cellsize in y-direction for each cell (in meter) $inpath_grid//$regio_grid regression_regions 1 # region grid if using multiple regression perameter files for meteorological data interpolation $inpath_grid//$slope_grid slope_angle 1 # grid with slope angle data $inpath_grid//$aspect_grid slope_aspect 1 # grid with slope aspect data $inpath_grid//$subcatchments zonegrid_soilmodel 1 # zone grid for the runoff generation model (and unsaturated zone model) $inpath_grid//$land_use_grid landuse fillcode = 1 # writecode = 8 readcode = 1 outname = $outpath//$land_use_grid # grid with land use data (will be written out after reading in for getting gthe filles values) $inpath_grid//$soil_types soil_types fillcode = 1 # writecode = 8 readcode = 1 outname = $outpath//$soil_types # soil types as codes for the soil table $inpath_grid//$flow_time_grid flow_times fillcode = 1 # writecode = 8 readcode = 1 outname = $outpath//$flow_time_grid # grid with flow times for surface runoff to the subbasin outlet # $inpath_grid//$ice_firn_grid ice_firn 0 # grid with firn or ice cells (code 0: nodata values should not be replaced by nearest neighbour) # $inpath_ini//$lai_grid leaf_area_index1 fillcode = 1 defaultValue = 3 writecode = 0 readcode = 0 outname = $outpath//$lai_grid statfile = $outpath//$lai_stat statcode = $once_per_interval # $inpath_ini//$z0_grid RoughnessLength1 fillcode = 1 defaultValue = 0.1 writecode = 0 readcode = 0 outname = $outpath//$z0_grid statfile = $outpath//$z0_stat statcode = $once_per_interval # $inpath_ini//$root_grid root_depth1 fillcode = 1 defaultValue = 1.0 writecode = 0 readcode = 0 outname = $outpath//$root_grid statfile = $outpath//$root_stat statcode = $once_per_interval # $inpath_ini//$vcf_grid vegetation_coverage_degree1 fillcode = 1 defaultValue = 0.9 writecode = 0 readcode = 0 outname = $outpath//$vcf_grid statfile = $outpath//$vcf_stat statcode = $once_per_interval # $inpath_ini//$rse_grid SurfaceEvaporationResistance fillcode = 1 defaultValue = 300 writecode = 0 readcode = 0 outname = $outpath//$rse_grid statfile = $outpath//$rse_stat statcode = $once_per_interval # $inpath_ini//$rsi_grid SurfaceIntercepResistance1 fillcode = 1 defaultValue = 5 writecode = 0 readcode = 0 outname = $outpath//$rsi_grid statfile = $outpath//$rsi_stat statcode = $once_per_interval # $inpath_ini//$rsc_grid SurfaceCanopyResistance1 fillcode = 1 defaultValue = 75 writecode = 0 readcode = 0 outname = $outpath//$rsc_grid statfile = $outpath//$rsc_stat statcode = $once_per_interval $inpath_ini//$albedo_grid albedo fillcode = 1 defaultValue = 0.2 writecode = 0 readcode = 0 outname = $outpath//$albedo_grid statfile = $outpath//$albedo_stat statcode = $once_per_interval # $inpath_grid//$lake_grid lake_codes fillcode = 0 # grid with a unique code for each lake # $inpath_grid//$max_pond_grid max_ponding_storage fillcode = 1 defaultValue = 0 # grid with height of small dams around the fields for water ponding (in m). 0 if no ponding occurs. For a call which is active in the lake grid, this value is the theoretical value when the pond overflows. $inpath_grid//$T_lower_boundary_grid T_Lower_Boundary_Condition fillcode = 0 # defaultValue = 5 writecode = 8 readcode = 0 $inpath_grid//$debris_on_glaciers debris_on_glaciers fillcode = 0 # $inpath_grid//$river_depth_grid river_depth 1 # grid with the depth of all streams in the stream network in m $inpath_grid//$river_width_grid river_width 1 # grid with the witdh of all streams in m $inpath_grid//$river_links_grid river_links 0 # grid with codes of tributaries, from which a channel was routed (only for real routing channels!!!) $inpath_grid//$kolmationsgrid kolmation 1 # grid with codes of tributaries, from which a channel was routed (only for real routing channels!!!) $inpath_grid//$aquiferthick1 aquifer_thickness_1 fillcode = 1 # writecode = 8 readcode = 1 outname = $outpath//$aquiferthick1 # grid with thickness of first aquifer (m from soil surface to the aquifer bottom) $inpath_grid//$gw_storage_coeff_1 gw_storage_coeff_1 fillcode = 1 # writecode = 8 readcode = 1 outname = $outpath//$gw_storage_coeff_1 # storage coefficients for 1. aquifer $inpath_grid//$gw_bound_h_1_grid gw_boundary_fix_h_1 0 # periodicity = 1 D 12 persistent = 0 # boundary conditions 1 constant head for layer 1 $inpath_grid//$gw_bound_q_1_grid gw_boundary_fix_q_1 0 # boundary conditions 2 (given flux perpendicular to the border) for layer 1 $inpath_grid//$gw_kx_1_grid gw_k_x_1 fillcode = 1 # writecode = 8 readcode = 1 outname = $outpath//$gw_kx_1_grid # lateral hydraulic conductivities for the 1. aquifer in x direction $inpath_grid//$gw_ky_1_grid gw_k_y_1 fillcode = 1 # writecode = 8 readcode = 1 outname = $outpath//$gw_ky_1_grid # lateral hydraulic conductivities for the 1. aquifer in y direction # $inpath_grid//$gw_kolmation_1 gw_kolmation_1 1 # kolmation (leakage factor) between 1st and 2nd aquifer # $inpath_grid//$aquiferthick2 aquifer_thickness_2 1 # grid with thickness of first aquifer (m from soil surface to the aquifer bottom) # $inpath_grid//$gw_storage_coeff_2 gw_storage_coeff_2 1 # storage coefficients for 1. aquifer # $inpath_grid//$gw_bound_h_2_grid gw_boundary_fix_h_2 0 # boundary conditions 1 constant head for layer 1 # $inpath_grid//$gw_bound_q_2_grid gw_boundary_fix_q_2 0 # boundary conditions 2 (given flux perpendicular to the border) for layer 1 # $inpath_grid//$gw_kx_2_grid gw_k_x_2 1 # lateral hydraulic conductivities for the 1. aquifer in x direction # $inpath_grid//$gw_ky_2_grid gw_k_y_2 1 # lateral hydraulic conductivities for the 1. aquifer in y direction # $inpath_grid//$gw_kolmation_2 gw_kolmation_2 1 # kolmation (leakage factor) between 2nd and 3rd aquifer # $inpath_grid//$drain_depth_grid drainage_depth 1 # grid with depth of drainage pipes in the soil # $inpath_grid//$drain_distance_grid drainage_distance 1 # grid with distances of the drainage pipes or hoses from each other # $inpath_grid//$clay_depth_grid clay_depth 1 # grid with the depth of an unpermeable layer (0 if no clay layer exists # $inpath_grid//$irrigationcodes irrigation_codes 1 # grid with codes according to the irrigation table # $inpath_grid//$taucrit_grid CriticalShearStress 1 # # $inpath_grid//$ThawCoeffPermaFrost ThawCoeffPermaFrost 0 # grid with coefficients for a simple permafrost thawing model (nodata if no permafrost soil is present, else a suiteable alpha value) # $inpath_grid//$T_lower_boundary_grid T_Lower_Boundary_Condition fillcode = 0 # defaultValue = 5 writecode = 8 readcode = 0 # $inpath_grid//$FlowDirection_grid FlowDirection 1 # grid with flow directions from tanalys # $inpath_grid//$hydr_cond_grid hydraulic_conductivity 1 # grid with hydraulic conductivity of the soil -> old soilmodel # $inpath_grid//$field_capacity_grid available_soil_moisture 1 # grid with available soil moisture at field capacity [mm] -> old soil model # $inpath_grid//$ATBgrid topographic_faktor 1 # soil-topograhic-factor ln(A/(T*tanb)) # $inpath_grid//$tracer_1 concflux_tracer_1_input 1 # $inpath_grid//$tracer_2 concflux_tracer_2_input 1 # $inpath_grid//$tracer_3 concflux_tracer_3_input 1 # $inpath_grid//$tracer_4 concflux_tracer_4_input 1 # $inpath_grid//$tracer_5 concflux_tracer_5_input 1 # $inpath_grid//$tracer_6 concflux_tracer_6_input 1 # $inpath_grid//$tracer_7 concflux_tracer_7_input 1 # $inpath_grid//$tracer_8 concflux_tracer_8_input 1 # $inpath_grid//$tracer_9 concflux_tracer_9_input 1 # variable grids are used by more than one module or can be changed (like albedo and soil storage) $set $SurfStorSiltingUp = sfstsu//$grid//.//$suffix $set $pondgridtopmodel = pond_top//$grid//.//$suffix $set $VegetationStart = vegstart//$grid//.//$suffix $set $VegetationStop = vegstop//$grid//.//$suffix $set $VegetationDuration = vegduration//$grid//.//$suffix [variable_grids] 0 # Number of variable grids to read $outpath//$etr2etpgrid ETR2ETP 1 1 # effectice for wasim-richards only: ETR/ETP fraction, used for dynamic irrigation amount modelling in irrigation method 4 $Writegrid # effectice for wasim-richards only 0 # effectice for wasim-richards only $outpath//$pondgridtopmodel ponding_storage_top 0 0 # effectice for wasimtop only: pond grid for lake modelling, nodata values grid must not be filled $Writegrid # effectice for wasimtop only: Writegrid for topmodel-ponds 0 # effectice for wasimtop only: 0, if ponds should be initialized in routing model by the volume-waterlevel relation, 1 if actual pond content should be read in from existing pond grid # $outpath//$albedo albedo 1 0 # albedo; for time without snow derived from land use data # $Writegrid # Writegrid for $albedo # $readgrids # 0, if albedo is derived from land use at model start time, 1, if albedo is read from file # $outpath//$glacierizedCells_grid GlacierizedCells 0 -9999 # glacierized fraction of each cell (0...1, -9999 for all-time non-glacierized cells) when using the dynamic glacier model; wasim will check if there are only nodata. If yes, the _ice_firn_ grid will be used for initialization of the glacier cells # $Writegrid # Writegrid for glacerized cells # $readgrids # should always be 1 since otherwise no glacier would be created # $outpath//$glacier_codes_grid GlacierCodes 0 -9999 # codes for each single glacier. This grid is required when using the dynamic glacier model. It separates multiple glaciers even in the same subbasin for a applying the V-A-relation correctly # $Writegrid # Writegrid for glacier codes # $readgrids # should always be 1 since otherwise no glacier zones could be created in the dynamic glacier model $outpath//$VegetationStart VegetationStart1 0 -1 # JD for start of vegetation period (is set to actual JD when Landusetable indicates the JD for start of vegetation is reached); $Writegrid # Writegrid for $VegetationStart $readgrids # 0, will only be read in when a simulation starts within the year somewhen $outpath//$VegetationStop VegetationStop1 0 -1 # JD for end of vegetation period (is set to actual JD when Landusetable indicates the JD for the end of vegetation is reached); $Writegrid # Writegrid for $VegetationStop $readgrids # 0, will only be read in when a simulation starts within the year somewhen $outpath//$VegetationDuration VegetationDuration1 0 -1 # Daycount for actual vegetation period; $Writegrid # Writegrid for $VegetationDuration $readgrids # 0, will only be read in when a simulation starts within the year somewhen $outpath//$soilstoragegrid soil_storage 1 0 # soil water storage $Writegrid # Writegrid for this grid $readgrids # 0, if soil_storage should be derived from soil types, 1, if it should be read from file $outpath//$SurfStorSiltingUp SurfStorSiltingUp 1 0 # storage for surface runoff which was routed into other grid cells but not into a cell with a river $Writegrid # Writegrid for this grid $readgrids # 0, if soil_storage should be derived from soil types, 1, if it should be read from file $outpath//$forcingunitsgrid1 SumOfForcingUnits1 0 -1 # Sum of forcing units until phenological cycle starts $Writegrid # Writegrid for this grid 0 # 0, if forcing units will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$TStartPhenoGrid1 Pheno_start1 0 -1 # actual starting day as calculated by forcing units sum $Writegrid # Writegrid for this grid 0 # 0, if TStart-Day units will be initialized to -1, otherwise it will be read in from a file (what for?) $outpath//$chillingunitsgrid1 SumOfChillingUnits1 0 -1 # Sum of chilling units until DP2_t1_dorm is reached -> FStar is calculated dependent on this values $Writegrid # Writegrid for this grid $DPreadgrids # 0, if chilling units will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$FStargrid1 FStar_ForcingThreshold1 0 -1 # FStar value to be reached by the sum of forcing untis until dynamic phenology starts (only used by Method 4 in Landuse) $Writegrid # Writegrid for this grid 0 # 0, if FStar will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$forcingunitsgrid2 SumOfForcingUnits2 0 -1 # Sum of forcing units until phenological cycle starts $Writegrid # Writegrid for this grid 0 # 0, if forcing units will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$TStartPhenoGrid2 Pheno_start2 0 -1 # actual starting day as calculated by forcing units sum $Writegrid # Writegrid for this grid 0 # 0, if TStart-Day units will be initialized to -1, otherwise it will be read in from a file (what for?) $outpath//$chillingunitsgrid2 SumOfChillingUnits2 0 -1 # Sum of chilling units until DP2_t1_dorm is reached -> FStar is calculated dependent on this values $Writegrid # Writegrid for this grid $DPreadgrids # 0, if chilling units will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$FStargrid2 FStar_ForcingThreshold2 0 -1 # FStar value to be reached by the sum of forcing untis until dynamic phenology starts (only used by Method 4 in Landuse) $Writegrid # Writegrid for this grid 0 # 0, if FStar will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$forcingunitsgrid3 SumOfForcingUnits3 0 -1 # Sum of forcing units until phenological cycle starts $Writegrid # Writegrid for this grid 0 # 0, if forcing units will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$TStartPhenoGrid3 Pheno_start3 0 -1 # actual starting day as calculated by forcing units sum $Writegrid # Writegrid for this grid 0 # 0, if TStart-Day units will be initialized to -1, otherwise it will be read in from a file (what for?) $outpath//$chillingunitsgrid3 SumOfChillingUnits3 0 -1 # Sum of chilling units until DP2_t1_dorm is reached -> FStar is calculated dependent on this values $Writegrid # Writegrid for this grid $DPreadgrids # 0, if chilling units will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$FStargrid3 FStar_ForcingThreshold3 0 -1 # FStar value to be reached by the sum of forcing untis until dynamic phenology starts (only used by Method 4 in Landuse) $Writegrid # Writegrid for this grid 0 # 0, if FStar will be initialized to 0, otherwise it will be read in from a file (what for?) $outpath//$tavg_grid TemperatureAVG 1 0 # Average Temperature for a day (will be updated each day at the last interval (if time step is smaller than 1d) $Writegrid # Writegrid for this grid 0 # readgrid (not necessary for single day average temperature $outpath//$ThawDepthGrid PermafrostThawDepth 0 1 # grid with depth of thawed soil for permafrost soils, initialised with 0 in this case (counts positive downwards) $Writegrid # Writegrid for this grid 0 # readgrid -> 0 if grid is not read in, 1 if grid will be read in $outpath//$SnowFreeDaysGrid SnowFreeDaysGrid 0 0 # grid with number of effective snow-free days for permafrost soils (even if there is snow, snow free days will be reset only after a certain number of snow cover days is reached) $Writegrid # Writegrid for this grid 0 # readgrid -> 0 if grid is not read in, 1 if grid will be read in $outpath//$SnowCoverDaysGrid SnowCoverDaysGrid 0 31 # grid with number of snow cover days $Writegrid # Writegrid for this grid 0 # readgrid -> 0 if grid is not read in, 1 if grid will be read in [model_time] $starthour # start hour $startday # start day $startmonth # start month $startyear # start year $endhour # end hour $endday # end day $endmonth # end month $endyear # end year [meteo_data_count] 8 [meteo_names] # the name of the temperature interpolation result is mandatory if dynamic phenology is used (calculating forcing units depends on a grid called "temperature") temperature_reg1 temperature_reg2 precipitation_reg1 precipitation_reg2 wind_speed #air_humidity vapor_pressure global_radiation sunshine_duration # methods: # 1 = idw # 2 = regress # 3 = idw+regress # 4 = thiessen # 5 = bilinear # 6 = bilinear gradients and residuals linarly combined # 7 = bicubic spline # 8 = bicubic splines of gradients and residuals linearly combined # 9 = read grids according to the name in a grid list file # 10 = regression from Stationdata instead from outputfile of regr.exe (similar to method 1, except that no station selection may be applied) # 11 = regression and IDW from station data (equivalent to method 3, except that no station selection may be applied) # 12 = Thiessen with given lapse rate (as single next line parameter or with multiple parameters: lower lapse rate, upper limit, upper lapse rate, type (P-type or T-type, important for continuous or discontinuous data modelling)) [temperature_reg1] 10 # methods, see comments above $inpath_meteo//t2m_reg1//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//t2m_reg1_//$year//.out # file name with regression data (if method = 2 or 3) 820 1400 200 1 300 # lower inversion [m asl], upper inversion [m asl], tolerance [m], overlap [0/1 for true/false], clusterlimit [m] $outpath//$tempegrid1 # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $outpath//$tempegrid1 # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) 5//$Writegrid # 0, if no grid-output is needed, else one of the codes described above 1.0 # correction faktor for results $outpath//t2m_reg1_//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.2 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell -65 # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) 0.8 # ratio of the short to the long axis of the anisotropy-ellipsis -40 # lower limit of interpolation results -40 # replace value for results below the lower limit 40 # upper limit for interpolation results 40 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 1 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 4 # number of scenario cells [temperature_reg2] 10 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//t2m_reg2_//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//t2m_reg2_//$year//.out # file name with regression data (if method = 2 or 3) 900 1500 300 1 150 # lower inversion [m asl], upper inversion [m asl], tolerance [m], overlap [0/1 for true/false], clusterlimit [m] $outpath//$tempegrid1 # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $outpath//$tempegrid2 # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) 5//$Writegrid # 0, if no grid-output is needed, else one of the codes described above 1.0 # correction faktor for results $outpath//t2m_reg2_//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.2 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell 25 # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) 0.6 # ratio of the short to the long axis of the anisotropy-ellipsis -40 # lower limit of interpolation results -40 # replace value for results below the lower limit 40 # upper limit for interpolation results 40 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 1 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 4 [wind_speed] 3 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined, 7=bicubic spline, 8=bicubic splines of gradients and residuals linearly combined, 9=read grids according to the name in a grid list file, 10=regression from Stationdata, 11=regression and IDW from station data $inpath_meteo//wind__//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) $inpath_meteo//wind__//$year//.out # file name with regression data (if method = 2 or 3) $outpath//$windgrid # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $Writegrid # 0, if no grid-output is needed, else one of the codes described above 0.1 # correction faktor for results $outpath//wind//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.3 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell $Anisoslope # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) $Anisotropie # ratio of the short to the long axis of the anisotropy-ellipsis 0 # lower limit of interpolation results 0 # replace value for results below the lower limit 90 # upper limit for interpolation results 90 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 3 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 4 # number of scenario cells [precipitation_reg1] 10 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined, 7=bicubic spline, 8=bicubic splines of gradients and residuals linearly combined, 9=read grids according to the name in a grid list file, 10=regression from Stationdata, 11=regression and IDW from station data $inpath_meteo//prec_reg1_//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1,3,4,5,6,7,8 or 9 else ignored) 820 1400 200 1 300 # lower inversion [m asl], upper inversion [m asl], tolerance [m], overlap [0/1 for true/false], clusterlimit [m] $outpath//$preci_grid1 # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above 0.1 # correction faktor for results $outpath//prec_reg1_//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.75 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell $Anisoslope # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) $Anisotropie # ratio of the short to the long axis of the anisotropy-ellipsis 0.1 # lower limit of interpolation results 0 # replace value for results below the lower limit 900 # upper limit for interpolation results 900 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 2 # 3 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 1 # 4 # number of scenario cells [precipitation_reg2] 4 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined, 7=bicubic spline, 8=bicubic splines of gradients and residuals linearly combined, 9=read grids according to the name in a grid list file, 10=regression from Stationdata, 11=regression and IDW from station data $inpath_meteo//prec_reg2_//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) # $inpath_meteo//prec_reg1.out # file name with regression data (if method = 2 or 3) 700 1400 400 1 300 # lower inversion [m asl], upper inversion [m asl], tolerance [m], overlap [0/1 for true/false], clusterlimit [m] $outpath//$tempegrid1 # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $outpath//$preci_grid2 # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above 0.1 # correction faktor for results $outpath//prec_reg2_//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.75 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell $Anisoslope # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) $Anisotropie # ratio of the short to the long axis of the anisotropy-ellipsis 0.1 # lower limit of interpolation results 0 # replace value for results below the lower limit 900 # upper limit for interpolation results 900 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 2 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 1 # number of scenario cells [sunshine_duration] 1 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//ssd___//$year//.rel AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) $inpath_meteo//ssd___//$year//.out # file name with regression data (if method = 2 or 3) $outpath//$sunshinegrid # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $Writegrid # 0, if no grid-output is needed, else one of the codes described above 1.0 # correction faktor for results $outpath//ssd_//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.5 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell $Anisoslope # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) $Anisotropie # ratio of the short to the long axis of the anisotropy-ellipsis 0 # lower limit of interpolation results 0 # replace value for results below the lower limit 1.0 # upper limit for interpolation results 1.0 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 3 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 1 # number of scenario cells [global_radiation] 2 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//glob__//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) $inpath_meteo//glob__//$year//.out # file name with regression data (if method = 2 or 3) $outpath//$radiationgrid # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $Writegrid # 0, if no grid-output is needed, else one of the codes described above 1.0 # correction faktor for results $outpath//rad_//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 9998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.5 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell $Anisoslope # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) $Anisotropie # ratio of the short to the long axis of the anisotropy-ellipsis 0 # lower limit of interpolation results 0 # replace value for results below the lower limit 1367 # upper limit for interpolation results 1367 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 1 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 4 # number of scenario cells [air_humidity] 2 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//humi__//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) $inpath_meteo//humi__//$year//.out # file name with regression data (if method = 2 or 3) $outpath//$humiditygrid # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $Writegrid # 0, if no grid-output is needed, else one of the codes described above 0.001 # correction faktor for results $outpath//humi//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 9998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.5 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell $Anisoslope # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) $Anisotropie # ratio of the short to the long axis of the anisotropy-ellipsis 0.01 # lower limit of interpolation results 0.01 # replace value for results below the lower limit 1.0 # upper limit for interpolation results 1.0 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 3 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 1 # number of scenario cells [vapor_pressure] 2 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//vapr__//$year//.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) $inpath_meteo//vapr__//$year//.out # file name with regression data (if method = 2 or 3) $outpath//$vaporgrid # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages) $Writegrid # 0, if no grid-output is needed, else one of the codes described above 0.1 # correction faktor for results $outpath//vapo//$grid//.//$code//$year $once_per_interval # file name for the statistic output (statially averaged values per time step and subcatchment...) 998 # error value: all data in the input file greater than this values or lesser the negative value are nodata $IDWweight # weighting of the reciprocal distance for IDW 0.5 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell $Anisoslope # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive) $Anisotropie # ratio of the short to the long axis of the anisotropy-ellipsis 0 # lower limit of interpolation results 0 # replace value for results below the lower limit 90 # upper limit for interpolation results 90 # replace value for results with larger values than the upper limit $SzenUse # 1=use scenario data for correction, 0=dont use scenarios 1 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 4 # number of scenario cells 699000 235000 0.994 1.187 1.021 1.035 1.201 1.021 .635 .566 .538 1.021 .800 1.007 699000 290000 1.021 1.201 1.035 1.049 1.201 1.035 .649 .566 .538 1.035 .787 1.021 737000 235000 0.980 1.173 1.035 1.076 1.173 1.021 .621 .635 .593 1.035 .800 1.007 737000 290000 1.007 1.187 1.049 1.090 1.173 1.035 .649 .649 .607 1.049 .800 1.007 # ---------- parameter for model components ----------------- # [RegionalSuperposition] 1 $time NumberOfEntities = 2; temperature { entityinputgrid = temperature_reg1 ; regions = 1 2 ; weights = 1.0 0.0 ; entityinputgrid = temperature_reg2 ; regions = 1 2 ; weights = 0.0 1.0 ; outputgrid = $outpath//$tempegrid ; writecode = 5//$Writegrid ; outputtable = $outpath//t2m_//$grid//.//$code//$year; statcode = $once_per_interval; } precipitation { entityinputgrid = precipitation_reg1 ; regions = 1 2 ; weights = 1.0 0.0 ; entityinputgrid = precipitation_reg2 ; regions = 1 2 ; weights = 0.0 1.0 ; outputgrid = $outpath//$preci_grid ; writecode = 1//$Writegrid ; outputtable = $outpath//prec//$grid//.//$code//$year; statcode = $once_per_interval; } # for precipitation correction the paragraphs "precipitation" "temperature" and # "wind_speed" are searched in the memory. If thea are not there (no definition in the control file for precipitation, wind or temperature), # the prec. corr. will not be calculated [precipitation_correction] 1 # 0=ignore this module, 1 = run the module 0.0 # Snow-rain-temperature 1.05 # liquid: b in: y = p(ax + b) 0.05 # liquid: a in: y = p(ax + b) = 1% more per m/s + 0.5% constant 1.20 # Snow: b in: y = p(ax + b) 0.25 # Snow: a in: y = p(ax + b) = 15% more per m/s + 45% constant # correction factors for direct radiation are calculated # if the cell is in the shadow of another cell, or if a cell is not in the sun (slope angle!) # then the factor is 0. # control_parameter: 1 = radiation correction WITH shadow WITHOUT temperature correction # 2 = radiation correction WITH shadow WITH temperature correction # 3 = radiation correction WITHOUT shadow WITHOUT temperature correction, # 4 = radiation correction WITHOUT shadow WITH Temperatur [radiation_correction] 1 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 2 # control parameter for radiation correction (see above) $outpath//$Tcorrgrid # name of the grids with the corrected temperatures $Writegrid # Writegrid for corrected temperatures 5 # factor x for temperature correction x * (-1.6 .... +1.6) $outpath//$ExpoCorrgrid # name of the grids with the correction factors for the direct radiation $Writegrid # Writegrid $outpath//$Shapegrid # name of the grids for codes 1 for theor. shadow, 0 for theor. no shadow (day; assumed: SSD=1.0) $Writegrid # Writegrid 1 # interval counter, after reaching this value, a new correction is calculated (3=all 3 hours a.s.o.) 1 # Splitting of the interval, usefull for time step=24 hours (then: split=24, -> each hour one correction calculation) [evapotranspiration] 1 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 1 # Method: 1=Penman-Monteith, 2=Hamon (only daily), 3=Wendling (only daily) 4= Haude (only daily) 0.2 0.2 0.35 0.4 0.4 0.4 0.4 0.4 0.35 0.2 0.2 0.2 # PEC correction factor for HAMON-evapotranspiration 0.20 0.20 0.21 0.29 0.29 0.28 0.26 0.25 0.22 0.22 0.20 0.20 # fh (only for method 4: Haude) monthly values (Jan ... Dec) (here: for Grass) 0.5 # fk -> factor for Wendling-evapotranspiration (only for Method = 3) $outpath//$ETPgrid # result grid for pot. evapotranspiration in mm/dt 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//etp_//$grid//.//$code//$year $once_per_interval # statisticfile for Teilgebiete of pot. evapo-Transpiration $outpath//$ETRgrid # result grid for real evapotranspiration in mm/dt 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//etr_//$grid//.//$code//$year $once_per_interval # statistic for subcatchments (zones) of the real evapotranspiration $outpath//$EVAPgrid # result grid for real evapotranspiration in mm/dt 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//evap//$grid//.//$code//$year $once_per_interval # statistic for subcatchments (zones) of the potential evaporation $outpath//$EVARgrid # result grid for real evapotranspiration in mm/dt 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//evar//$grid//.//$code//$year $once_per_interval # statistic for subcatchments (zones) of the real evaporation $outpath//$ETRSgrid # result grid for real snow evapotranspiration in mm/dt 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//etrs//$grid//.//$code//$year $once_per_interval # statistic for subcatchments (zones) of the real snow evaporation $outpath//$EIPgrid # result grid for pot. interception evaporation in mm/dt 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//eip_//$grid//.//$code//$year $once_per_interval # statisticfile for zones of pot. interception evaporation $outpath//rgex//$grid//.//$code//$year $once_per_interval # statistic for subcatchments (zones) of the corrected radiation +0.23 +1.77 -2.28 +1.28 # coefficients c for Polynom of order 3 RG = c1 + c2*SSD + c3*SSD^2 + c4*SSD^3 +0.072 -0.808 +2.112 -0.239 # coefficients x for Polynom of order 3 SSD = x1 + x2*RG + x3*RG^2 + x4*RG^3 0.88 0.05 # Extinktion coefficient for RG-modeling (Phi and dPhi) (summer phi = phi-dphi, winter phi=phi+dphi) 1654.0 # recession constant (e-function for recession of the daily temperature amplitude with altitude [m] 3.3 4.4 6.1 7.9 9.4 10.0 9.9 9.0 7.8 6.0 4.2 3.2 # monthly values of the max. daily T-amplitudes (for 0 m.a.s.l) 0.62 0.1 # part of the temperature amplitude (dt), that is added to the mean day-temperature # (followed by the range of changing within a year ddt) to get the mean temperature of light day # in the night: mean night temperature is mean day temperature minus (1-dt)*(temp. amplitude) [snow_model] 1 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 2 # method 1=T-index, 2=t-u-index, 3=Anderson comb., 4=extended com. 1.0 # transient zone for rain-snow (T0R +- this range) 0.0 # T0R temperature limit for rain (Grad Celsius) 0.0 # T0 temperature limit snow melt 0.1 # CWH storage capacity of the snow for water (relative part) 1.0 # CRFR coefficient for refreezing 1.8 # C0 degree-day-factor mm/d/C 0.8 # C1 degree-day-factor without wind consideration mm/(d*C) 0.17 # C2 degree-day-factor considering wind mm/(d*C*m/s) 0.07 # z0 roughness length cm for energy bilance methods (not used) 1.0 # RMFMIN minimum radiation melt factor mm/d/C comb. method 2.5 # RMFMAX maximum radiation melt factor mm/d/C comb. method 0.45 # Albedo for snow (Min) 0.90 # Albedo for snow (Max) $outpath//$rain_rate # rain rate 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//rain//$grid//.//$code//$year $once_per_interval # rain rate $outpath//$snow_rate # snow rate 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//snow//$grid//.//$code//$year $once_per_interval # snow rate $outpath//$days_snow # days with snow (SWE > 5mm) $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//sday//$grid//.//$code//$year $once_per_interval # days with snow (SWE > 5mm) $outpath//$snow_age # snow age (days without new snow) $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//sage//$grid//.//$code//$year $once_per_interval # days since last snowfall $outpath//albe//$grid//.//$code//$year $once_per_interval # Albedo $outpath//$snowcover_outflow # discharge from snow, input (precipitation) for following modules $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//qsch//$grid//.//$code//$year $once_per_interval # melt flow (or rain, if there is no snow cover) in mm/dt $outpath//$melt_from_snowcover # discharge from snow, input (precipitation) for following modules $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//qsme//$grid//.//$code//$year $once_per_interval # melt flow in mm/dt $outpath//$SSNOgrid # name of the grids with the snow storage solid in mm $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//$SLIQgrid # name of the grids with the snow storage liquid in mm $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//ssto//$grid//.//$code//$year $once_per_interval # total snow storage, in mm, (liquid and solid fraction) $outpath//$SSTOgrid # name of the grids with the total snow storage solid AND liquid in mm $Writegrid # 0, if no grid-output is needed, else one of the codes described above $readgrids # 1=read snow storage solid, liquid grids from disk, 0=generate new grids [ice_firn] 2 # method for glacier melt: 1=classical t-index, 2=t-index with correction by radiation, 11 = dynamic glacier model with classical t-index, 12 = dynamic glacier model with radiation correction 6.0 # t-index factor for ice 5.0 # t-index factor for firn 4.0 # t-index factor for snow 1.8 # melt factor -0.00010 # radiation coefficient for ice_min (for method 2/12) +0.00070 # radiation coefficient for ice_max (for method 2/12) +0.00006 # radiation coefficient for snow_min (for method 2/12) +0.00030 # radiation coefficient for snow_max (for method 2/12) 12 # els-konstante for ice 120 # els-konstante for firn 24 # els-konstante for snow 0.01 # initial reservoir content for ice discharge (single linear storage approach) 0.01 # initial reservoir content for firn discharge (single linear storage approach) 0.01 # initial reservoir content for snow discharge (single linear storage approach) $outpath//$firn_melt # melt from firn $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//qfir//$grid//.//$code//$year $once_per_interval # melt from firn as statistic file $outpath//$ice_melt # melt from ice $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//qice//$grid//.//$code//$year $once_per_interval # melt from ice as statistic file $outpath//qglc//$grid//.//$code//$year $once_per_interval # melt from ice and firn as statistic file # ----------------------------------------------------------------------------- # now some new parameters for the new dynamic glacier model (methods 11 and 12) $outpath//qsgl//$grid//.//$code//$year $once_per_interval # melt from snow from glacier only as statistic file (but still with respect tothe subbasins areas!) --> new in version 8.07.00 $readgrids # 1=read grids and stacks from disk, 0=generate new grids and stacks (using the parameters in the following line for WE_Firn stack) 7 2500 1.8 # number of layers for the firn stack, followed by two initialization parameters: average Equilibrium line elevation in m (e.g. 2500) and change rate of WE per m in mm (e.g. 2) -> every 100m the WE of firn in each layer will grow by 200mm 09 30 # month and day (hour is set automatically to 24) for which the Volume-Area-Relation will be applied newly (and temporary (i.e. internal) Balances are reset to 0) 28.5 1.36 10 4 # VAscaling and VAexponent for Volume-Area-Relation of glaciers and number of iterations (elevation belts) and extraWeightFactorBand0 (elevation band 0 will be processed in each iteration this given number of times more than once. Default = 0) $outpath//$firn_WE_stack # water equivalent for firn (given as stack, number of layers taken from the parameter given before); layer 0 will contain the total WE for all firn layers $Writestack # 0, if no grid-output is needed, else one of the codes described above $outpath//glfirn//$grid//.//$code//$year $once_per_interval # water equivalent for firn as statistics file (sum over all firn layers) $outpath//$GlacierMassBalance # output grid with mass balance of the glacier 3 # 3: write at end of simulation (important to start another model run with correct initialization values) $outpath//$OldGlacierMassBalance # output grid with mass balance of the glacier 3 # 3: write at end of simulation (important to start another model run with correct initialization values) $outpath//glmb//$grid//.//$code//$year $once_per_interval # mass balance for the glaciers as statistics file (mass balance for each time step with respect to the entire subbasin the glaciers are located in) $outpath//glmb2//$grid//.//$code//$year $once_per_interval # mass balance for the glaciers as statistics file (mass balance for each time step with respect to the glaciers only!) 1.0 # additional parameter when using a debris grid: this value is used to globally scale the values of the debris grid. Only values > nodata are regarded, i.e. when a cells melt coefficient should not be altered, the debris grid should contain -9999 at this location # permafrost parameter # note: # - parameter alpha must be read in as a grid with valid cells marked by an alpha value > 0 (all other cells must be nodata, NOT 0) # - two grids are used within the mode: SnowCoverDaysGrid and SnowFreeDaysGrid. If these grids should be initialized, they must be read in as variable grid # otherwise they will be generated internally (and cannot be written) # - parameters are then: minimum number of days with snowcover, after which the soild will fereeze (happens suddenly - this is NOT # a refreezing model, only a state change in order to initialize the next thawing period # - minimum SWE (snow water equivalent) to be counted as snow cover days [permafrost] 1 # method: 1=simple Alpha*sqrt(snow-free-days) approach to estimate thawdepth 30 # number of days with snow cover after which the soil is assumed to be froozen again 5 # maximum snow water equivalent for the interval to be counted as snow covered (then, the snow-cover-days grid will be incremented by the length of an interval [interception_model] 1 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 1 # method: 1 = use ETP for calculating EI; 2 = use EIP for calculating EI (only effective for method 1 in evapotranspiration model -> for other methods, ETP = EIP) $outpath//$throughfall # result grid : = outflow from the interception storage $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//qi__//$grid//.//$code//$year $once_per_interval # statistic file interception storage outflow $outpath//$EIgrid # Interzeption evaporation, grid 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//ei__//$grid//.//$code//$year $once_per_interval # zonal statistic $outpath//$SIgrid # storage content of the interception storage 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//si__//$grid//.//$code//$year $once_per_interval # zonal statistic For interception storage content 0.35 # layer thickness of the waters on the leaves (multiplied with LAI -> storage capacity) $readgrids # 1=read grids from disk, else generate internal [infiltration_model] 0 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes $outpath//$INFEXgrid # grid with infiltration excess in mm (surface runoff) $Writegrid # Writegrid for surface discharge (fraction 1) $outpath//infx//$grid//.//$code//$year $once_per_interval # statistic file for the infiltration excess $outpath//$SATTgrid # grid with code 1=saturation at interval start, 0 =no saturation. $Writegrid # Writegrid for saturation code grids 0.1 # fraction of reinfitrating water (of the infiltration excess) $set $SDISPgrid = sdis//$grid//.//$suffix $set $RPAUSgrid = paus//$grid//.//$suffix $set $EKIN_grid = ekin//$grid//.//$suffix $set $TSBB_grid = tsbb//$grid//.//$suffix $set $QDSU_grid = qdsu//$grid//.//$suffix [SiltingUpModel] 0 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 1 # method: 1=traditional (default if this line is missing), 2=read eight regresion parameters for individual control over i0, ie and Cv, 3=use free expressions $outpath//sdis//$grid//.//$code//$year $once_per_interval # statistics for silting up disposition (Verschlämmungsneigung) $outpath//qdsu//$grid//.//$code//$year $once_per_interval # direct discharge from silting up module $outpath//$SDISPgrid # grid with actual silting up disposition $Writegrid # writegrid for this grid $outpath//$RPAUSgrid # grid with actual rain pause length (for getting ekin for events longer than a time step and for regeneration of soil) $Writegrid # writegrid for this grid $outpath//$EKIN_grid # grid with actual kinetic energy of the event $Writegrid # writegrid for this grid $outpath//$TSBB_grid # grid with actual time since last soil tillage $Writegrid # writegrid for this grid $outpath//$QDSU_grid # grid with direct runoff from silting up model (will be used in unsatzonmodel!) $Writegrid # writegrid for this grid 1 2 3 4 5 6 7 8 9 10 11 12 13 # range for subbasin codes 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 # minimum rainy break to separate two precipitation events (in days) $readgrids # readgrid code 0 do not read, 1 = read grids 65.1 # for method 2: parameter A in I_0 = A (initial infiltration capacity, in method 0 defined as 65 mm/h) 12.21 # for method 2: parameter B in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 12.2) 0.521 # for method 2: parameter C in I_end = B*(dg^C)*(fd^D) (in method 0 defined as 0.52) -0.641 # for method 2: parameter D in I_end = B*(dg^C)*(fd^D) (in method 0 defined as -0.64) 0.0131 # for method 2: parameter E in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.013) -1.031 # for method 2: parameter F in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -1.03) 0.71 # for method 2: parameter G in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as 0.7) -0.191 # for method 2: parameter H in Cv = E*(fd^F)*(dg^G)*(t_cult^H) (in method 0 defined as -0.19) SiltingUpExpressions { # please read the short documentation on the expression parser below the expression list definition W = ((P>0.05)&(P<76.2))*(11.89+8.73*log10(Abs(P+0.001))) + (P>=76.2)*28.33; # this is the energy for the actual rain intensity P; the two terms are valid for 0.05=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 only for compatibility here. numfiles = 2; outputfile { header = glacierdata; filename = $outpath//special_output_glaciers.//$code//$year; entity { ID = GlacierMassBalance; Symbol = GMB; Xcoords = 771371, 801115, 771211; Ycoords = 214666, 194848, 164323; } entity { ID = melt_from_firn; Symbol = Mfirn; Xcoords = 771371, 801115, 771211; Ycoords = 214666, 194848, 164323; } entity { ID = melt_from_ice; Symbol = Mice; Xcoords = 771371, 801115, 771211; Ycoords = 214666, 194848, 164323; } } outputfile { header = soildata; filename = $outpath//special_output_soilmodel.//$code//$year; entity { ID = theta; Symbol = th; Xcoords = 748503, 748503, 748503, 748503, 748503, 748503, 770572, 770572, 770572, 770572, 770572, 770572; Ycoords = 196127, 196127, 196127, 196127, 196127, 196127, 256698, 256698, 256698, 256698, 256698, 256698; Level = 1, 2, 3, 4, 5, 6, 1, 2, 3, 4, 5, 6; } entity { ID = hydraulic_heads; Symbol = hh; Xcoords = 748503, 748503, 748503, 748503, 748503, 748503, 770572, 770572, 770572, 770572, 770572, 770572; Ycoords = 196127, 196127, 196127, 196127, 196127, 196127, 256698, 256698, 256698, 256698, 256698, 256698; Level = 1, 2, 3, 4, 5, 6, 1, 2, 3, 4, 5, 6; } }