%PIC_DETIM_v4: NOTE: YOU NEED TO ADJUST THE INPUT AND OUTPUT FOLDER (2. AND 3. ROW BELOW), AND MAKE AN OUTPUT FOLDER
%************************************************************************************
100   %output to screen every X day.  daysscreenoutput.
/reginefolder/MODEL/studies/NicoleSchaefer2016/indata/
/reginefolder/MODEL/studies/NicoleSchaefer2016/output/
260  2008  %first julian day to be calculated.  jdbeg  yearbeg
259  2010   %last julian day to be calculated.  jdend  yearend
0    %discharge to be calculated: 1=yes,0=no,2=yes, but no discharge data.  disyes
1    %1=whole grid computed,  2=only grid cell of weather station.  calcgridyes
%********* 1.) MODEL OUTPUT PARAMETERS: DEFINE WHAT IS WRITTEN TO FILE ***************
2    %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          0          0          0         0             1                 1      1
0    %surface conditions to grid file (2=for specified JDs).  surfyes
0    %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
%******************************** 2.) MASS BALANCE OUTPUT ********************************
1       %output grid of winter mass balance     yes=1, no=0.  winterbalyes
1       %output grid of summer mass balance, yes=1, no=0.  summerbalyes
260     %day of year when winter starts. Summer ends day before. winterjdbeg 
118     %day of year when winter ends. Summer starts next day.   winterjdend
0       %0=fixed dates above, 1=dates for each year read from file.  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 **colu/y-coord ** glob and net data included from input data
pointoutputAWS.txt 717  1267  0    %outgridname[][], stn_xcoordinate[i], stn_ycoordinate[i], &outglobnet[i]
%******** 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.25    %gamma in V-A scaling.  gamma
0.6    %coefficient in V-A scaling.  c_coefficient
%********* 4.) NAMES OF INPUT FILES *********************************************
026_dem_utm_cut.bin    %name of Digital Terrain Model.  namedgm
026dm2014v22_watershed_cut.bin    %name of DTM with drainage basin.  namedgmdrain
026dm2014v22_cut.bin    %name of DTM glacier.  namedgmglac
026_dem_slope_cut.bin    %name of DTM slope.  namedgmslope
026_dem_aspect_cut.bin    %name of DTM aspect.  namedgmaspec
026_sky_dummy.bin    %name of DTM sky view factor.  namedgmskyview
026dem1646v22_firn_cut.bin    %name of DTM firnarea.  namedgmfirn
026dm2014v22_SWE_cut.bin   %name of DTM initial snow cover.  nameinitialsnow 
026_thickness_dummy.bin  %name of DTM ice thickness.  namedgmthickness 
RACMO2_climate_m0714.txt    %name of climate data file.  nameklima
%********* 5.) GRID INFORMATION**************************************************
-65.855    %geographical longitude [degree].  laenge
67.254    %latitude.  breite
-60    %longitude time refers to.  reflongitude
717    %row in DTM where climate station is located.   rowclim  ...this is the grid cell for the AWS
1267    %column of climate station.  colclim...this is the grid cell for the AWS
0  1817    %take this elevation for AWS yes/no.  climoutsideyes  heightclim
60.90279    %gridsize in m.  gridsize
24    %time step in hours.  timestep
%********* 6.) CLIMATE DATA *****************************************************
2    %1=midnight time is 0, 2=time is 24, 3=24 but previous day.  formatclimdata
5    %number of columns in climate file.  maxcol
4    %columns in climate input file: temperature.  coltemp
30    %column containing relative humidity.  colhum
30    %column wind speed [m/s].  colwind
30    %global radiation.  colglob
30  %reflected shortwave radiation.  colref
30   %net radiation.  colnet
30    %longwave incoming radiation.  collongin
30    %longwave outgoing radiation.  collongout
5    %precipitation.  colprec
30    %cloud cover (number of eigths).  colcloud
30    %time variant lapse rate (neg=decrease).  coltempgradvarying
30    %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.447    %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
1      -0.394    -0.399  -0.343   -0.452 -0.514 -0.423 -0.381 -0.426 -0.568 -0.491 -0.51 -0.468    %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
0.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
1    %direct radiation read from file = 1.  directfromfile
/cygdrive/c/meltmodel/Ipot_cut3/    %Path for direct files.  pathdirectfile
1    %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
3.6   %precipitation change with elevation [%/100m].  precgrad
0.0 1900    %precipitation change with elevation beyond certain elevation.  precgradhigh precgradelev
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.0    %threshold temperature rain/snow precipitation.  threshtemp
1    %0=if massbal calculations for whole drainage basin, 1=only glacier.  onlyglacieryes
100   %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.  rockkons
%************* 14.) 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.) 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
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
%========================================================
%   16.) 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
678  1269
724  1249
