$set $startyear           = 1985
$set $endyear             = 1985

#
# comment lines are indicated by an "#". They are allowed after the entries or at the beginning of the lines
#
# paragraphs of the control file:
#
# [output_list]
# [output_interval]
# [coordinates]            (geogr. coordinates)
# [region_transition_distance] for multiple regression regions
# [elevation_model]        (name of the elevation model)
# [zonengrid]              (name of the zone grid)
# [standardgrids]          (name of other static grids like slope angle, slope aspect, topogr. factor)
# [variable_grids]         (names of albedo and soil storage - used by more than one modules)
# [model_time]             (start end end-dates of model time)
# [meteo_data_count]       (number of mete data to interpolate)
# [meteo_names]            (names of meteo data to interpolate - each name is later the headline of a paragraph for interpolation
# [precipitation_correction] (paragraph for parameters of the prec.-correction)
# [radiation_correction]   (paragraph with parameters for radiation correction)
# [evapotranspiration]     (paragraph with parameters for evapotranspiration)
# [snow_model]             (paragraph with parameters for the snow model)
# [ice_firn]               (paragraph with parameters for the glacier model)
# [interzeption_model]     (paragraph with parameters for the interzeption modell)
# [infiltration_model]     (paragraph with parameters for the infiltration modell)
# [soil_model]             (paragraph with parameters for the soilmodel)
# [SiltingUpModel]         (paragraph with parameters for the silting up model)
# [unsatzon_model]         (paragraph with parameters for the unsaturated zone model)
# [SurfaceRoutingModel]    (paragraph with parameters for the surface routing model)
# [irrigation]             (paragraph with parameters for irrigation model)
# [groundwater_flow]       (paragraph with parameters for groundwater model)
# [ExternalCoupling]       (paragraph with parameters for external coupling)
# [routing_model]          (paragraph with Parametern for discharge routing)
# [landuse_table]          (paragraph with land use (vegetation) parameters)
# [soil_table]             (paragraph with soil properties)
# [substance_transport]    (paragraph with tracer properties)
# [abstraction_rule_reservoir_n] rules for routing submodel, each reservoir has its own rule
# [multilayer_landuse]      new multilayer landuse table definitio
#
# symbol definitions begin with the set command:
# max. 200 symbols can be defined

$set $mainpath	     = d:\Wasim\TestData\

# use an empty value here will cause grids to be read in from $InitialStateDirectory and to be written to $DefaultOutputDirectory 
# which is much more convenient that having the input grids overwritten each time wasim runs at the end.
$set $outpath        =  

$set $inpath         =  $mainpath//Input\
$set $exchngpath     =  $mainpath//Exchng\
$set $InitialStateDirectory = $mainpath//StateIni\
$set $DefaultOutputDirectory = $mainpath//Output\


$set $time           =  1440.0

$set $grid           =  col500
$set $stack          =  cols500
$set $suffix         =  grd
$set $code           =  s

# variables for standardgrids
# first section: grids, which differ for different subdivisions of the basin
$set $zone_grid             =  $grid//.ezg57d
$set $subcatchments         =  $grid//.ezg57d
$set $flow_time_grid        =  $grid//.fzs30
$set $river_links_grid      =  $grid//.lnk57
$set $regio_grid            =  $grid//.reg

#second section: grids, which doesn't depend on subdivision (only pixel-values are of interest)
$set $elevation_model       =  $grid//.dhk
$set $RelCellArea_grid      =  $grid//.rca
$set $CellSizeX_grid        =  $grid//.csx
$set $CellSizeY_grid        =  $grid//.csy
$set $slope_grid            =  $grid//.slp
$set $FlowDirection_grid    =  $grid//.fld
$set $aspect_grid           =  $grid//.exp
$set $land_use_grid         =  $grid//.use
$set $ice_firn_grid         =  $grid//.ice
$set $field_capacity_grid   =  $grid//.nfk
$set $ATBgrid               =  $grid//.atb
$set $hydr_cond_grid        =  $grid//.k
$set $soil_types            =  $grid//.art
$set $sky_view_factor_grid  =  $grid//.hor
$set $drain_depth_grid      =  $grid//.drn
$set $drain_distance_grid   =  $grid//.dis
$set $irrigationcodes       =  $grid//.irr19
$set $max_pond_grid         =  $grid//.maxpond
#$set $mean_pond_grid        =  $grid//.meanpond
$set $clay_depth_grid       =  $grid//.cly
$set $river_depth_grid      =  $grid//.dep
$set $river_width_grid      =  $grid//.wit
$set $tracer_1              =  $grid//.c1
$set $tracer_2              =  $grid//.c2
$set $tracer_3              =  $grid//.c3
$set $tracer_4              =  $grid//.c4
$set $tracer_5              =  $grid//.c5
$set $tracer_6              =  $grid//.c6
$set $tracer_7              =  $grid//.c7
$set $tracer_8              =  $grid//.c8
$set $tracer_9              =  $grid//.c9
$set $kolmationsgrid        =  $grid//.kol
$set $gw_kx_1_grid          =  $grid//.kx005
$set $gw_kx_2_grid          =  $grid//.kx2
$set $gw_kx_3_grid          =  $grid//.kx3
$set $gw_ky_1_grid          =  $grid//.ky005
$set $gw_ky_2_grid          =  $grid//.ky2
$set $gw_ky_3_grid          =  $grid//.ky3
$set $gw_bound_h_1_grid     =  $grid//.bh1
$set $gw_bound_h_2_grid     =  $grid//.bh2
$set $gw_bound_h_3_grid     =  $grid//.bh3
$set $gw_bound_q_1_grid     =  $grid//.bq1
$set $gw_bound_q_2_grid     =  $grid//.bq2
$set $gw_bound_q_3_grid     =  $grid//.bq3
$set $aquiferthick1         =  $grid//.aq1
$set $aquiferthick2         =  $grid//.aq2
$set $aquiferthick3         =  $grid//.aq3
$set $gw_storage_coeff_1    =  $grid//.s01
$set $gw_storage_coeff_2    =  $grid//.s02
$set $gw_storage_coeff_3    =  $grid//.s03
$set $gw_kolmation_1        =  $grid//.gk1
$set $gw_kolmation_2        =  $grid//.gk2
$set $gw_kolmation_3        =  $grid//.gk3
$set $lake_grid             =  $grid//.lak
$set $taucrit_grid          =  $grid//.tau
$set $ThawCoeffPermaFrost   =  $grid//.alpha
$set $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_12_grid        =  prec//$grid//.//$suffix
$set $preci_idw_grid       =  prec//$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   = glac//$grid//.//$suffix
$set $glacier_codes_grid      = gid_//$grid//.//$suffix

# result-stacks for Unsatzonmodel
$set $Thetastack              = thet//$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
$set $SEstack	      	= SE__//$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

# 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

# result stackd for layered snow model
$set $SWESstack		= swes_stack//$grid//.//$suffix		# stack with SWE solid values per snow layer
$set $SWELStack		= swel_stack//$grid//.//$suffix		# stack with SWE liquid values per snow layer
$set $temperSnowStack		= tempersnow_stack//$grid//.//$suffix	# stack with snow temperatures per snow layer
$set $thetaSnowStack		= thetasnow_stack//$grid//.//$suffix	# stack with theta values per snow layer
$set $hhSnowStack		= hhydsnow_stack//$grid//.//$suffix	# stack with hydraulic heads per snow layer
$set $hgeoSnowStack		= hgeosnow_stack//$grid//.//$suffix	# stack with geodetic elevation per snow layer (bottom of layer)
$set $quSnowStack		= qusnow_stack//$grid//.//$suffix	# stack with (bottom) water fluxes per snow layer
$set $thickSnowStack		= thicksnow_stack//$grid//.//$suffix	# stack with the thickness of each snow layer
$set $densSnowStack		= denssnow_stack//$grid//.//$suffix	# stack with density of each snow layer
$set $upperIndexSnow		= upperIndexSnow//$grid//.//$suffix	# grid with the index of the first (uppermost) layer that constains snow

# parameters for interpolation of meteorological input data
$set $SzenUse        =  0
$set $IDWmaxdist     =  100000
$set $IDWweight      =  2
$set $Anisoslope     =  35.0
$set $Anisotropie    =  0.9

# explanation of writegrid and outputcode some lines below
$set $Writegrid      =  3
$set $Writestack     =  3

[output_list]
0  # numer of following rows indicating each a subbasin for which output should be generated
1  # ascending order is required !
2
10
12

#$set $outputcode     =  5001
#$set $output_meteo   =  5001
#$set $day_sum        =  6024
#$set $day_mean       =  5024
#$set $hour_mean      =  5001
#$set $routing_code   =  6001

$set $outputcode     =  2001
$set $output_meteo   =  2001
$set $day_sum        =  4024
$set $day_mean       =  2024
$set $hour_mean      =  2001
$set $routing_code   =  4001

# readgrids : 1 = read storage grids (as SI, SSNOW,SLIQ...) from hard disk, 0=generate and initialize with 0
$set $readgrids     =  0
# 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 

#
#  Wrigegrid : max. 3 digits (nnn)
#
#     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<nnn>  : spatial mean values for the entire basin, averaged in time over <nnn> intervals (timesteps)
# 2<nnn>  : spatial mean values for all zones (subbasin) and for the entire basin, averaged in time over <nnn> intervals (timesteps)
# 3<nnn>  : spatial means for the entire basin, added up in time over <nnn> intervals (timesteps)
# 4<nnn>  : spatial means for all zones (subbasin) and for the entire basin, added up in time over <nnn> intervals (timesteps)
# 5<nnn>  : spatial means for the entire basin and for those subbasins which are specified in the output-list, averaged in time over <nnn> intervals
# 6<nnn>  : spatial means for the entire basin and for those subbasins which are specified in the output-list, added up in time over <nnn> 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]
0	# doSpinUp: 0=no, 1=yes
1	# doTemperaturePreSpinUp: 0=no, 1=yes
1	# repeatCnt: how many times should the model run for the defined time period (in [model_time]), e.g. 5 for 5 tiumes
3	# 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]
18     # 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

[output_interval]
10             # 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              # output of runoff results in 0=mm/time step, 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 = 0 # 0 if grids should be written in WaSiM native format, 1 if in ESRI ASCII format
WriteDoubleGrids = 0 # 0 if grids should be written in WaSiM native format, 1 if in ESRI ASCII format
InitialStateDirectory = $InitialStateDirectory # if using this parameter, all state grids  as well as the storage_richards.ftz file will be expected in that directory for reading
DefaultOutputDirectory = $DefaultOutputDirectory # this is the default output directory, all output is written to unless the given filename contains an absolute path (starting with either / or ~ for UNIX or a drive letter and :\ for Windows
# there are some exceptions, though: for external coupling no default output path is used
# relative pathnames may be used as well.
# for compatibility reasins with older control files and WaSiM versions, both directories will only be used if the given filename has no absolute path,
# so in order to use the new features, all $outpath uses should be reviewed and removed if necessary (or the variable should be set to an empty string)

[coordinates]
51.0           # geogr. latitude (center of the basin -> for radiation calculations)
-116.6         # geogr. longitude (center of the basin)
-120.0           # meridian according to the official time (middle europe: 15)(east: 0 ... +180 degree, west: 0 ... -180 (or 360 ... 180)
0              # 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//$elevation_model    # grid with the digital elevation data

[zone_grid]
$inpath//$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//$suffix
$set $z0_stat = z0_//$grid//.//$code//$suffix
$set $root_stat = root_//$grid//.//$code//$suffix
$set $rse_stat = rse_//$grid//.//$code//$suffix
$set $rsi_stat = rsi_//$grid//.//$code//$suffix
$set $rsc_stat = rsc_//$grid//.//$code//$suffix
$set $albedo_stat = albedo_//$grid//.//$code//$suffix 
$set $vcf_stat = vcf_//$grid//.//$code//$suffix

$set $T_Lower_Boundary = col500.tlow

# 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]
17                     # number of standard grids
# path                         # identification         # fillcode 0=no, 1=yes (fill missing values with values of nearest neighbor)
#$inpath//$regio_grid			regression_regions        1    # region grid if using multiple regression perameter files for meteorological data interpolation
$inpath//$land_use_grid         landuse                   1    # grid with land use data
#$inpath//$ice_firn_grid         ice_firn                   0   # grid with firn or ice cells (code 0: nodata values should not be replaced by nearest neighbour)
$inpath//$subcatchments         zonegrid_soilmodel        1    # zone grid for the runoff generation model (and unstaurated zone model)
$inpath//$soil_types            soil_types                1    # soil types as codes for the soil table
$inpath//$flow_time_grid        flow_times                1    # grid with flow times for surface runoff to the subbasin outlet
$inpath//$river_depth_grid      river_depth               0    		# grid with the depth of all streams in the stream network in m
$inpath//$river_width_grid      river_width               0    		# grid with the witdh of all streams in m
$inpath//$river_links_grid      river_links               0    		# grid with codes of tributaries, from which a channel was routed (only for real routing channels!!!)
$inpath//$kolmationsgrid        kolmation                 0    		# grid with codes of tributaries, from which a channel was routed (only for real routing channels!!!)
$inpath//$aquiferthick1         aquifer_thickness_1       1    		# grid with thickness of first aquifer (m from soil surface to the aquifer bottom)
$inpath//$gw_storage_coeff_1    gw_storage_coeff_1        1    		# storage coefficients for 1. aquifer
$inpath//$gw_bound_h_1_grid     gw_boundary_fix_h_1       0   # periodicity = 1 D 12 persistent = 0 2 # boundary conditions 1 constant head for layer 1, periodicity = <period> <type of period H=hour, D=day, M=month, Y=year> <subintervals> persistent = 0 means that the grid is valid for the next number of intervals as given as last parameter only (then, the default grid is valid again, 1 = new read in grid is persistent)
$inpath//$gw_bound_q_1_grid     gw_boundary_fix_q_1       0    		# boundary conditions 2 (given flux perpendicular to the border) for layer 1
$inpath//$gw_kx_1_grid          gw_k_x_1                  1    		# lateral hydraulic conductivities for the 1. aquifer in x direction
$inpath//$gw_ky_1_grid          gw_k_y_1                  1    		# lateral hydraulic conductivities for the 1. aquifer in y direction
$inpath//$gw_kolmation_1        gw_kolmation_1            1    		# kolmation (leakage factor) between 1st and 2nd aquifer
$inpath//$ATBgrid		topographic_faktor	  1             # topographic factor, used for topmodel version only
$inpath//$T_Lower_Boundary      _T_Lower_Boundary_Condition_ 1 # initial grid for lower boundary condition of the temperatuere in the soil

# more possible options and grids that may be read in if the algorithms require these grids (like slide fractions for snow model methods woth slides etc.)
$inpath//$river_depth_grid      	river_depth               0    		# grid with the depth of all streams in the stream network in m
$inpath//$river_width_grid      	river_width               0    		# grid with the witdh of all streams in m
$inpath//$river_links_grid      	river_links               0    		# grid with codes of tributaries, from which a channel was routed (only for real routing channels!!!)
$inpath//$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//$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//$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//$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//$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//$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//$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//$gw_kolmation_1        	gw_kolmation_1            fillcode = 1 # defaultValue = 1e-5   writecode = 8 readcode = 0 # kolmation (leakage factor) between 1st and 2nd aquifer
$inpath//$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//$slidefraction1		slidefraction1		  fillcode = 0 # new snow model extensions: fraction of the area flowing to the North (slides)
$inpath//$slidefraction2		slidefraction2		  fillcode = 0 # new snow model extensions: fraction of the area flowing to the West (slides)
$inpath//$slidefraction3		slidefraction3		  fillcode = 0 # new snow model extensions: fraction of the area flowing to the East (slides)
$inpath//$slidefraction4		slidefraction4		  fillcode = 0 # new snow model extensions: fraction of the area flowing to the South (slides)
$inpath//$depositionindex          deposition_index          fillcode = 1 # new snow model extensions: correction factor for wind impact on snow fall
$inpath//$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]
4                                                     # 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//$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//$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
7 #$Writegrid                                             # Writegrid for glacerized cells
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
7 #$Writegrid                                             # Writegrid for glacier codes
1 #$readgrids           # should always be 1 since otherwise no glacier zones could be created in the dynamic glacier model
$outpath//$pondgridtopmodel  ponding_storage_top  0  0 # effectice for wasimtop only: pond grid for lake modelling, nodata values grid must not be filled 
$Writegrid                                             # effectice for wasimtop only: Writegrid for topmodel-ponds
0 #$readgrids                                          # effectice for wasimtop only: 0, if ponds should be initialized in routing model by the volume-waterlevel relation, 1 if actual pond content should be read in from existing pond grid

$outpath//$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?)

[model_time]
24                     #  start hour
1                     #  start day
1                     #  start month
$startyear        #  start year
24                    #  end hour
1                    #  end day
1                    #  end month
$endyear             #  end year

[meteo_data_count]
3

[meteo_names]
temperature	  # the name of the temperature interpolation result is mandatory if dynamic phenology is used (calculating forcing units depends on a grid called "temperature")
precipitation_12
precipitation_idw
air_humidity
vapor_pressure
global_radiation
wind_speed
sunshine_duration
air_humidity
# 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]
12                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined, 7=bicubic spline, 8=bicubic splines of gradients and residuals linearly combined, 9=read grids according to the name in a grid list file
$inpath//tmean.txt AdditionalColumns=0 # file name with station data (if method = 1, 3 or 4, else ignored)
-0.0066 # method 12 parameter: temperature gradient
#1200 1400 200 1 300 # Method 10 parameter: lower inversion [m asl], upper inversion [m asl], tolerance [m], overlap [0/1 for true/false], clusterlimit [m] 
#$inpath//tmean.reg # Method 1 or 3 paremeter: file name with regression data (if method = 2 or 3)
$outpath//$tempegrid          # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages)
$Writegrid                 # 0, if no grid-output is needed, else one of the codes described above
1.0                           # correction faktor for results
$outpath//temp//$grid//.//$code//$startyear $hour_mean # 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
1.0                           #  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
-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

[wind_speed]
2                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined
$inpath//wind//$startyear//a.dat AdditionalColumns=0  # file name with station data (if method = 1, 3 or 4, else ignored)
$inpath//wind//$startyear//.out    # file name with regression data (if method = 2 or 3)
$outpath//$windgrid           # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages)
$Writegrid                    # 0, if no grid-output is needed, else one of the codes described above
0.1                           # correction faktor for results
$outpath//wind//$grid//.//$code//$startyear $hour_mean  # file name for the statistic output (statially averaged values per time step and subcatchment...)
998                           # error value: all data in the input file greater than this values or lesser the negative value are nodata
$IDWweight                    # weighting of the reciprocal distance for IDW
0.3                           # for interpolation method 3: relative weight of IDW-interpolation in the result
$IDWmaxdist                   # max. distance of stations to the actual interpolation cell
$Anisoslope                   # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive)
$Anisotropie                  #  ratio of the short to the long axis of the anisotropy-ellipsis
0                             # lower limit of interpolation results
0                             # replace value for results below the lower limit
90                            # upper limit for interpolation results
90                            # replace value for results with larger values than the upper limit
$SzenUse                      # 1=use scenario data for correction, 0=dont use scenarios
3                             # 1=add scenarios, 2=multiply scenarios, 3=percentual change
4                             # number of scenario cells

[precipitation_12]
12                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined 7=spline 8=spline+regress
$inpath//precip.txt AdditionalColumns=0   # file name with station data (if method = 1,3,4,5,6,7,8 or 9 else ignored)
0.007 2070 -0.004 P-type # method 12 parameter: altitudinal gradient in (mm/d)/m --> 1mm per 100m 
#4000 4100 200 1 300 # regression parameters for precipitation
#$inpath//precip.reg   # file name with regression data (if method = 2 or 3)
$outpath//$preci_12_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_12//$grid//.//$code//$startyear $hour_mean  # 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.80                          # for interpolation method 3: relative weight of IDW-interpolation in the result
$IDWmaxdist                   # max. distance of stations to the actual interpolation cell
$Anisoslope                   #  slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive)
$Anisotropie                  #  ratio of the short to the long axis of the anisotropy-ellipsis
0.1                           # lower limit of interpolation results
0                             # replace value for results below the lower limit
900                           # upper limit for interpolation results
900                           # replace value for results with larger values than the upper limit
$SzenUse                      # 1=use scenario data for correction, 0=dont use scenarios
3                         # 1=add scenarios, 2=multiply scenarios, 3=percentual change
1 # 4                         # number of scenario cells
699000 235000    0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5     # koordinates of the cells, then one value for each month of a year

[precipitation_idw]
1                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined 7=spline 8=spline+regress
$inpath//precip.txt AdditionalColumns=0   # file name with station data (if method = 1,3,4,5,6,7,8 or 9 else ignored)
#0.01 # method 12 parameter: altitudinal gradient in (mm/d)/m --> 1mm per 100m 
4000 4100 200 1 300 # regression parameters for precipitation
#$inpath//precip.reg   # file name with regression data (if method = 2 or 3)
$outpath//$preci_idw_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.25                           # correction faktor for results
$outpath//prec_idw//$grid//.//$code//$startyear $hour_mean  # 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.80                          # for interpolation method 3: relative weight of IDW-interpolation in the result
$IDWmaxdist                   # max. distance of stations to the actual interpolation cell
$Anisoslope                   #  slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive)
$Anisotropie                  #  ratio of the short to the long axis of the anisotropy-ellipsis
0.1                           # lower limit of interpolation results
0                             # replace value for results below the lower limit
900                           # upper limit for interpolation results
900                           # replace value for results with larger values than the upper limit
$SzenUse                      # 1=use scenario data for correction, 0=dont use scenarios
3                         # 1=add scenarios, 2=multiply scenarios, 3=percentual change
1 # 4                         # number of scenario cells
699000 235000    0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5     # koordinates of the cells, then one value for each month of a year

[sunshine_duration]
1                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined
$inpath//sonne//$startyear//a.rel AdditionalColumns=0  # file name with station data (if method = 1, 3 or 4, else ignored)
$inpath//sonne//$startyear//.out   # file name with regression data (if method = 2 or 3)
$outpath//$sunshinegrid       # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages)
$Writegrid                    # 0, if no grid-output is needed, else one of the codes described above
1.0                           # correction faktor for results
$outpath//ssd_//$grid//.//$code//$startyear  $hour_mean  # file name for the statistic output (statially averaged values per time step and subcatchment...)
998                           # error value: all data in the input file greater than this values or lesser the negative value are nodata
$IDWweight                    # weighting of the reciprocal distance for IDW
0.5                           # for interpolation method 3: relative weight of IDW-interpolation in the result
$IDWmaxdist                   # max. distance of stations to the actual interpolation cell
$Anisoslope                   # slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive)
$Anisotropie                  # ratio of the short to the long axis of the anisotropy-ellipsis
0                             # lower limit of interpolation results
0                             # replace value for results below the lower limit
1.0                           # upper limit for interpolation results
1.0                           # replace value for results with larger values than the upper limit
$SzenUse                      # 1=use scenario data for correction, 0=dont use scenarios
3                             # 1=add scenarios, 2=multiply scenarios, 3=percentual change
1                             # number of scenario cells

[global_radiation]
2                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined
$inpath//Kglob//$startyear//a.dat AdditionalColumns=0  # file name with station data (if method = 1, 3 or 4, else ignored)
$inpath//global//$startyear//.out  # file name with regression data (if method = 2 or 3)
$outpath//$radiationgrid      # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages)
$Writegrid                    # 0, if no grid-output is needed, else one of the codes described above
1.0                           # correction faktor for results
$outpath//rad_//$grid//.//$code//$startyear $hour_mean  # file name for the statistic output (statially averaged values per time step and subcatchment...)
13000                         # error value: all data in the input file greater than this values or lesser the negative value are nodata
$IDWweight                    # weighting of the reciprocal distance for IDW
0.5                           # for interpolation method 3: relative weight of IDW-interpolation in the result
$IDWmaxdist                   # max. distance of stations to the actual interpolation cell
$Anisoslope                   #  slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive)
$Anisotropie                  #  ratio of the short to the long axis of the anisotropy-ellipsis
0                             # lower limit of interpolation results
0                             # replace value for results below the lower limit
1367                          # upper limit for interpolation results
1367                          # replace value for results with larger values than the upper limit
$SzenUse                      # 1=use scenario data for correction, 0=dont use scenarios
1                             # 1=add scenarios, 2=multiply scenarios, 3=percentual change
4                             # number of scenario cells

[air_humidity]
2                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined
$inpath//RH.txt AdditionalColumns=0  # file name with station data (if method = 1, 3 or 4, else ignored)
$inpath//RH.reg  # file name with regression data (if method = 2 or 3)
$outpath//$humiditygrid       # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages)
$Writegrid                    # 0, if no grid-output is needed, else one of the codes described above
1.00                          # correction faktor for results
$outpath//humi//$grid//.//$code//$startyear $hour_mean # file name for the statistic output (statially averaged values per time step and subcatchment...)
9998                          # error value: all data in the input file greater than this values or lesser the negative value are nodata
$IDWweight                    # weighting of the reciprocal distance for IDW
0.5                           # for interpolation method 3: relative weight of IDW-interpolation in the result
$IDWmaxdist                   # max. distance of stations to the actual interpolation cell
$Anisoslope                   #  slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive)
$Anisotropie                  #  ratio of the short to the long axis of the anisotropy-ellipsis
0.01                          # lower limit of interpolation results
0.01                          # replace value for results below the lower limit
1.0                           # upper limit for interpolation results
1.0                           # replace value for results with larger values than the upper limit
$SzenUse                      # 1=use scenario data for correction, 0=dont use scenarios
3                             # 1=add scenarios, 2=multiply scenarios, 3=percentual change
1                             # number of scenario cells
699000 235000    0 0 0 0 0 0 0 0 0 0 0 0


[vapor_pressure]
2                             # method: 1=idw 2=regress 3=idw+regress 4=thiessen 5=bilinear 6=bilinear gradients and residuals linarly combined
$inpath//dampf//$startyear//a.dat AdditionalColumns=0  # file name with station data (if method = 1, 3 or 4, else ignored)
$inpath//dampf//$startyear//.out   # file name with regression data (if method = 2 or 3)
$outpath//$vaporgrid          # name of the output grid (is also used for deriving names of daily, monthly, yearly sums or averages)
$Writegrid                    # 0, if no grid-output is needed, else one of the codes described above
0.1                           # correction faktor for results
$outpath//vapo//$grid//.//$code//$startyear $hour_mean # file name for the statistic output (statially averaged values per time step and subcatchment...)
998                           # error value: all data in the input file greater than this values or lesser the negative value are nodata
$IDWweight                    # weighting of the reciprocal distance for IDW
0.5                           # for interpolation method 3: relative weight of IDW-interpolation in the result
$IDWmaxdist                   # max. distance of stations to the actual interpolation cell
$Anisoslope                   #  slope of the mean axis of the anisotropy-ellipsis (-90 ... +90 degree, mathem. positive)
$Anisotropie                  #  ratio of the short to the long axis of the anisotropy-ellipsis
0                             # lower limit of interpolation results
0                             # replace value for results below the lower limit
90                            # upper limit for interpolation results
90                            # replace value for results with larger values than the upper limit
$SzenUse                      # 1=use scenario data for correction, 0=dont use scenarios
1                             # 1=add scenarios, 2=multiply scenarios, 3=percentual change
4                             # number of scenario cells
699000 235000  0.994   1.187   1.021   1.035   1.201   1.021    .635    .566    .538   1.021    .800   1.007 # 4 scenario cells for vapor pressure
699000 290000  1.021   1.201   1.035   1.049   1.201   1.035    .649    .566    .538   1.035    .787   1.021
737000 235000  0.980   1.173   1.035   1.076   1.173   1.021    .621    .635    .593   1.035    .800   1.007
737000 290000  1.007   1.187   1.049   1.090   1.173   1.035    .649    .649    .607   1.049    .800   1.007


# ----------  parameter for model components  -----------------
#

# for precipitation correction the paragraphs "precipitation" "temperature" and
# "wind_speed" are searched in the memory. If thea are not there (no definition in the control file for precipitation, wind or temperature),
# the prec. corr. will not be calculated

[RegionalSuperposition]
1
$time
NumberOfEntities = 1;
precipitation {
	entityinputgrid = precipitation_12 ;
		regions = 1   ;
		weights = 0.3 ;
	entityinputgrid = precipitation_idw ;
		regions = 1   ;
		weights = 0.7 ;
	outputgrid  = $outpath//$preci_grid ; 
		writecode = $Writegrid ;        
	outputtable = $outpath//prec//$grid//.//$code//$startyear;
		statcode  = $hour_mean;
}
temperature {
	entityinputgrid = temperature_reg;
		regions = 1   2   ;
		weights = 1.0 1.0 ;
	entityinputgrid = temperature_idw;
		regions = 1   2   ;
		weights = 0.0 0.0 ;
	outputgrid  = $outpath//$tempegrid ; 
		writecode = $Writegrid ;        
	outputtable = $outpath//temp//$grid//.//$code//$startyear ;
		statcode  = $hour_mean ;
}

# for precipitation correction the paragraphs "precipitation" "temperature" and
# "wind_speed" are searched in the memory. If thea are not there (no definition in the control file for precipitation, wind or temperature),
# the prec. corr. will not be calculated

[precipitation_correction]
1                     # 0=ignore this module, 1 = run the module
1.9                   # Snow-rain-temperature
1.01                  # liquid:   b in:  y = p(ax + b)
0.00                  # 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.00                  # 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
$Writegrid        # Writegrid for corrected temperatures
2                    # factor x for temperature correction x * (-1.6 .... +1.6)
$outpath//$ExpoCorrgrid   # name of the grids with the correction factors for the direct radiation
$Writegrid           # Writegrid
$outpath//$Shapegrid # name of the grids for codes 1 for theor. shadow, 0 for theor. no shadow (day; assumed: SSD=1.0)
$Writegrid           # Writegrid
1                    # interval counter, after reaching this value, a new correction is calculated (3=all 3 hours a.s.o.)
24                   # 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
2                        # Method: 1=Penman-Monteith, 2=Hamon (only daily), 3=Wendling (only daily) 4= Haude (only daily)
0.65  0.75  0.85  1.2  1.3 1.4 1.5 1.4 1.3 1.1 0.9 0.74  # 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
$Writegrid            # 0, if no grid-output is needed, else one of the codes described above
$outpath//etp_//$grid//.//$code//$startyear $hour_mean # statisticfile for Teilgebiete of pot. evapo-Transpiration
$outpath//$ETRgrid       # result grid for real evapotranspiration in mm/dt
$Writegrid            # 0, if no grid-output is needed, else one of the codes described above
$outpath//etr_//$grid//.//$code//$startyear $hour_mean # statistic for subcatchments (zones) of the real evapotranspiration
$outpath//$EVAPgrid       # result grid for real evapotranspiration in mm/dt
$Writegrid            # 0, if no grid-output is needed, else one of the codes described above
$outpath//evap//$grid//.//$code//$startyear $hour_mean # statistic for subcatchments (zones) of the potential evaporation
$outpath//$EVARgrid       # result grid for real evapotranspiration in mm/dt
$Writegrid            # 0, if no grid-output is needed, else one of the codes described above
$outpath//evar//$grid//.//$code//$startyear $hour_mean # statistic for subcatchments (zones) of the real evaporation
$outpath//$ETRSgrid	 # result grid for real snow evapotranspiration in mm/dt
$Writegrid        	 # 0, if no grid-output is needed, else one of the codes described above
$outpath//etrs//$grid//.//$code//$startyear $hour_mean # statistic for subcatchments (zones) of the real snow evaporation
$outpath//$EIPgrid       # result grid for pot. interception evaporation in mm/dt
$Writegrid            # 0, if no grid-output is needed, else one of the codes described above
$outpath//eip_//$grid//.//$code//$startyear $hour_mean # statisticfile for zones of pot. interception evaporation
$outpath//rgex//$grid//.//$code//$startyear $hour_mean # statistic for subcatchments (zones) of the corrected radiation
+0.23   +1.77    -2.28    +1.28    # coefficients c for Polynom of order 3 RG = c1 + c2*SSD + c3*SSD^2 + c4*SSD^3
+0.072  -0.808   +2.112   -0.239   # coefficients x for Polynom of order 3 SSD = x1 + x2*RG + x3*RG^2 + x4*RG^3
0.88 0.05                          # Extinktion coefficient for RG-modeling (Phi and dPhi) (summer phi = phi-dphi, winter phi=phi+dphi)
1654.0                             # recession constant (e-function for recession of the daily temperature amplitude with altitude [m]
3.3  4.4  6.1  7.9  9.4  10.0  9.9  9.0  7.8  6.0  4.2  3.2  # monthly values of the max. daily T-amplitudes (for 0 m.a.s.l)
0.62  0.1                          # part of the temperature amplitude (dt), that is added to the mean day-temperature # (followed by the range of changing within a year ddt) to get the mean temperature of light day in the night: mean night temperature is mean day temperature minus (1-dt)*(temp. amplitude)
$outpath//cloud//$grid//.//$code//$suffix $hour_mean 	# 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
1                    # method 1=T-index, 2=t-u-index, 3=Anderson comb., 4=extended com.
3.0                  # transient zone for rain-snow (T0R +- this range)
1.9                  # T0R    temperature limit for rain (Grad Celsius)
1.9                  # T0     temperature limit snow melt
0.02                 # CWH    storage capacity of the snow for water (relative part)
0.2                  # CRFR   coefficient for refreezing
1.65                  # C0     degree-day-factor mm/d/C
1.1                  # C1     degree-day-factor without wind consideration  mm/(d*C)
1.0                  # C2     degree-day-factor considering wind mm/(d*C*m/s)
0.07                 # z0     roughness length cm for energy balance methods (not used)
2.3                  # RMFMIN minimum radiation melt factor      mm/d/C comb. method
3.6                  # RMFMAX maximum radiation melt factor      mm/d/C comb. method
0.40                 # Albedo for snow (Min)
0.85                 # Albedo for snow (Max)
$outpath//$rain_rate         # rain rate
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//rain//$grid//.//$code//$startyear $hour_mean # rain rate
$outpath//$snow_rate         # snow rate
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//snow//$grid//.//$code//$startyear $hour_mean # snow rate
$outpath//$days_snow         # days with snow (SWE > 5mm)
$Writegrid                # 0, if no grid-output is needed, else one of the codes described above
$outpath//sday//$grid//.//$code//$startyear $hour_mean # days with snow (SWE > 5mm)
$outpath//$snow_age          # snow age (days without new snow)
$Writegrid                # 0, if no grid-output is needed, else one of the codes described above
$outpath//sage//$grid//.//$code//$startyear $hour_mean # days since last snowfall
$outpath//albe//$grid//.//$code//$startyear $hour_mean # Albedo
$outpath//$snowcover_outflow # discharge from snow, input (precipitation) for following modules
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//qsch//$grid//.//$code//$startyear $hour_mean # melt flow (or rain, if there is no snow cover) in mm/dt
$outpath//$melt_from_snowcover # discharge from snow, input (precipitation) for following modules
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//qsme//$grid//.//$code//$startyear $hour_mean # melt flow in mm/dt
$outpath//$SSNOgrid        # name of the grids with the snow storage solid in mm
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$SLIQgrid          # name of the grids with the snow storage liquid in mm
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//ssto//$grid//.//$code//$startyear $hour_mean # total snow storage, in mm, (liquid and solid fraction)
$outpath//$SSTOgrid          # name of the grids with the total snow storage solid AND liquid in mm
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$readgrids                   # 1=read snow storage solid, liquid grids from disk, 0=generate new grids
# 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//$suffix $hour_mean 	# 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//$suffix $hour_mean 	# 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//$suffix $hour_mean 	# 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//$suffix $hour_mean 	# 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//$suffix $hour_mean 	# 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
0.5	# scaling factor for energy balance (taking into account that the complete snow pack is modelled as a single bulk mass instead of as several layers)

# 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 = like 12 but melt methods are taken from snow model (for snow only)
[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)
5.5   # t-index factor for ice
3.3   # t-index factor for firn
1.65  # t-index factor for snow
1.5   # melt factor (used for snow, firn and ice when using radiation correction)
0.0001   # radiation coefficient for ice_min  (for method 2)
0.0007   # radiation coefficient for ice_max  (for method 2)
0.0001   # radiation coefficient for snow_min (for method 2)
0.00055  # radiation coefficient for snow_max (for method 2)
8   # els-konstante for ice
60  # els-konstante for firn
25  # els-konstante for snow
0.1  # initial reservoir content for ice discharge (single linear storage approach)
0.1  # initial reservoir content for firn discharge (single linear storage approach)
0.1  # initial reservoir content for snow discharge (single linear storage approach)
$outpath//$firn_melt # melt from firn
$Writegrid           # 0, if no grid-output is needed, else one of the codes described above
$outpath//qfir//$grid//.//$code//$startyear $hour_mean # melt from firn as statistic file
$outpath//$ice_melt  # melt from ice
$Writegrid           # 0, if no grid-output is needed, else one of the codes described above
$outpath//qice//$grid//.//$code//$startyear $hour_mean # melt from ice as statistic file
$outpath//qglc//$grid//.//$code//$startyear $hour_mean # melt from ice and firn as statistic file
# -----------------------------------------------------------------------------
# now some new parameters for the new dynamic glacier model (methods 11 and 12)
$outpath//qsgl//$grid//.//$code//$startyear $hour_mean # melt from snow from glacier only as statistic file (but still with respect tothe subbasins areas!) --> new in version 8.07.00
$readgrids                # 1=read grids and stacks from disk, 0=generate new grids and stacks (using the parameters in the following line for WE_Firn stack)
7 2500 1.8 # number of layers for the firn stack, followed by two initialization parameters: average Equilibrium line elevation in m (e.g. 2500) and change rate of WE per m in mm (e.g. 2) -> every 100m the WE of firn in each layer will grow by 200mm
09 30 # month and day (hour is set automatically to 24) for which the Volume-Area-Relation will be applied newly (and temporary (i.e. internal) Balances are reset to 0)
28.5 1.36 8 7 # 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//$startyear $hour_mean # water equivalent for firn as statistics file (sum over all firn layers)
$outpath//$GlacierMassBalance # output grid with mass balance of the glacier
3                         # 3: write at end of simulation (important to start another model run with correct initialization values)
$outpath//$OldGlacierMassBalance # output grid with mass balance of the glacier (last year; used for calculating the next years balance correctly when starting somewhen inbetween two balancing days)
3                         # 3: write at end of simulation (important to start another model run with correct initialization values)
$outpath//glmb//$grid//.//$code//$startyear $hour_mean # 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//$startyear $hour_mean # 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//$suffix $hour_mean 	# 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
500 0.25 # additional parameter for versions >= 9.09.08 that describes the SWE on a cell that must be reached before the exceeding SWE is either out to an existing glacier on that cell or a new glacier cell is created (if glid is nodata for that cell). Followed by the fraction of the cell that should be turned into glacier (>0...1.0)

# Layered Snow model is a component that can be used within unsaturated zone model. It handles the snow like soil layers, so heat transfer and water movement will be computed as in a soil.
# therefor it is required to have a section [LayeredSoilProperties] in the control which defines the thermal and hydraulic prperties of snow (like a soil type)
[LayeredSnowModel]
1               # do_LayeredSnow: 0=do not integrate snow pack into the soil model; 1=integrate snow pack into soil model, using several layers of snow (supports heat transfer and Richards-approach of soil/snow water movement)
1		# LayeredSnowMethod: method of dealing with snow layers: 1 = constant thickness, 2=constant density (still under construction)
0.1		# --> not yet active: LayeredSnowMinDensity: minimum density of freshly fallen dry snow (will later be modeled in more detail, i.e. depending on temperatur and wind or so)
0.85		# --> not yet active: LayeredSnowMaxDensity: maximum density of old snow (will later be modeled dynamically depending on snow age, temperature and water content etc.)
-2.0		# recommended: -2.0 ... 0.0 temperature difference of freshly fallen snow to Tair (because snow comes from higher atmospheric layers, the actual snow temperature may be lower than the air temperature in 2m elevation -> otherwise snow fall at 2°C weren't possible). This parameter leads to a cooler snow pack and less melt when the snow pack starts to accumulate (and soils are still warm as well)
0.0		# --> not yet active: scaling parameter for snow internal runoff (like interflow, in analogy to drainage density in interflow module) --> use with care, only effective if k_recession < 1
0.5		# maximum water storage capacity of snow (this parameter exists because snow is looked at as a soil type, so the soil matrix may change into water due to warming up and become "soil moisture" itself). To make the numerical solution more stable, this parameter should be used (0.1...0.2)
1.0		# --> set always to 1.0, otherwise unpredictable behavior! SE effectiveness for snow: effectiveness of SE on reduction of hydraulic conductivity and suction for snow. When set to 0, hydraulic properties of SNOW only are like without snow model (only heat transfer uses SE then). When set to 1, the effectiveness is 100% like in older WaSiM versions. Range 0...1. Recommended: 1
0.0		# SE effectiveness for soil: effectiveness of SE on reduction of hydraulic conductivity and suction for soil. When set to 0, hydraulic properties of SOIL only are like without snow model (only heat transfer uses SE then). When set to 1, the effectiveness is 100% like in older WaSiM versions. Range 0...1. Recommended for first using layered snow model when the calibrated model didn't use heat transfer yet: 0, else 1
# some stacks grids that may be of interest 
$outpath//$SWESstack		# stack with SWE solid values per snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$SWELStack		# stack with SWE liquid values per snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$temperSnowStack	# stack with snow temperatures per snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$thetaSnowStack	# stack with theta values per snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$hhSnowStack		# stack with hydraulic heads per snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$hgeoSnowStack	# stack with geodetic elevation per snow layer (bottom of layer)
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$quSnowStack		# stack with (bottom) water fluxes per snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$thickSnowStack	# stack with thickness of each snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$densSnowStack	# stack with density of each snow layer
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//$upperIndexSnow	# grid that holds the index of the first (uppermost) layer that has a positive SWES value, i.e. the first layer with snow in it
1//$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above

# snow layer properties: this is a short soil table, since the snow layers are handled like a soil in the model
# all entries have the same syntax as for the soil table, However, some values may not be used internally or even changed dynamically, see comments
# actually, ONLY ONE layer is really used, though.
[LayeredSnowProperties]
1 # number of entries --> ATTENTION: do not change until the model can use multiple layers
1 Snow { method = MultipleHorizons; 
	 ThicknessScaling	= 1.0 ; # use with care - effects on SWE may be unpredictable when changing this parameter value
	 PMacroThresh	= 0 ; # better to use no macropores for snow (not implemented anyway)
	 MacroCapacity	= 0 ; # better to use no macropores for snow (not implemented anyway)
	 CapacityRedu	= 0 ; # better to use no macropores for snow (not implemented anyway)
	 MacroDepth		= 0 ; # better to use no macropores for snow (not implemented anyway)
	 EvapMaxDepth	= 0 ; # there is no evaporation from snow that counts as soil evaporation. 
	 horizon		= 1 ; # ATTENTION: do not use more than one horizon for snow! The layers will be filled up dynamically from bottom to top and also some of the properties will be changed by the model itself. So one horizon is sufficient for now
	 Name		= Snow ; 
	 DryHeatCapacity	= 2000 ; # dry heat capacity in J/(Kg*K) (This needs to be checked, I used Cvol/density)
	 DryDensity 	= 400  ; # dry density in kg/m^3
	 DryThermalConduct	= 2.33 ; # dry thermal conductivity in W/(m*K) (or J/(m*s*K) --> 1J = 1Ws)
	 KMinFrozenSoil	= 1e-6;  # minimum hydraulic conductivity in m/s when the soil is coimpletely frozen -> should be quite high for snow
	 ksat		= 1e-4; # like a highly porous sand or small gravel
	 k_recession	= 1.0  ; # although the saturated hydraulic conductivity decreases, there will be no interflow, only increaed storage in the lower layers
	 theta_sat		= 0.7  ; # (should be 1.0) this will be dynamically adjusted (at least later, wehn the dynamic snow layers are implemented completely)
	 theta_res		= 0.1  ; # (should be 0.0) 
	 alpha		= 18   ; 12  # like ?
	 Par_n		= 4    ; 3   # like a highly porous sand or small gravel
	 Par_tau		= 0.5  ; # like a highly porous sand or small gravel (should always be 0.5 anyway)
	 thickness		= 0.3  ; # this can be changed if required (measure in meters). Do not use too thin layers, since the performance will go down dramatically then (recommended minimum 0.1m, better use 0.2 m)
	 layers		= 30   ; # this can be changed if required (1...inf)
}

# 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
$Writegrid                   # 0, if no grid-output is needed, else one of the codes described above
$outpath//qi__//$grid//.//$code//$startyear $hour_mean # 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//$startyear $hour_mean # 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//$startyear $hour_mean # zonal statistic For interception storage content
0.35                         # layer thickness of the waters on the leaves (multiplied with LAI -> storage capacity)
$readgrids                   # 1=read grids from disk, else generate internal

[infiltration_model]
0                            # 0=ignore this module, 1 = run the module
$time                        # duration of a time step in minutes
$outpath//$INFEXgrid         # grid with infiltration excess in mm (surface runoff)
$Writegrid                # Writegrid for surface discharge (fraction 1)
$outpath//infx//$grid//.//$code//$startyear $hour_mean # 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
$outpath//sdis//$grid//.//$code//$startyear $hour_mean # statistics for silting up disposition (Verschlämmungsneigung)
$outpath//qdsu//$grid//.//$code//$startyear $hour_mean # 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<P<76.2 and P>=76.2, respectively
	X = A; # start infiltration rate; here, A was taken from the soil table, since parameters SU_PAR01 ... SU_PAR10 are mapped to internal variables A to J, see description below. Also possible: X = 65.2; but then no variation for different soil types would be possible
	Y = B*K^C*(L*100)^D; # end infiltration rate
	C1 = (100*L)^F; 	# F = SU_PAR06  C1 will be stored in a new internal variable
	C2 = K^G; 			# G = SU_PAR07 C2 will be stored in a new internal variable
	C3 = (O+0.001)^H; 	# H = SU_PAR08  C3 will be stored in a new internal variable
	Z = (O<=0) + (O>0)*(E*C1*C2*C3); # E = SU_PAR05, O = time since last soil tillage, see below
	# attention: only starting with V= and following following expressions will be called after internal update of Q. V must be set only after this internal update, but any other expression may be placed herunder for preparation of the V-call. However, they will be called before the internal update of Q, so they should not touch any of the variables needed for EKIN update
	V = ((X-Y)*exp(-Z*Q)+Y)*R/60; # potential infiltration, will be limited internally by available precipitation. 
}
# Short description of the expression parser and the expression list syntax for method 3
# -	Expressions can be defined following algebraic rules:
# -	Each line contains a single expression which must be closed with a semi colon.
# -	Each assignment (e.g. A = 15) results in creating or updating a value in the internal variable list.
# -	A number of values is already defined by WaSiM (as interface from the calling module), and WaSiM expects some other values to be defined after all expressions were called
# The expression parser  is based on the source code of the expression parser used in SpeQ Mathematics (http://www.speqmath.com/tutorials/expression_parser_cpp/index.html), 
# written by Jos de Jong, 2007. It was adopted to the usage in WaSiM by simplifying the error handling (exceptions are to be handled by WaSiM), 
# extracting the variable list as an external class (to be handled by WaSiM) and some other minor technical changes
#     Operators (ascending precedence per line, no precedence within a line):
#        & | << >>   			(AND, OR, BITSHIFTLEFT, BITSHIFTRIGHT)
#        = <> < > <= >=			(EQUAL, UNEQUAL, SMALLER, LARGER, SMALLEREQ, LARGEREQ)
#        + -				(PLUS, MINUS)
#        * / % ||			(MULTIPLY, DIVIDE, MODULUS, XOR)
#        ^				(POW)
#        !				(FACTORIAL)
#    Functions (must be used with brackets):
#        Abs(arg), Exp(arg), Sign(arg), Sqrt(arg), Log(arg), Log10(arg)
#        Sin(arg), Cos(arg), Tan(arg), ASin(arg), ACos(arg), ATan(arg)
#        Factorial(arg)
#    Variables:
#        Pi, Euler (not only e, e is a predefined variable used by WaSiM to deliver a value to the expression parser interface)
#        you can define your own variables, even with with more than one significant character length, e.g. Inf0 or Help etc.
#	there is no distinction between upper and lower case in function names and variables.
#    Other:
#        Scientific notation supported
#
# ====> what values WaSiM defines forinput (can be used in any expression)
# A to J: values as used in soiltable with names SU_PAR01 to SU_PAR10
# K: grain size distribution Dg, internally calculated after 
#	double FClay   = log004+log2;
#	double FSilt   = 0.3326 * (log2+log6_3) +  0.3348 *(log6_3+log20)  + 0.1704 * (log20+log36)  + 0.1622 * (log36+log63) ;
#	double FSand   = 0.1336 * (log63+log100) + 0.2005 *(log100+log200) + 0.3318 *(log200+log630) + 0.3341 *(log630+log2000);
#	double FStones1 = (log2000+log6300);
#	double FStones2 = (log6300+log20000);
#	double FStones3 = (log20000+log63000);
#	double FStones4 = (log63000+log200000);
#	double dg = (FClay*dFractionClay + FSilt*dFractionSilt + FSand*dFractionSand + FStones1*dFractionStones1 + FStones2*dFractionStones2 + FStones3*dFractionStones3 + FStones4*dFractionStones4)/2.0;
# 	with fractions of each grain size class taken from the soil table
# L: 	fraction of sand
# M: 	fraction of clay
# N: 	fraction of silt
# O: 	t_cult, time since last soil cultivation (in days)
# P: 	rain intensity in mm/h, taken from precipitation input
# Q: 	e_kin: accumulated cinetic energy: for all expressions resulting in W, X, Y or Z: result value of the last time time step; for V: value of the actual time step
# R:	internal time step in minutes
# ====> What WaSiM expects for output: (ranging from Z downwards, will be used by WaSiM when going ahead)
# Z: 	silting up disposition SDISP
# Y: 	end infiltration rate i_inf
# X: 	start infiltration rate i0
# W:	actual cinetic energy
# V: 	potential infiltration rate inf_pot, depending on energy, siting up disposition, inf_start and inf_infinite
# order of expressions evanulated by WaSiM: 
# 	expressions returning W, X, Y and Z are independently of each other.
# 	expression V must be called as last call in any case, since WaSiM will update EKIN internally using the energy-result (in W) and V depends on all the other results W to Z
# 	other expressions for storing intermediate results may be defined at any position in the expression list before the results will be used in another expression

[SurfaceRoutingModel]
0                            # 0=ignore this module, 1 = run the module
$time                        # duration of a time step in minutes
2                            # method: 1=MultipleFlowPaths for diverging areas, 2=single flowpaths (nearest direction as given by aspect)
$outpath//qdsr//$grid//.//$code//$startyear $hour_mean # direct discharge from surface routing module 
$outpath//qisr//$grid//.//$code//$startyear $hour_mean # interflow from surface routing module
$outpath//qbsr//$grid//.//$code//$startyear $hour_mean # baseflow from surface routing module 
$outpath//qgsr//$grid//.//$code//$startyear $hour_mean # total discharge from surface routing module 
$outpath//$surfspeed_grid    # grid with actual flow velocity of surface flow in m/s
$Writegrid                   # writegrid for this grid
$outpath//$surfflux_grid     # grid with actual flow amounts of surface flow in m^3/s
$Writegrid                   # writegrid for this grid
0.001                        # maximum wake lenght iteration difference (if Delta_A_nl < this value, iteration for a_NL stops)
40                           # maximum number of iterations for a_NL
0.0001                       # maximum flow velocity iteration difference (if Delta v is less than this value, iteration stops)
40                           # maximum number of iterations for v
30                           # shortest sub-time step in seconds
3600                         #longest allowed sub time step (even if flow travel times are longer, the time step is subdivided into sub timesteps of this lenght) be careful: tracers are mixed much faster when multiple sub time steps are applied
0.02                         # minimum water depth for regarding roughenss of crops in m (shallower sheet flow: only roughness of bare soil will be regarded)
2.0                          # ConcentrationFactor takes into account the micro scale concentration of flow pathes, flow will take place on a fraction of the cell only, so the amount flowing per meter width will be multiplied by this factor (1..n)
$readgrids                   # readgrid code 0 do not read, 1 = read grids 
$outpath//sfstsr//$grid//.//$code//$suffix $hour_mean # 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

$set $kelsqd1 	= 25
$set $kelsqd2 	= 20
$set $kelsqd3 	= 25
$set $kelsqd4 	= 20
$set $kelsqd5 	= 25
$set $kelsqd6 	= 25
$set $kelsqd7 	= 25

$set $kelsqi1 	= 60
$set $kelsqi2 	= 50
$set $kelsqi3 	= 60
$set $kelsqi4 	= 50
$set $kelsqi5 	= 60
$set $kelsqi6 	= 60
$set $kelsqi7 	= 60

$set $dd1 		= 20
$set $dd2 		= 22
$set $dd3 		= 20
$set $dd4 		= 22
$set $dd5 		= 20
$set $dd6 		= 20
$set $dd7 		= 20

$set $surfRO1 	= 0.07
$set $surfRO2 	= 0.07
$set $surfRO3 	= 0.07
$set $surfRO4 	= 0.07
$set $surfRO5 	= 0.07
$set $surfRO6 	= 0.07
$set $surfRO7 	= 0.07

[unsatzon_model]
1                            # 0=ignore this module, 1 = run the module
$time                        # duration of a time step in minutes
3                            # method, 1=simple method (will not work anymore from version 7.x), 2 = FDM-Method 3 = FDM-Method with dynamic time step down to 1 secound
2  # controlling interaction with surface water: 0 = no interaction, 1 = exfitration possible 2 = infiltration and exfiltration possible
0  # controlling surface storage in ponds:       0 = no ponds,       1 = using ponds for surface storage (pond depth as standard grid needed -> height of dams oround fields)
0  # controlling artificial drainage:            0 = no artificial drainage 1 = using drainage (drainage depth and horizontal pipe distances as standard grids needed!)
0  # controlling clay layer:                     0 = no clay layer,  1 = assuming a clay layer in a depth, specified within a clay-grid (declared as a standard grid)
5e-8                         # permeability of the clay layer (is used for the clay layer only)
4  # parameter for the initialization of the gw_level (range between 1..levels (standard: 4))
$outpath//qdra//$grid//.//$code//$startyear $hour_mean # results drainage discharge in mm per zone
$outpath//gwst//$grid//.//$code//$startyear $hour_mean # results groundwater depth
$outpath//gwn_//$grid//.//$code//$startyear $hour_mean # results mean groundwater recharge per zone
$outpath//sb05//$grid//.//$code//$startyear $hour_mean # results rel. soil moisture within the root zone per zone
$outpath//sb1_//$grid//.//$code//$startyear $hour_mean # results rel. soil moisture within the unsat. zone (0m..GW table) per zone
$outpath//wurz//$grid//.//$code//$startyear $hour_mean # results statistic of the root depth per zone
$outpath//infx//$grid//.//$code//$startyear $hour_mean # results statistic of the infiltration excess
$outpath//pond//$grid//.//$code//$startyear $hour_mean # results statistic of the ponding water storage content
$outpath//qdir//$grid//.//$code//$startyear $hour_mean # results statistic of the direct discharge
$outpath//qifl//$grid//.//$code//$startyear $hour_mean # results statistic of the interflow
$outpath//qbas//$grid//.//$code//$startyear $hour_mean # results statistic of the baseflow
$outpath//qges//$grid//.//$code//$startyear $hour_mean # results statistic of the total discharge
$outpath//gwin//$grid//.//$code//$startyear $hour_mean # statistic of the infiltration from surface water into groundwater (from rivers and lakes)
$outpath//gwex//$grid//.//$code//$startyear $hour_mean # statistic of the exfiltration from groundwater into surface water (into rivers and lakes)
$outpath//macr//$grid//.//$code//$startyear $hour_mean # statistic of infiltration into macropores
$outpath//qinf//$grid//.//$code//$startyear $hour_mean # statistic of total infiltration into the first soil layer
$outpath//$SB_1_grid         # grid with actual soil water content for the root zone
$Writegrid                # writegrid for this grid
$outpath//$SB_2_grid         # grid with actual soil water content for the entire unsaturated zone
$Writegrid                # writegrid for this grid
$outpath//$ROOTgrid          # grid with root depth 
$Writegrid                # writegrid for this grid
$outpath//$Thetastack       # stack, actual soil water content for all soil levels
$Writestack                 # Writecode for this stack
$outpath//$hydraulic_heads_stack   # stack, contaiing hydraulic heads
$Writestack                 # Writecode for this stack
$outpath//$geodetic_altitude_stack   # stack, containig geodaetic altitudes of the soil levels (lower boudaries)
$Writestack                 # Writecode for this stack
$outpath//$flowstack        # stack, containing the outflows from the soil levels
$Writestack                 # Writecode for this stack
$outpath//$GWdepthgrid      # grid with groudwaterdepth
$Writegrid                  # writegrid for this grid
$outpath//$GWthetagrid      # grid with theta in GWLEVEL
$Writegrid                  # writegrid for this grid
$outpath//$GWNgrid          # grid with groundwater recharge
$Writegrid               # writegrid for this grid
$outpath//$GWLEVELgrid      # grid with level index of groundwater surface (Index der Schicht)
$Writegrid                  # writegrid for this grid
$outpath//$QDRAINgrid       # grid with the drainage flows
$Writegrid                  # writegrid for this grid
$outpath//$SATTgrid         # grid with code 1=saturation at interval start, 0 no sat.
$Writegrid                  # writegrid for this grid
$outpath//$INFEXgrid        # grid with infiltration excess in mm (surface discharge)
$Writegrid               # writegrid for this grid
$outpath//$QDgrid           # grid with direct discharge
$Writegrid               # writegrid for this grid
$outpath//$QIgrid           # grid with Interflow
$Writegrid               # writegrid for this grid
$outpath//$QBgrid           # grid with baseflow
$Writegrid               # write code for baseflow
$outpath//$GWINgrid         # grid with infiltration from rivers into the soil (groundwater)
$Writegrid               # writegrid for re-infiltration
$outpath//$GWEXgrid         # grid with exfiltration (baseflow) from groundwater (is only generated, if groundwater module is active, else baseflow is in QBgrid)
$Writegrid               # writegrid for exfiltration
$outpath//$act_pond_grid    # grid with content of ponding storge
$Writegrid                  #
$outpath//$UPRISEgrid       # grid with amount of capillary uprise (mm)
$Writegrid                  #
$outpath//$PERCOLgrid       # grid with amount of percolation (mm)
$Writegrid                  #
$outpath//$MACROINFgrid     # grid with amount of infiltration into macropores (mm)
$Writegrid                  #
$outpath//$irrig_grid        # grid with irrigation amount (will be written when irrigation is used, only)
$Writegrid                  # writegrid for this grid (however: will be written when irrigation is used, only)
140 140 # row/col 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//$startyear  # name of a file containing the flows between the layers of the control point
$outpath//thet//$grid//.//$code//$startyear  # name of a file containing the soil moisture as theta values of the layers of the control point
$outpath//hhyd//$grid//.//$code//$startyear  # name of a file containing the hydraulic head of the layers of the control point
$outpath//otherdata//$grid//.//$code//$startyear  # name of a file containing some other water balance data of the control point (non layer data)
$outpath//etrd//$grid//.//$code//$startyear  # name of a file containing the hydraulic head of the layers of the control point
$outpath//intd//$grid//.//$code//$startyear  # name of a file containing the hydraulic head of the layers of the control point
1			2			3			4			5			6			7			8			9			16			21			27			28		54		55			57 			# range for subbasin codes		
$kelsqd1	$kelsqd2	$kelsqd3	$kelsqd4	$kelsqd5	$kelsqd6	$kelsqd7	$kelsqd2	$kelsqd4	$kelsqd6	$kelsqd5	$kelsqd7	3		3		$kelsqd7	$kelsqd7	# kelsqd		
$kelsqi1	$kelsqi2	$kelsqi3	$kelsqi4	$kelsqi5	$kelsqi6	$kelsqi7	$kelsqi2	$kelsqi4	$kelsqi6	$kelsqi5	$kelsqi7	15		15		$kelsqi7	$kelsqi7	# kelsqi		
$dd1		$dd2		$dd3		$dd4		$dd5		$dd6		$dd7		$dd2		$dd4		$dd6		$dd5		$dd7		18		18		$dd7		$dd7		# drainage density		
0.43		0.43		0.43		0.43		0.43		0.43		0.43		0.43		0.43		0.43		0.43		0.43		0.43	0.43	0.43		0.43		# k in qb = Q0 * exp(-k/z) with z = depth to groundwater		
1			1			1			1			1			1			1			1			1			1			1			1			1		1		1			1			# Q0 in the above formula		
$surfRO1	$surfRO2	$surfRO3	$surfRO4	$surfRO5	$surfRO6	$surfRO7	$surfRO2	$surfRO4	$surfRO6	$surfRO5	$surfRO7	0.08	0.08	$surfRO7	$surfRO7	# fraction of surface runoff on snow melt		
$readgrids               # meanings are extended now! read the follwing comments
$outpath//storage_richards.ftz    # if readgrids = 1, then this file contains the contents of the flow travel time zones for interflow and surface flow and for the tracers
300 # minimum dynamic time step in secounds. the smaller this number, the longer the model runs but the results will be more accurate due to a maintained Courant condition
$outpath//step//$grid//.//$code//$startyear $hour_mean # results statistic of the number of substeps 
$outpath//$SUBSTEPSgrid     # grid with number of substeps --> a good idea is to use writecode 5x (e.g. 53) to get the average number of substeps per cell for the model run
$Writegrid               # for substeps, the areal distribution is of interest for the annual average value. This is code 6 as first digit in 2-digit codes. Or use 5 for the entire model run

[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
12.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
120	# 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//$startyear    # results soil temperature for control point
$outpath//ts_avg//$grid//.//$code//$startyear $hour_mean 	# results soil temperature thaw depth or active layer thickness as average value for subbasins
$outpath//$Temperaturestack        		# stack, actual soil water content for all soil levels
$Writegrid                  			# Writecode for this stack
$outpath//$ThawDepthGridTMod        		# grid containing the active layer thickness relative to the soil surface (deepest thawing front in the soil profile)
$Writegrid                  			# Writecode for this grid
$outpath//$SEstack        				# stack, SE values (relative ice saturation)
$Writegrid                  			# Writecode for this stack
# is not used but taken from unsatzon_model: $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
$exchngpath//wasim.inf  # name of the semaphore file to inform wasim that all grids written by the groundwater model are available now
50   # wait interval in ms between scanning the directory for the new semaphore file (wasim for windows will use a second thread to minimize CPU time, whereas wasim as 
# console application will use 100% CPU time while waiting for the output file of the groundwater model. The wait time is then used to minimize disk access
# A follwing version will use a DLL with a memory pointer to the required grid and a flag, which is used by both programs to couple the models. 
# But this is music for the future yet...
H     # Coupling mode: I=each interval, H=each hour, D=each Day, M=each month, Y=each year
60  # time interval in minutes, the external model uses. This is important to convert changes in groundwater level into fluxes as used by WaSiM
1     # number of following grid names which must be available once the semaphore file was written. Each following row (1..n) will contain a symbolic name 
# (grid names from modconst.h) Thus, any grid may be read in, even for other sub models like the boundary conditions as gw_boundary_fix_h_1 for the first aquifer or a changed landuse grid a.s.o.
$exchngpath//gwtable38.grd GWTableExtern 1 0 # the first value is the file name, the second the internal grid name (see English-section in modconst.h ), the third parameter is the fillMissings-parameter (0=no fill, 1=fill with nearest neighbors value), the last ine is the rename(1)/delete(0) parameter
#$exchngpath//bh.grd      gw_boundary_fix_h_1 0 0 # the first value is the file name, the second the internal grid name (see English-section in modconst.h ), the third parameter is the fillMissings-parameter (0=no fill, 1=fill with nearest neighbors value), the last ine is the rename(1)/delete(0) parameter
3     # number of grids (each matching one of the following rows) which should be written when the next synchronisation is due
$exchngpath//gwn.grd      groundwater_recharge  Mean   # hier als Beispiel das Grid mit der Grundwasserneubildung
$exchngpath//gwstand.grd  groundwater_distance  Last   # hier als Beispiel das Grid mit der Grundwasserneubildung
$exchngpath//balance.grd  Balance  Sum  # hier als Beispiel das Grid mit der Bilanz aller Wasserinhaltsänderungen durch die Kopplung (sollte 0 sein)
2     # number of subbasin correlated statistics (mean values) which should be written as table (in ASCII-Format) (this is actually limited to directflow and interflow)
$exchngpath//qdir.table direct_discharge Sum # direct flow per subbasin/zone in mm
$exchngpath//difl.table Interflow Sum # interflow per subbasin/zone in mm
$exchngpath//geofim.inf # name of the semaphore file wasim will write after all of the output above was written
geofim # content of the semaphore file 

[irrigation]
0                                                 # 0=ignore this module, 1 = run the module
$time                                             # duration of a time step in minutes
$outpath//irgw//$grid//.//$code//$startyear $hour_mean # statistic of the irrigation water from groundwater
$outpath//irsw//$grid//.//$code//$startyear $hour_mean # statistic of the irrigation water from surface water

[groundwater_flow]
1                        # 0=ignore the module, 1 = run the module
$time                    # duration of a time step in minutes; doen't change the value unless you have strong reasons to do so!!
1                        # solving method: 1=Gauss-Seidel-iteration (using alpha for control wether it is explicite, partly or fully implicite), 2=PCG (not yet implemented
1000                     # if iterative solving method (1): max.numberof iterations
0.000001                  # if iterative solving method (1): max. changes between two iterations
0.0                      # Alpha for estimation of central differences 0.5 = Crank-Nicholson Method, 0 = fully explicite, 1 = fully implicite
-1.20                    # factor for relaxing the iteration if using iterativemethod (successive over[/under] relaxation) 
$readgrids               # 1=read grids for heads from disk, 0=do not read but initialize with gw-level from unsaturated zone
1                        # number of layers 
140 140 #384                  # r/c coordinates of a control point for all fluxes and for each layer : q0..q4, leakage up and down
$outpath//glog//$grid//.//$code//$startyear  # name of a file containing the flows between of the control point
1                        # use Pond Grid -> this enables the model to use the hydraulic head of a pond in addition to the groundwater itself 0=use traditional method without pond (default), 1=use ponds
$outpath//$head1grid     # (new) grid for hydraulic heads for layer 1
$Writegrid               # writecode for hydraulic heads for layer 1
$outpath//$flowx1grid    # (new) grid for fluxes in x direction for layer 1
$Writegrid               # writecode for flux-x-grid in layer 1
$outpath//$flowy1grid    # (new) grid for fluxes in y direction for layer 1
$Writegrid               # writecode for flux-y-grid in layer 1
$outpath//$GWbalance1grid # (new) grid for balance (difference of storage change vs. balance of fluxes -> should be 0 or the amount of in-/outflows by boundary conditions
13                       # writecode for balance control grid in layer 1 (should be at least one sum grid per year --> Code = 20 or 23 (if old grids must be read in)
$outpath//$head2grid     # (new) grid for hydraulic heads for layer 2
$Writegrid               # writecode for hydraulic heads for layer 2
$outpath//$flowx2grid    # (new) grid for fluxes in x direction for layer 2
$Writegrid               # writecode for flux-x-grid in layer 2
$outpath//$flowy2grid    # (new) grid for fluxes in y direction for layer 2
$Writegrid               # writecode for flux-y-grid in layer 2
$outpath//$GWbalance2grid # (new) grid for balance (difference of storage change vs. balance of fluxes -> should be 0 or the amount of in-/outflows by boundary conditions
13                       # writecode for balance control grid in layer 2 (should be at least one sum grid per year --> Code = 20 or 23 (if old grids must be read in)
$outpath//$head3grid     # (new) grid for hydraulic heads for layer 3
$Writegrid               # writecode for hydraulic heads for layer 3
$outpath//$flowx3grid    # (new) grid for fluxes in x direction for layer 3
$Writegrid               # writecode for flux-x-grid in layer 3
$outpath//$flowy3grid    # (new) grid for fluxes in y direction for layer 3
$Writegrid               # writecode for flux-y-grid in layer 3
$outpath//$GWbalance3grid # (new) grid for balance (difference of storage change vs. balance of fluxes -> should be 0 or the amount of in-/outflows by boundary conditions
13                       # writecode for balance control grid in layer 3 (should be at least one sum grid per year --> Code = 20 or 23 (if old grids must be read in)

[groundwater_tuning]
frozenDepthScaling = 0.0 # if frozen soils impact on groundwater flows should be reduced, this parameter can be set in the range of 0...1 (1 is fully effectiv as without section groundwater_tuning at all)

# this paragraph is not needed for WaSiM-uzr but for the WaSiM-version with the variable saturated area approach (after Topmodel)
[soil_model]
1                        # 0=ignore this module, 1 = run the module
$time                    # duration of a time step in minutes
1                        # method, 1 = without slow baseflow, 2 = with slow baseflow (not recommended)
$outpath//$sat_def_grid  # (new) saturation deficite-grid (in mm)
$Writegrid               # writegrid for this grid
$outpath//$SUZgrid       # (new) storage grid for unsat. zone
$Writegrid               # writegrid for this grid
$outpath//$SIFgrid       # (new) storage grid for interflow storage
$Writegrid               # writegrid for this grid
$outpath//$SBiagrid      # (new) grid for soil moisture in the inaktive soil storage
$Writegrid               # Writegrid for inaktive soil moisture
$outpath//$fcia_grid     # (new) grid for plant available field capacity in the  inaktiven soil storage
$Writegrid               # writegrid for this grid
$outpath//$SSPgrid       # (new) grid for the relative fraction of the soil storages, which is in contact with ground water
$Writegrid               # writegrid for this grid
$outpath//$QDgrid        # (new) grid for surface runoff
$Writegrid            # writegrid for this grid
$outpath//$QIgrid        # (new) grid for Interflow
$Writegrid            # writegrid for this grid
$outpath//$Peakgrid      # (new) grid for Peakflow (maximum peakflow for the entire model time)
$outpath//qdir//$grid//.//$code//$startyear $hour_mean # statistic of the surfeca discharge
$outpath//qifl//$grid//.//$code//$startyear $hour_mean # statistic of the Interflows
$outpath//qbas//$grid//.//$code//$startyear $hour_mean # statistic of the base flow
$outpath//qbav//$grid//.//$code//$startyear $hour_mean # statistic of the slow base flow
$outpath//qges//$grid//.//$code//$startyear $hour_mean # statistic of the total discharge
$outpath//sb__//$grid//.//$code//$startyear $hour_mean # soil storage in mm per zone
$outpath//suz_//$grid//.//$code//$startyear $hour_mean # drainage storage in mm per zone
$outpath//sifl//$grid//.//$code//$startyear $hour_mean # interflow storage in mm per zone
$outpath//sd__//$grid//.//$code//$startyear $hour_mean # saturation deficite per zone in mm
1	2	3	4	5	6	7	8	9	16	21	27	28	54	55	57 	# range for subbasin codes		
0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	0.01	# Rezessionsparameterter m fuer Saettigungsflaechenmodell in Metern
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     # Korrekturfaktor fuer Transmissivitaeten
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     # Korrekturfaktor fuer K-Wert (vertikale Versickerung), Modell erwartet k in m/s
6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     6.0     # Speicherrueckgangskonstante Direktabflus ELS in h
25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    25.0    # Saettigungsdefizit, bei dessen Unterschreitung lokaler Interflow gebildet wird
80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    80.0    # Speicherrueckgangskonstante Interflow ELS in h
3600    3600    3600    3600    3600    3600    3600    3600    3600    3600    3600    3600    3600    3600    3600    3600    # Rueckgangskonstante verzoegerter Basisabfluss in h
0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    0.03    # maximale Tiefenversickerungsrate bei Saettigung in mm/h
0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    0.01    # Anfangswert QBB
0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     0.0     # Anfangsfuellung des SUZ-Speichers in n*nFK
0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    0.75    # Anfangssaettigungsdefizit in n*nFK, beeinflusst den ersten Basisabfluss
3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     3.0     # Anspringpunkt fuer Makroporenabfluss (in mm/h!, bezogen auf Stundenniederschlag!), alles darueber geht direkt in den Drainspeicher!
0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     0.9     # Reduktionsfaktor fuer Auffuellung von Verdunstungsverlusten aus dem Grundwasser und aus dem Interflowspeicher
0.1     0.1     0.1     0.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     # Anteil an der effektiven Schneeschmelze, der bei geschlossener Schneedecke direkt abfliesst und nicht in den Boden gelangen kann
$readgrids          # 1=read grids from disk, else generate internal
$outpath//storage_topmodel.ftz    # if readgrids = 1, then this file contains the contents of the flow travel time zones for interflow and surface flow and for the tracers

#
# now the routing descriptions follow
# TG (Teilgebiet) is the code for subbasin, OL (Oberlieger) is the code for a upper subbasin (tributary) , flowing into the actual subbasin TG
# AE and AErel (in TG (...)) are the area of the subbasin (TG) and the relative fraction of the modelled
# compared to the real subbasin area (normally 1.0, but if inflows have to be considered or only parts of a subbasin are modelled, AErel can be smaller)
#
# example:
# TG 9  (AE= 484.0, AErel=0.8125) # TG 11 is external inflow ZL 1:
# from OL 10 (kh=0.4, kv=0.4, Bh=3.5,  Bv=20.0, Th=0.5, Mh=15.0, Mv=15.0, I=0.0128, L=10408.3, AE=3.75  )
#  and ZL  1 (modus = extern $outpath//abstract11.dat 4 5 , kh=0.4, kv=0.4, Bh=7.0,  Bv=50.0, Th=1.0, Mh=15.0, Mv=15.0, I=0.0066, L=10838.5, AE=90.75 )
#  and OL 12 (kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=36339.6, AE=81.5 )
# TG 9  (AE= 484.0, AErel=0.6435)
# from OL 10 (kh=0.4, kv=0.4, Bh=3.5,  Bv=20.0, Th=0.5, Mh=15.0, Mv=15.0, I=0.0128, L=10408.3, AE=3.75  ) # TG 9 with external inflow from OL 11 and 12 
#  and ZL  2 (modus = extern $outpath//abstract11.dat 4 5, kh=0.4, kv=0.4, Bh=7.0,  Bv=50.0, Th=1.0, Mh=15.0, Mv=15.0, I=0.0066, L=10838.5, AE=90.75 )
#  and ZL  1 (modus = extern $outpath//abstract12.dat 4 5, kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=36339.6, AE=81.5  )
#
# kh, kv are storage coefficients for mean channel and flood plain flow; measure in h
# Mh, Mv are Manning roughness for mean channel (Mh) and flood plains (Mv); measure in m^(1/3)/s
# Bh, Bv are width of mean channel resp. floodplains; measure in m
# Th is the water depth in the mean channel, measure in m
# I is the slope angle in flow direction in m/m, L is flow length in m, AE is the area of the sub-basin in km^2
#
# for abstractions: modus = intern or extern (if extern, a file name must follow), 
# then 3 numeric parameters: min. threshould, fraction of abstracted water, max. capacity of abstraction channel, 
# then the measure (standard is mm/h, other poss.: m3/s or m^3/s; all other measures are interpreted as mm/h)
# example:
#  and AL 1 (modus = intern 0.0 1.0 1 m^3/s )  -> there must be an corresponding internal inflow ZL 1 ! 
#  and AL 2 (modus = extern $outpath//abstrac2.dat 0 1 15 m^3/s )
#
# for inflows: like normal routings except that the first parameters are modus = intern / extern,
# if modus = extern, a filename and two integer parameters indicating the start row/start column for the given file must follow
# example:
#     from  ZL 1 (modus = intern , kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=10, AE=81.5 )
# or: from  ZL 1 (modus = extern $outpath//abstrac1.dat 4 5, kh=0.4, kv=0.4, Bh=10.0, Bv=60.0, Th=1.4, Mh=27.0, Mv=15.0, I=0.0084, L=36339.6, AE=81.5 )
#
# for reservoirs: a filename for the output of reservoir content/tracer concentrations, an initial filling and initial tracer concentrations must be given
# additionally for each reservoir a paragraph containing the abstraction rule has to be found within the control file
# example: 
# and SP 1 ( file = $outpath//spv___01.//$code//$startyear , V0 = 1E8, C0 = 1.0  0.1  3.0  0.0  0.0  0.0  0.0  0.0  0.0 )  # C0: initial concentration for all of max. 9 tracers in the reservoir water
# ..
# [abstraction_rule_reservoir_1] 
# 2                              # number of following reservoir-content - abstraction points
# 0        20.0                  # first value: reservoir content (storage), second value: abstraction(outflow) at this filling level
# 1e8      20.0                  
#
#

[routing_model]
1                   # 0=ignore this module, 1 = run the module, 2=run the module with observed inflows into the routing channels (from discharge files)
$time               # duration of a time step in minutes
1 1500 120 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)
# 5 1200 90 24 means: the flow velocities are calculated befor the model start for a range of 5 to 1200 l/s/km^2 with a logarithmic splitting into 90 flow range intervals. each timestep is also splitted into 24 time sub-intervals
$outpath//qgko//$grid//.//$code//$startyear $routing_code  # name of the statistic file with routed discharges
$inpath//SpecificDischargeSorted.txt  # name of the file with observed discharges (l/s/km^2)
7 # 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
#13 0   # this is not required, if there is no column in the file for a subbasin, then this is indicated by an absent entry for this subbasin or by setting column index to 0
0 #  timeoffset (for r-square calculation. intervals up to this parameter are not evaluated in r-square calculation. e.g. 12: first 12 intervals are neglected )
TG 8 (AE=283.280, AErel=1.0)
 from OL 35 (kh=0.1, kv=0.4, Bh=10.1, Bv= 40.3, Th= 1.29, Mh=15.0, Mv=10.0, I=0.0165, L=11076.9, AE=40.0)
TG 3 (AE=79.800, AErel=1.0)
 from OL 30 (kh=0.1, kv=0.4, Bh= 0.8, Bv=  3.2, Th= 0.23, Mh=15.0, Mv=10.0, I=0.0819, L=13536.7, AE=1.0)
TG 2 (AE=585.720, AErel=1.0)
 from OL  8 (kh=0.1, kv=0.4, Bh=17.8, Bv= 71.1, Th= 1.78, Mh=15.0, Mv=10.0, I=0.0069, L=32426.6, AE=283.280)
  and OL  3 (kh=0.1, kv=0.4, Bh= 9.3, Bv= 37.2, Th= 1.03, Mh=15.0, Mv=10.0, I=0.0124, L=10479.9, AE=79.800)
TG 11 (AE=267.850, AErel=1.0)
 from OL 38 (kh=0.1, kv=0.4, Bh= 1.7, Bv=  6.8, Th= 0.30, Mh=15.0, Mv=10.0, I=0.0302, L=20000, AE=4.250)
TG 9 (AE=198.320, AErel=1.0)
 from OL 36 (kh=0.1, kv=0.4, Bh= 6.3, Bv= 25.2, Th= 0.92, Mh=15.0, Mv=10.0, I=0.0213, L=20000, AE=60.0)
TG 10 (AE=123.760, AErel=1.0)
 from OL 37 (kh=0.1, kv=0.4, Bh= 3.4, Bv= 13.8, Th= 0.68, Mh=15.0, Mv=10.0, I=0.0388, L=20000, AE=4.250)
TG 13 (AE=679.630, AErel=1.0)
 from OL 11 (kh=0.5, kv=1.4, Bh=19.6, Bv= 78.5, Th= 1.96, Mh=15.0, Mv=10.0, I=0.0036, L=10000, AE=267.850)
  and SUMTRIB 9&10 (kh=0.5, kv=1.4, Bh=14.6, Bv= 58.5, Th= 1.46, Mh=15.0, Mv=10.0, I=0.0095, L=30000, AE=322.08)
  and OL 40 (kh=0.1, kv=0.4, Bh= 1.2, Bv=  4.8, Th= 0.24, Mh=15.0, Mv=10.0, I=0.0405, L=30000, AE=0.5)
TG 12 (AE=234.570, AErel=1.0)
 from OL 39 (kh=0.1, kv=0.4, Bh= 2.1, Bv=  8.5, Th= 0.30, Mh=15.0, Mv=10.0, I=0.0193, L=30000, AE=4.250)
TG 15 (AE=1181.220, AErel=1.0)
 from SUMTRIB 12&13 (kh=0.5, kv=1.4, Bh=29.2, Bv=117.0, Th= 2.92, Mh=15.0, Mv=10.0, I=0.0028, L=30000, AE=914.2)
  and OL 42 (kh=0.1, kv=0.4, Bh= 0.6, Bv=  2.5, Th= 0.23, Mh=15.0, Mv=10.0, I=0.1361, L=20000, AE=2.00)
TG 14 (AE=350.660, AErel=1.0)
 from OL 41 (kh=0.1, kv=0.4, Bh= 1.9, Bv=  7.7, Th= 0.35, Mh=15.0, Mv=10.0, I=0.0336, L=50000, AE=4.0)
TG 4 (AE=1839.420, AErel=1.0)
 from SUMTRIB 14&15 (kh=0.5, kv=1.4, Bh=30.0, Bv=120.1, Th= 3.00, Mh=15.0, Mv=10.0, I=0.0073, L=60000, AE=1531.88)
TG 19 (AE=569.870, AErel=1.0)
 from OL 46 (kh=0.1, kv=0.4, Bh= 4.6, Bv= 18.4, Th= 0.50, Mh=15.0, Mv=10.0, I=0.0119, L=66554.4, AE=4.00)
TG 17 (AE=174.540, AErel=1.0)
 from OL 44 (kh=0.1, kv=0.4, Bh= 2.2, Bv=  9.0, Th= 0.40, Mh=15.0, Mv=10.0, I=0.0318, L=25056.2, AE=10.25)
TG 16 (AE=313.390, AErel=1.0)
 from OL 43 (kh=0.1, kv=0.4, Bh= 3.9, Bv= 15.5, Th= 0.69, Mh=15.0, Mv=10.0, I=0.0321, L=13525.4, AE=47.25)
TG 18 (AE=404.980, AErel=1.0)
 from OL 45 (kh=0.1, kv=0.4, Bh= 2.7, Bv= 10.7, Th= 0.49, Mh=15.0, Mv=10.0, I=0.0333, L=14843.8, AE=0.25)
  and OL 16 (kh=0.1, kv=0.4, Bh=17.0, Bv= 67.9, Th= 1.72, Mh=15.0, Mv=10.0, I=0.0103, L=13246.7, AE=313.390)
TG 20 (AE=784.630, AErel=1.0)
 from SUMTRIB 17&18 (kh=0.1, kv=0.4, Bh=14.3, Bv= 57.3, Th= 1.43, Mh=15.0, Mv=10.0, I=0.0082, L=29547.3, AE=579.52)
  and OL 47 (kh=0.1, kv=0.4, Bh= 2.6, Bv= 10.3, Th= 0.49, Mh=15.0, Mv=10.0, I=0.0364, L=27860.3, AE=0.750)
TG 6 (AE=1458.250, AErel=1.0)
 from SUMTRIB 19&20 (kh=0.1, kv=0.4, Bh=20.1, Bv= 80.4, Th= 2.16, Mh=15.0, Mv=10.0, I=0.0115, L=13333.7, AE=1354.5)
TG 23 (AE=516.470, AErel=1.0)
 from OL 50 (kh=0.1, kv=0.4, Bh= 4.7, Bv= 18.7, Th= 0.76, Mh=15.0, Mv=10.0, I=0.0261, L=41068.6, AE=14.50)
TG 25 (AE=612.140, AErel=1.0)
 from OL 52 (kh=0.1, kv=0.4, Bh= 3.4, Bv= 13.5, Th= 0.49, Mh=15.0, Mv=10.0, I=0.0213, L=61916.5, AE=50.750)
TG 26 (AE=673.030, AErel=1.0)
 from OL 53 (kh=0.1, kv=0.4, Bh= 4.8, Bv= 19.3, Th= 0.58, Mh=15.0, Mv=10.0, I=0.0142, L=56533.7, AE=22.50)
TG 27 (AE=691.550, AErel=1.0)
 from OL 54 (kh=0.1, kv=0.4, Bh= 3.2, Bv= 12.9, Th= 0.41, Mh=15.0, Mv=10.0, I=0.0165, L=59139.6, AE=0.750)
TG 7 (AE=897.150, AErel=1.0)
 from OL 27 (kh=0.1, kv=0.4, Bh=30.5, Bv=121.9, Th= 1.25, Mh=15.0, Mv=10.0, I=0.0023, L=4338.5, AE=691.550)
 and SP 1 ( file = $outpath//resv01.//$code//$startyear , V0 = 7.0E+7, C0 = 1.0  0.1  3.0  0.0  0.0  0.0  0.0  0.0  0.0 , dTmin=60 )  # C0: initial concentration for all of max. 9 tracers in the reservoir water
TG 24 (AE=1370.680, AErel=1.0)
 from OL  7 (kh=0.1, kv=0.4, Bh=90, Bv=300, Th= 1.5, Mh=15.0, Mv=10.0, I=0.0010, L=50000, AE=897.150)
  and SP 2 ( file = $outpath//resv02.//$code//$startyear , V0 = 1.1E8, C0 = 1.0  0.1  3.0  0.0  0.0  0.0  0.0  0.0  0.0 , dTmin=60 )  # C0: initial concentration for all of max. 9 tracers in the reservoir water
TG 55 (AE=3309.730, AErel=1.0)
 from OL 23 (kh=0.1, kv=0.4, Bh=31.9, Bv=127.8, Th= 3.19, Mh=25.0, Mv=10.0, I=0.0010, L=3000.0, AE=516.470)
  and OL 25 (kh=0.1, kv=0.4, Bh=30.4, Bv=121.5, Th= 3.04, Mh=25.0, Mv=10.0, I=0.0018, L=3000.0, AE=612.140)
  and SUMTRIB 24&26 (kh=0.5, kv=2, Bh=150, Bv=300.3, Th= 1.5, Mh=15.0, Mv=10.0, I=0.0010, L=15000, AE=2043.71)
TG 22 (AE=3842.900, AErel=1.0)
 from OL 55 (kh=0.5, kv=2, Bh=200, Bv=500, Th= 2.5, Mh=15.0, Mv=10.0, I=0.0010, L=60000, AE=3309.730)
TG 21 (AE=376.470, AErel=1.0)
 from OL 48 (kh=0.1, kv=0.4, Bh= 5.9, Bv= 23.7, Th= 0.75, Mh=25.0, Mv=10.0, I=0.0161, L=45995.8, AE=33.440)
TG 56 (AE=112.310, AErel=1.0)
 from OL 49 (kh=0.1, kv=0.4, Bh= 2.4, Bv=  9.7, Th= 0.52, Mh=25.0, Mv=10.0, I=0.0458, L=31149.1, AE=11.290)
TG 57 (AE=4457.950, AErel=1.0)
 from OL 21 (kh=0.1, kv=0.4, Bh=28.4, Bv=113.5, Th= 2.84, Mh=25.0, Mv=10.0, I=0.0010, L=100.0, AE=376.470)
  and SUMTRIB 22&56 (kh=0.5, kv=2, Bh=200, Bv=500, Th= 2.6, Mh=15.0, Mv=10.0, I=0.0010, L=12000, AE=3955.21)
TG 5 (AE=6524.440, AErel=1.0)
 from SUMTRIB 6&57 (kh=0.5, kv=2, Bh=250, Bv=600, Th= 2.7, Mh=15.0, Mv=10.0, I=0.0010, L=100000 AE=5916.2)
TG 1 (AE=9715.940, AErel=1.0)
 from OL  2 (kh=0.1, kv=0.4, Bh=28.1, Bv=112.2, Th= 1.0, Mh=15.0, Mv=10.0, I=0.0026, L=27997.0, AE=585.720)
  and SUMTRIB 4&5  (kh=0.5, kv=2, Bh=150, Bv=600, Th= 2.50, Mh=15.0, Mv=10.0, I=0.0010, L=50000, AE=8363.86)
  and OL 28 (kh=0.5, kv=2, Bh= 4.9, Bv= 19.7, Th= 0.61, Mh=15.0, Mv=10.0, I=0.0154, L=69359.2, AE=0.250)



[abstraction_rule_reservoir_1] #AE ca. 900 km^2, = ca 10 m^3/s pro 1 mm/d
14
0	0
5.00E+07	1
1.00E+08	3
1.30E+08	6
1.50E+08	12
1.70E+08	25
1.90E+08	40
2.20E+08	70
2.40E+08	100
2.60E+08	150
2.80E+08	220
3.00E+08	320
3.30E+08	470
3.50E+08	600

[abstraction_rule_reservoir_2] #AE ca. 2800 km^2, = ca 30 m^3/s pro 1 mm/d
14
0	0
5.00E+07	1
1.00E+08	3
1.30E+08	6
1.50E+08	12
1.70E+08	25
1.90E+08	40
2.20E+08	70
2.40E+08	100
2.60E+08	150
2.80E+08	220
3.00E+08	320
3.30E+08	470
3.50E+08	600

# the following section defines combinations of single landuse types to combinations of them.
# e.g. a landuse type deciduous forest may contain of oaks, bushes, and herbs, so each of those three components
# must be parameterised in the traditional landuse table. Example: oaks = code 1, bushes = code 2, herbs = code 3
# here, the combination of oaks, bushes and herbs will be parameterised like: 1 deciduous_forest { layers = 1, 2, 3;}
# The VCF (vegetation covered fraction) of each landuse will define the amount of water and radiation (except diffuse 
# radiation which will go through the canopy layer) reaching the next layer. The uppermost layer must be listed first, 
# the next layer follows then a.s.o.
# All multilayer-landuses must have an equal number of layers. Missing layers can be filled up from the end of the 
# list using landuse code 9999, e.g. grassland would be defined in a 3-layer configuration by "2 grass {layers = 4, 9999, 9999;}
# When the multilayer_landuse table is used, the codes of the LANDUSE-Grid are referring no longer to the landuse_table 
# anymore but to the multilayer_table following. The codes in the old landuse table are reffering to the entries in the 
# multilayer_landuse table

# data taken from BCeast.e00, downloaded from http://geogratis.cgdi.gc.ca/geogratis/en/option/select.do?id=99 AVHRR Land Cover Data BC, Eastern Part
# 1 Mixed forest 
# 2 Deciduous trees 
# 3 Water 
# 	4 Transitional forest 
# 5 Coniferous forest 
# 	6 Tundra 
# 7 Barren land 
# 8 Ice, snow 
# 	9 Agr - cropland 
# 	10 Agr - rangeland 
# 	11 Built-up area 


[multilayer_landuse] 
6 # count of multilayer landuses
1  MixedForest	{ Landuse_Layers = 1 ; 	# references to the landuse table: the first land use code is assumed to be the uppermosts layer, the next one comes as scond a.s.o.
 		          k_extinct = 0.3; 				# extinction coefficient of d'Lambert-Beer's law for reducing radiation in its way through the canopies (after COUPmodel, Jansson and Karlberg)
 		          LAI_scale = 20;   			# Scaling factor for calculating the aerodynamic resistencies of layer 2..n dependent on the cumulated leaf area index (after COUPmodel, Jansson and Karlberg)
		        }
2  DeciduousForest 	{ Landuse_Layers = 2;	k_extinct = 0.3; LAI_scale = 20;}
3  Water			{ Landuse_Layers = 3;	k_extinct = 0.3; LAI_scale = 20;}
5  ConiferousForest	{ Landuse_Layers = 5;	k_extinct = 0.3; LAI_scale = 20;}
7  BarrenLand		{ Landuse_Layers = 7;	k_extinct = 0.3; LAI_scale = 20;}
8  IceSnow		 	{ Landuse_Layers = 8;	k_extinct = 0.3; LAI_scale = 20;}


# declaring some common variables for vegetation period dependent grid-writing
# default (if not used in land use table at all) is JDVegReset = 1 and JDVegWrite = 365
$set $JDVegReset = 3
$set $JDVegWrite = 365

[landuse_table]
6   # number of following land use codes, then aequences for land uses will follow
1 mixed_forest  {  method    	       = VariableDayCount; 	# valid methods: "VariableDayCount" with variable number of fix points, DynamicPhenology_1 which estimates the begin of the diurnal cycle in spring dependent on daily temperature sums only (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
	         RootDistr         = 1.0;		# parameter for root density distribution
	         TReduWet          = 0.95;		# relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
	         LimitReduWet      = 0.5;		# minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
	         HReduDry          = 3.45;		# hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
	         IntercepCap       = 0.35; 		# optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
	         JulDays           = 15    46    74    105   135   166    196    227    258    288    319    349 ;  # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
	         Albedo            = 0.15  0.15  0.15  0.15  0.15  0.15   0.15   0.15   0.15   0.15   0.15   0.15;  # Albedo (snow free)
	         rsc               = 90    90    85    70    60    60     60     60     60     80     90     90;    # leaf surface resistance in s/m
	         rs_interception   = 90    90    85    70    60    60     60     60     60     80     90     90;    # INTERCEPTION surface resistance in s/m
	         rs_evaporation    = 230   230   230   230   230   230	  230	 230	230	   230	  230	 230;    # SOIL surface resistance in s/m (for evaporation only)
	         LAI               = 2     2     2     4     8     10     10     10     8      5      3      3;     # Leaf Area Index (1/1)
	         Z0                = 3.0   3.0   3.0   5.0   8.0   10.0   10.0   10.0   9.0    5.0    3.0    3.0;   # Roughness length in m
	         VCF               = 0.8   0.8   0.8   0.9   0.92  0.92   0.92   0.92   0.9    0.8    0.8    0.8;   # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
	         RootDepth         = 1.3   1.3   1.3   1.3   1.3   1.3    1.3    1.3    1.3    1.3    1.3    1.3;   # Root depth in m
	         AltDep            = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
	       }	
2 decidous_forest { method   	      = VariableDayCount; 	# valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
                   RootDistr        = 1.0;		# parameter for root density distribution
	         TReduWet         = 0.95;		# relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
	         LimitReduWet     = 0.5;		# minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
	         HReduDry         = 3.45;		# hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
	         IntercepCap      = 0.3; 		# optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
	         JulDays          = 15    46    74    105   135   166    196    227    258    288    319    349 ;  # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
	         Albedo           = 0.17  0.17  0.17  0.17  0.17  0.17   0.17   0.17   0.17   0.17   0.17   0.17;  # Albedo (snow free)
	         rsc              = 100   100   95    75    65    65     65     65     65     85     100    100;   # leaf surface resistance in s/m
	         rs_interception  = 100   100   95    75    65    65     65     65     65     85     100    100;   # INTERCEPTION surface resistance in s/m
	         rs_evaporation   = 220   220   220  220    220   220    220    220    220    220    220    220;   # SOIL surface resistance in s/m (for evaporation only)
	         LAI              = 0.5   0.5   0.5   3     8     8      8      8      8      3      0.5    0.5;   # Leaf Area Index (1/1)
	         Z0               = 0.3   0.3   0.3   3.0   8.00  10.0   10.0   10.0   10.0   3.0    0.5    0.3;   # Roughness length in m
	         VCF              = 0.7   0.7   0.7   0.8   0.95  0.95   0.95   0.95   0.9    0.8    0.7    0.7;   # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
	         RootDepth        = 1.4   1.4   1.4   1.4   1.4   1.4    1.4    1.4    1.4    1.4    1.4    1.4;   # Root depth in m
	         AltDep           = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
	       }	
3 water    { method           = VariableDayCount;	# valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
	         RootDistr        = 1.0;		# parameter for root density distribution
	         TReduWet         = 0.95;		# relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
	         LimitReduWet     = 0.5;		# minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
	         HReduDry         = 3.45;		# hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
	         IntercepCap      = 0; 		# optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
	         JulDays          = 15    46    74    105   135   166    196    227    258    288    319    349 ;  # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter. The count of days doesn't matter.
	         Albedo           = 0.05  0.05  0.05  0.05  0.05  0.05   0.05   0.05   0.05   0.05   0.05   0.05;  # Albedo (snow free)
	         rsc              = 20    20    20    20    20    20     20     20     20     20     20     20;    # leaf surface resistance in s/m
	         rs_interception  = 20    20    20    20    20    20     20     20     20     20     20     20;    # INTERCEPTION surface resistance in s/m
	         LAI              = 1.0   1.0   1.0   1.0   1.0   1.0    1.0    1.0    1.0    1.0    1.0    1.0;   # Leaf Area Index (1/1)
	         Z0               = 0.01  0.01  0.01  0.01  0.01  0.01   0.01   0.01   0.01   0.01   0.01   0.01;  # Roughness length in m
	         VCF              = 0.1   0.1   0.1   0.1   0.1   0.1    0.1    0.1    0.1    0.1    0.1    0.1;   # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
	         RootDepth        = 0.01  0.01  0.01  0.01  0.01  0.01   0.01   0.01   0.01   0.01   0.01   0.01;  # Root depth in m
	         AltDep           = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
	       }
5 coniferous_forest    { method           = VariableDayCount;	# valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
	         SoilTillage      = 90 240;     # optional set of 1..n Julian days, depicting days with soil tillage. Important for silting up model
	         RootDistr        = 1.0;		# parameter for root density distribution
	         TReduWet         = 0.95;		# relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
	         LimitReduWet     = 0.5;		# minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
	         HReduDry         = 3.45;		# hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
	         IntercepCap      = 0.4; 		# optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
	         JulDays          = 15    46    74    105   135   166    196    227    258    288    319    349 ;  # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
	         Albedo           = 0.12  0.12  0.12  0.12  0.12  0.12   0.12   0.12   0.12   0.12   0.12   0.12;  # Albedo (snow free)
	         rsc              = 80    80    75    65    55    55     55     55     55     75     80     80;    # leaf surface resistance in s/m
	         rs_interception  = 80    80    75    65    55    55     55     55     55     75     80     80;    # INTERCEPTION surface resistance in s/m
	         rs_evaporation   = 220   220   220  220    220   220    220    220    220    220    220    220;   # SOIL surface resistance in s/m (for evaporation only)
	         LAI              = 8.0   8.0   8.0   9.0   12.0  12.0   12.0   12.0   11.0   10.0   8.0    8.0;   # Leaf Area Index (1/1)
	         Z0               = 10.0  10.0  10.0  10.0  10.0  10.0   10.0   10.0   10.0   10.0   10.0   10.0;  # Roughness length in m
	         VCF              = 0.9   0.9   0.9   0.9   0.9   0.9    0.9    0.9    0.9    0.9    0.9    0.9;   # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
	         RootDepth        = 1.2   1.2   1.2   1.2   1.2   1.2    1.2    1.2    1.2    1.2    1.2    1.2;   # Root depth in m
	         AltDep           = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
	       }	              
7 BarrenLand    {method           = VariableDayCount; 					# valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
	        RootDistr        = 1.0;		# parameter for root density distribution
	        TReduWet         = 0.95;	# relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
	        LimitReduWet     = 0.5;		# minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
	        HReduDry         = 0.5;  	# hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
	        IntercepCap      = 0.1; 	# optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
	        JulDays          = 15    46    74    105   135   166    196    227    258    288    319    349 ;  # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
	        Albedo           = 0.12  0.12  0.12  0.12  0.12  0.12   0.12   0.12   0.12   0.12   0.12   0.12;  # Albedo (snow free)
	        rsc              = 250   250   200   150   150   120    100    90     100    150    200    250;   # leaf surface resistance in s/m
            rs_interception  = 80    80    75    65    55    55     55     55     55     75     80     80;    # INTERCEPTION surface resistance in s/m
	        rs_evaporation   = 220   220   220   220   220   220    220    220    220    220    220    220;   # SOIL surface resistance in s/m (for evaporation only)
	        LAI              = 1.0   1.0   1.0   1.2   1.5   1.7    1.8    1.8    1.4    1.0    1.0    1.0;   # Leaf Area Index (1/1)
	        Z0               = 0.05  0.05  0.05  0.10  0.15  0.15   0.15   0.15   0.15   0.10   0.05   0.05;  # Roughness length in m
	        VCF              = 0.8   0.8   0.8   0.8   0.8   0.8    0.8    0.8    0.8    0.8    0.8    0.8;   # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
	        RootDepth        = 0.1   0.1   0.1   0.15  0.15  0.2    0.2    0.2    0.2    0.15   0.1    0.1;   # Root depth in m
	        AltDep           = 0.025 0.025 0.025 0.025 0.025 0.025 -0.025 -0.025 -0.025 -0.025 -0.025 -0.025; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
	        }	
8 IceSnow    {method           = VariableDayCount; 					# valid methods: "VariableDayCount" with variable number of fix points (other methods will follow) --> old method: if the table is structured like the old ones, they are still valid
	        RootDistr        = 1.0;		# parameter for root density distribution
	        TReduWet         = 0.95;	# relative Theta value for beginning water stress (under wet conditions -> set >= 1 for crop which doesn't depend on an aeral zone
	        LimitReduWet     = 0.5;		# minimum relative reduction factor of real transpiration when water content reaches saturation. The reduction factor value will go down linearly starting at 1.0 when relative Theta equals TReduWet (e.g. 0.95) to LimitReduWet when the soil is saturated (Theta rel = 1.0)
	        HReduDry         = 0.5;  	# hydraulic head (suction) for beginning dryness stress (for water content resulting in higher suctions, ETR will be reduced down to 0 at suction=150m)
	        IntercepCap      = 0.1; 	# optional: specific thickness of the water layer on the leafes in mm. if omitted here, the dedfault parameter from interception_model is used
	        JulDays          = 15    46    74    105   135   166    196    227    258    288    319    349 ;  # Julian days for all following rows. Each parameter must match the number of julian days given here! The count of days doesn't matter.
	        Albedo           = 0.82  0.82  0.82  0.82  0.72  0.62   0.52   0.52   0.52   0.62   0.72   0.82;  # Albedo (snow free)
	        rsc              = 50    50    50    50    50    50     50     50     50     50     50     50;   # leaf surface resistance in s/m
            rs_interception  = 50    50    50    50    50    50     50     50     50     50     50     50;    # INTERCEPTION surface resistance in s/m
	        rs_evaporation   = 50    50    50    50    50    50     50     50     50     50     50     50;   # SOIL surface resistance in s/m (for evaporation only)
	        LAI              = 1.0   1.0   1.0   1.0   1.0   1.0    1.0    1.0    1.0    1.0    1.0    1.0;   # Leaf Area Index (1/1)
	        Z0               = 0.02  0.02  0.02  0.02  0.02  0.02   0.02   0.02   0.02   0.02   0.02   0.02;  # Roughness length in m
	        VCF              = 0.8   0.8   0.8   0.8   0.8   0.8    0.8    0.8    0.8    0.8    0.8    0.8;   # Vegetation covered fraction ("Vegetationsbedeckungsgrad")
	        RootDepth        = 0.1   0.1   0.1   0.1   0.1   0.1    0.1    0.1    0.1    0.1    0.1    0.1;   # Root depth in m
	        AltDep           = 0	 0	   0 	 0	   0	 0	    0	   0	  0	     0	    0 	   0; # Verschiebung des Juldays pro Meter (positiv: wird nach hinten geschoben, negativ: wird nach vorne geschoben -> Limit: Wenn zwei Punkte aufeinandertreffen, dann wird nicht weiter verschoben)
	        }	

#original landuse table up to WaSiM 6.4
 1  water           0.05   20  20  20  20  20  20  20  20  20  20  20  20   110 150 250 280  1   1   1   1     .01 .01 .01 .01   .1  .1  .1  .1  0.01 0.01 0.01 0.01     1.0         3.45
 2  settlements     0.1   100 100 100 100 100 100 100 100 100 100 100 100   110 150 250 280  1   1   1   1     10  10  10  10    .5  .5  .5  .5  0.2  0.2  0.2  0.2      1.0         3.45
 3  pine_forest     0.12   80  80  75  65  55  55  55  55  55  75  80  80   110 150 250 280  8  12  12   8     10  10  10  10    .9  .9  .9  .9  1.2  1.2  1.2  1.2      1.0         3.45
 4  decidous_forest 0.17  100 100  95  75  65  65  65  65  65  85 100 100   110 150 250 280  .5  8   8  .5     .3  10  10  .3    .7  .95 .95 .7  1.4  1.4  1.4  1.4      1.0         3.45
 5  mixed_forest    0.15   90  90  85  70  60  60  60  60  60  80  90  90   110 150 250 280  2  10  10   2     3   10  10  3     .8  .92 .92 .8  1.3  1.3  1.3  1.3      1.0         3.45
 6  agriculture     0.25   80  80  75  75  65  65  65  65  65  75  90  90   110 150 250 255  1   5   3   1     .05 .5  .2  .05   .3  .8  .7  .3  0.05 0.4  0.3  0.05     1.0         3.45
 7  grass           0.25   90  90  75  65  50  55  55  55  60  70  90  90   110 150 250 280  2   4   4   2     .15 .4  .3  .15   .95 .95 .95 .95 0.4  0.4  0.4  0.4      1.0         3.45
 8  bushes          0.20   80  80  70  70  50  50  50  55  55  70  80  80   110 150 250 280  3   5   5   3     1.5 2.5 2.5 1.5   .9  .95 .95 .9  0.5  0.5  0.5  0.5      1.0         3.45
15  rock            0.12  250 250 250 250 250 250 250 250 250 250 250 250   110 150 250 280  1   1   1   1     .05 .05 .05 .05   .8  .8  .8  .8  0.1  0.1  0.1  0.1      1.0         0.5
19  horticulture    0.25  100 100  90  70  60  60  60  60  60  80 100 100   110 150 250 280  .5  5   5  .5     .4  3   3   .4    .75 .75 .75 .75 0.8  0.8  0.8  0.8      1.0         3.45
	        

#[soil_table]
#14                  # number of following entries
#Code  name                       FC(Vol.%) mSB(Vol.%) ksat(m/s)  suction.  parameter     #k   theta-values                                               rel.k-values for the given thetas                              #h    theta-values for the following h-values                                              h-values for the left Thetas [m]                  #+thickns. maxratio      for each layer (downwards) one entry for k_sat [m/s] has to follow																										
#                                                                   [mm]   1=par 2=table                                                                                                                                                                                                                                                                              of layers	 ko_rel/ku_rel																																										
#---- ----------------           ----------  ---------  --------  -------- ------------- ----- ---------------------------------------------------------  ------------------------------------------------------------- -----  ---------------------------------------------------------------------------------   -------------------------------------------------  ---------  ------------- -----------------------------------------------------------------------------------------------------------------------------------------------------------------------
  1   Sand_(S)                        6.21     38.5     8.25E-5     385         2        10   .43 .4286 .3537 .2143 .1071 .0650 .0513 .0470 .0455 .0450  1 .9221 .2543 .0212 4.8E-04 7.1E-6 9.9E-8 1.3E-9 1.1E-11 0.0    14    0.43 .4286 .3537 .2143 .1071 .0650 .0513 .0470 .0455 .0452 .0450 .0450 .0450 .0450   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      8.16E-5 7.86E-5 7.38E-5 6.85E-5 6.36E-5 5.90E-5 5.48E-5 5.09E-5 4.72E-5 4.39E-5 4.07E-5 3.78E-5 3.51E-5 3.26E-5 3.03E-5 2.81E-5 2.61E-5 2.42E-5 2.25E-5 2.09E-5 1.94E-5		  
  2   loamy_sand_(LS)                10.91     37.3     4.05E-5     373         2        10   .43 .4282 .3740 .2736 .1661 .1043 .0767 .0651 .0600 .0570  1 .8651 .2679 .0422 .0022   7.3E-5 2.1E-6 5.6E-8 1.0E-09 0.0    14    0.43 .4282 .3740 .2736 .1661 .1043 .0767 .0651 .0600 .0583 .0573 .0571 .0570 .0570   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      4.01E-5 3.86E-5 3.62E-5 3.37E-5 3.12E-5 2.90E-5 2.69E-5 2.50E-5 2.32E-5 2.15E-5 2.00E-5 1.86E-5 1.72E-5 1.60E-5 1.49E-5 1.38E-5 1.28E-5 1.19E-5 1.10E-5 1.02E-5 9.52E-6
  3   sandy_loam_(SL)                12.28     34.5     1.23E-5     345         2        10   .41 .4088 .3871 .3431 .2659 .1878 .1339 .1026 .0841 .0657  1 .8097 .3595 .1269 .0207   1.7E-3 1.1E-4 6.0E-6 2.4E-07 0.0    14    0.41 .4088 .3871 .3431 .2659 .1878 .1339 .1026 .0841 .0753 .0690 .0668 .0660 .0657   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      1.21E-5 1.17E-5 1.10E-5 1.02E-5 9.47E-6 8.79E-6 8.16E-6 7.57E-6 7.03E-6 6.53E-6 6.06E-6 5.63E-6 5.22E-6 4.85E-6 4.50E-6 4.18E-6 3.88E-6 3.60E-6 3.34E-6 3.10E-6 2.88E-6
  4   silty_loam_(SIL)               22.58     38.3     1.25E-6     383         2        10   .45 .4496 .4458 .4392 .4239 .3936 .3469 .2928 .2373 .1039  1 .6382 .3764 .2441 .1251   4.5E-2 1.1E-2 1.9E-3 2.2E-04 0.0    14    0.45 .4496 .4458 .4392 .4239 .3936 .3469 .2928 .2373 .1966 .1513 .1249 .1106 .1039   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      1.24E-6 1.19E-6 1.12E-6 1.04E-6 9.64E-7 8.95E-7 8.31E-7 7.71E-7 7.16E-7 6.65E-7 6.17E-7 5.73E-7 5.32E-7 4.94E-7 4.58E-7 4.26E-7 3.95E-7 3.67E-7 3.41E-7 3.16E-7 2.93E-7
  5   loam_(L)                        12.9     35.2     2.89E-6     352         2        10   .43 .4293 .4217 .4074 .3754 .3223 .2608 .2071 .1633 .0884  1 .7131 .3875 .2154 .0811   1.8E-2 2.7E-3 3.0E-4 2.4E-05 0.0    14    0.43 .4293 .4217 .4074 .3754 .3223 .2608 .2071 .1633 .1361 .1101 .0972 .0910 .0884   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      2.86E-6 2.75E-6 2.58E-6 2.40E-6 2.23E-6 2.07E-6 1.92E-6 1.78E-6 1.65E-6 1.54E-6 1.43E-6 1.32E-6 1.23E-6 1.14E-6 1.06E-6 9.83E-7 9.13E-7 8.48E-7 7.87E-7 7.31E-7 6.78E-7
  6   sandy_clay_(SC)                19.43     28.0     3.33E-7     280         2        10   .38 .3794 .3758 .3706 .3606 .3437 .3206 .2943 .2656 .1704  1 .3191 .1423 .0797 .0357   1.2E-2 3.1E-3 6.6E-4 1.0E-04 0.0    14    0.38 .3794 .3758 .3706 .3606 .3437 .3206 .2943 .2656 .2422 .2117 .1906 .1772 .1704   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      3.30E-7 3.18E-7 2.98E-7 2.77E-7 2.57E-7 2.39E-7 2.21E-7 2.06E-7 1.91E-7 1.77E-7 1.65E-7 1.53E-7 1.42E-7 1.32E-7 1.22E-7 1.13E-7 1.05E-7 9.78E-8 9.08E-8 8.43E-8 7.83E-8
  7   silty_clay_(SIC)               27.65     29.0     5.56E-8     290         2        10   .36 .3599 .3596 .3591 .3581 .3562 .3526 .3465 .3363 .2665  1 .1439 .0803 .0569 .0369   2.1E-2 1.0E-2 4.2E-3 1.2E-03 0.0    14    0.36 .3599 .3596 .3591 .3581 .3562 .3526 .3465 .3363 .3245 .3042 .2865 .2737 .2665   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      5.49E-8 5.29E-8 4.97E-8 4.61E-8 4.28E-8 3.98E-8 3.69E-8 3.43E-8 3.18E-8 2.95E-8 2.74E-8 2.55E-8 2.36E-8 2.19E-8 2.04E-8 1.89E-8 1.76E-8 1.63E-8 1.51E-8 1.40E-8 1.30E-8
  8   clay_(C)                       29.12     31.2     5.56E-7     312         2        10   .38 .3799 .3792 .3784 .3767 .3735 .3679 .3592 .3461 .2707  1 .1244 .0641 .0429 .0258   1.3E-2 5.8E-3 2.0E-3 5.1E-04 0.0    14    0.38 .3799 .3792 .3784 .3767 .3735 .3679 .3592 .3461 .3324 .3101 .2915 .2782 .2707   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      5.49E-7 5.29E-7 4.97E-7 4.61E-7 4.28E-7 3.98E-7 3.69E-7 3.43E-7 3.18E-7 2.95E-7 2.74E-7 2.55E-7 2.36E-7 2.19E-7 2.04E-7 1.89E-7 1.76E-7 1.63E-7 1.51E-7 1.40E-7 1.30E-7
  9   moor_(M)                       47.31     75.0     8.25E-5     750         2        10   .80 .7995 .7965 .7920 .7821 .7617 .7248 .6694 .5937 .2868  1 .4269 .2497 .1733 .1049   5.2E-2 2.0E-2 5.8E-3 1.1E-03 0.0    14    0.80 .7995 .7965 .7920 .7821 .7617 .7248 .6694 .5937 .5231 .4248 .3545 .3099 .2868   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      8.16E-5 7.86E-5 7.38E-5 6.85E-5 6.36E-5 5.90E-5 5.48E-5 5.09E-5 4.72E-5 4.39E-5 4.07E-5 3.78E-5 3.51E-5 3.26E-5 3.03E-5 2.81E-5 2.61E-5 2.42E-5 2.25E-5 2.09E-5 1.94E-5
 10   settlements_rock_(R)           14.00     15.0     5.56E-8     150         2        10   .20 .1999 .1996 .1992 .1984 .1969 .1942 .1900 .1837 .1474  1 .1244 .0641 .0429 .0258   1.3E-2 5.8E-3 2.0E-3 5.1E-04 0.0    14    0.20 .1999 .1996 .1992 .1984 .1969 .1942 .1900 .1837 .1771 .1664 .1575 .1510 .1474   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      5.49E-8 5.29E-8 4.97E-8 4.61E-8 4.28E-8 3.98E-8 3.69E-8 3.43E-8 3.18E-8 2.95E-8 2.74E-8 2.55E-8 2.36E-8 2.19E-8 2.04E-8 1.89E-8 1.76E-8 1.63E-8 1.51E-8 1.40E-8 1.30E-8
 11   clay_loam_(CL)                 21.24     31.5     7.22E-7     315         2        10   .41 .4096 .4067 .4021 .3920 .3729 .3434 .3074 .2675 .1496  1 .5005 .2721 .1720 .0882   3.4E-2 9.3E-3 1.9E-3 2.7E-04 0.0    14    0.41 .4096 .4067 .4021 .3920 .3729 .3434 .3074 .2675 .2357 .1968 .1717 .1569 .1496   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      7.14E-7 6.88E-7 6.46E-7 6.00E-7 5.57E-7 5.17E-7 4.80E-7 4.46E-7 4.14E-7 3.84E-7 3.56E-7 3.31E-7 3.07E-7 2.85E-7 2.65E-7 2.46E-7 2.28E-7 2.12E-7 1.97E-7 1.83E-7 1.70E-7
 12   silt_(SI)                      28.17     42.6     6.94E-7     426         2        10   .46 .4596 .4565 .4511 .4386 .4129 .3702 .3157 .2549 .0901  1 .6139 .3711 .2503 .1384   5.7E-2 1.6E-2 3.3E-3 4.3E-04 0.0    14    0.46 .4596 .4565 .4511 .4386 .4129 .3702 .3157 .2549 .2075 .1519 .1181 .0991 .0901   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      6.87E-7 6.62E-7 6.21E-7 5.77E-7 5.35E-7 4.97E-7 4.61E-7 4.28E-7 3.98E-7 3.69E-7 3.43E-7 3.18E-7 2.95E-7 2.74E-7 2.55E-7 2.36E-7 2.19E-7 2.04E-7 1.89E-7 1.76E-7 1.63E-7
 13   silty_clay_loam_(SICL)         28.16     34.1     1.94E-7     341         2        10   .43 .4298 .4284 .4264 .4218 .4126 .3958 .3706 .3362 .1967  1 .4269 .2497 .1733 .1049   5.2E-2 2.0E-2 5.8E-3 1.1E-03 0.0    14    0.43 .4298 .4284 .4264 .4218 .4126 .3958 .3706 .3362 .3041 .2594 .2275 .2072 .1967   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      1.92E-7 1.85E-7 1.74E-7 1.61E-7 1.50E-7 1.39E-7 1.29E-7 1.20E-7 1.11E-7 1.03E-7 9.60E-8 8.91E-8 8.27E-8 7.68E-8 7.13E-8 6.62E-8 6.15E-8 5.71E-8 5.30E-8 4.92E-8 4.57E-8
 14   sandy_clay_loam_(SCL)          13.35     29.0     3.64E-6     290         2        10   .39 .3886 .3760 .3566 .3221 .2772 .2335 .1976 .1680 .1112  1 .5525 .2157 .0922 .0256   4.6E-3 6.4E-4 7.6E-5 6.7E-06 0.0    14    0.39 .3886 .3760 .3566 .3221 .2772 .2335 .1976 .1680 .1489 .1294 .1189 .1136 .1112   0 .01 .05 .1 .2 .4 .8 1.6 3.45 6.9 20 50 100 150   21 .3333       200      3.60E-6 3.47E-6 3.25E-6 3.02E-6 2.81E-6 2.60E-6 2.42E-6 2.24E-6 2.08E-6 1.93E-6 1.80E-6 1.67E-6 1.55E-6 1.44E-6 1.33E-6 1.24E-6 1.15E-6 1.07E-6 9.91E-7 9.20E-7 8.54E-7

[soil_table]
12                  # number of following entries
#co- name of the 
#de  soil profile
#-- ---------------
030303 SL_SL_SL {method = MultipleHorizons; 		
				 FCap = 12.28; mSB = 34.5; ksat_topmodel = 1.23E-5; suction = 345; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.4 0.4 0.2 0.0 0.0 0.0 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
				 PMacroThresh   = 20  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 3   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.5 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.0 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3       ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = SL01m   SL02m   SL10m   ;  # short descriptions
        		 ksat      	    = 1.23e-5 1.03e-5 1.03e-5 ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0     ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.40    0.40    0.40    ;	 # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.065   0.065   0.065   ;	 # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 7.50    7.50    7.50    ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.89    1.89    1.89    ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	  ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333  ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        } 
040303 L_SL_SL {method = MultipleHorizons; 		
				 FCap = 10.91;  mSB = 37.3; ksat_topmodel = 4.05E-5; suction = 373; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.55 0.15 0.2 0.05 0.03 0.02 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 5   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.9 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3       ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = L1m    SL2m    SL10m   ;  # short descriptions
        		 ksat      	    = 2.89e-6 1.03e-5 1.03e-5 ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0     ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.40    0.40    0.40    ;	 # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.078   0.065   0.065   ;	 # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 3.6     7.50    7.50    ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.56    1.89    1.89    ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	  ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333  ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
040403 L_L_SL {method = MultipleHorizons; 		
				 FCap = 10.91;  mSB = 37.3; ksat_topmodel = 4.05E-5; suction = 373; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.55 0.15 0.2 0.05 0.03 0.02 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 5   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.9 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3       ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = L01m    L02m    SL10m   ;  # short descriptions
        		 ksat      	    = 2.89e-6 2.89e-6 1.03e-5 ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0     ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.40    0.40    0.40    ;	 # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.078   0.078   0.065   ;	 # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 3.6     3.60    7.50    ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.56    1.56    1.89    ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	  ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333  ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
040404 L_L_L {method = MultipleHorizons; 		
				 FCap = 10.91;  mSB = 37.3; ksat_topmodel = 4.05E-5; suction = 373; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.55 0.15 0.2 0.05 0.03 0.02 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 5   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.9 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3       ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = L01m    L02m    L10m   ;  # short descriptions
        		 ksat      	    = 2.89e-6 2.89e-5 2.89e-5 ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0     ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.40    0.40    0.40    ;	 # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.078   0.078   0.078   ;	 # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 3.6     3.60    3.60    ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.56    1.56    1.56    ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	  ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333  ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
040406 L_L_SCL {method = MultipleHorizons; 		
				 FCap = 10.91;  mSB = 37.3; ksat_topmodel = 4.05E-5; suction = 373; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.55 0.15 0.2 0.05 0.03 0.02 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 5   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.9 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3       ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = L01m    L02m    SCL10m  ;  # short descriptions
        		 ksat      	    = 2.89e-6 2.89e-5 1.5e-7 ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0     ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.40    0.40    0.43    ;	 # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.078   0.078   0.089   ;	 # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 3.6     3.60    1.00    ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.56    1.56    1.23    ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	  ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333  ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
				 }
040606_L_SCL_SCL {method = MultipleHorizons; 		
				 FCap = 10.91;  mSB = 37.3; ksat_topmodel = 4.05E-5; suction = 373; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.55 0.15 0.2 0.05 0.03 0.02 0.0; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 5   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.9 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3       ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = L01m    SCL02m    SCL10m  ;  # short descriptions
        		 ksat      	    = 2.89e-6 1.5e-7 1.5e-7 ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0     ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.40    0.43    0.43    ;	 # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.078   0.089   0.089   ;	 # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 3.6     1.00    1.00    ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.56    1.23    1.23    ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	  ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333  ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
060303 SCL_SL_SL {method = MultipleHorizons; 		
				 FCap = 12.28; mSB = 34.5; ksat_topmodel = 1.23E-5; suction = 345; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.39 0.20 0.32 0.06 0.02 0.01 0.00; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 4   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 1.0 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3        ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = SCL01m  SL02m   SL10m   ;  # short descriptions
        		 ksat      	    = 1.5e-7  1.23e-5 1.03e-6  ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0      ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.43    0.40    0.40     ;  # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.089   0.065   0.065    ;  # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 1.00    7.50    7.50     ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.23    1.89    1.89     ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	   ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333   ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59       ;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
060603 SCL_SCL_SL {method = MultipleHorizons; 		
				 FCap = 12.28; mSB = 34.5; ksat_topmodel = 1.23E-5; suction = 345; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.39 0.20 0.32 0.06 0.02 0.01 0.00; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 4   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 1.0 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3        ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = SCL01m  SCL02m  SL10m   ;  # short descriptions
        		 ksat      	    = 1.5e-7  1.5e-7  1.03e-6  ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0      ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.43    0.43    0.40     ;  # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.089   0.089   0.065    ;  # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 1.00    1.00    7.50     ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.23    1.23    1.89     ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	   ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333   ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59       ;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
060606 SCL_SCL_SCL {method = MultipleHorizons; 		
				 FCap = 12.28; mSB = 34.5; ksat_topmodel = 1.23E-5; suction = 345; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.39 0.20 0.32 0.06 0.02 0.01 0.00; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
				 PMacroThresh   = 15  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 4   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 1.0 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.5 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3        ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = SCL01m  SCL02m  SCL10m    ;  # short descriptions
        		 ksat      	    = 1.5e-7  1.5e-7  1.23e-7  ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0      ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.43    0.43    0.43     ;  # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.089   0.089   0.089    ;  # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 1.00    1.00    1.00     ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.23    1.23    1.23     ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	   ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333   ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59       ;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
070606 CL_SCL_SCL {method = MultipleHorizons; 		
				 FCap = 21.24; mSB = 31.5; ksat_topmodel = 7.22E-7; suction = 315; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.1 0.5 0.3 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
				 PMacroThresh   = 12  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 3   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.5 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.2 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3        ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = CL01m   SCL02m  SCL10m    ;  # short descriptions
        		 ksat      	    = 7.22e-7 1.5e-7  1.23e-7  ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0      ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.41    0.43    0.43     ;  # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.095   0.089   0.089    ;  # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 1.9     1.00    1.00     ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.31    1.23    1.23     ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	   ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333   ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59       ;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
070706 CL_CL_SCL {method = MultipleHorizons; 		
				 FCap = 21.24; mSB = 31.5; ksat_topmodel = 7.22E-7; suction = 315; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.1 0.5 0.3 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
				 PMacroThresh   = 12  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 3   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.5 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.2 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3        ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = CL01m   CL02m   SCL10m    ;  # short descriptions
        		 ksat      	    = 7.22e-7 7.22e-7 1.23e-7  ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0      ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.41    0.41    0.43     ;  # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.095   0.095   0.089    ;  # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 1.9     1.90    1.00     ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.31    1.31    1.23     ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	   ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333   ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59       ;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }
070707 CL_CL_CL {method = MultipleHorizons; 		
				 FCap = 21.24; mSB = 31.5; ksat_topmodel = 7.22E-7; suction = 315; # optional parameters which are needed for Topmodel only
				 GrainSizeDist  = 0.1 0.5 0.3 0.05 0.03 0.01 0.01; # optional: when using silting up model, the grain size fractions for sand, silt, clay, and Stones1..4  must be given here. Stones1 = 2-6.3mm, Stones2=6.3-20mm, Stones3=20-63mm, Stones4=63-200mm.; # optional: when using silting up model, the grain size fractions for sand, silt and clay must be given here
				 PMacroThresh   = 12  ;  # precipitation capacity thresholding macropore runoff in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 MacroCapacity  = 3   ;  # capacity of the macropores in mm per hour (not in m/s, because it's more convenient than to write it down in m/s, e.g. 5mm/h = 1.38e-6)
				 CapacityRedu   = 0.5 ;  # reduction of the macropore capacity with depth -> pores become less dense. This Factor describes the reduction ratio per meter 
				 MacroDepth     = 1.2 ;  # maximum depth of the macropores
		         horizon   	    = 1       2       3        ;  # ID of the horizon (must be ascendent) it's recommended to name the horizons shortly in the following row
				 Name		    = CL01m   CL02m   CL10m    ;  # short descriptions
        		 ksat      	    = 7.22e-7 7.22e-7 6.00e-7  ;	 # saturated hydraulic conductivity	in m/s
	       		 k_recession    = 0.8     0.9     1.0      ;  # k sat recession with depth (could also be controlled by different layers if no k decrease is wanted (set this parameter to 1.0
		         theta_sat 	    = 0.41    0.41    0.41     ;  # saturated water content (fillable porosity in 1/1)
		         theta_res 	    = 0.095   0.095   0.095    ;  # residual water content (in 1/1, water content which cannot be poured by transpiration, only by evaporation)
		         alpha     	    = 1.9     1.90    1.90     ;	 # van Genuchten Parameter Alpha
		         Par_n     	    = 1.31    1.31    1.31     ;	 # van Genuchten Parameter n
		         Par_tau		= 0.5	  0.5	  0.5	   ;  # sog. Mualem-Parameter tau in der van-Genuchten-Gleichung (dort normalerweise 0.5)
		         thickness 	    = 0.1     0.1     0.3333   ;	 # thickness of each single numerical layer in this horizon in m
		         layers    	    = 1       2       59       ;	 # numerical number of layers in this horizon. The thickness of the layer is given by layers x thickness. All profiles must have an identical number of layers (for memory handling reasons only)
		        }


# allowed keywords for substance transport (without ""-chars):
# "radioactive" resp. "non_radioactive"
# "evaporating" resp. "non_evaporating"
# "half_time" with its unit "d"
# "min_conc" and "max_conc"
# measures: "mg/l", "g/l", "kg/kg", "Kg/Kg"; all other units will be interpreted as kg/kg (relative concentration)

[substance_transport]
0    # number of tracers to be considered (max. 9)
#
# name  radioact. or not half time in days   evapor. or not  minim. concentr.  max.conc. with unit    initial        initial        output code           writecode             output path      output extension       output extension
#3chars                                                      if no: -9999      mg/l g/l kg/kg         conc. in soil  conc. in gr.w  statfiles             for grids             with closing "\" for stat-files         for grid files
#------ ---------------- ------------------- --------------- ----------------- ---------------------- -------------  -------------  --------------------- --------------------- ---------------- ---------------------- ----------------
  3H        radioactive  half_time =  4500 d     evaporating min_conc = 0     max_conc = 3500  kg/kg  soilini = 1.0   gwini = 2.0   statcode = $hour_mean gridcode = $Writegrid path = $outpath  statext = $code//$startyear gridext = $suffix
  18O   non_radioactive  half_time = -9999 d     evaporating min_conc = -9999 max_conc = -9999 kg/kg  soilini = 1.0   gwini = 1.0   statcode = $hour_mean gridcode = $Writegrid path = $outpath  statext = $code//$startyear gridext = $suffix
  NACL  non_radioactive  half_time = -9999 d non_evaporating min_conc = 0     max_conc = 0.35  kg/kg  soilini = 0.01  gwini = 0.01  statcode = $hour_mean gridcode = $Writegrid path = $outpath  statext = $code//$startyear gridext = $suffix

# irrigation descriptions
# method 1: count  MM1 DD1 amount1 MM2 DD2 amount2 MM3 DD3 amount3 MM4 DD4 amount4 MM5 DD5 amount5 MM6 DD6 amount6 MM7 DD7 amount7 MM8 DD8 amount8 MM9 DD9 amount9 MM10 DD10 amount10  
# method 2a: "starting from MM DD with XX mm to MM DD with YY mm every ZZ days" here, the start end end date are explicitly given 
# method 2b: "starting from MM DD with XX mm YY times every ZZ days" Here, the number of irrigation events is given explicitly */
# method 3: by demand: without additional parameters
[irrigation_table]
0                  # 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
# e.g.
#  6   agriculture     1        1           3  0.3       5     5  15    80     6   1     80    6  15    80     7   1    80     7  15    80     8   1    80     8  15    80     9   1    80     9  15    80     10   15    80
#  6   agriculture     2        1           3  0.3       starting from 05 01 with 20 mm to 07 31 with 20 mm every 15 days
#  6   agriculture     2        1           3  0.3       starting from 05 01 with 20 mm 7 times every 14 days
#  6   agriculture     3        1           3  0.3       
#
#---- ------------ ---------- --------  --------------  ------ --------------- --------------- --------------- --------------- --------------- --------------- --------------- --------------- --------------- ------------------  
 19   horticulture     3        2           0.3 0.1
  1   water            3        2           0.3 0.1
  2   settlements      3        2           0.3 0.1
  3   pine_forest      3        2           0.3 0.1
  4   decidous_forest  3        2           0.3 0.1
  5   mixed_forest     3        2           0.3 0.1
  6   agriculture      3        2           0.3 0.1
  7   grass            3        2           0.3 0.1
  8   bushes           3        2           0.3 0.1
 15   rock             3        2           0.3 0.1
 
                
# the actual configuration of special output is only an example for a model configuration using layered snow model (which is switvhed off for this test data set)
# if you would like to use special output, you need to adopt the fields ID, Xcoords, Ycoords and Level (use as many  levels your stacks have, some stacks have useful information in level 0)
# other entities names are described in the manual at the very end.
[special_output]
1 # 0=do not run this module, 1=run this module
$time # duration of a time step in minutes --> only for compatibility here.
numfiles = 10;
outputfile { header = LayeredSnowData_SWES;
	filename = $outpath//special_output_SWES.//$code//$startyear;
	entity {
		ID 	= SnowWaterEquivalentSolid;
		Symbol 	= swes;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = LayeredSnowData_SWEL;
	filename = $outpath//special_output_SWEL.//$code//$startyear;
	entity {
		ID 	= SnowWaterEquivalentLiquid;
		Symbol 	= swel;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = SnowDensity;
	filename = $outpath//special_output_densSnow.//$code//$startyear;
	entity {
		ID 	= SnowDensity;
		Symbol 	= dens;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = SnowLayerThickness;
	filename = $outpath//special_output_thickSnow.//$code//$startyear;
	entity {
		ID 	= SnowLayerThickness;
		Symbol 	= thick;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = WaterFluxesInSnowLayers;
	filename = $outpath//special_output_qu_down_Snow.//$code//$startyear;
	entity {
		ID 	= q_down_Snow;
		Symbol 	= qu_s;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = Temperature_Stack_Snow;
	filename = $outpath//special_output_Temperature_Stack_Snow.//$code//$startyear;
	entity {
		ID 	= Temperature_Stack_Snow;
		Symbol 	= Tsnow;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = theta_Snow;
	filename = $outpath//special_output_theta_Snow.//$code//$startyear;
	entity {
		ID 	= theta_Snow;
		Symbol 	= thetaS;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = temperature_Soil;
	filename = $outpath//special_output_Temperature_Stack_Soil.//$code//$startyear;
	entity {
		ID 	= Soil_Temperature_Stack;
		Symbol 	= Tsoil;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = theta_Soil;
	filename = $outpath//special_output_theta_Soil.//$code//$startyear;
	entity {
		ID 	= theta;
		Symbol 	= theta;
		Xcoords	= -93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ,-93900 ;
		Ycoords	= 1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250,1252250;
		Level	= 0	 ,1	 ,2	 ,3	 ,4	 ,5	 ,6	 ,7	 ,8	 ,9	 ,10	 ,11	 ,12	 ,13	 ,14	 ,15	 ,16	 ,17	 ,18	 ,19	 ,20	 ,21	 ,22	 ,23	 ,24	 ,25	 ,26	 ,27	 ,28	 ,29	 ,30	;
	}
}
outputfile { header = other_snow_data;
	filename = $outpath//special_output_otherSnowData.//$code//$startyear;
	entity {
		ID 	= snowstorage_solid;
		Symbol 	= ssno;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
	}
	entity {
		ID 	= snowstorage_liquid;
		Symbol 	= sliq;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
	}
	entity {
		ID 	= snow_storage_total;
		Symbol 	= ssto;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
	}
	entity {
		ID 	= snow_storage_total;
		Symbol 	= ssto;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
	}
	entity {
		ID 	= SnowWaterEquivalentSolid;
		Symbol 	= swes0;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
		Level	= 0	 ,0	 ;
	}
	entity {
		ID 	= SnowWaterEquivalentLiquid;
		Symbol 	= swel0;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
		Level	= 0	 ,0	 ;
	}
	entity {
		ID 	= directflow;
		Symbol 	= Qdir;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
	}
	entity {
		ID 	= infiltration_excess;
		Symbol 	= INFEX;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
	}
	entity {
		ID 	= interflow;
		Symbol 	= Qifl;
		Xcoords	= -93900 ,-93900 ;
		Ycoords	= 1252250,1252250;
	}
}
