-
Notifications
You must be signed in to change notification settings - Fork 10
/
settings_CWatM_template_30min.ini
860 lines (662 loc) · 29.5 KB
/
settings_CWatM_template_30min.ini
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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
# ------------------------------------------------
######## ## ## #### ###### ## ##
## ## ## ## ## ## #### ####
## ## ## ## ## ## ## #### ##
## ## ## ## ######## ## ## ## ##
## ## #### ## ## ## ## ## ##
## #### #### ## ## ## ## ##
########## ## ## ## ## ## ## ##
# Community Water Model Version 0.99
# SETTINGS FILE
# ------------------------------------------------
[OPTIONS]
#-------------------------------------------------------
# OPTION - to switch on/off
#-------------------------------------------------------
# Data otions
# if temperature is stored in Kelvin instead Celsius
TemperatureInKelvin = True
# if lat/lon the area has to be user defined = precalculated
gridSizeUserDefined = True
#-----------------------------------------------
# Evaporation: calculate pot. evaporation (True) or use precalculated pot.evaporation map stacks (False)
calc_evaporation = True
#-----------------------------------------------
# Irrigation and water demand
# if irrigation is included, otherwise paddy and non paddy is put into 'grassland'
includeIrrigation = True
# if water demand from irrigation, industry and domestic is included
includeWaterDemand = True
# Water allocation
# if water demand and availability is calculated for region to compare demand vs. avail
usingAllocSegments = False
# limit abstraction to available groundwater (True) include fossil groundwater (False)
limitAbstraction = False
# fossil water that is not consumed is not lost, but treated normally
fossil_water_treated_normally = True
#Sector- source-abstraction fractions are specified, in WATERDEMAND
sectorSourceAbstractionFractions = True
# Environmental Flow
calc_environflow = False
#-----------------------------------------------
# Soil
# use preferential flow, that bypasses the soil matrix and drains directly to the groundwater (not for irrPaddy)
preferentialFlow = True
# Capillar rise
CapillarRise = True
#-----------------------------------------------
# Routing
# if runoff concentration to the edge of a cell is included
includeRunoffConcentration = True
# Waterbodies like lakes and reservoirs
includeWaterBodies = True
# kinematic wave routing, if False no routing is calculated
includeRouting = True
#-----------------------------------------------
# Inflow from outside of the modelled area
inflow = False
# --- Reporting & Debugging ------------------- ----------------------------------------------
# Reporting options
writeNetcdfStack = True
reportMap = True
reportTss = True
# Checking water balance (for debugging)
calcWaterBalance = False
sumWaterBalance = False
# use additional PCRaster GIS commands
PCRaster = False
includeCrops = False
use_GeneralCropnonIrr = 0
use_GeneralCropIrr = 0
activate_fallow = 1
automaticFallowingIrr = True
leftoverIrrigatedCropIsRainfed = 0
moveIrrFallowToNonIrr = False
#-------------------------------------------------------
# DEFINITIONS OF PARAMETERS
#-------------------------------------------------------
#-------------------------------------------------------
[FILE_PATHS]
#-------------------------------------------------------
PathRoot = C:\GitHub\CWatM-Input\CWatM-Earth-30min
PathOut = C:\nonGithub\Output2
PathMaps = $(PathRoot)
#PathMeteo = P:\watxene\ISIMIP\ISIMIP3a\InputData\climate\obsclim_updated\GSWP3-W5E5
PathMeteo = C:\Users\smilovic\OneDrive - IIASA\cwatm_input\Climate\GSWP3-W5E5
PathWaterdemand = $(PathMaps)/landsurface/waterDemand
Excel_settings_file = $(PathMaps)\settings\cwatm_settings.xlsx
#-------------------------------------------------------
[NETCDF_ATTRIBUTES]
institution = IIASA
title = Rhine Water Model - WATCH WDFEI
metaNetcdfFile = $(FILE_PATHS:PathMaps)/metaNetcdf.xml
#-------------------------------------------------------
# AREA AND OUTLETS
#-------------------------------------------------------
[MASK_OUTLET]
# Area mask
# A pcraster map, tif or netcdf map e.g. $(FILE_PATHS:PathRoot)/data/areamaps/area_indus.map
# or a retancle: Number of Cols, Number of rows, cellsize, upper left corner X, upper left corner Y
#MaskMap = $(FILE_PATHS:PathRoot)/areamaps/rhine30min.tif
# Nile 30 arcmin
#MaskMap = 32.7082 12.4218
# Danube 30 arcmin
#MaskMap = 29.618 45.2613
# LV 30 arcmin
#MaskMap = 31.209 2.7517
# Mekong 30 arcmin
#MaskMap = 106.61 10.23
# Weser 30 arcmin
MaskMap = 8.234 53.324
# Rhine 30 arcmin
#MaskMap = 90.742 22.701
#Saudi Arabia
#MaskMap = 25 25 0.5 34.0 31.0
#Indus
#MaskMap = 68.1986 24.7312
# Yellow River
#MaskMap = 117.7316 37.287
#-------------------------------------------------------
# Station data
# either a map e.g. $(FILE_PATHS:PathRoot)/data/areamaps/area3.map
# or a location coordinates (X,Y) e.g. 5.75 52.25 9.25 49.75 )
# Lobith/Rhine
# Longitude Latitude -180 to 180, same as Earth Engine
#Gauges = 32.7082 12.4218
#Gauges = 29.618 45.2613
#Gauges = 31.209 2.7517
#Gauges = 106.61 10.23
#Gauges = 8.234 53.324
#Gauges = 90.742 22.701
#Gauges = 35 30
#Gauges = 68.1986 24.7312
#Gauges = 117.7316 37.287
Gauges = 8.234 53.324
# if .tif file for gauges, this is a flag if the file is global or local
# e.g. Gauges = $(FILE_PATHS:PathRoot)/data/areamaps/gaugesRhine.tif
GaugesLocal = True
#-------------------------------------------------------
[TIME-RELATED_CONSTANTS]
#-------------------------------------------------------
# StepStart has to be a date e.g. 01/06/1990
# SpinUp or StepEnd either date or numbers
# SpinUp: from this date output is generated (up to this day: warm up)
StepStart = 01/01/2000
SpinUp = None
StepEnd = 02/01/2001
#-------------------------------------------------------
[INITITIAL CONDITIONS]
#-------------------------------------------------------
# for a warm start initial variables a loaded
# e.g for a start on 01/01/2010 load variable from 31/12/2009
load_initial = False
initLoad = $(FILE_PATHS:PathRoot)/init/Yellow_20000201.nc
#initLoad = $(FILE_PATHS:PathRoot)/init/Rhine_19891231.nc
# saving variables from this run, to initiate a warm start next run
# StepInit = saving date, can be more than one: 10/01/1973 20/01/1973
save_initial = False
initSave = $(FILE_PATHS:PathRoot)/init/Yellow
StepInit = 01/02/2000
#-------------------------------------------------------
# CALIBARTION PARAMETERS
#-------------------------------------------------------
[CALIBRATION]
# These are parameter which are used for calibration
# could be any parameter, but for an easier overview, tehey are collected here
# in the calibration template a placeholder (e.g. %arnoBeta) instead of value
# Snow
SnowMeltCoef = 0.0027
# Cropf factor correction
crop_correct = 1.11
crop_correct_forest = 1
crop_correct_grassland = 1
crop_correct_irrpaddy = 1
crop_correct_irrnonpaddy = 1
#Soil
soildepth_factor = 1.28
#Soil preferentialFlowConstant = 4.0, arnoBeta_add = 0.1
preferentialFlowConstant = 4.5
arnoBeta_add = 0.19
# interflow part of recharge factor = 1.0
factor_interflow = 2.8
# groundwater recessionCoeff_factor = 1.0
recessionCoeff_factor = 5.278
# runoff concentration factor runoffConc_factor = 1.0
runoffConc_factor = 0.1
#Routing manningsN Factor to Manning's roughness = 1.0 [0.1-10.]
manningsN = 1.86
# reservoir normal storage limit (fraction of total storage, [-]) [0.15 - 0.85] default 0.5
normalStorageLimit = 0.44
# lake parameter - factor to alpha: parameter of of channel width and weir coefficient [0.33 - 3.] dafault 1.
lakeAFactor = 0.33
# lake parameter - factor for wind evaporation
lakeEvaFactor = 1.52
#-------------------------------------------------------
# TOPOGRAPHY MAPS
#-------------------------------------------------------
[TOPOP]
# local drain direction map (1-9)
Ldd = $(FILE_PATHS:PathMaps)/routing/ldd.map
# Elevation standard deviation [m], i.e. altitude difference elevation within pixel.
# Used for sub-pixel modelling of snow accumulation and melt
ElevationStD = $(FILE_PATHS:PathMaps)/landsurface/topo/elvstd.map
# Area of pixel [m2] (for lat/lon every cell has a different area)
CellArea = $(FILE_PATHS:PathMaps)/routing/cellarea.map
#-------------------------------------------------------
# INPUT METEOROLOGICAL TIMESERIES AS MAPS
#-------------------------------------------------------
[METEO]
# precipitation [kg m-2 s-1]
PrecipitationMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_pr*
# average daily temperature [K]
TavgMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_tas*
# ----------------------------------------
# This is used if calc_evaporation = False
# daily reference evaporation (free water)
E0Maps = $(FILE_PATHS:PathMeteo)/EWRef_rhine.nc
#E0Maps = $(FILE_PATHS:PathMeteo)/EWRef_daily*
# daily reference evapotranspiration (crop)
ETMaps = $(FILE_PATHS:PathMeteo)/ETRef_rhine.nc
#ETMaps = $(FILE_PATHS:PathMeteo)/ETRef_daily*
# --------------------------------
# from kg m-2s-1 to m : 86.4
precipitation_coversion = 86.4
# from MM to m : 0.001
#precipitation_coversion = 0.001
evaporation_coversion = 1.00
# OUTPUT maps and timeseries
#OUT_Dir = $(FILE_PATHS:PathOut)
#OUT_MAP_Daily = Precipitation, prec1
#-------------------------------------------------------
# CALCULATE EVAPORATION - PENMAN - MONTEITH
#-------------------------------------------------------
[EVAPORATION]
# This is used if calc_evaporation = True
# use albedo maps
albedo = True
albedoMaps = $(FILE_PATHS:PathMaps)/landsurface/albedo/albedo.nc
# if not albedo maps use fixed albedo
# Albedo of bare soil surface (Supit et. al.)
AlbedoSoil = 0.15
# Albedo of water surface (Supit et. al.)
AlbedoWater = 0.05
# Albedo of vegetation canopy (FAO,1998)
AlbedoCanopy = 0.23
# use specific humidity (TRUE) QAir, or relative humidity (FALSE) - rhs
useHuss = True
# map stacks Temperature [K}]
TminMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_tasmin*
TmaxMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_tasmax*
# Instantaneous surface pressure[Pa]
PSurfMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_ps*
# 2 m istantaneous specific humidity[kg /kg] (QAir) or relative humidity [%] (rhs)
RhsMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_hurs*
QAirMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_huss*
# wind speed maps at 10m [m/s]
WindMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_sfcwind*
# radiation surface downwelling shortwave maps [W/m2]
RSDSMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_rsds*
# radiation surface downwelling longwave maps [W/m2] [W/m2]
RSDLMaps = $(FILE_PATHS:PathMeteo)/gswp3-w5e5_obsclim_rlds*
#-------------------------------------------------------
[SNOW]
#-------------------------------------------------------
# Number of vertical Snow layers
NumberSnowLayers = 7
# up to which layer the ice melt is calculated with the middle temperature
GlacierTransportZone = 3
# Temperature lapse rate with altitude [deg C / m]
TemperatureLapseRate = 0.0065
# Multiplier applied to precipitation that falls as snow
SnowFactor = 1.0
# Range [m C-1 d-1] of the seasonal variation, SnowMeltCoef is the average value
SnowSeasonAdj = 0.001
# Average temperature at which snow melts
TempMelt =1.0
# Average temperature below which precipitation is snow
TempSnow = 1.0
# Snow melt coefficient: default: 4.0
# SRM: 0.0045 m/C/day ( = 4.50 mm/C/day), Kwadijk: 18 mm/C/month (= 0.59 mm/C/day)
# See also Martinec et al., 1998.
# use in CALIBRATION -> copied to CALIBRATION
#SnowMeltCoef = 0.004
IceMeltCoef = 0.007
#--------------------------------------------------------------
# INITIAL CONDITIONS - Initial snow depth in snow zone 1-7 [mm] - SnowCoverIni
[FROST]
# Snow water equivalent, (based on snow density of 450 kg/m3) (e.g. Tarboton and Luce, 1996)
SnowWaterEquivalent = 0.45
# Daily decay coefficient, (Handbook of Hydrology, p. 7.28)
Afrost = 0.97
# Snow depth reduction coefficient, [cm-1], (HH, p. 7.28)
Kfrost = 0.57
# Degree Days Frost Threshold (stops infiltration, percolation and capillary rise)
# Molnau and Bissel found a value 56-85 for NW USA.
FrostIndexThreshold = 56
#--------------------------------------------------------------
# INITIAL CONDITIONS: FrostIndexIni
[VEGETATION]
cropgroupnumber = $(FILE_PATHS:PathMaps)/soil/cropgrp.nc
# soil water depletion fraction, Van Diepen et al., 1988: WOFOST 6.0, p.86, Doorenbos et. al 1978
#-------------------------------------------------------
[SOIL]
#-------------------------------------------------------
PathTopo = $(FILE_PATHS:PathMaps)/landsurface/topo
PathSoil = $(FILE_PATHS:PathMaps)/soil
# Topography mapsNetcdf/maps
# tangent slope, slope length and orographybeta
tanslope = $(PathTopo)/tanslope.map
slopeLength = $(PathTopo)/slopeLength.map
# maps of relative elevation above flood plains
relativeElevation = $(PathTopo)/dzRel_hydro1k.nc
# Soil hydraulic properties
# soil (Hypres pedotransfer function - http://esdac.jrc.ec.europa.eu/ESDB_Archive/ESDBv2/popup/hy_param.htm)
KSat1 = $(PathSoil)/ksat1.map
KSat2 = $(PathSoil)/ksat2.map
KSat3 = $(PathSoil)/ksat3.map
# Alpha: an Genuchten’s shape parameter
alpha1 = $(PathSoil)/alpha1.map
alpha2 = $(PathSoil)/alpha2.map
alpha3 = $(PathSoil)/alpha3.map
#Lambda: an Genuchten’s shape parameter = n-1-> n = lamda+1, m = 1 - (1/n)
lambda1 = $(PathSoil)/lambda1.map
lambda2 = $(PathSoil)/lambda2.map
lambda3 = $(PathSoil)/lambda3.map
# thetas is the volumetric water content θ saturated
thetas1 = $(PathSoil)/thetas1.map
thetas2 = $(PathSoil)/thetas2.map
thetas3 = $(PathSoil)/thetas3.map
# thetar is the volumetric water content θ residual
thetar1 = $(PathSoil)/thetar1.map
thetar2 = $(PathSoil)/thetar2.map
thetar3 = $(PathSoil)/thetar3.map
percolationImp = $(PathSoil)/percolationImp.map
maxGWCapRise = 5.0
minCropKC = 0.2
minTopWaterLayer = 0.0
# Soil depth
StorDepth1 = $(PathSoil)/storageDepth1.map
StorDepth2 = $(PathSoil)/storageDepth2.map
# preferential flow (between 1.0 and 8.0)
# used in CALIBRATION -> copied to CALIBRATION
#preferentialFlowConstant = 4.0
#-------------------------------------------------------
[LANDCOVER]
PathLandcover = $(FILE_PATHS:PathMaps)/landsurface
coverTypes = forest, grassland, irrPaddy, irrNonPaddy, sealed, water
coverTypesShort = f, g, i, n, s, w
fractionLandcover = $(PathLandcover)/fractionLandcover.nc
# Landcover can vary from year to year
dynamicLandcover = True
# if landcover cannot vary, which year should be taken as fixed year
fixLandcoverYear = 1961
#-------------------------------------------------------
[__forest]
PathForest = $(FILE_PATHS:PathMaps)/landcover/forest
PathSoil = $(FILE_PATHS:PathMaps)/soil
# Parameters for the Arno's scheme
# arnoBeta is defined by orographic,+ land cover add + calibration add, the soil water capacity distribution is based on this
# range [0.01 - 1.2]
forest_arnoBeta = 0.2
#forest_soil
forest_KSat1 = $(PathSoil)/forest_ksat1.map
forest_KSat2 = $(PathSoil)/forest_ksat2.map
forest_KSat3 = $(PathSoil)/ksat3.map
forest_alpha1 = $(PathSoil)/forest_alpha1.map
forest_alpha2 = $(PathSoil)/forest_alpha2.map
forest_alpha3 = $(PathSoil)/alpha3.map
forest_lambda1 = $(PathSoil)/forest_lambda1.map
forest_lambda2 = $(PathSoil)/forest_lambda2.map
forest_lambda3 = $(PathSoil)/lambda3.map
forest_thetas1 = $(PathSoil)/forest_thetas1.map
forest_thetas2 = $(PathSoil)/forest_thetas2.map
forest_thetas3 = $(PathSoil)/thetas3.map
forest_thetar1 = $(PathSoil)/forest_thetar1.map
forest_thetar2 = $(PathSoil)/forest_thetar2.map
forest_thetar3 = $(PathSoil)/thetar3.map
# other paramater values
forest_minInterceptCap = 0.001
forest_cropDeplFactor = 0.0
forest_fracVegCover = $(PathForest)/fracVegCover.map
forest_rootFraction1 = $(PathForest)/rootFraction1.map
forest_rootFraction2 = $(PathForest)/rootFraction2.map
#forest_maxRootDepth = 2.0
forest_maxRootDepth = $(PathForest)/maxRootDepth.map
forest_minSoilDepthFrac = $(PathForest)/minSoilDepthFrac.map
forest_cropCoefficientNC = $(PathForest)/cropCoefficientForest_10days.nc
forest_interceptCapNC = $(PathForest)/interceptCapForest_10days.nc
# initial conditions: forest_interceptStor, forest_w1, forest_w2, forest_w3,
[__grassland]
PathGrassland = $(FILE_PATHS:PathMaps)/landcover/grassland
# Parameters for the Arno's scheme:
grassland_arnoBeta = 0.0
# arnoBeta is defined by orographic,+ land cover add + calibration add, the soil water capacity distribution is based on this
# range [0.01 - 1.2]
# other paramater values
grassland_minInterceptCap = 0.001
grassland_cropDeplFactor = 0.0
grassland_fracVegCover = $(PathGrassland)/fracVegCover.map
grassland_rootFraction1 = $(PathGrassland)/rootFraction1.map
grassland_rootFraction2 = $(PathGrassland)/rootFraction2.map
grassland_maxRootDepth = $(PathGrassland)/maxRootDepth.map
grassland_minSoilDepthFrac = $(PathGrassland)/minSoilDepthFrac.map
grassland_cropCoefficientNC = $(PathGrassland)/cropCoefficientGrassland_10days.nc
grassland_interceptCapNC = $(PathGrassland)/interceptCapGrassland_10days.nc
# initial conditions: grassland_interceptSto, grassland_w1, grassland_w2, grassland_w3
[__irrPaddy]
PathIrrPaddy = $(FILE_PATHS:PathMaps)/landcover/irrPaddy
# Parameters for the Arno's scheme:
irrPaddy_arnoBeta = 0.2
# arnoBeta is defined by orographic,+ land cover add + calibration add, the soil water capacity distribution is based on this
# range [0.01 - 1.2]
# other paramater values
irrPaddy_minInterceptCap = 0.001
irrPaddy_cropDeplFactor = 0.0
irrPaddy_fracVegCover = $(PathIrrPaddy)/fracVegCover.map
irrPaddy_rootFraction1 = $(PathIrrPaddy)/rootFraction1.map
irrPaddy_rootFraction2 = $(PathIrrPaddy)/rootFraction2.map
irrPaddy_maxRootDepth = $(PathIrrPaddy)/maxRootDepth.map
irrPaddy_minSoilDepthFrac = $(PathIrrPaddy)/minSoilDepthFrac.map
irrPaddy_cropCoefficientNC = $(PathIrrPaddy)/cropCoefficientirrPaddy_10days.nc
# maximum flooding depth for paddy
irrPaddy_maxtopwater = 0.05
# initial conditions: irrPaddy_interceptStor, irrPaddy_w1, irrPaddy_w2, irrPaddy_w3
[__irrNonPaddy]
PathIrrNonPaddy = $(FILE_PATHS:PathMaps)/landcover/irrNonPaddy
# Parameters for the Arno's scheme:
irrNonPaddy_arnoBeta = 0.2
# arnoBeta is defined by orographic,+ land cover add + calibration add, the soil water capacity distribution is based on this
# range [0.01 - 1.2]
# other paramater values
irrNonPaddy_minInterceptCap = 0.001
irrNonPaddy_cropDeplFactor = 0.0
irrNonPaddy_fracVegCover = $(PathIrrNonPaddy)/fracVegCover.map
irrNonPaddy_rootFraction1 = $(PathIrrNonPaddy)/rootFraction1.map
irrNonPaddy_rootFraction2 = $(PathIrrNonPaddy)/rootFraction2.map
irrNonPaddy_maxRootDepth = $(PathIrrNonPaddy)/maxRootDepth.map
irrNonPaddy_minSoilDepthFrac = $(PathIrrNonPaddy)/minSoilDepthFrac.map
irrNonPaddy_cropCoefficientNC = $(PathIrrNonPaddy)/cropCoefficientirrNonPaddy_10days.nc
# initial conditions: irrNonPaddy_interceptStor, irrNonPaddy_w1, irrNonPaddy_w2, irrNonPaddy_w3
[__sealed]
PathSealed = $(FILE_PATHS:PathMaps)/landcover/sealed
sealed_minInterceptCap = 0.001
# initial conditions: sealed_interceptStor
[__open_water]
PathWater = $(FILE_PATHS:PathMaps)/landcover/water
water_minInterceptCap = 0.0
#-------------------------------------------------------
[GROUNDWATER]
#-------------------------------------------------------
PathGroundwater = $(FILE_PATHS:PathMaps)/groundwater
recessionCoeff = $(PathGroundwater)/recessionCoeff.map
# baseflow = recessionCoeff * storage groundwater
specificYield = $(PathGroundwater)/specificYield.map
kSatAquifer = $(PathGroundwater)/kSatAquifer.map
# both not used at the moment in groundwater modul, but already loaded
#--------------------------------------------------------------
# INITIAL CONDITIONS: storGroundwater
#-------------------------------------------------------
[WATERDEMAND]
#-------------------------------------------------------
# For water demand vs. availability: areas have to be aggregated
# Allocation map
allocSegments = $(FILE_PATHS:PathWaterdemand)/catchx.nc
# Demand in m/day [True] (default) or as mio m3 per year or month [False]
demand_unit = False
historicalIrrigationArea = $(FILE_PATHS:PathWaterdemand)/irrigationArea.nc
# Demand data monthly or yearly, as m day-1 or mio m3 per month - here m3/year
domesticWaterDemandFile = $(FILE_PATHS:PathWaterdemand)/dom_isimip_year_millionm3_1850_2018.nc
# time = True for monthly, False for yearly
domesticTimeMonthly = False
# name of the variable Withrawal = Gross, consumption = Netto
domesticWithdrawalvarname = domww
domesticConsuptionvarname = domwc
industryWaterDemandFile = $(FILE_PATHS:PathWaterdemand)/ind_isimip_year_millionm3_1850_2018.nc
industryTimeMonthly = False
industryWithdrawalvarname = indww
industryConsuptionvarname = indwc
uselivestock = True
livestockWaterDemandFile = $(FILE_PATHS:PathWaterdemand)/liv_isimip_year_millionm3_1850_2018.nc
livestockTimeMonthly = False
livestockvarname = livestockConsumption
irrNonPaddy_efficiency = $(FILE_PATHS:PathMaps)/landsurface/waterDemand/efficiency.nc
irrPaddy_efficiency = $(FILE_PATHS:PathMaps)/landsurface/waterDemand/efficiency.nc
# using environmental flow (EF) (per month) as input value
# EF will be treated as one part of overall water demand
use_environflow = False
EnvironmentalFlowFile = $(FILE_PATHS:PathOut)/MQ90_12month.nc
#irrNonPaddy_efficiency = 0.8
#irrPaddy_efficiency = 0.8
irrigation_returnfraction = 0.5
swAbstractionFraction_Channel_Domestic = 1
swAbstractionFraction_Channel_Livestock = 1
swAbstractionFraction_Channel_Industry = 1
swAbstractionFraction_Channel_Irrigation = 1
swAbstractionFraction_Lift_Domestic = 0
swAbstractionFraction_Lift_Livestock = 0
swAbstractionFraction_Lift_Industry = 0
swAbstractionFraction_Lift_Irrigation = 0
swAbstractionFraction_Lake_Domestic = 1
swAbstractionFraction_Lake_Livestock = 1
swAbstractionFraction_Lake_Industry = 1
swAbstractionFraction_Lake_Irrigation = 1
swAbstractionFraction_Res_Domestic = 1
swAbstractionFraction_Res_Livestock = 1
swAbstractionFraction_Res_Industry = 1
swAbstractionFraction_Res_Irrigation = 1
gwAbstractionFraction_Domestic = 1
gwAbstractionFraction_Livestock = 1
gwAbstractionFraction_Industry = 1
gwAbstractionFraction_Irrigation = 1
# -----------------------------------------------------------
# Estimate of fractions of groundwater and surface water abstractions
# Either a fixed fraction for surface water abstration
# based on fraction of average baseflow and upstream average discharge
# if swAbstractionFrac < 0: fraction is taken from baseflow / discharge
# if swAbstractionFrac > 0 this value is taken as a fixed value
swAbstractionFrac = 0.9
averageDischarge = $(FILE_PATHS:PathOut)/discharge_totalavg_rhine30min.nc
# in [m3/s]
averageBaseflow = $(FILE_PATHS:PathOut)/baseflow_totalavg_rhine30min.nc
# in [m3/s]
baseflowInM = True
# if baseflow is in [m] instead of [m3/s] it will be converted
#-------------------------------------------------------
# RUNOFF CONCENTRATION
#-------------------------------------------------------
[RUNOFF_CONCENTRATION]
# using triagular weigthning method
# the bigger the factor, more lag time
forest_runoff_peaktime = 1.0
grassland_runoff_peaktime = 0.5
irrPaddy_runoff_peaktime = 0.5
irrNonPaddy_runoff_peaktime = 0.5
sealed_runoff_peaktime = 0.15
water_runoff_peaktime = 0.01
interflow_runoff_peaktime =1.0
baseflow_runoff_peaktime = 2.0
# initial conditions:
# here only 1 layer is shown, but there are up to 10: runoff_concIni
#-------------------------------------------------------
# ROUTING MAPS and PARAMETERSD
#-------------------------------------------------------
[ROUTING]
PathRouting = $(FILE_PATHS:PathMaps)/routing
# Number of substep per day
# should be 10 for 0.5 deg but 24 for 0.1 deg
NoRoutingSteps = 10
#kinematic wave parameter: 0.6 is for broad sheet flow
chanBeta = 0.6
# Channel gradient (fraction, dy/dx)
chanGrad = $(PathRouting)/kinematic/changrad.nc
# Minimum channel gradient (for kin. wave: slope cannot be 0)
chanGradMin = 0.0001
#Channel Manning's n
chanMan = $(PathRouting)/kinematic/chanman.nc
#Channel length [meters]
chanLength = $(PathRouting)/kinematic/chanleng.nc
#Channel bottom width [meters]
chanWidth = $(PathRouting)/kinematic/chanbw.nc
#Bankfull channel depth [meters]
chanDepth = $(PathRouting)/kinematic/chanbnkf.nc
# initial conditions: channelStorageIni, riverbedExchangeIni, dischargeIni
#-------------------------------------------------------
# LAKES AND RESERVOIRS
#-------------------------------------------------------
[LAKES_RESERVOIRS]
PathLakesRes = $(FILE_PATHS:PathMaps)/routing/lakesreservoirs
# Use reservoirs and lakes (otherwise use only lakes Lake ID=1 and 3 => natural conditions)
useResAndLakes = True
# Reservoirs do have a year of implementation
dynamicLakesRes = True
# if Reservoirs does not have a year of implemtation, which year should be taken as fixed year
fixLakesResYear = 1950
#----------------------------
#Big lakes and Reservoirs
# ID of every lake, reservoir from HydroLakes database
waterBodyID = $(PathLakesRes)/lakesResID.nc
# 1 for lake, 2 for reservoir, 3 for lake and reservoir
waterBodyTyp = $(PathLakesRes)/lakesResType.nc
# Avergae discharge from HydroLakes Database
waterBodyDis = $(PathLakesRes)/lakesResDis.nc
# Lakes surface area from HydroLakes Database
waterBodyArea = $(PathLakesRes)/lakesResArea.nc
# a factor to scale the outlet of a lake
#lakeAFactor = 1.0 -> calibration
# Create a buffer around water bodies as command areas for lakes and reservoirs
buffer_waterbodies = 2
#----------------------------
# Small lakes and reservoirs
useSmallLakes = False
smallLakesRes = $(PathLakesRes)/smallLakesRes.nc
smallwaterBodyDis = $(PathLakesRes)/smalllakesresDis.nc
# averageRunoff in [m] (if not given smallwaterBodyDis is taken instead)
#averageRunoff = $(FILE_PATHS:PathOut)/runoff_totalavg_cali.nc
# for water demand
#min storage in [m3] (if not give it is calculated)
#minStorage = $(FILE_PATHS:PathOut)/minsmalllakeStorage_cali.nc
# initial conditions: lakeInflowIni, lakeStorageIni, outLakeIni, lakeOutflowIni, reservoirStorageIni
#------------------
# Reservoirs
# reservoir volume from HydroLakes database
waterBodyVolRes = $(PathLakesRes)/lakesResVolRes.nc
# reservoir starting year from HydroLakes database
waterBodyYear = $(PathLakesRes)/lakesResYear.nc
# Conservative, normal and flood storage limit (fraction of total storage, [-])
conservativeStorageLimit = 0.1
#normalStorageLimit = 0.5 # --> put into calibration
floodStorageLimit = 0.9
# adjusting the balance between normal and flood storage
# [0 ..1] 0: NormalstorageLimit 1: (= closer to flood) results in keeping the normal qoutflow longer constant
adjust_Normal_Flood = 0.5
# Minimum, Normal and Non-damaging reservoir outflow (fraction of average discharge, [-])
MinOutflowQ = 0.2
NormalOutflowQ = 1.0
NonDamagingOutflowQ = 4.0
# initial conditions: lakeInflowIni, lakeStorageIni, outLakeIni, lakeOutflowIni, reservoirStorageIni
#-------------------------------------------------------
[INFLOW]
#-------------------------------------------------------
# if option inflow = true
# the inflow from outside is added at inflowpoints
In_Dir = $(FILE_PATHS:PathRoot)/in
# nominal map with locations of (measured)inflow hydrographs [cu m / s]
InflowPoints = $(In_Dir)/in.map
#InflowPoints = 8.25 49.75 7.75 50.25
# if InflowPoints is a map, this flag is to identify if it is global (False) or local (True)
# observed or simulated input hydrographs as time series [cu m / s]
# Note: that identifiers in time series have to correspond to InflowPoints
# can be several timeseries in one file or different files e.g. main.tss mosel.tss
#QInTS = main1.tss mosel1.tss
QInTS = mm.tss
#-------------------------------------------------------
[ENVIRONMENTALFLOW]
#-------------------------------------------------------
# Either calculate without run with predone discharge (set calc_ef_after = False)
calc_ef_after = True
# Or calculate after run (set calc_ef_after = False) and defining the file to be used
EFDis = $(FILE_PATHS:PathOut)/discharge_rhine.nc
# if predone discharge, do the maps need to be cut to fit to the mask?
cut_ef_map = False
# MAF: Mean, Q90: percentile 90, MMF: monthly average, MQ90: monthly Q90 9averagwed over al Jan, Feb..
# EF_VMF: Environmental flow - variable monthly flow, EF_VMF_LIH - EF- variable monthly flow, high intermediate, low class
OUT_Dir = $(FILE_PATHS:PathOut)
#OUT_MAP_Once = MAF, Q90
#OUT_MAP_12month = MMF, MQ90, EF_VMF, EF_VMF_LIH
#OUT_MAP_12month = MQ90, EF_VMF
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++
[OUTPUT]
# OUTPUT maps and timeseries
OUT_Dir = $(FILE_PATHS:PathOut)
OUT_TSS_Daily = discharge
#OUT_TSS_MonthAvg = discharge
#OUT_TSS_AnnualAvg = discharge
#OUT_Map_Daily = discharge
#OUT_Map_MonthAvg = discharge, precipitation, runoff
#OUT_Map_AnnualAvg = discharge
#OUT_MAP_TotalAvg = discharge, baseflow
OUT_MAP_Daily = Rain, Snow, IceMelt, SnowMelt, snowEvap, iceEvap, act_nonpaddyConsumption, act_paddyConsumption, act_livConsumption, returnflowIrr, returnflowNonIrr, GW_Irrigation, Res_Irrigation, Lake_Irrigation, Channel_Irrigation, GW_Industry, GW_Livestock, GW_Domestic, Res_Industry, Res_Livestock, Res_Domestic, Lake_Industry, Lake_Livestock, Lake_Domestic, Channel_Domestic, Channel_Livestock, Channel_Industry, prefFlow_GW, perc3toGW_GW, lakeResInflowM, lakeResOutflowM, act_indConsumption, act_domConsumption, act_livConsumption, act_irrWithdrawal, act_nonIrrWithdrawal, act_domWithdrawal, act_indWithdrawal, act_livWithdrawal, snowEvap, capillar, baseflow, actTransTotal_forest, actTransTotal_grasslands, actTransTotal_paddy, actTransTotal_nonpaddy, unmet_lost, unmetDemand, pot_GroundwaterAbstract, storGroundwater, discharge, nonFossilGroundwaterAbs, Precipitation, totalET, EvapoChannel, EvapWaterBodyM, act_nonIrrConsumption, channelStorage, lakeResStorage, totalSto, sum_actTransTotal, sum_actBareSoilEvap, sum_interceptEvap, sum_openWaterEvap, addtoevapotrans, lakeResInflowM, act_bigLakeResAbst, sum_gwRecharge, sum_capRiseFromGW, act_totalIrrConsumption, sum_runoff, returnFlow, act_SurfaceWaterAbstract, dis_outlet, act_irrNonpaddyWithdrawal
OUT_MAP_MonthAvg = tws, storGroundwater
OUT_MAP_MonthTot = unmetDemand
OUT_MAP_TotalEnd = cellArea