%Example input.txt: NOTE, YOU NEED TO ADJUST THE INPUT AND OUTPUT FOLDER (2. AND 3. ROW BELOW), AND MAKE AN OUTPUT FOLDER, THEN THE MODEL SHOULD WORK !!!
%***********************************************************
1    %output to screen every X day.  daysscreenoutput.
/home/fcovi/meltmodel_MODTEST/SiteJ_io/input/     %Path for Inputfiles.   inpath
/home/fcovi/meltmodel_MODTEST/SiteJ_io/output/    %Path for Outputfiles.  outpath
125  2017   %first julian day to be calculated.  jdbeg  yearbeg
9    2018   %last julian day to be calculated.   jdend  yearend
0    %discharge to be calculated: 1=yes,0=no,2=yes, but no discharge data.  disyes
2    %1=whole grid computed,  2=only grid cell of weather station.  calcgridyes
%********* 1.) MODEL OUTPUT PARAMETERS: DEFINE WHAT IS WRITTEN TO FILE ***************
1    %number of stakes for melt output.  maxmeltstakes
1    %cum massbal multiplied by this factor in melting.dat.  plusminus
3    %0=no output 1=every step, 2=daily, 3=whole period 4=2and3.  do_out
%shayes exkyes solyes diryes dir2yes difyes gloyes albyes  swbyes linyes loutyes
   0     0      0      0      0       0      0       0       0      0      0
%netyes senyes latyes rainyes enbyes melyes  ablyes surftempyes  posyes ddfyes
0       0      0      1       1      1       1      1            1      1
1    %surface conditions to grid file (2=for specified JDs).  surfyes
1    %snow cover to grid file at midnight.  snowyes
1    %snow or surface written to file if jd dividable by this value.  daysnow
0    %number of jd for output of surface type/snow cover.  numbersnowdaysout
190 210
%******************************** 2.) MASS BALANCE OUTPUT ********************************
0    %gridout winter mass balance yes=1, no=0.   winterbalyes
244  135    %julian day winter starts and ends.  winterjdbeg  winterjdend
1    %gridout summer mass balance yes=1, no=0.  summerbalyes
136  243    %julian day summer starts and ends.  summerjdbeg  summerjdend
0    %1=dates for MB meas read from file, 0=fixed dates.  datesfromfileyes
dummydates.dat    %file containing the dates of massbal meas.  namedatesmassbal
20    % vertical extent of elevation bands (m) for mass balance profiles.  beltwidth
1    %set snow cover to 0 at start of each massbal year.  snow2zeroeachyearyes
1   %times series file with number of pixels snowfree written to file.  snowfreeyes 
-----------------------------------------
0   %gridoutput of melt cumulated=1 or mean=0.  cumulmeltyes
1000    %if cumulated, output in cm=10 or m=1000.  cm_or_m
1    %time series of spatial mean to output (yes=1 no=0).  do_out_area
1    %number of individual grid points for which model result output.  outgridnumber
%***====read if number > 0============================***
%***Outputfilename ** row/x-coord ***column/y-coord *** glob and net data included from input data
pointoutput.txt        50                50       0
%******** 3.) METHODS ENERGY BALANCE COMPONENTS ********************************
2    %1=surface maps 2=start with initial snow cover.  methodinisnow
1    %1=constant for surface types, >2=alb generated.  methodsnowalbedo, 6=zongo
1    %1=direct and diffuse not separated, 2=separated.  methodglobal
2    %1=longin from net,glob,ref (Tsurf=0),2=meas,3-6=from paramet. methodlonginstation
1    %1=longin constant, 2=spatially variable. methodlongin
1    %1=surftemp=0, 2=iteration 3=measurement+(decrease height), 4=snowmodel.  methodsurftempglac
%********************* TURBULENCE OPTION
1    %1=turbulence accord. to Escher-Vetter, 2=Ambach 3=stabil.  methodturbul
1    %1=z0T/z0w fixed ratio 2=according to Andreas (1987).  method_z0Te
1    %1=no ice heat flux 2=ice heat flux.  methodiceheat
1    %1=neglected 2=neg energy balance stored to retard melt.  methodnegbal
%********* SCALING ************************************************
0    %V-A scaling yes=1, no=0.  scalingyes
1.375    %gamma in V-A scaling.  gamma
0.6    %coefficient in V-A scaling.  c_coefficient
%********* 4.) NAMES OF INPUT FILES *********************************************
dav_dem.bin    %name of Digital Terrain Model.  namedgm
dav_watershed.bin    %name of DTM with drainage basin.  namedgmdrain
dav_glacier.bin    %name of DTM glacier.  namedgmglac
dav_slope.bin    %name of DTM slope.  namedgmslope
dav_aspect.bin    %name of DTM aspect.  namedgmaspec
dummy.bin    %name of DTM sky view factor.  namedgmskyview
dav_firn.bin    %name of DTM firnarea.  namedgmfirn
dav_snow0.bin   %name of DTM initial snow cover.  nameinitialsnow 
dav_climate_data.txt    %name of climate data file.  nameklima
%********* 5.) GRID INFORMATION**************************************************
-134.50479    %geographical longitude [degree].  laenge
58.51169    %latitude.  breite
-120.0    %longitude time refers to.  reflongitude
100    %row in DTM where climate station is located.   rowclim
100    %column of climate station.  colclim
1  7    %take this elevation for AWS yes/no.  climoutsideyes  heightclim
100    %gridsize in m.  gridsize
1    %time step in hours.  timestep
%********* 6.) CLIMATE DATA *****************************************************
1    %1=midnight time is 0, 2=time is 24, 3=24 but previous day.  formatclimdata
14    %number of columns in climate file.  maxcol
4    %columns in climate input file: temperature.  coltemp
5    %column containing relative humidity.  colhum
6    %column wind speed [m/s].  colwind
8    %global radiation.  colglob
9    %reflected shortwave radiation.  colref
10    %net radiation.  colnet
11    %longwave incoming radiation.  collongin
12    %longwave outgoing radiation.  collongout
13    %precipitation.  colprec
30    %cloud cover (number of eigths).  colcloud
30    %time variant lapse rate (neg=decrease).  coltempgradvarying
14    %column of discharge data.  coldis
%********** (BIAS) CORRECTIONS TO INPUT DATA (TEMP, WIND) ********************************
0    %add this to air temp (e.g., to get to elevationstation).  ERAtempshift
0    %add this to wind speed (e.g., to get to elevationstation).  ERAwindshift
%********** 7.) LAPSE RATE / SCENARIOS ********************************
1    %1=const lapse rate 2=variable 2AWS 3=grid.  methodtempinterpol
-0.6    %temperature change with elevation [degree/100m].  tempgrad
0    %climate perturbation: temp + this amount.  tempscenario
0    %climate perturbation: precip + this amount in percent.  precscenario
%on/off Jan   Feb   Mar   Apr   May   Jun   Jul   Aug   Sep   Oct   Nov   Dec
0      0.0    0.0  0.0   0.0 -0.432 -0.308 -0.388 -0.322 -0.362 -0.362 0.0 0.0    %monthtempgrad(yes)
0      0.5    1.0  1.0   1.0  1.0     1.0   1.0   1.0   1.0   1.0  1.0   1.0    %monthtempscen(yes)
0     10.0   10.0  10.0  10.0 10.0   10.0   10.0  10.0  10.0  10.0 10.0 10.0    %monthprecipscen(yes)
%******** 8.) SURFACE TYPE / ALBEDO ***********************************
0    %number of transient surface type files
0.875    %albedo for snow and firn (fixed value).  albsnow
0.55    %albedo for slush.  albslush
0.3    %albedo for ice.  albice
0.6    %albedo for firn (method 2).  albfirn
0.1    %albedo for rock outside glacier.  albrock
0.5    %minimum albedo for snow if generated.  albmin
8.0    %increase snowalb/100m elevation for 1. time. snowalbincrease
0.0    %decrease of ice albedo with elevation.  albiceproz
0    %number of days since snow fall at start.  ndstart
%******** 9.) RADIATION *******************************************
6    %number of shade calculation per time step.  split
15.0    %percent diffuse radiation of global radiation.  prozdiffuse
0.75    %transmissivity.  trans
1.0    %first ratio of global radiation and direct rad.  ratio
0.95    %first ratio of direct and clear-sky direct rad.  ratiodir2dir
-0.5    %decrease in surftemp with height if longout meas.  surftemplapserate
0    %direct radiation read from file = 1.  directfromfile
/Users/luser/Documents/example_data/direct/    %Path for direct files.  pathdirectfile
0    %files only exist every number of days defined here.  daysdirect
1.0    %0=slope at climate station is set to 0.  slopestation
%******** 10.) TURBULENCE ***********************************************
0.25    %step for surface temp lowering for iteration.  iterstep
0.001    %roughness length for wind for ice in m   .0027.  z0wice
100.0    %z0Temp is zow divided by this value.  dividerz0T
1.0    %z0snow is z0wice divided by this value.  dividerz0snow
0.0    %increase of z0 with decreasing elevation.  z0proz
0.00000001    %min z0w ice.  icez0min
0.02    %max z0w ice.  icez0max
%********** 11.) PRECIPITATION *******************************************
1    %1=lapse rate 2=scale AWS precip with index map 3=read grids.  methodprecipinterpol
10.0    %precipitation change with elevation [%/100m].  precgrad
0.0 1900    %precipitation change with elevation beyond certain elevation.  precgradhigh precgradelev
50.0    %35 precipitation correction, caused by losses.  preccorr
1.0    %snow precip is multiplied by this factor.  snowmultiplierglacier
1.0    %snow precip is multiplied by this factor.  snowmultiplierrock
1.5    %threshold temperature rain/snow precipitation.  threshtemp
0    %0=if massbal calculations for whole drainage basin, 1=only glacier.  onlyglacieryes
50.0    %percentage of glacierization.  glacierpart
%************* 13.) DISCHARGE STARTING VALUES**********************
discharge.txt     %name of discharge output file.  nameqcalc
-9999    %nodata value of discharge file.  nodis
350    %storage constant k for firn.  firnkons
50    %storage constant k for snow.  snowkons
5    %storage constant k for ice.  icekons
30    %storage constant k for rock(outside glacier non-snowcovered.  icekons
%************* 13.) DISCHARGE STARTING VALUES**********************
5.0    %start value for firn discharge (previous time step).  qfirnstart
0.0    %start value for snow discharge (m3/s).  qsnowstart
5.0    %start value for ice discharge.  qicestart
5.0    %start value for rock discharge (outside glacier non-snowcovered).  qicestart
0.0    %groundwater discharge[m3].  qground
0.0    %difference between start of calculation and start r2.  jdstartr2diff
%%%%%%%%%%%% 15.) OPTIMIZATION %%%%%%%%%%%%%%%%%%%%%
0    %optimization run for k-values = 1; simulation=0.  disyesopt
radfactorsnow    %1. parameter to optimize.  optkA
0.4    %startvalue of 1. parameter to optimize.  startopt1
0.2    %step length.  stepopt1
7    %number of steps no.1  anzahlopt1
radfactorice    %2. parameter to optimize.  optkB
0.4    %startvalue of 2. parameter.  startopt2
0.2    %steplength no.2 of optimal r2.  stepopt2
8    %number of steps no.2.  anzahlopt2
matriz.dat    %name of r2-outputfile.  namematrix
%============ 16.) SNOW MODEL by C. Tijm-Reijmer 2/2005 ========================
1    %0=no percolation, 1=percolation+refreezing in snowlayer.  percolationyes
1    %0=no slush, 1=meltwater accumulation in snowlayer.  slushformationyes
1    %0=no densification, 1=densific. of dry snow due to aging.  densificationyes
0    %0=dry start, 1=wet start.  wetstartyes
1    %1=gradient method, 2=read snow model initialization profiles from file.  methodinitialverticalgrid
30    %maximum number of vertical layers.  ndepths
8    %number of subtimesteps for interpolation per main timestep 8.  factinter
%----------------------------
0.04    %layer thickness of first layer (m snow).  thicknessfirst
5.0    %layer thickness at deepest layer (m snow).  thicknessdeep
30.0    %maximum depth model (m).  depthdeep
200    %density of fresh snowfall kg/m3.  denssnow
1    %0=constant irreducible water content, 1=density dep. Schneider, 2= density dep. Coleou.  irrwatercontyes
0.02    %fraction of space irreducible filled with water.  irrwatercont
%---- Output ----------
1    %factor for subsurf output, 1=every hour, 24=once per day at midnight.  factsubsurfout
12    %offfsetfactor for subsurf output to make print at noon possible.  offsetsubsurfout
%runoffyes superyes   wateryes   surfwateryes slushyes   coldsnowyes   coldtotyes.      for grid output
     1       1	 	 1		 1		     1		   1		    1
%========================================================
%   17.) TEMPERATURE INDEX METHOD
%========================================================
1    %which temp index method (1,2 or 3).  ddmethod
8.0    %degree day factor for ice (only simple DDF method 1).  DDFice
5.5    %degree day factor for snow (only simple DDF method 1).  DDFsnow
%---------------------------------------------------
2.7    %meltfactor (only for modified temp index method 2 or 3).  meltfactor
0.9    %radiation melt factor for ice.  radfactorice
0.6    %radiation melt factor for snow.  radfactorsnow
1    %factor to reduce melt over debris.  debrisfactor
%******* 18.) OUTPUT STAKES **************************************
3    %1=stake locations given in local coordinates (x,y), 2=as 1 but center, 3=row/col.  coordinatesyes
100  100
102  103
102  99
102  104
103  105
