# == protected part to be changed by the module adapter == # == all manual changes in this section will be overwritten by the module adapter == $set $runID = 0073_1991 $set $startyear = 1991 $set $endyear = 1991 $set $mainpath = d:\Daten\Calibration\ $set $InitialStateDirectory = $mainpath//StateIni_$runID//\ $set $DefaultOutputDirectory = $mainpath//Output_$runID//\ $set $inpath_grid = $mainpath//input\grids\ $set $inpath_meteo = $mainpath//input\meteo\ $set $inpath_hydro = $mainpath//input\hydro\ $set $inpath_ini = $mainpath//input\ini\ $set $time = 60 $set $year = 96 $set $starthour = 01 $set $startday = 01 $set $startmonth = 01 $set $endhour = 24 $set $endday = 31 $set $endmonth = 12 # 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 # == end of protected part == $set $time = 60.0 # it is important to set $outpath to an empty string in order to activate $DefaultOutputDirectory $set $outpath = # variables for parameters in unsatzon model: since several subbasins will be parameterized with identical parameters # it is very convenient to have them here defined as variables, so we have only 18 parameter sets instead of 194 (for each subbasin one set) # kd --> recession constant for single linear reservoir for direct runoff $set $kd1 = 3 $set $kd2 = 3 $set $kd3 = 3 $set $kd4 = 3 $set $kd5 = 3 $set $kd6 = 3 $set $kd7 = 3 $set $kd8 = 3 $set $kd9 = 3 $set $kd10 = 3 $set $kd11 = 3 $set $kd12 = 3 $set $kd13 = 3 $set $kd14 = 3 $set $kd15 = 3 $set $kd16 = 3 $set $kd17 = 3 $set $kd18 = 3 $set $kd19 = 3 $set $kd20 = 3 $set $kd21 = 3 $set $kd22 = 3 $set $kd23 = 3 $set $kd14gl = 3 # 30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46, (70,71) $set $kd4gl = 3 # 47,48 $set $kd8gl = 3 # 49,50 $set $kd9gl = 3 # 51,52,53,54,55, (68) $set $kd5gl = 3 # 56 $set $kd16gl = 3 # 57,58,59,60,61,62 $set $kd23gl = 3 # 63 $set $kd11gl = 3 # 64,65,66 $set $kd20gl = 3 # 67 $set $kd68gl = 3 # wie dr9gl $set $kd70gl = 3 # wie dr14gl # ki --> recession constant for single linear reservoir for interflow $set $ki1 = 6 $set $ki2 = 6 $set $ki3 = 6 $set $ki4 = 6 $set $ki5 = 6 $set $ki6 = 6 $set $ki7 = 6 $set $ki8 = 6 $set $ki9 = 6 $set $ki10 = 6 $set $ki11 = 6 $set $ki12 = 6 $set $ki13 = 6 $set $ki14 = 6 $set $ki15 = 6 $set $ki16 = 6 $set $ki17 = 6 $set $ki18 = 6 $set $ki19 = 6 $set $ki20 = 6 $set $ki21 = 6 $set $ki22 = 6 $set $ki23 = 6 $set $ki14gl = 6 # 30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46, (70,71) $set $ki4gl = 6 # 47,48 $set $ki8gl = 6 # 49,50 $set $ki9gl = 6 # 51,52,53,54,55, (68) $set $ki5gl = 6 # 56 $set $ki16gl = 6 # 57,58,59,60,61,62 $set $ki23gl = 6 # 63 $set $ki11gl = 6 # 64,65,66 $set $ki20gl = 6 # 67 $set $ki68gl = 6 # wie dr9gl $set $ki70gl = 6 # wie dr14gl # dr --> drainage density (interflow generation parameter) $set $dr1 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr2 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr3 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr4 = 1.2 # 4.8 # 2.4 # 1.2 $set $dr5 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr6 = 3.6 # 14.4 # 7.2 # 3.6 $set $dr7 = 3.6 # 14.4 # 7.2 # 3.6 $set $dr8 = 3.6 # 14.4 # 7.2 # 3.6 $set $dr9 = 3.6 # 14.4 # 7.2 # 3.6 $set $dr10 = 2.8 # 11.2 # 5.6 # 2.8 $set $dr11 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr12 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr13 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr14 = 1.2 # 4.8 # 2.4 # 1.2 $set $dr15 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr16 = 4.0 # 16.0 # 8 # 4.0 $set $dr17 = 0.6 # 2.4 # 1.2 # 0.6 $set $dr18 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr19 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr20 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr21 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr22 = 3.6 # 14.4 # 7.2 # 3.6 $set $dr23 = 1.8 # 7.2 # 3.6 # 1.8 $set $dr14gl = 1.2 # 4.8 # 2.4 # 1.2 # 30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46, (70,71) $set $dr4gl = 1.2 # 4.8 # 2.4 # 1.2 # 47,48 $set $dr8gl = 3.6 # 14.4 # 7.2 # 3.6 # 49,50 $set $dr9gl = 3.6 # 14.4 # 7.2 # 3.6 # 51,52,53,54,55, (68) $set $dr5gl = 1.8 # 7.2 # 3.6 # 1.8 # 56 $set $dr16gl = 2.8 # 11.2 # 5.6 # 2.8 # 57,58,59,60,61,62 $set $dr23gl = 1.8 # 7.2 # 3.6 # 1.8 # 63 $set $dr11gl = 3.6 # 14.4 # 7.2 # 3.6 # 64,65,66 $set $dr20gl = 3.6 # 14.4 # 7.2 # 3.6 # 67 $set $dr68gl = 1.8 # 7.2 # 3.6 # 1.8 # wie dr9gl $set $dr70gl = 1.8 # 7.2 # 3.6 # 1.8 # wie dr14gl # sdf --> Snow melt: Direct Runoff fraction $set $sdf1 = 0.05 $set $sdf2 = 0.05 $set $sdf3 = 0.05 $set $sdf4 = 0.05 $set $sdf5 = 0.05 $set $sdf6 = 0.05 $set $sdf7 = 0.05 $set $sdf8 = 0.05 $set $sdf9 = 0.05 $set $sdf10 = 0.05 $set $sdf11 = 0.05 $set $sdf12 = 0.05 $set $sdf13 = 0.05 $set $sdf14 = 0.05 $set $sdf15 = 0.05 $set $sdf16 = 0.05 $set $sdf17 = 0.05 $set $sdf18 = 0.05 $set $sdf19 = 0.05 $set $sdf20 = 0.05 $set $sdf21 = 0.05 $set $sdf22 = 0.05 $set $sdf23 = 0.05 $set $sdf14gl = 0.05 # 30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46, (70,71) $set $sdf4gl = 0.05 # 47,48 $set $sdf8gl = 0.05 # 49,50 $set $sdf9gl = 0.05 # 51,52,53,54,55, (68) $set $sdf5gl = 0.05 # 56 $set $sdf16gl = 0.05 # 57,58,59,60,61,62 $set $sdf23gl = 0.05 # 63 $set $sdf11gl = 0.05 # 64,65,66 $set $sdf20gl = 0.05 # 67 $set $sdf68gl = 0.05 # wie dr9gl $set $sdf70gl = 0.05 # wie dr14gl $set $grid = beo500 $set $stack = beo500 $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_LakesConst $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 # use_c stands for resampling using central value instead of most frequent value, which would be use_f $set $ice_firn_grid = $grid//.glc $set $field_capacity_grid = $grid//.nfk $set $ATBgrid = $grid//.atb $set $hydr_cond_grid = $grid//.k $set $soil_types = $grid//.soil_origCodes $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//.dep05all $set $river_width_grid = $grid//.wit3all $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//.kol2e-6 $set $gw_kx_1_grid = $grid//.kx1-4 $set $gw_kx_2_grid = $grid//.kx2 $set $gw_kx_3_grid = $grid//.kx3 $set $gw_ky_1_grid = $grid//.ky1-4 $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_OrigCodes $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//.lakes $set $taucrit_grid = $grid//.tau $set $ThawCoeffPermaFrost = $grid//.alpha $set $debris_on_glaciers = $grid//.debris $set $slidefraction1 = $grid//.sd1 $set $slidefraction2 = $grid//.sd2 $set $slidefraction3 = $grid//.sd3 $set $slidefraction4 = $grid//.sd4 $set $depositionindex = $grid//.svfdir_2 $set $elevationordercols = $grid//.eoc $set $elevationorderrows = $grid//.eor # output 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 $SSNOOnGlacgrid = ssno_onGlac_//$grid//.//$suffix $set $SLIQOnGlacgrid = sliq_onGlac_//$grid//.//$suffix $set $SSTOOnGlacgrid = ssto_onGlac_//$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 $set $ts_avg = ts_avg//$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 $set $Temperaturestack = tsoil//$stack//.//$suffix # result-grids for sow model $set $snowtemperaturgrid = snowtemp//$grid//.//$suffix $set $snowsurftemperaturgrid = ssrftemp//$grid//.//$suffix $set $inputmassgrid = inpmass//$grid//.//$suffix $set $mobilemassgrid = mobmass//$grid//.//$suffix $set $depositiongrid = deposit//$grid//.//$suffix $set $snowtemperaturgridGL = snowtempgl//$grid//.//$suffix $set $snowsurftemperaturgridGL = ssrftempgl//$grid//.//$suffix $set $ETRSgridGl = etrsgl//$grid//.//$suffix $set $QSNOWOnGlacgrid = qsmegl//$grid//.//$suffix # parameters for interpolation of meteorological input data $set $SzenUse = 0 $set $IDWmaxdist = 200000 $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 = 2001 $set $sum_per_24Invs = 4001 $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 # SpinUp is a new run mode that tries to initialize many grids and stacks like soil moisture, soil temperatures, firn on glaciers and snow during a multi-year spinup period. # to be as fast as possible, the basin is divided into HRUs (hydrologic respond units). So no lateral connections will be regarded here (i.e. no groudnwater model and surface routing can be done). # during SpinUp, lakes, groundwater and glacier resizing is not computed. # After [SpinUp] 1 # doSpinUp: 0=no, 1=yes 1 # doTemperaturePreSpinUp: 0=no, 1=yes 0 # repeatCnt: how many times should the model run for the defined time period (in [model_time]), e.g. 5 for 5 tiumes 5 # numGrids : numbe rof following entries for grids to use for HRU creation. 0 means, that a (standard) grid with identifier HydrologicRespondUnits should be read in, otherwise the follwoing lines will be evaluated elevation_model 100 lin # 1st parameter: internal grid identifier | 2nd parameter: class width for reclassifying, e.g. 50 means that 0...50 in the input will be mapped to 1 in the outpout, (50-100] will be mapped to 2 in the output etc | 3rd parameter lin or lg2: how input is transformed first before reclassifying soil_types 1 lin # each landuse class is used for HRU creation landuse 1 lin # each soil type is used for HRU creation max_ponding_storage 1 lin # each lake will have at least one specific HRU slope_angle 2 sqrt # sqrt means, that reclassifying is done with quadratic growing class widths, i.e. the input slopes in degree will be classified as 0-1, 1-2, 2-4, 4-8, 8-16, 16-32, 32-64, > 64 slope_aspect 90 lin # each soil type is used for HRU creation [output_list] 23 # number of subbasins which are scheduled for output (is only of interest, if the code for the statistic files are >5000) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 70 71 [output_interval] 25 # 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 # 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] 46.1 # geogr. latitude (center of the basin -> for radiation calculations) 7.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 = alb//$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] 24 # number of standard grids # path # identification # fillcode 0=no, 1=yes (fill missing values with values of nearest neighbor) $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//$regio_grid regression_regions fillcode = 1 # defaultValue = 1 writecode = 8 readcode = 0 outname = $outpath//$regio_grid # region grid if using multiple regression perameter files for meteorological data interpolation #$inpath_ini//$lai_grid leaf_area_index1 fillcode = 2 defaultValue = 3 writecode = 8 readcode = 0 outname = $outpath//$lai_grid statfile = $outpath//$lai_stat statcode = $sum_per_24Invs #$inpath_ini//$z0_grid RoughnessLength1 fillcode = 2 defaultValue = 0.1 writecode = 8 readcode = 0 outname = $outpath//$z0_grid statfile = $outpath//$z0_stat statcode = $sum_per_24Invs #$inpath_ini//$root_grid root_depth1 fillcode = 2 defaultValue = 1.0 writecode = 8 readcode = 0 outname = $outpath//$root_grid statfile = $outpath//$root_stat statcode = $sum_per_24Invs #$inpath_ini//$vcf_grid vegetation_coverage_degree1 fillcode = 2 defaultValue = 0.9 writecode = 8 readcode = 0 outname = $outpath//$vcf_grid statfile = $outpath//$vcf_stat statcode = $sum_per_24Invs #$inpath_ini//$rse_grid SurfaceEvaporationResistance fillcode = 2 defaultValue = 300 writecode = 8 readcode = 0 outname = $outpath//$rse_grid statfile = $outpath//$rse_stat statcode = $sum_per_24Invs #$inpath_ini//$rsi_grid SurfaceIntercepResistance1 fillcode = 2 defaultValue = 5 writecode = 8 readcode = 0 outname = $outpath//$rsi_grid statfile = $outpath//$rsi_stat statcode = $sum_per_24Invs #$inpath_ini//$rsc_grid SurfaceCanopyResistance1 fillcode = 2 defaultValue = 75 writecode = 8 readcode = 0 outname = $outpath//$rsc_grid statfile = $outpath//$rsc_stat statcode = $sum_per_24Invs # $inpath_grid//$albedo_grid albedo fillcode = 1 defaultValue = 0.2 writecode = 8 readcode = 0 outname = $outpath//$albedo_grid # statfile = $outpath//$albedo_stat statcode = $sum_per_24Invs $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//$ice_firn_grid ice_firn 0 # grid with firn or ice cells (code 0: nodata values should not be replaced by nearest neighbour) $inpath_grid//$subcatchments zonegrid_soilmodel 1 # zone grid for the runoff generation model (and unstaurated zone model) $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//$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//$lake_grid lake_codes fillcode = 0 # grid with a unique code for each lake $inpath_grid//$river_depth_grid river_depth 0 # grid with the depth of all streams in the stream network in m $inpath_grid//$river_width_grid river_width 0 # 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 fillcode = 1 # defaultValue = 1e-5 writecode = 8 readcode = 0 #outname = $outpath//$kolmationsgrid # grid with codes of tributaries, from which a channel was routed (only for real routing channels!!!) $inpath_grid//$aquiferthick1 aquifer_thickness_1 fillcode = 1 # defaultValue = 15 writecode = 8 readcode = 0 #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 # defaultValue = 0.1 writecode = 8 readcode = 0 #outname = $outpath//$gw_storage_coeff_1 # storage coefficients for 1. aquifer $inpath_grid//$gw_bound_h_1_grid gw_boundary_fix_h_1 fillcode = 0 # defaultValue = -9999 writecode = 8 readcode = 0 #outname = $outpath//$gw_bound_h_1_grid # 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 fillcode = 0 # defaultValue = -9999 writecode = 8 readcode = 0 #outname = $outpath//$gw_bound_q_1_grid # boundary conditions 2 (given flux perpendicular to the border) for layer 1 $inpath_grid//$gw_kx_1_grid gw_k_x_1 fillcode = 1 # defaultValue = 1e-5 writecode = 8 readcode = 0 #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 # defaultValue = 1e-5 writecode = 8 readcode = 0 #outname = $outpath//$gw_ky_1_grid # lateral hydraulic conductivities for the 1. aquifer in y direction $inpath_grid//$gw_kolmation_1 gw_kolmation_1 fillcode = 1 # defaultValue = 1e-5 writecode = 8 readcode = 0 # kolmation (leakage factor) between 1st and 2nd aquifer $inpath_grid//$soil_types soil_types fillcode = 1 # defaultValue = 1 writecode = 8 readcode = 0 outname = $outpath//$soil_types # soil types as codes for the soil table $inpath_grid//$slidefraction1 slidefraction1 fillcode = 0 # new snow model extensions: fraction of the area flowing to the North (slides) $inpath_grid//$slidefraction2 slidefraction2 fillcode = 0 # new snow model extensions: fraction of the area flowing to the West (slides) $inpath_grid//$slidefraction3 slidefraction3 fillcode = 0 # new snow model extensions: fraction of the area flowing to the East (slides) $inpath_grid//$slidefraction4 slidefraction4 fillcode = 0 # new snow model extensions: fraction of the area flowing to the South (slides) $inpath_grid//$depositionindex deposition_index fillcode = 1 # new snow model extensions: correction factor for wind impact on snow fall $inpath_grid//$debris_on_glaciers debris_on_glaciers fillcode = 0 # # 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] 3 # Number of variable grids to read $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; 7=am Jahresende in der Form glc_a500_20051231.024 1 #$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 7=am Jahresende in der Form glida500_20051231.024 1 #$readgrids # should always be 1 since otherwise no glacier zones could be created in the dynamic glacier model $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//$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//$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) 3 # 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 3 # 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] 6 [meteo_names] temperature precipitation wind_speed air_humidity global_radiation sunshine_duration vapor_pressure # 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] 11 # methods, see comments above $inpath_meteo//T_1990-2013.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//t2m_reg1.out # file name with regression data (if method = 2 or 3) 1200 1500 300 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.3 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell -35 # 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 [precipitation] 11 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//P_all_SMN1+KTN+MM+NIME_1990-2013.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) #0.0003 2070 0.0001 P-type # method 12 parameter: altitudinal gradient in (mm/d)/m --> 1mm per 100m with thrshold elevation to an upper gradient and type of data (P-type for intermittend data. T-type for continuous data) 500 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_grid # 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 1.0 # correction faktor for results $outpath//prec//$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.8 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell -35 # 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 0.01 # 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 [precipitation_reg1] 11 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//P_all_1990-2013.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) #0.0003 2070 0.0001 P-type # method 12 parameter: altitudinal gradient in (mm/d)/m --> 1mm per 100m with thrshold elevation to an upper gradient and type of data (P-type for intermittend data. T-type for continuous data) 1200 1500 300 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_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 1.0 # 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.8 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell -35 # 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 0.01 # 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 [precipitation_reg2] 11 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//P_all_1990-2013.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) #0.0003 2070 0.0001 P-type # method 12 parameter: altitudinal gradient in (mm/d)/m --> 1mm per 100m with thrshold elevation to an upper gradient and type of data (P-type for intermittend data. T-type for continuous data) 1200 1500 300 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 1.0 # 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.8 # for interpolation method 3: relative weight of IDW-interpolation in the result $IDWmaxdist # max. distance of stations to the actual interpolation cell 45 # 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 0.01 # 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 [wind_speed] 11 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//u_1990-2013.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//wind9610.out # file name with regression data (if method = 2 or 3) 1200 1500 300 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//$windgrid # 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//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.50 # 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 [global_radiation] 11 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//RG_1990-2013.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//glob8009.out # file name with regression data (if method = 2 or 3) 1200 1500 300 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//$radiationgrid # 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//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 [sunshine_duration] 11 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//SSD_1990-2013.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//ssd_9610.out # file name with regression data (if method = 2 or 3) 1200 1500 300 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//$sunshinegrid # 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//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.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 # 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 [air_humidity] 11 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//RH_1990-2013.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//humi8009.out # file name with regression data (if method = 2 or 3) 1200 1500 300 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//$humiditygrid # 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 0.01 # 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.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.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] 10 # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined $inpath_meteo//e_1990-2013.dat AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored) #$inpath_meteo//vapo9610.out # file name with regression data (if method = 2 or 3) 1200 1500 300 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//$vaporgrid # 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//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.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 1 # 1=add scenarios, 2=multiply scenarios, 3=percentual change 4 # number of scenario cells # ---------- parameter for model components ----------------- [RegionalSuperposition] 0 $time NumberOfEntities = 1; 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; } 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_correction] 1 # 0=ignore this module, 1 = run the module 0.7 # Snow-rain-temperature 1.04 # liquid: b in: y = p(ax + b) 0.04 # liquid: a in: y = p(ax + b) = 1% more per m/s + 0.5% constant 1.15 # Snow: b in: y = p(ax + b) 0.20 # Snow: a in: y = p(ax + b) = 15% more per m/s + 45% constant # corretion 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 5//$Writegrid # Writegrid for corrected temperatures 8 # scaling factor for temperature correction $outpath//$ExpoCorrgrid # name of the grids with the correction factors for the direct radiation 5//$Writegrid # Writegrid $outpath//$Shapegrid # name of the grids for codes 1 for theor. shadow, 0 for theor. no shadow (day; assumed: SSD=1.0) 5//$Writegrid # Writegrid 1 # interval counter, after reaching this value, a new correction is calculated (3=all 3 hours a.s.o.) 1 # Spitting 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.5 0.6 0.8 1.0 1.1 1.1 1.2 1.1 1.0 0.9 0.7 0.5 # 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 $sum_per_24Invs # 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 $sum_per_24Invs # 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 $sum_per_24Invs # 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 $sum_per_24Invs # 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 $sum_per_24Invs # 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 $sum_per_24Invs # statisticfile for zones of pot. interception evaporation $outpath//rgex//$grid//.//$code//$year $sum_per_24Invs # 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 $outpath//cloud//$grid//.//$code//$year $once_per_interval # statistic for subcatchments (zones) of the corrected radiation # snow model methods: # conventional: # 1 = T-index # 2 = T-u-Index # 3 = energy balance approach after Anderson # 4 = extended energy balance approach after Braun (based on Anderson) with observed vapor pressure # extended methods (introduced by Michael Warscher in 2013, implemented into WaSiMs main distribution by J. Schulla, 2014): # 5 = Enhanced energy balance approach (enhEnbal) (with more sophisticated parameter estimations) # 6 = enhEnbal + gravitational slides # 7 = enhEnbal + gravitational slides + wind redistribution # 8 = T-Index + gravitational slides # 9 = T-Index + gravitational slides + wind redistribution # 10 = T-Index + wind redistribution # 11 = enhEnbal + wind redistribution # 12 = T-u-Index + gravitational slides # 13 = T-u-Index + wind redistribution # 14 = T-u-Index + gravitational slides + wind redistribution [snow_model] 1 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 14 # 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.75 #0.7 # T0R temperature limit for rain (Grad Celsius) 0.75 #0.7 # T0 temperature limit snow melt 0.15 # CWH storage capacity of the snow for water (relative part) 0.5 # CRFR coefficient for refreezing 1.8 # C0 degree-day-factor mm/d/C 1.0 # C1 degree-day-factor without wind consideration mm/(d*C) 0.7 # C2 degree-day-factor considering wind mm/(d*C*m/s) 0.07 # z0 roughness length cm for energy bilance methods (not used) 3.5 # RMFMIN minimum radiation melt factor mm/d/C comb. method 5.0 # RMFMAX maximum radiation melt factor mm/d/C comb. method 0.50 #0.45 # Albedo for snow (Min) 0.90 #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) 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//sday//$grid//.//$code//$year $sum_per_24Invs # days with snow (SWE > 5mm) $outpath//$snow_age # snow age (days without new snow) 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//sage//$grid//.//$code//$year $sum_per_24Invs # days since last snowfall $outpath//albe//$grid//.//$code//$year $sum_per_24Invs # Albedo $outpath//$snowcover_outflow # discharge from snow, input (precipitation) for following modules 1//$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 1//$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 1//$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 1//$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 1//$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 # now some new parameters and file names follow for the snow model extensions implemented in 2014 (originally done by M. Warscher and adopted and a little bit extended to the publicly available release in 2014 by J. Schulla) $outpath//$snowtemperaturgrid # result grid with snow pack temperature 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//snowtemp//$grid//.//$code//$year $once_per_interval # temperature of the snow pack (used for new energy balance model) $outpath//$snowsurftemperaturgrid # result grid with 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//ssrftemp//$grid//.//$code//$year $once_per_interval # temperature of the snow surface (used for new energy balance model) $outpath//$inputmassgrid # result grid with snow surface temperature 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//inpmass//$grid//.//$code//$year $once_per_interval # input mass for snow redistribution (gravitational snow slides) $outpath//$mobilemassgrid # result grid with input mass for gravitational snow slides 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//mobmass//$grid//.//$code//$year $once_per_interval # moving/mobile mass for snow redistribution (gravitational snow slides) $outpath//$depositiongrid # result grid with deposition amount for each cell (gravitational snow slides) 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//deposit//$grid//.//$code//$year $once_per_interval # deposition amount (mm) for gravitational slides 55 # maximum deposition slope (0...90) 3 # scaling for maximum deposition (0..inf) 30 # minimum slope for creating slides (0...90) (scale dependent! The coarser the resolution, the smaller this value must be since the average slope decreases). 0.007 # fraction of snow pack that forms the slide (0...1] 0.52 # LWINcorr: correction factor for incoming long wave radiation for fine tuning the energy balance (accouting together with LWOUTcorr for errors in cloudiness and albedo); recommended Values: 0.8...1.2 0.40 # 0.85 # LWOUTcorr: correction factor for outgoing long wave radiation for fine tuning the energy balance (accouting together with LWINcorr for errors in cloudiness and albedo): recommended Values: 0.8...1.2 [ice_firn] 13 # 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, 13 = melt methods are taken from snow model (for snow only) 6.0 # t-index factor for ice 5.0 # t-index factor for firn 4.0 # t-index factor for snow 1.7 # melt factor (when using radiation correction in methods 11...13) -0.00010 # -0.00010 # radiation coefficient for ice_min (for method 2/12/13) +0.00060 # +0.00070 # radiation coefficient for ice_max (for method 2/12/13) +0.00006 # +0.00006 # radiation coefficient for snow and firn min (for method 2/12/13) +0.00030 # +0.00030 # radiation coefficient for snow and firn max (for method 2/12/13) 12 # els-konstante for ice 240 # 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 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//qfir//$grid//.//$code//$year 2001 #$once_per_interval # melt from firn as statistic file $outpath//$ice_melt # melt from ice 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//qice//$grid//.//$code//$year 2001 #$once_per_interval # melt from ice as statistic file $outpath//qglc//$grid//.//$code//$year 2001 #$once_per_interval # discharge from snow, 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 2001 #$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 2800 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 3.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 $Writegrid # 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 $Writegrid # 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 # some new parameters for better statistic output (for balance verfication) and also separate handling of snow on glaciers and beside glaciers on the same cell $outpath//ssto_OffGlac_//$grid//.//$code//$year $once_per_interval # total snow storage, in mm, (liquid and solid fraction) for the unglacierized part of a cell (usefull for annual balances of all inputs/outputs/storages, since snow on glaciers is handled in the glacier mass balance already) $outpath//$SSNOOnGlacgrid # name of the grids with the snow storage solid in mm valid for the glacierized part of a cell 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//$SLIQOnGlacgrid # name of the grids with the snow storage liquid in mm valid for the glacierized part of a cell 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//$SSTOOnGlacgrid # name of the grids with the total snow storage solid AND liquid in mm valid for the glacierized part of a cell 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//$QSNOWOnGlacgrid # name of the grids with the total snow outflow in mm valid for the glacierized part of a cell 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//$snowsurftemperaturgridGL # name of the grids with snow surface temperature, valid for the glacierized part of a cell 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//$snowtemperaturgridGL # name of the grids with snow pack temperature, valid for the glacierized part of a cell 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//$ETRSgridGl # name of the grids with snow evaporation in mm, valid for the glacierized part of a cell 1//$Writegrid # 0, if no grid-output is needed, else one of the codes described above # 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 model: 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 1//$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 $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//ei__//$grid//.//$code//$year $sum_per_24Invs # zonal statistic $outpath//$SIgrid # storage content of the interception storage $Writegrid # 0, if no grid-output is needed, else one of the codes described above $outpath//si__//$grid//.//$code//$year $avrg_per_24Invs # 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 # 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 Rainpause 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 evanulated by WaSiM: # expressions returning W, X, Y and Z are independently of each other. # expression V must be called as last call in any case, since WaSiM will update EKIN internally using the energy-result (in W) and V depends on all the other results W to Z # other expressions for storing intermediate results may be defined at any position in the expression list before the results will be used in another expression [SurfaceRoutingModel] 0 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 2 # method: 1=MultipleFlowPaths for diverging areas, 2=single flowpaths (nearest direction as given by aspect), 3=multiple flow paths dynamically based on water level, 4 = same as 3 but with single flow paths, 5 = simple v calculation without iteration (suited for back water conditions) $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] 1 # 0=ignore this module, 1 = run the module 2 # method for recalculating DHM, # 1 = do not change the DHM, it refects already the ground surface of the lakes, # 2 = use mean_pond_grid to calculate dhm corrections # max_pond_grid will be used for mapping the cells pond content to a lake during model runs - so the lake level may well rise above the normal surface 0.1 # Albedo_OpenWater (will be used only, when the pond is filled with water when calculating potential evaporation -> otherwise, the normal landuse for this cell is referenced for this parameter) # 10 # rsc for water (usage as above) 0.4 # z0 for water (usage as above) # 10.0 # LAI_OpenWater (usage as above) # 1.0 # VCF_OpenWater (usage as above) $readgrids # readgrid code 0 do not read, 1 = read grids --> # if 0, the initial valte for the POND-grid as Volume of Lakes and Reservoirs is set by V0 from the routing description, # if readgrids=1, no initialization in done (POND-Grid is read in) but the Vakt-Value is set by the various grids [unsatzon_model] 1 # 0=ignore this module, 1 = run the module $time # duration of a time step in minutes 3 # method, 1=simple method (will not work anymore from version 7.x), 2 = FDM-Method 3 = FDM-Method with dynamic time step down to 1 secound 1 # controlling interaction with surface water: 0 = no interaction, 1 = exfiltration possible 2 = infiltration and exfiltration possible 1 # 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) 3 # parameter for the initialization of the gw_level (range between 1..levels (standard: 4)) $outpath//qdra//$grid//.//$code//$year $sum_per_24Invs # results drainage discharge in mm per zone $outpath//gwst//$grid//.//$code//$year $once_per_interval # results groundwater depth $outpath//gwn_//$grid//.//$code//$year $sum_per_24Invs # results mean groundwater recharge per zone $outpath//sb05//$grid//.//$code//$year $avrg_per_24Invs # 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 $avrg_per_24Invs # 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 $avrg_per_24Invs # results statistic of the ponding water storage content $outpath//qdir//$grid//.//$code//$year 2001 #$once_per_interval # results statistic of the direct discharge $outpath//qifl//$grid//.//$code//$year 2001 #$once_per_interval # results statistic of the interflow $outpath//qbas//$grid//.//$code//$year 2001 #$once_per_interval # results statistic of the baseflow $outpath//qges//$grid//.//$code//$year 2001 #$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 1//$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) 1//$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) 80 80 # 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 l $outpath//thet//$grid//.//$code//$year # name of a file containing the soil moisture as theta values of the layers of the control point l $outpath//hhyd//$grid//.//$code//$year # name of a file containing the hydraulic head of the layers of the control point l $outpath//otherdata//$grid//.//$code//$year # name of a file containing some other water balance data of the control point (non layer data) l $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) l $outpath//intd//$grid//.//$code//$year # name of a file containing the interflow for the soil layers of the control point 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 30 47 49 51 56 57 63 64 67 68 70 71 #Code $kd1 $kd2 $kd3 $kd4 $kd5 $kd6 $kd7 $kd8 $kd9 $kd10 $kd11 $kd12 $kd13 $kd14 $kd15 $kd16 $kd17 $kd18 $kd19 $kd20 $kd21 $kd22 $kd23 $kd14gl $kd4gl $kd8gl $kd9gl $kd5gl $kd16gl $kd23gl $kd11gl $kd20gl $kd68gl $kd70gl $kd70gl # kelskd $ki1 $ki2 $ki3 $ki4 $ki5 $ki6 $ki7 $ki8 $ki9 $ki10 $ki11 $ki12 $ki13 $ki14 $ki15 $ki16 $ki17 $ki18 $ki19 $ki20 $ki21 $ki22 $ki23 $ki14gl $ki4gl $ki8gl $ki9gl $ki5gl $ki16gl $ki23gl $ki11gl $ki20gl $ki68gl $ki70gl $ki70gl # kelski $dr1 $dr2 $dr3 $dr4 $dr5 $dr6 $dr7 $dr8 $dr9 $dr10 $dr11 $dr12 $dr13 $dr14 $dr15 $dr16 $dr17 $dr18 $dr19 $dr20 $dr21 $dr22 $dr23 $dr14gl $dr4gl $dr8gl $dr9gl $dr5gl $dr16gl $dr23gl $dr11gl $dr20gl $dr68gl $dr70gl $dr70gl # dr 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 # k in qb = Q0 * exp(-k/z) with z = depth to groundwater 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.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.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 # Q0 in the above formula $sdf1 $sdf2 $sdf3 $sdf4 $sdf5 $sdf5 $sdf7 $sdf8 $sdf9 $sdf10 $sdf11 $sdf12 $sdf13 $sdf14 $sdf15 $sdf16 $sdf17 $sdf18 $sdf19 $sdf20 $sdf21 $sdf22 $sdf23 $sdf14gl $sdf4gl $sdf8gl $sdf9gl $sdf5gl $sdf16gl $sdf23gl $sdf11gl $sdf20gl $sdf68gl $sdf70gl $sdf70gl # sdf $readgrids # meanings are extended now! read the follwing comments $outpath//storage_richards.ftz # if readgrids = 1, then this file contains the contents of the flow travel time zones for interflow and surface flow and for the tracers 300 # minimum dynamic time step in secounds. the smaller this number, the longer the model runs but the results will be more accurate due to a maintained Courant condition $outpath//step//$grid//.//$code//$year $avrg_per_24Invs # results statistic of the number of substeps $outpath//$SUBSTEPSgrid # grid with number of substeps --> a good idea is to use writecode 5x (e.g. 53) to get the average number of substeps per cell for the model run $Writegrid # for substeps, the areal distribution is of interest for the annual average value. This is code 6 as first digit in 2-digit codes. Or use 5 for the entire model run [unsatzon_tuning] maxerror = 1e-5 max_iterations = 100 # 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)), 11=heat diffusion 1D with implicit solution method, 12=like 11 including advection (recommended are 11 or 12) 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 mean annual surface temperature and the lapse rate as defined in the next two lines -9.0 # used when no grid "_T_Lower_Boundary_Condition_" was read in only: mean annual surface 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 frozen 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 literature 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)) 0.98 # SE value which must be underrun to evaluate the soil layer as frozen for the thaw depth output grid and statistics 6 # 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) 6//$Writegrid # Writecode for this stack $readgrids # like in all other models: 1= read stack for temperature, 0 = create new stack according to boundary condition )linear interpolation between upper and lower boundary condition) [ExternalCoupling] 0 # 0 = no coupling, 1=coupling [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] 3 # 3 # 0=ignore the module, 1 = run the module, 3=GW-flow follows hydrological subbasins striktly, no border crossing allowed $time # duration of a time step in minutes; doen't change the value unless you have strong reasons to do so!! 1 # solving method: 1=Gauss-Seidel-iteration (using alpha for control wether it is explicite, partly or fully implicite), 2=PCG (not yet implemented 100 # if iterative solving method (1): max.numberof iterations 0.0005 # if iterative solving method (1): max. changes between two iterations 0.5 # Alpha for estimation of central differences 0.5 = Crank-Nicholson Method, 0 = fully explicite, 1 = fully implicite 1.0 # factor for relaxing the iteration if using iterativemethod (successive over[/under] relaxation) $readgrids # 1=read grids for heads from disk, 0=do not read but initialize with gw-level from unsaturated zone 1 # number of layers 112 39 # 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 1//$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 $set $dtmin1 = 10 # Sekundenzeitschritt, in dem Stauseen aktualisiert werden (insbesindere kleine ABs müssen hier kleiner 60 haben um numerisch stabil zu bleiben $set $dtmin2 = 10 # Sekundenzeitschritt, in dem Stauseen aktualisiert werden (insbesindere kleine ABs müssen hier kleiner 60 haben um numerisch stabil zu bleiben [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 1 4800 150 48 # 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 $inpath_hydro//specific_runoff_1990-2013.dat 23 # number of following column descripotr (which column in the spec. disch. file corresponds to which subbasin 1 1 # first number: subbasin, second: column index 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23 72 # timeoffset # internal abstractions with target TG 40 (SP 1 Curnera) TG 12 (AE=118.020, AErel=1.0) from OL 13 (kh=0.1, kv=0.4, Bh= 9.2, Bv= 37.0, Th= 0.92, Mh=25.0, Mv=10.0, I=0.0036, L=13274.0, AE=53.610) TG 33 (AE=4.000, AErel=1.0) from OL 34 (kh=0.1, kv=0.4, Bh= 0.6, Bv= 2.2, Th= 0.20, Mh=25.0, Mv=10.0, I=0.1226, L=1065.7, AE=1.140) TG 32 (AE=5.900, AErel=1.0) from OL 33 (kh=0.1, kv=0.4, Bh= 1.5, Bv= 5.9, Th= 0.30, Mh=25.0, Mv=10.0, I=0.0417, L=965.7, AE=4.000) TG 31 (AE=15.930, AErel=1.0) from OL 32 (kh=0.1, kv=0.4, Bh= 1.5, Bv= 5.9, Th= 0.35, Mh=25.0, Mv=10.0, I=0.0557, L=1931.4, AE=5.900) TG 37 (AE=33.430, AErel=1.0) from OL 35 (kh=0.1, kv=0.4, Bh= 0.9, Bv= 3.7, Th= 0.29, Mh=25.0, Mv=10.0, I=0.0968, L=882.8, AE=3.370) and OL 36 (kh=0.1, kv=0.4, Bh= 0.8, Bv= 3.2, Th= 0.30, Mh=25.0, Mv=10.0, I=0.1410, L=1989.9, AE=3.410) and SP 4 ( file = $outpath//spv_04_sub37.//$code//$year , V0 = 8.6E6, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 ) # Triftsee TG 71 (AE=35.720, AErel=1.0) from OL 45 (kh=0.1, kv=0.4, Bh= 2.2, Bv= 8.7, Th= 0.63, Mh=25.0, Mv=10.0, I=0.0837, L=1789.9, AE=26.610) and SP 5 ( file = $outpath//spv_05_sub71.//$code//$year , V0 = 5.7E6, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 ) # Gruebensee TG 46 (AE=44.500, AErel=1.0) from OL 71 (kh=0.1, kv=0.4, Bh= 1.7, Bv= 7.0, Th= 0.73, Mh=25.0, Mv=10.0, I=0.1749, L=2938.5, AE=35.720) TG 70 (AE=154.50, AErel=1.0) from OL 38 (kh=0.1, kv=0.4, Bh= 0.8, Bv= 3.3, Th= 0.36, Mh=25.0, Mv=10.0, I=0.1876, L=5662.7, AE=5.290) and OL 39 (kh=0.1, kv=0.4, Bh= 1.0, Bv= 4.0, Th= 0.40, Mh=25.0, Mv=10.0, I=0.1585, L=3255.6, AE=7.130) and OL 41 (kh=0.1, kv=0.4, Bh= 1.0, Bv= 4.0, Th= 0.49, Mh=25.0, Mv=10.0, I=0.2448, L=1065.7, AE=12.030) and OL 44 (kh=0.1, kv=0.4, Bh= 4.0, Bv= 16.2, Th= 0.82, Mh=25.0, Mv=10.0, I=0.0415, L=13982.2, AE=59.250) and OL 40 (kh=0.1, kv=0.4, Bh= 1.0, Bv= 4.0, Th= 0.41, Mh=25.0, Mv=10.0, I=0.1689, L=4828.4, AE=7.710) and SP 6 ( file = $outpath//spv_06_sub70.//$code//$year , V0 = 1.92E8, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 ) # KW-Stauseen (5 Stück zusammen) and AL 1 (modus = intern_with_rule_from_reservoir ) # used for hydro power generation weekdays and AL 2 (modus = intern_with_rule_from_reservoir ) # used for hydro power generation weekends TG 43 (AE=7.850, AErel=1.0) from SP 3 ( file = $outpath//spv_03_sub43.//$code//$year , V0 = 10.7E6, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 ) # Engstlensee. TG 14 (AE=454.060, AErel=1.0) from OL 43 (kh=0.1, kv=0.4, Bh= 1.4, Bv= 5.8, Th= 0.39, Mh=25.0, Mv=10.0, I=0.0744, L=15209.5, AE=7.850) and OL 30 (kh=0.1, kv=0.4, Bh= 2.6, Bv= 10.2, Th= 0.48, Mh=25.0, Mv=10.0, I=0.0351, L=15753.8, AE=14.240) and OL 31 (kh=0.1, kv=0.4, Bh= 2.1, Bv= 8.4, Th= 0.51, Mh=25.0, Mv=10.0, I=0.0592, L=20575.1, AE=15.930) and OL 37 (kh=0.1, kv=0.4, Bh= 2.7, Bv= 10.6, Th= 0.68, Mh=25.0, Mv=10.0, I=0.0648, L=14098.2, AE=33.430) and OL 46 (kh=0.1, kv=0.4, Bh= 3.1, Bv= 12.3, Th= 0.75, Mh=25.0, Mv=10.0, I=0.0591, L=10218.3, AE=44.500) and OL 42 (kh=0.1, kv=0.4, Bh= 1.4, Bv= 5.8, Th= 0.28, Mh=25.0, Mv=10.0, I=0.0364, L=14398.2, AE=3.270) and OL 70 (kh=0.1, kv=0.4, Bh= 6.1, Bv= 24.2, Th= 1.17, Mh=25.0, Mv=10.0, I=0.0375, L=16305.3, AE=154.50) and ZL 1 (modus = intern , kh=0.1, kv=0.1, Bh=10.0, Bv=10.0, Th=1.4, Mh=40.0, Mv=40.0, I=0.01, L=100, AE=154.5 ) # TG 70 and ZL 2 (modus = intern , kh=0.1, kv=0.1, Bh=10.0, Bv=10.0, Th=1.4, Mh=40.0, Mv=40.0, I=0.01, L=100, AE=154.5 )# TG 70 TG 4 (AE=556.310, AErel=1.0) from OL 14 (kh=0.1, kv=0.4, Bh=16.4, Bv= 65.4, Th= 1.70, Mh=25.0, Mv=10.0, I=0.0108, L=11666.9, AE=454.060) and OL 48 (kh=0.1, kv=0.4, Bh= 2.2, Bv= 8.8, Th= 0.50, Mh=25.0, Mv=10.0, I=0.0513, L=17385.2, AE=15.200) and OL 47 (kh=0.1, kv=0.4, Bh= 2.0, Bv= 7.8, Th= 0.43, Mh=25.0, Mv=10.0, I=0.0492, L=15371.0, AE=10.560) TG 52 (AE=22.940, AErel=1.0) from OL 51 (kh=0.1, kv=0.4, Bh= 0.9, Bv= 3.4, Th= 0.42, Mh=25.0, Mv=10.0, I=0.2448, L=3272.8, AE=8.020) TG 68 (AE=35.670, AErel=1.0) from OL 55 (kh=0.1, kv=0.4, Bh= 1.4, Bv= 5.7, Th= 0.54, Mh=25.0, Mv=10.0, I=0.1399, L=2838.5, AE=16.300) and OL 54 (kh=0.1, kv=0.4, Bh= 1.0, Bv= 4.2, Th= 0.42, Mh=25.0, Mv=10.0, I=0.1625, L=2979.9, AE=8.510) TG 9 (AE=165.030, AErel=1.0) from OL 52 (kh=0.1, kv=0.4, Bh= 3.9, Bv= 15.7, Th= 0.56, Mh=25.0, Mv=10.0, I=0.0207, L=8194.1, AE=22.940) and OL 53 (kh=0.1, kv=0.4, Bh= 2.3, Bv= 9.2, Th= 0.45, Mh=25.0, Mv=10.0, I=0.0374, L=13439.6, AE=11.660) and OL 68 (kh=0.1, kv=0.4, Bh= 3.7, Bv= 14.6, Th= 0.68, Mh=25.0, Mv=10.0, I=0.0343, L=13198.2, AE=35.670) TG 8 (AE=380.430, AErel=1.0) from OL 49 (kh=0.1, kv=0.4, Bh= 3.6, Bv= 14.2, Th= 0.57, Mh=25.0, Mv=10.0, I=0.0257, L=19992.3, AE=23.140) and OL 9 (kh=0.1, kv=0.4, Bh=10.3, Bv= 41.4, Th= 1.17, Mh=25.0, Mv=10.0, I=0.0127, L=5028.4, AE=165.030) and OL 50 (kh=0.1, kv=0.4, Bh= 4.9, Bv= 19.6, Th= 0.73, Mh=25.0, Mv=10.0, I=0.0220, L=17743.8, AE=45.130) TG 3 (AE=1140.840, AErel=1.0) from OL 4 (kh=0.1, kv=0.4, Bh=28.2, Bv=112.8, Th= 2.82, Mh=25.0, Mv=10.0, I=0.0010, L=23124.3, AE=556.310) and OL 8 (kh=0.1, kv=0.4, Bh=17.8, Bv= 71.3, Th= 1.78, Mh=25.0, Mv=10.0, I=0.0054, L=2989.9, AE=380.430) and SP 2 ( file = $outpath//spv_02_sub03.//$code//$year , V0 = 5.2E9, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 ) # Brienzersee. TG 11 (AE=34.520, AErel=1.0) from OL 66 (kh=0.1, kv=0.4, Bh= 0.9, Bv= 3.7, Th= 0.35, Mh=25.0, Mv=10.0, I=0.1477, L=3014.2, AE=5.260) and OL 65 (kh=0.1, kv=0.4, Bh= 0.5, Bv= 2.1, Th= 0.26, Mh=25.0, Mv=10.0, I=0.2699, L=2997.0, AE=2.350) and OL 64 (kh=0.1, kv=0.4, Bh= 1.1, Bv= 4.2, Th= 0.48, Mh=25.0, Mv=10.0, I=0.2071, L=4528.4, AE=11.500) TG 20 (AE=200.250, AErel=1.0) from OL 11 (kh=0.1, kv=0.4, Bh= 6.6, Bv= 26.5, Th= 0.66, Mh=25.0, Mv=10.0, I=0.0088, L=18219.5, AE=34.520) and OL 67 (kh=0.1, kv=0.4, Bh= 1.5, Bv= 6.0, Th= 0.32, Mh=25.0, Mv=10.0, I=0.0450, L=24913.6, AE=4.730) TG 7 (AE=346.040, AErel=1.0) from OL 20 (kh=0.1, kv=0.4, Bh=12.7, Bv= 50.7, Th= 1.27, Mh=25.0, Mv=10.0, I=0.0092, L=16109.5, AE=200.250) TG 6 (AE=562.650, AErel=1.0) from OL 7 (kh=0.1, kv=0.4, Bh=15.6, Bv= 62.5, Th= 1.56, Mh=25.0, Mv=10.0, I=0.0091, L=13219.5, AE=346.040) and AL 3 (modus = intern 2.5 1.0 30.0 m^3/s ) # parameters are: threshold, fraction, capacity -> there must be an corresponding internal inflow ZL 1 (TG 22)! TG 22 (AE=585.920, AErel=1.0) from OL 6 (kh=0.1, kv=0.4, Bh=18.6, Bv= 74.5, Th= 1.86, Mh=25.0, Mv=10.0, I=0.0094, L=4062.7, AE=562.650) and ZL 3 (modus = intern , kh=0.1, kv=0.1, Bh=8.0, Bv=60.0, Th=2.4, Mh=40.0, Mv=40.0, I=0.01, L=100, AE=562 )# TG 6 and AL 4 (modus = intern 26.0 1.0 150.0 m^3/s ) # parameters are: threshold, fraction, capacity -> there must be an corresponding internal inflow ZL 2 (TG 2)! TG 23 (AE=144.640, AErel=1.0) from OL 10 (kh=0.1, kv=0.4, Bh= 3.4, Bv= 13.6, Th= 0.62, Mh=25.0, Mv=10.0, I=0.0335, L=16199.4, AE=28.780) and OL 63 (kh=0.1, kv=0.4, Bh= 2.0, Bv= 8.2, Th= 0.49, Mh=25.0, Mv=10.0, I=0.0568, L=20147.9, AE=14.160) TG 59 (AE=21.360, AErel=1.0) from OL 58 (kh=0.1, kv=0.4, Bh= 5.9, Bv= 23.6, Th= 0.59, Mh=25.0, Mv=10.0, I=0.0010, L=1600.0, AE=8.590) and SP 11 ( file = $outpath//spv_11_sub59.//$code//$year , V0 = 33.0E6, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 ) # Oeschinensee TG 60 (AE=49.310, AErel=1.0) from OL 57 (kh=0.1, kv=0.4, Bh= 3.1, Bv= 12.3, Th= 0.66, Mh=25.0, Mv=10.0, I=0.0457, L=3738.5, AE=32.050) TG 61 (AE=83.670, AErel=1.0) from OL 60 (kh=0.1, kv=0.4, Bh= 7.2, Bv= 28.7, Th= 0.74, Mh=25.0, Mv=10.0, I=0.0106, L=3148.5, AE=49.310) TG 16 (AE=147.880, AErel=1.0) from OL 59 (kh=0.1, kv=0.4, Bh= 1.9, Bv= 7.8, Th= 0.58, Mh=25.0, Mv=10.0, I=0.0886, L=4772.8, AE=21.360) and OL 61 (kh=0.1, kv=0.4, Bh= 4.8, Bv= 19.0, Th= 0.93, Mh=25.0, Mv=10.0, I=0.0387, L=5087.0, AE=83.670) and OL 62 (kh=0.1, kv=0.4, Bh= 0.7, Bv= 2.9, Th= 0.23, Mh=25.0, Mv=10.0, I=0.1051, L=12078.1, AE=1.800) TG 21 (AE=188.000, AErel=1.0) from OL 16 (kh=0.1, kv=0.4, Bh= 6.0, Bv= 23.8, Th= 1.16, Mh=25.0, Mv=10.0, I=0.0377, L=10191.1, AE=147.880) TG 5 (AE=498.290, AErel=1.0) from OL 23 (kh=0.1, kv=0.4, Bh=10.7, Bv= 42.8, Th= 1.11, Mh=25.0, Mv=10.0, I=0.0107, L=9932.5, AE=144.640) and OL 21 (kh=0.1, kv=0.4, Bh=12.0, Bv= 48.1, Th= 1.22, Mh=25.0, Mv=10.0, I=0.0103, L=11774.0, AE=188.000) and OL 56 (kh=0.1, kv=0.4, Bh= 2.2, Bv= 8.8, Th= 0.49, Mh=25.0, Mv=10.0, I=0.0486, L=19355.0, AE=14.290) TG 2 (AE=2474.240, AErel=1.0) from OL 3 (kh=0.1, kv=0.4, Bh=36.9, Bv=147.7, Th= 3.69, Mh=25.0, Mv=10.0, I=0.0010, L=31862.7, AE=1140.840) and OL 15 (kh=0.1, kv=0.4, Bh= 9.4, Bv= 37.6, Th= 0.94, Mh=25.0, Mv=10.0, I=0.0022, L=29821.3, AE=43.990) and OL 22 (kh=0.1, kv=0.4, Bh=21.1, Bv= 84.3, Th= 2.11, Mh=25.0, Mv=10.0, I=0.0053, L=13794.1, AE=585.920) and OL 5 (kh=0.1, kv=0.4, Bh=19.3, Bv= 77.2, Th= 1.93, Mh=25.0, Mv=10.0, I=0.0061, L=16791.1, AE=498.290) and ZL 4 (modus = intern , kh=0.1, kv=0.1, Bh=8.0, Bv=60.0, Th=2.4, Mh=40.0, Mv=40.0, I=0.01, L=100, AE=585 )# TG 22 and SP 1 ( file = $outpath//spv_01_sub02.//$code//$year , V0 = 6.5E9, C0 = 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0, dTmin = $dtmin1 ) # Thunersee. TG 1 (AE=2957.950, AErel=1.0) from OL 17 (kh=0.1, kv=0.4, Bh= 8.5, Bv= 34.0, Th= 0.85, Mh=25.0, Mv=10.0, I=0.0048, L=28369.2, AE=50.000) and OL 12 (kh=0.1, kv=0.4, Bh=14.3, Bv= 57.3, Th= 1.43, Mh=25.0, Mv=10.0, I=0.0017, L=9476.9, AE=118.020) and OL 18 (kh=0.1, kv=0.4, Bh= 8.2, Bv= 32.9, Th= 0.82, Mh=25.0, Mv=10.0, I=0.0024, L=22899.4, AE=32.180) and OL 19 (kh=0.1, kv=0.4, Bh=12.8, Bv= 51.4, Th= 1.28, Mh=25.0, Mv=10.0, I=0.0017, L=28227.8, AE=89.430) and OL 2 (kh=0.1, kv=0.4, Bh=44.5, Bv=178.0, Th= 4.45, Mh=25.0, Mv=10.0, I=0.0017, L=28869.2, AE=2474.240) # 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) # or reservoir volume in m^3, followed by the abstraction in m^3/s --> to differentiate between discharge # in m^3/s and reservoir content in m^3, the keyword "modus = intern_with_rule" must be extended by the # keyword "_from_reservoir", i.e. intern_with_rule_from_reservoir [abstraction_rule_reservoir_1] # Thunersee #V 01.01. 20.01. 20.02. 01.04. 20.08. 01.10. 15.11. 01.01. # Datum 71 1 20 51 90 231 273 318 366 # Julianischer Tag 0 0 0 0 0 0 0 0 0 # 423.98 (=558 - mittlere Tiefe von 138.02m) 6402970000 0 0 0 0 0 0 0 0 # 556.00 6451470000 5 10 10 5 5 5 5 5 # 557.00 6456320000 8.6 30 30 7.8 7.8 6.8 6.8 8.6 # 557.10 6457290000 9.3 34.0 34.0 8.3 8.3 7.1 7.1 9.3 # 557.12 6458260000 10 38.0 38.0 8.9 8.9 7.5 7.5 10 # 557.14 6459230000 13.3 42.0 42.0 9.4 9.4 7.9 7.9 13.3 # 557.16 6460200000 16.7 46.0 46.0 10 10 8.2 8.2 16.7 # 557.18 6461170000 20.0 50 50 11.3 11.3 8.6 8.6 20.0 # 557.20 6462140000 23.3 53.3 53.3 12.5 12.5 8.9 8.9 23.3 # 557.22 6463110000 26.7 56.7 56.7 13.8 13.8 9.3 9.3 26.7 # 557.24 6464080000 30.0 60.0 60.0 15.0 15.0 9.6 9.6 30.0 # 557.26 6465050000 33.3 63.3 63.3 16.3 16.3 10 10 33.3 # 557.28 6466020000 36.7 66.7 66.7 17.5 17.5 11.2 11.2 36.7 # 557.30 6466990000 40.0 70.0 70.0 18.8 18.8 12.4 12.4 40.0 # 557.32 6467960000 43.3 73.3 73.3 20 20 13.5 13.5 43.3 # 557.34 6468445000 45.0 75 75 20.9 20.9 14.1 14.1 45.0 # 557.35 6468930000 46.7 76.7 76.7 21.8 21.8 14.7 14.7 46.7 # 557.36 6469900000 50.0 80.0 80.0 23.6 23.6 15.9 15.9 50.0 # 557.38 6470870000 54.5 83.3 83.3 25.5 25.5 17.1 17.1 54.5 # 557.40 6471840000 59.1 86.7 86.7 27.3 27.3 18.2 18.2 59.1 # 557.42 6472810000 63.6 90.0 90.0 29.1 29.1 19.4 19.4 63.6 # 557.44 6473295000 65.9 91.7 91.7 30 30 20 20 65.9 # 557.45 6473780000 68.2 93.3 93.3 34.0 34.0 20.6 20.6 68.2 # 557.46 6474750000 72.7 96.7 96.7 42.0 42.0 21.8 21.8 72.7 # 557.48 6475235000 75.0 98.3 98.3 46.0 46.0 22.4 22.4 75.0 # 557.49 6475720000 77.1 100 100 50 50 22.9 22.9 77.1 # 557.50 6476690000 81.3 105.0 105.0 57.1 57.1 24.1 24.1 81.3 # 557.52 6477660000 85.4 110.0 110.0 64.3 64.3 25.3 25.3 85.4 # 557.54 6478630000 89.6 115.0 115.0 71.4 71.4 26.5 26.5 89.6 # 557.56 6479115000 91.7 117.5 117.5 75 75 27.1 27.1 91.7 # 557.57 6479600000 93.8 120.0 120.0 80.0 80.0 27.6 27.6 93.8 # 557.58 6480570000 97.9 125.0 125.0 90.0 90.0 28.8 28.8 97.9 # 557.60 6481055000 100.0 127.5 127.5 95.0 95.0 29.4 29.4 100.0 # 557.61 6481540000 103.3 130.0 130.0 100 100 30 30 103.3 # 557.62 6482510000 110.0 135.0 135.0 107.7 107.7 32.7 32.7 110.0 # 557.64 6483480000 116.7 140.0 140.0 115.4 115.4 35.3 35.3 116.7 # 557.66 6484450000 123.3 145.0 145.0 123.1 123.1 38.0 38.0 123.3 # 557.68 6485420000 130.0 150 150 130.8 130.8 40.7 40.7 130.0 # 557.70 6486390000 136.7 160.0 160.0 138.5 138.5 43.3 43.3 136.7 # 557.72 6487360000 143.3 170.0 170.0 146.2 146.2 46.0 46.0 143.3 # 557.74 6487845000 146.7 175.0 175.0 150 150 47.3 47.3 146.7 # 557.75 6488330000 150.0 180.0 180.0 160.0 160.0 48.7 48.7 150 # 557.76 6488815000 162.5 185.0 185.0 170.0 170.0 50 50 162.5 # 557.77 6489300000 175.0 190.0 190.0 180.0 180.0 55.0 55.0 175.0 # 557.78 6490270000 200.0 200 200 200 200 65.0 65.0 200 # 557.80 6491240000 220.0 220.0 220.0 220.0 220.0 75 75 220.0 # 557.82 6492210000 240.0 240.0 240.0 240.0 240.0 91.7 91.7 240.0 # 557.84 6492695000 250.0 250 250 250 250 100 100 250 # 557.85 6493180000 260.0 260.0 260.0 260.0 260.0 116.7 116.7 260.0 # 557.86 6494150000 280.0 280 280 280 280 150 150 280 # 557.88 6495120000 282.5 282.5 282.5 282.5 282.5 183.3 183.3 282.5 # 557.90 6495605000 283.8 283.8 283.8 283.8 283.8 200 200 283.8 # 557.91 6496090000 285.0 285.0 285.0 285.0 285.0 216.7 216.7 285.0 # 557.92 6497060000 287.5 287.5 287.5 287.5 287.5 250 250 287.5 # 557.94 6497545000 288.8 288.8 288.8 288.8 288.8 280 280 288.8 # 557.95 6498030000 290.0 290.0 290.0 290.0 290.0 283.0 283.0 290.0 # 557.96 6499000000 292.5 292.5 292.5 292.5 292.5 289.0 289.0 292.5 # 557.98 6499970000 295 295 295 295 295 295 295 295 # 558.00 6505020000 314 314 314 314 314 314 314 314 # 558.10 6510070000 336 336 336 336 336 336 336 336 # 558.20 6515120000 357 357 357 357 357 357 357 357 # 558.30 6520170000 380 380 380 380 380 380 380 380 # 558.40 6525220000 402 402 402 402 402 402 402 402 # 558.50 6530270000 423 423 423 423 423 423 423 423 # 558.60 6535320000 442 442 442 442 442 442 442 442 # 558.70 6540370000 468 468 468 468 468 468 468 468 # 558.80 6550470000 509 509 509 509 509 509 509 509 # 559.00 6555650000 535 535 535 535 535 535 535 535 # 559.10 6559270000 564 564 564 564 564 564 564 564 # 559.17 6602220000 1000 1000 1000 1000 1000 1000 1000 1000 # 560.00 # Brienzersee, etwas stärkere Abgaben 5153590000 13.0 13.0 # 562.4 5168380000 13.0 13.0 # 562.91 5168670000 13.7 13.7 # 562.92 5169250000 15.0 15.0 # 562.94 5169830000 16.3 16.3 # 562.96 5170410000 17.7 17.7 # 562.98 5170990000 19.0 19.0 # 563 5171570000 19.1 19.1 # 563.02 5172150000 19.2 19.2 # 563.04 5172730000 19.2 19.2 # 563.06 5173310000 19.3 19.3 # 563.08 5173890000 19.4 19.4 # 563.1 5174470000 19.5 19.5 # 563.12 5175050000 20.7 20.7 # 563.14 5175630000 21.8 21.8 # 563.16 5176210000 23.0 23.0 # 563.18 5176790000 24.2 24.2 # 563.2 5177370000 25.3 25.3 # 563.22 5177950000 26.5 26.5 # 563.24 5178530000 27.1 27.1 # 563.26 5179110000 27.2 27.2 # 563.28 5180270000 27.3 27.3 # 563.32 5181430000 27.4 27.4 # 563.36 5182590000 27.5 27.5 # 563.4 5183750000 27.6 27.6 # 563.44 5184910000 27.7 27.7 # 563.48 5186650000 27.8 27.8 # 563.54 5187230000 27.8 27.8 # 563.56 5188970000 28.0 28.0 # 563.62 5189550000 28.4 28.4 # 563.64 5190130000 29.1 29.1 # 563.66 5190710000 29.9 29.9 # 563.68 5191290000 30.6 30.6 # 563.7 [abstraction_rule_reservoir_2] # Brienzersee 39 1 366 # Seepegel müM 0 0 0 # 384.7 (564-mittlere Tiefe von 179.3m) 5083990000 0 0 # 560 5141990000 4 4 # 562 5153590000 13.0 13.0 # 13.0 13.0 # 562.4 5168380000 13.0 13.0 # 13.0 13.0 # 562.91 5175050000 20.7 20.7 # 20.7 20.7 # 563.14 5190130000 25.1 25.1 # 25.1 25.1 # 563.66 5190710000 29.9 29.9 # 29.9 29.9 # 563.68 5191290000 30.6 30.6 # 30.6 30.6 # 563.7 5194190000 40.4 40.4 # 40.4 40.4 # 563.8 5195350000 46.7 46.7 # 46.7 46.7 # 563.84 5196510000 53.0 53.0 # 53.0 53.0 # 563.88 5197090000 55.3 55.3 # 56.3 56.3 # 563.9 5197670000 57.6 57.6 # 59.7 59.7 # 563.92 5198250000 59.9 59.9 # 63.0 63.0 # 563.94 5198830000 62.1 62.1 # 66.3 66.3 # 563.96 5199410000 64.4 64.4 # 69.7 69.7 # 563.98 5199990000 66.7 66.7 # 73.0 73.0 # 564 5203070000 78.8 78.8 # 90.7 90.7 # 564.1 5206140000 90.9 90.9 # 109.0 109.0 # 564.2 5209210000 103.0 103.0 # 128.0 128.0 # 564.3 5212290000 115.2 115.2 # 146.0 146.0 # 564.4 5215370000 125.0 125.0 # 166.0 166.0 # 564.5 5218440000 136.6 136.6 # 176.0 176.0 # 564.6 5221520000 148.3 148.3 # 187.0 187.0 # 564.7 5224590000 160.0 160.0 # 197.0 197.0 # 564.8 5227660000 171.6 171.6 # 207.0 207.0 # 564.9 5230740000 183.3 183.3 # 220.0 220.0 # 565 5233820000 195.0 195.0 # 231.0 231.0 # 565.1 5236890000 212.3 212.3 # 242.0 242.0 # 565.2 5246120000 264.2 264.2 # 280 280 # 565.5 5255340000 316.1 316.1 # 320 320 # 565.8 5264570000 368.1 368.1 # 370 370 # 566.1 5273790000 420.0 420.0 # 420 420 # 566.4 5283020000 470.0 470.0 # 470 470 # 566.7 5292240000 530.0 530.0 # 530 530 # 567 5301460000 600.0 600.0 # 600 600 # 567.3 5310690000 700.0 700.0 # 700 700 # 567.6 5322990000 1000.0 1000.0 # 1000 1000 # 568 [abstraction_rule_reservoir_3] # Engstlensee, AE ca. 8.25km^2, mittlerer Abfluss ca. 0.2-0.4 m3/s, 10.7mio m3, TG 43 5 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # = leer 4.5e06 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # 9.0e06 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 # 10.7e06 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 # 11.6e06 30 30 30 30 30 30 30 30 30 30 30 30 # maximale HW-Entlastung [abstraction_rule_reservoir_4] # Triftsee, AE ca. 33km^2, mittlerer Abfluss ca. 1-1.5 m^3/s, 8.6e6 m^3, TG 37 5 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # = leer 4.0e06 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # 8.0e06 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 # 8.6e06 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 # 9.5e06 30 30 30 30 30 30 30 30 30 30 30 30 # maximale HW-Entlastung [abstraction_rule_reservoir_5] # Gruebensee, ca. 35km^2, mittlerer Abfluss ca. 1-1.5m^3/s, 5.7e6 m^3, TG 71+45 5 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # = leer 3.0e06 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # 5.3e06 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 # 5.7e06 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 # 6.2e06 30 30 30 30 30 30 30 30 30 30 30 30 # maximale HW-Entlastung [abstraction_rule_reservoir_6] # KWO, 5 Stauseen (Gelmersee, Tàterichsbodensee, Grimselsee, Totensee, Oberaarsee), werden als ein See betrachtet, AE ~ 153 km^2, 5-10 m^3 mittl. Abfluss 4 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # = leer 192e06 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # 193e06 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 1.5 # 200e06 300 300 300 300 300 300 300 300 300 300 300 300 # maximale HW-Entlastung [abstraction_rule_reservoir_11] # Oeschinensee, AE ca. 12.5km^2, mittlerer Abfluss ca. 0.3-0.6m^3/s, 33mio m^3 Inhalt, TG 59 5 32 60 91 121 152 182 213 244 274 305 335 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # = leer 17e06 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # = leer 30.0e06 0.2 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 # darueber hinaus: outflow = inflow 33.0e06 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 # 36.3e06 30 30 30 30 30 30 30 30 30 30 30 30 # maximale HW-Entlastung [abstraction_rule_abstraction_1] # power generation water from KWO Grimsel et al., SP 1, TG 70, Monday to Friday 4 10 70 91 121 152 182 213 244 274 320 356 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) 0 0 0 0 0 0 0 0 0 0 0 0 0 # empty -> no abstraction possible 20.00e06 20 20 20 20 20 20 20 20 20 20 20 20 # 192.00e06 20 20 20 20 20 20 20 20 20 20 20 20 # 193.00e06 100 100 100 100 100 100 100 100 100 100 100 100 # TargetCap = 120 120 120 120 120 120 120 120 120 120 120 120 # maximum capacity of the target river: abstraction is reduced automatically if the resulting total runoff would be higher (works best with single infliws into a river) WeekDays = 12345 12345 12345 12345 12345 12345 12345 12345 12345 12345 12345 12345 # defines the days of week (ranging from 1 for Monday to 7 for Sunday) the rule is active. On other days, no abstraction is abstracted start_hour = 6 6 6 6 6 6 6 6 6 6 6 6 # at which time should the abstraction start (0=at midnight, 12=at noon, 21=in the evening at 21:00 etc.) stop_hour = 18 18 18 18 18 18 18 18 18 18 18 18 # at which time should the abstraction stop (time format as for start time) [abstraction_rule_abstraction_2] # power generation water from KWO Grimsel et al., SP 1, TG 70, weekend 4 10 70 91 121 152 182 213 244 274 320 356 366 # Julian Days; here: end of the months (rules are valid for the period BEFORE the given JD) 0 0 0 0 0 0 0 0 0 0 0 0 0 # empty -> no abstraction possible 20.00e06 5 5 5 5 5 5 5 5 5 5 5 5 # 192.00e06 5 5 5 5 5 5 5 5 5 5 5 5 # 193.00e06 100 100 100 100 100 100 100 100 100 100 100 100 # TargetCap = 120 120 120 120 120 120 120 120 120 120 120 120 # maximum capacity of the target river: abstraction is reduced automatically if the resulting total runoff would be higher (works best with single infliws into a river) WeekDays = 67 67 67 67 67 67 67 67 67 67 67 67 # defines the days of week (ranging from 1 for Monday to 7 for Sunday) the rule is active. On other days, no abstraction is abstracted start_hour = 7 7 7 7 7 7 7 7 7 7 7 7 # at which time should the abstraction start (0=at midnight, 12=at noon, 21=in the evening at 21:00 etc.) stop_hour = 19 19 19 19 19 19 19 19 19 19 19 19 # at which time should the abstraction stop (time format as for start time) # 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 [multilayer_landuse] 17 # count of multilayer landuses 1 Wasserflächen { Landuse_Layers = 16, -9999; k_extinct = 0.3; LAI_scale = 20;} 2 Bebauung { Landuse_Layers = 2, -9999; k_extinct = 0.3; LAI_scale = 20;} 3 Teilbebauung { Landuse_Layers = 1, -9999; k_extinct = 0.3; LAI_scale = 20;} 4 Wald { Landuse_Layers = 13, 4; k_extinct = 0.3; LAI_scale = 20;} 5 offener_Wald { Landuse_Layers = 14, 6; k_extinct = 0.3; LAI_scale = 20;} 6 Gebueschwald { Landuse_Layers = 10, 4; k_extinct = 0.3; LAI_scale = 20;} 7 Ackerland { Landuse_Layers = 7, -9999; k_extinct = 0.3; LAI_scale = 20;} 8 Gruenland { Landuse_Layers = 5, -9999; k_extinct = 0.3; LAI_scale = 20;} 9 Obstbauflaechen { Landuse_Layers = 17, 4; k_extinct = 0.3; LAI_scale = 20;} 10 Gebuesch { Landuse_Layers = 10, -9999; k_extinct = 0.3; LAI_scale = 20;} 11 verbuschtes_Weidel { Landuse_Layers = 9, -9999; k_extinct = 0.3; LAI_scale = 20;} 12 uebrige_Gruenfl { Landuse_Layers = 6, -9999; k_extinct = 0.3; LAI_scale = 20;} 13 Moore_Suempfe { Landuse_Layers = 15, -9999; k_extinct = 0.3; LAI_scale = 20;} 14 versteinte_wiesen { Landuse_Layers = 4, -9999; k_extinct = 0.3; LAI_scale = 20;} 15 Felsflaechen_Oedland { Landuse_Layers = 3, -9999; k_extinct = 0.3; LAI_scale = 20;} 16 Gletscher { Landuse_Layers = 22, -9999; k_extinct = 0.3; LAI_scale = 20;} 21 Deponien { Landuse_Layers = 3, -9999; k_extinct = 0.3; LAI_scale = 20;} [landuse_table] 19 1 teilversiegelte_Flaechen {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.3; 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 80 70 70 70 70 70 70 100 100 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 150 150 150 150 150 150 150 150 150 150 150 150 ; 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.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; 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; } 2 versiegelte_Flaechen {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 = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.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.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; } 3 vegetationslose_Flaechen {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 = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 100 100 100 100 100 100 100 100 100 100 100 100 ; LAI = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; Z0 = 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.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.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 ; 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 spaerliche_Vegetation {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 = 90 90 80 60 50 45 45 45 50 60 90 90 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 100 100 100 100 100 100 100 100 100 100 100 100 ; LAI = 1 1 1 1.5 1.5 2 2 2 1.5 1.5 1 1 ; 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.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.3 0.2 0.2 0.2 ; 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; } 5 Intensiv-Gruenland {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.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 60 50 45 40 45 50 60 90 90 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 2 2 2 3 4 4 4 4 3 2 2 2 ; 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.95 0.95 0.95 0.95 1.0 1.0 1.0 1.0 0.95 0.95 0.95 0.95 ; RootDepth = 0.3 0.3 0.4 0.5 0.5 0.6 0.6 0.5 0.5 0.4 0.3 0.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 Extensiv-Gruenland {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 60 50 45 40 45 50 60 90 90 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 2 2 2 2 3 3 3 3 3 2 2 2 ; 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.95 1.0 1.0 1.0 1.0 0.95 0.9 0.9 0.9 ; 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; } 7 Intensiv-Ackerland_unbewaessert {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 65 55 45 45 45 55 65 90 90 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.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.1 0.1 0.1 0.15 0.2 0.2 0.2 0.2 0.2 0.1 0.1 0.1 ; 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; } 8 Extensiv-Ackerland {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 65 55 45 45 45 55 65 90 90 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 1 1 2 3 4 4 4 4 3 2 1 1 ; Z0 = 0.1 0.1 0.1 0.15 0.2 0.2 0.2 0.2 0.2 0.1 0.1 0.1 ; VCF = 0.3 0.3 0.3 0.5 0.8 0.8 0.8 0.7 0.6 0.3 0.3 0.3 ; 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; } 9 Heidevegetation {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 60 50 40 40 50 50 60 70 80 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 3 3 3 4 4 4 4 4 3 3 3 3 ; Z0 = 0.1 0.1 0.1 0.15 0.2 0.2 0.2 0.2 0.2 0.1 0.1 0.1 ; VCF = 1 1 1 1 1 1 1 1 1 1 1 1 ; 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; } 10 Busch-Kraut-Vegetation_mitteldicht {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 60 50 40 40 50 50 60 70 80 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; 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; } 11 Laubwald {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 65 45 45 45 45 50 75 100 100 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 1 1 4 4 6 8 8 7 6 4 1 1 ; Z0 = 2 2 2 3 5 5 5 5 4 3 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; } 12 Nadelwald {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 55 40 40 40 40 45 65 80 80 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; 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; } 13 Mischwald {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 60 45 45 45 45 50 70 90 90 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 3 3 3 6 8 8 8 8 8 6 3 3 ; Z0 = 2.5 2.5 2.5 3.0 3.5 4.5 4.5 4.5 4.0 3.0 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; } 14 locker_baumbestanden {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.5; 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 60 45 45 45 45 50 70 90 90 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 200 200 200 200 200 200 200 200 200 200 200 200 ; LAI = 2.5 2.5 2.5 3.5 4.0 4.0 4.0 4.0 3.5 2.5 2.5 2.5 ; 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.4 0.4 0.4 0.4 0.8 0.9 0.9 0.9 0.7 0.4 0.4 0.4 ; 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; } 15 Moore_Suempfe {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.5; IntercepCap = 0.0; 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 = 50 50 40 40 30 20 20 30 40 40 50 50 ; rs_interception = 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ; rs_evaporation = 100 100 100 100 100 100 100 100 100 100 100 100 ; LAI = 3 3 3 3 4 4 4 4 3 3 3 3 ; Z0 = 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 ; 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.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; } 16 Wasserflaechen {method = VariableDayCount; RootDistr = 1; TReduWet = 1; LimitReduWet = 1; HReduDry = 150; IntercepCap = 0; JulDays = 365; Albedo = 0.15; rsc = 1; rs_interception = 0; rs_evaporation = 0; LAI = 1; Z0 = 0.1; VCF = 1; RootDepth = 1; AltDep = 0; } 17 horticulture {method = VariableDayCount; RootDistr = 1.0; TReduWet = 0.95; LimitReduWet = 0.5; HReduDry = 3.45; IntercepCap = 0.35; 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 = 100 100 90 60 45 45 45 45 45 70 100 100; rs_interception = 100 100 90 60 50 50 50 50 50 70 100 100; rs_evaporation = 280 280 280 280 280 280 280 280 280 280 280 280; LAI = 0.5 0.5 0.5 2 4 5 5 5 4 3 0.5 0.5; Z0 = 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0; VCF = 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75; RootDepth = 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8; 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; } 22 Eisflaechen {method = VariableDayCount; RootDistr = 1; TReduWet = 1; LimitReduWet = 1; HReduDry = 150; IntercepCap = 0; JulDays = 365; Albedo = 0.38; rsc = 100; rs_interception = 0; rs_evaporation = 400; LAI = 1.0; Z0 = 0.05; VCF = 1.0; RootDepth = 0.1; AltDep = 0; } 23 Firnflaechen {method = VariableDayCount; RootDistr = 1; TReduWet = 1; LimitReduWet = 1; HReduDry = 150; IntercepCap = 0; JulDays = 365; Albedo = 0.60; rsc = 100; rs_interception = 0; rs_evaporation = 400; LAI = 1.0; Z0 = 0.05; VCF = 1.0; RootDepth = 0.1; AltDep = 0; } #---------------------------------- #0.5x #$set $9e4 = 4.50e-4 #$set $8e4 = 4.00e-4 #$set $7e4 = 3.50e-4 #$set $6e4 = 3.00e-4 #$set $5e4 = 2.50e-4 #$set $4e4 = 2.00e-4 #$set $3e4 = 1.50e-4 #$set $2e4 = 1.00e-4 #$set $1e4 = 0.50e-4 # #$set $9e5 = 4.50e-5 #$set $8e5 = 4.00e-5 #$set $7e5 = 3.50e-5 #$set $6e5 = 3.00e-5 #$set $5e5 = 2.50e-5 #$set $4e5 = 2.00e-5 #$set $3e5 = 1.50e-5 #$set $2e5 = 1.00e-5 #$set $1e5 = 0.50e-5 # #$set $9e6 = 5.00e-6 #$set $8e6 = 4.00e-6 #$set $7e6 = 3.50e-6 #$set $6e6 = 3.00e-6 #$set $5e6 = 2.50e-6 #$set $4e6 = 2.00e-6 #$set $3e6 = 1.50e-6 #$set $2e6 = 1.00e-6 #$set $1e6 = 0.50e-6 # #$set $9e7 = 4.50e-7 #$set $8e7 = 4.00e-7 #$set $7e7 = 3.50e-7 #$set $6e7 = 3.00e-7 #$set $5e7 = 2.50e-7 #$set $4e7 = 2.00e-7 #$set $3e7 = 1.50e-7 #$set $2e7 = 1.00e-7 #$set $1e7 = 0.50e-7 # #$set $9e8 = 4.50e-8 #$set $8e8 = 4.00e-8 #$set $7e8 = 3.50e-8 #$set $6e8 = 3.00e-8 #$set $5e8 = 2.50e-8 #$set $4e8 = 2.00e-8 #$set $3e8 = 1.50e-8 #$set $2e8 = 1.00e-8 #$set $1e8 = 0.50e-8 # #$set $9e9 = 4.50e-9 #$set $8e9 = 4.00e-9 #$set $7e9 = 3.50e-9 #$set $6e9 = 3.00e-9 #$set $5e9 = 2.50e-9 #$set $4e9 = 2.00e-9 #$set $3e9 = 1.50e-9 #$set $2e9 = 1.00e-9 #$set $1e9 = 0.50e-9 #---------------------------------- # 1x $set $9e4 = 9.00e-4 $set $8e4 = 8.00e-4 $set $7e4 = 7.00e-4 $set $6e4 = 6.00e-4 $set $5e4 = 5.00e-4 $set $4e4 = 4.00e-4 $set $3e4 = 3.00e-4 $set $2e4 = 2.00e-4 $set $1e4 = 1.00e-4 $set $9e5 = 9.00e-5 $set $8e5 = 8.00e-5 $set $7e5 = 7.00e-5 $set $6e5 = 6.00e-5 $set $5e5 = 5.00e-5 $set $4e5 = 4.00e-5 $set $3e5 = 3.00e-5 $set $2e5 = 2.00e-5 $set $1e5 = 1.00e-5 $set $9e6 = 9.00e-6 $set $8e6 = 8.00e-6 $set $7e6 = 7.00e-6 $set $6e6 = 6.00e-6 $set $5e6 = 5.00e-6 $set $4e6 = 4.00e-6 $set $3e6 = 3.00e-6 $set $2e6 = 2.00e-6 $set $1e6 = 1.00e-6 $set $9e7 = 9.00e-7 $set $8e7 = 8.00e-7 $set $7e7 = 7.00e-7 $set $6e7 = 6.00e-7 $set $5e7 = 5.00e-7 $set $4e7 = 4.00e-7 $set $3e7 = 3.00e-7 $set $2e7 = 2.00e-7 $set $1e7 = 1.00e-7 $set $9e8 = 9.00e-8 $set $8e8 = 8.00e-8 $set $7e8 = 7.00e-8 $set $6e8 = 6.00e-8 $set $5e8 = 5.00e-8 $set $4e8 = 4.00e-8 $set $3e8 = 3.00e-8 $set $2e8 = 2.00e-8 $set $1e8 = 1.00e-8 $set $9e9 = 9.00e-9 $set $8e9 = 8.00e-9 $set $7e9 = 7.00e-9 $set $6e9 = 6.00e-9 $set $5e9 = 5.00e-9 $set $4e9 = 4.00e-9 $set $3e9 = 3.00e-9 $set $2e9 = 2.00e-9 $set $1e9 = 1.00e-9 #---------------------------------- # 2x #$set $9e4 = 18.00e-4 #$set $8e4 = 16.00e-4 #$set $7e4 = 14.00e-4 #$set $6e4 = 12.00e-4 #$set $5e4 = 10.00e-4 #$set $4e4 = 8.00e-4 #$set $3e4 = 6.00e-4 #$set $2e4 = 4.00e-4 #$set $1e4 = 2.00e-4 # #$set $9e5 = 18.00e-5 #$set $8e5 = 16.00e-5 #$set $7e5 = 14.00e-5 #$set $6e5 = 12.00e-5 #$set $5e5 = 10.00e-5 #$set $4e5 = 8.00e-5 #$set $3e5 = 6.00e-5 #$set $2e5 = 4.00e-5 #$set $1e5 = 2.00e-5 # #$set $9e6 = 18.00e-6 #$set $8e6 = 16.00e-6 #$set $7e6 = 14.00e-6 #$set $6e6 = 12.00e-6 #$set $5e6 = 10.00e-6 #$set $4e6 = 8.00e-6 #$set $3e6 = 6.00e-6 #$set $2e6 = 4.00e-6 #$set $1e6 = 2.00e-6 # #$set $9e7 = 18.00e-7 #$set $8e7 = 16.00e-7 #$set $7e7 = 14.00e-7 #$set $6e7 = 12.00e-7 #$set $5e7 = 10.00e-7 #$set $4e7 = 8.00e-7 #$set $3e7 = 6.00e-7 #$set $2e7 = 4.00e-7 #$set $1e7 = 2.00e-7 # #$set $9e8 = 18.00e-8 #$set $8e8 = 16.00e-8 #$set $7e8 = 14.00e-8 #$set $6e8 = 12.00e-8 #$set $5e8 = 10.00e-8 #$set $4e8 = 8.00e-8 #$set $3e8 = 6.00e-8 #$set $2e8 = 4.00e-8 #$set $1e8 = 2.00e-8 # #$set $9e9 = 18.00e-9 #$set $8e9 = 16.00e-9 #$set $7e9 = 14.00e-9 #$set $6e9 = 12.00e-9 #$set $5e9 = 10.00e-9 #$set $4e9 = 8.00e-9 #$set $3e9 = 6.00e-9 #$set $2e9 = 4.00e-9 #$set $1e9 = 2.00e-9 #---------------------------------- #5x #$set $9e4 = 4.50e-3 #$set $8e4 = 4.00e-3 #$set $7e4 = 3.50e-3 #$set $6e4 = 3.00e-3 #$set $5e4 = 2.50e-3 #$set $4e4 = 2.00e-3 #$set $3e4 = 1.50e-3 #$set $2e4 = 1.00e-3 #$set $1e4 = 0.50e-3 # #$set $9e5 = 4.50e-4 #$set $8e5 = 4.00e-4 #$set $7e5 = 3.50e-4 #$set $6e5 = 3.00e-4 #$set $5e5 = 2.50e-4 #$set $4e5 = 2.00e-4 #$set $3e5 = 1.50e-4 #$set $2e5 = 1.00e-4 #$set $1e5 = 0.50e-4 # #$set $9e6 = 5.00e-5 #$set $8e6 = 4.00e-5 #$set $7e6 = 3.50e-5 #$set $6e6 = 3.00e-5 #$set $5e6 = 2.50e-5 #$set $4e6 = 2.00e-5 #$set $3e6 = 1.50e-5 #$set $2e6 = 1.00e-5 #$set $1e6 = 0.50e-5 # #$set $9e7 = 4.50e-6 #$set $8e7 = 4.00e-6 #$set $7e7 = 3.50e-6 #$set $6e7 = 3.00e-6 #$set $5e7 = 2.50e-6 #$set $4e7 = 2.00e-6 #$set $3e7 = 1.50e-6 #$set $2e7 = 1.00e-6 #$set $1e7 = 0.50e-6 # #$set $9e8 = 4.50e-7 #$set $8e8 = 4.00e-7 #$set $7e8 = 3.50e-7 #$set $6e8 = 3.00e-7 #$set $5e8 = 2.50e-7 #$set $4e8 = 2.00e-7 #$set $3e8 = 1.50e-7 #$set $2e8 = 1.00e-7 #$set $1e8 = 0.50e-7 # #$set $9e9 = 4.50e-8 #$set $8e9 = 4.00e-8 #$set $7e9 = 3.50e-8 #$set $6e9 = 3.00e-8 #$set $5e9 = 2.50e-8 #$set $4e9 = 2.00e-8 #$set $3e9 = 1.50e-8 #$set $2e9 = 1.00e-8 #$set $1e9 = 0.50e-8 #---------------------------------- # 10x #$set $9e4 = 9.00e-3 #$set $8e4 = 8.00e-3 #$set $7e4 = 7.00e-3 #$set $6e4 = 6.00e-3 #$set $5e4 = 5.00e-3 #$set $4e4 = 4.00e-3 #$set $3e4 = 3.00e-3 #$set $2e4 = 2.00e-3 #$set $1e4 = 1.00e-3 # #$set $9e5 = 9.00e-4 #$set $8e5 = 8.00e-4 #$set $7e5 = 7.00e-4 #$set $6e5 = 6.00e-4 #$set $5e5 = 5.00e-4 #$set $4e5 = 4.00e-4 #$set $3e5 = 3.00e-4 #$set $2e5 = 2.00e-4 #$set $1e5 = 1.00e-4 # #$set $9e6 = 9.00e-5 #$set $8e6 = 8.00e-5 #$set $7e6 = 7.00e-5 #$set $6e6 = 6.00e-5 #$set $5e6 = 5.00e-5 #$set $4e6 = 4.00e-5 #$set $3e6 = 3.00e-5 #$set $2e6 = 2.00e-5 #$set $1e6 = 1.00e-5 # #$set $9e7 = 9.00e-6 #$set $8e7 = 8.00e-6 #$set $7e7 = 7.00e-6 #$set $6e7 = 6.00e-6 #$set $5e7 = 5.00e-6 #$set $4e7 = 4.00e-6 #$set $3e7 = 3.00e-6 #$set $2e7 = 2.00e-6 #$set $1e7 = 1.00e-6 # #$set $9e8 = 9.00e-7 #$set $8e8 = 8.00e-7 #$set $7e8 = 7.00e-7 #$set $6e8 = 6.00e-7 #$set $5e8 = 5.00e-7 #$set $4e8 = 4.00e-7 #$set $3e8 = 3.00e-7 #$set $2e8 = 2.00e-7 #$set $1e8 = 1.00e-7 # #$set $9e9 = 9.00e-8 #$set $8e9 = 8.00e-8 #$set $7e9 = 7.00e-8 #$set $6e9 = 6.00e-8 #$set $5e9 = 5.00e-8 #$set $4e9 = 4.00e-8 #$set $3e9 = 3.00e-8 #$set $2e9 = 2.00e-8 #$set $1e9 = 1.00e-8 #---------------------------------- [soil_table] 143 #co- name of the # van Genuchten Parameter nach Carsel & Parrish (1988) #de soil profile #-- --------------- 23 1_Seen_subhydrische_Boeden {method = MultipleHorizons; # AE = 0.28 % PMacroThresh = 0 ; MacroCapacity = 0 ; CapacityRedu = 0 ; MacroDepth = 0 ; horizon = 1 2 3 ; Name = SIL SIL CL ; ksat = $3e7 $1e7 $1e7 ; # SIL: $1e6 k_recession = 0.9 0.9 0.9 ; theta_sat = 0.45 0.45 0.40 ; theta_res = 0.067 0.067 0.095 ; alpha = 2.00 2.0 1.90 ; Par_n = 1.41 1.41 1.31 ; Par_tau = 0.5 0.5 0.5 ; thickness = 0.10 0.30 2.00 ; # 19.0 m layers = 4 2 9 ; } 26 6_Moorboeden {method = MultipleHorizons; #vgl. F1 # AE = 0.01 % PMacroThresh = 1.0 ; MacroCapacity = 20.0 ; CapacityRedu = 1.0 ; MacroDepth = 2.0 ; horizon = 1 2 3 ; Name = HMT NMT hHr_SL3 ; ksat = $1e5 $3e6 $1e6 ; k_recession = 0.9 0.9 0.9 ; theta_sat = 0.65 0.55 0.40 ; theta_res = 0.30 0.30 0.065 ; alpha = 4.0 4.0 7.50 ; Par_n = 1.2 1.2 1.89 ; Par_tau = 0.5 0.5 0.5 ; thickness = 0.10 0.30 2.00 ; # 19.0 m layers = 4 2 9 ; } 13 7_Siedlungsflaechen {method = MultipleHorizons; #vgl. 8 # AE = 0.42 % PMacroThresh = 4.0 ; MacroCapacity = 4.0 ; CapacityRedu = 1.0 ; MacroDepth = 2.0 ; horizon = 1 2 3 ; Name = SIL SIL R ; ksat = $1e6 $3e7 $5e8 ; # SIL: $1e6 k_recession = 0.9 0.9 0.9 ; theta_sat = 0.45 0.45 0.20 ; theta_res = 0.067 0.067 0.04 ; alpha = 2.0 2.0 8.0 ; Par_n = 1.41 1.41 1.80 ; Par_tau = 0.5 0.5 0.5 ; thickness = 0.10 0.35 0.35 ; # 4.0 m layers = 5 2 8 ; } 91 8_Felsflaechen {method = MultipleHorizons; # AE = 37 % PMacroThresh = 4.0 ; MacroCapacity = 4.0 ; CapacityRedu = 1.0 ; MacroDepth = 2.0 ; horizon = 1 2 3 4 ; Name = SIL CL R R ; # ksat = $1e6 $3e7 $1e8 $1e8 ; # SIL: $1e6 ksat = $1e5 $3e6 $1e6 $1e6 ; # SIL: $1e6 k_recession = 0.9 0.9 0.9 0.9 ; theta_sat = 0.45 0.41 0.20 0.20 ; theta_res = 0.067 0.095 0.04 0.04 ; alpha = 2.00 1.90 8.00 8.00 ; Par_n = 1.41 1.31 1.80 1.80 ; Par_tau = 0.5 0.5 0.5 0.5 ; thickness = 0.20 0.20 0.30 0.45 ; # = 4.0 m layers = 4 4 5 2 ; # thickness = 0.10 0.10 0.15 0.2 ; # = 4x10cm + 4x10cm + 4x14cm + 3x20cm = 40+40+60+60 = 200cm # layers = 4 4 4 3 ; } #...etc (put your own land uses in here... [substance_transport] 0 [irrigation_table] 2 # 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; } }