-
Notifications
You must be signed in to change notification settings - Fork 13
/
Copy pathmhm.nml
293 lines (275 loc) · 11.8 KB
/
mhm.nml
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
!> \file mhm.nml
!> \brief General Namelists of mHM, MPR, mRM for test domain 1.
!> \details This files provides all namelists for mHM, MPR, mRM.
!******************************************************************************************
! PROJECT DESCRIPTION (mandatory)
!******************************************************************************************
!-----------------------------------------------------------------------------
!> Provide details on the model simulations, to appear in the netcdf output attributes
!-----------------------------------------------------------------------------
&project_description
!> project name
project_details="mHM test domain project"
!> any specific description of simulation
setup_description="model run for the Mosel domain, forced with the E-OBS meteorologic data"
!> e.g. hindcast simulation, seasonal forecast, climate projection
simulation_type="historical simulation"
!> convention used for dataset
Conventions="XXX"
!> contact details, incl. PI name, modellers
contact="mHM developers (email:[email protected])"
!> developing institution, specific mHM revision, latest release version (automatically included)
mHM_details="Helmholtz Center for Environmental Research - UFZ, Department Computational Hydrosystems, Stochastic Hydrology Group"
!> some details on data/model run version (creation date is included automatically)
history="model run version 1"
/
!> Main namelist
&mainconfig
!> 0 -> regular X & Y coordinate system (e.g., GK-4 or Lambert equal area system)
iFlag_cordinate_sys = 0
!> Number of domains to be modeled.
nDomains = 1
!> resolution of Level-1 hydrological simulations in mHM [m] per domain
resolution_Hydrology(1) = 24000
!> specify same index for domains to share L0_data to save memory
L0Domain(1) = 1
!> flag for writing restart output
write_restart = .FALSE.
!> read domain specific optional data (the program decides)
read_opt_domain_data(1) = 0
/
!> main config for mHM and mRM (mHM and mRM-related)
&mainconfig_mhm_mrm
!> directory where restart input is located
mhm_file_RestartIn(1) = "restart/mHM_restart_001.nc"
mrm_file_RestartIn(1) = "restart/mRM_restart_001.nc"
!> resolution of Level-11 discharge routing [m] per domain
resolution_Routing(1) = 24000
!> model run timestep [h] either 1 or 24
timestep = 1
!> disabled restart and optimization (but needs to be given)
read_restart = .FALSE.
optimize = .FALSE.
optimize_restart = .FALSE.
opti_method = 1
opti_function = 10
/
!******************************************************************************************
! main config for mRM (mRM-related)
!******************************************************************************************
&mainconfig_mrm
! all defaults (does this rais an error?!)
/
!> DIRECTORIES
&directories_general
!> config run out file common to all modeled domains should be written to directory
dirConfigOut = "./"
!
!> directory where common input files should be located for all modeled domains
!> (only for *_classdefinition files)
dirCommonFiles = "input/morph/"
!> directory where morphological files are located
dir_Morpho(1) = "input/morph/"
!> directory where land cover files are located
dir_LCover(1) = "input/luse/"
!> directory where restart output should be written
mhm_file_RestartOut(1) = "restart/mHM_restart_001.nc"
mrm_file_RestartOut(1) = "restart/mRM_restart_001.nc"
!> directory where output should be written
dir_Out(1) = "output_b1/"
!> file containing latitude and longitude on the resolution_Hydrology
file_LatLon(1) = "input/latlon/latlon_1.nc"
/
!******************************************************************************************
! directories (mHM-related)
!******************************************************************************************
&directories_mHM
!> input format specification for the meteorological forcings: 'nc' only possible
inputFormat_meteo_forcings = "nc"
!> directory where meteorological input is located
dir_Precipitation(1) = "input/meteo/pre/"
dir_Temperature(1) = "input/meteo/tavg/"
!> if processCase(5) == 0 input directory of pet has to be specified
dir_ReferenceET(1) = "input/meteo/pet/"
!> switch to control read input frequency of the gridded meteo input: only at beginning of the run
time_step_model_inputs(1) = 0
/
!> directories (mRM-related)
&directories_mRM
!> directory where discharge files are located
dir_Gauges(1) = "input/gauge/"
!> directory where simulated runoff can be found (only required if coupling mode equals 0)
dir_Total_Runoff(1) = 'output_b1/'
!> directory where runoff at bankfull conditions can be found (only for coupling to groundwater model)
dir_Bankfull_Runoff(1) = 'input/optional_data/'
/
!> PROCESSES (mandatory)
&processSelection
!> interception
!> 1 - maximum Interception
processCase(1) = 1
!> snow
!> 1 - degree-day approach
processCase(2) = 1
!> soil moisture
!> 1 - Feddes equation for ET reduction, multi-layer infiltration capacity approach, Brooks-Corey like
processCase(3) = 1
!> directRunoff
!> 1 - linear reservoir exceedance approach
processCase(4) = 1
!> interflow
!> 1 - storage reservoir with one outflow threshold and nonlinear response
processCase(6) = 1
!> percolation
!> 1 - GW assumed as linear reservoir
processCase(7) = 1
!> routing
!> 3 - adaptive timestep, spatially varying celerity
processCase(8) = 3
!> baseflow
!> 1 - recession parameters (not regionalized yet)
processCase(9) = 1
/
!> LAND COVER (mandatory)
&LCover
!> number of land cover scenes to be used
nLCoverScene = 2
!> starting year of land cover scene 1
LCoverYearStart(1) = 1981
!> ending year of land cover scnene 1
LCoverYearEnd(1) = 1990
!> name of land cover file for scnene 1
LCoverfName(1) = 'lc_1981.asc'
!> starting year of land cover scene 2
LCoverYearStart(2) = 1991
!> ending year of land cover scnene 2
LCoverYearEnd(2) = 2000
!> name of land cover file for scnene 2
LCoverfName(2) = 'lc_1991.asc'
/
!******************************************************************************************
! Time periods (mHM and mRM-related)
!******************************************************************************************
&time_periods
!> specification of number of warming days [d] and the simulation period.
warming_Days(1) = 0
!> first year of wanted simulation period
eval_Per(1)%yStart = 1990
!> first month of wanted simulation period
eval_Per(1)%mStart = 01
!> first day of wanted simulation period
eval_Per(1)%dStart = 01
!> last year of wanted simulation period
eval_Per(1)%yEnd = 1993
!> last month of wanted simulation period
eval_Per(1)%mEnd = 12
!> last day of wanted simulation period
eval_Per(1)%dEnd = 31
/
!> INPUT SOIL DATABASE AND mHM LAYERING (MPR-related)
&soildata
!> iFlag_soilDB:
!> flag to handle multiple types of soil databases and their processing within the mHM.
!> This flag is unique and valid for all domains.
!> Depending on the choice of this flag you need to process your soil database differently.
!
!> iFlag_soilDB = 0:
!> Read and process the soil database in a classical mHM format which requires:
!> i) a single gridded ASCII file of soil-id (soil_class.asc - hard coded file name)
!> ii) a single soil look-up-table file (soil_classdefinition.txt) with information of
!> soil textural properties for every horizon.
!
!> Here mHM is quite flexible to handle multiple soil layers as specified in "nSoilHorizons_mHM"
!> and depths provided in "soil_Depth(:)".
!
!> The tillage depth is flexible in this case.
!
!> The depth of last mHM modeling layer is determined according the information given in the
!> input soil database, which could vary spatially depending on the soil type. Therefore the
!> user should not provide the depth of the last modeling layer. For example if you choose
!> nSoilHorizons_mHM = 3, then soil_Depth should be given for only soil_Depth(1) and soil_Depth(2).
!> Any additional depth related information would be discarded
iFlag_soilDB = 0
!> [mm] soil depth down to which organic matter is possible
tillageDepth = 200
!> No. of soil horizons to be modeled
nSoilHorizons_mHM = 2
!> Provide below the soil_Depth() for 1,2,..,*n-1* soil horizons. Depth of the last layer(n) is taken from the soil LUT
!> Soil_Horizon Depth[mm] ! bottom depth of soil horizons w.r.t. ground surface (positive downwards)
soil_Depth(1) = 200
/
!> INFORMATION RELATED TO LAI DATA (MPR-related)
&LAI_data_information
!> timeStep_LAI_input
!> 0: read LAI from long term monthly mean lookup table (related to land cover file).
!> The filename (LAI_classdefinition.txt) for the LUT is hard coded in mo_file.f90
!> Information regarding long-term monthly mean LAI for land cover classes
!> appearing in all modeled domains should be included in this LUT file.
!> This is an unique file applicable to all domains to be modeled.
!> The respective plant functional type is in LAI_class.asc, which must be also given
!> and should be located in each domain's morph directory.
timeStep_LAI_input = 0
!> nc - assume one file with name lai.nc with 12 monthly grids of mean LAI estimates
inputFormat_gridded_LAI = "nc"
/
!> LCover information (MPR-related)
&LCover_MPR
!>fraction of area within city assumed to be fully sealed [0.0-1.0]
fracSealed_cityArea = 0.6
/
!******************************************************************************************
! Specifcation of evaluation and inflow gauges (mRM-related)
!******************************************************************************************
!> namelist controlling the gauging station information
!> The ID has to correspond to the ID's given in the 'gaugelocation.asc' and
!> to the filename containing the time series
&evaluation_gauges
!> Gauges for model evaluation
!
!> Total number of gauges (sum of all gauges in all domains)
nGaugesTotal = 1
!> structure of gauge_id(i,j) & gauge_filename(i,j):
!> 1st dimension is the number of the subdomain i
!> 2nd dimension is the number of the gauge j within the subdomain i
!> numbering has to be consecutive
!
!> domain 1
!> number of gauges for subdomain (1)
NoGauges_domain(1) = 1
!> in subdomain(1), this is the id of gauge(1) --> (1,1)
Gauge_id(1,1) = 398
!> name of file with timeseries for subdomain(1) at gauge(1) --> (1,1)
gauge_filename(1,1) = "00398.txt"
/
&inflow_gauges
!> consider flows from upstream/headwater cells of inflow gauge(1) for subdomain(1) --> (1,1)
InflowGauge_Headwater(1,1) = .FALSE.
/
!******************************************************************************************
! ANNUAL CYCLE PAN EVAPORATION (mHM-related)
!******************************************************************************************
&panEvapo
! MONTH Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
!> monthly free pan evaporation
evap_coeff = 1.30, 1.20, 0.72, 0.75, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.50
/
!******************************************************************************************
! ANNUAL CYCLE METEOROLOGICAL FORCINGS (mHM-related)
!******************************************************************************************
&nightDayRatio
!> If read_meteo_weights is False than night fractions below are used
read_meteo_weights = .FALSE.
!> night ratio for precipitation
!> only night values required because day values are the opposite
fnight_prec = 0.46, 0.50, 0.52, 0.51, 0.48, 0.50, 0.49, 0.48, 0.52, 0.56, 0.50, 0.47
!> night ratio for PET
fnight_pet = 0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10, 0.10
!> night correction factor for temperature
fnight_temp = -0.76, -1.30, -1.88, -2.38, -2.72, -2.75, -2.74, -3.04, -2.44, -1.60, -0.94, -0.53
/
&Optimization
!> needed values (although optimization is turned of)
nIterations = 7
seed = 1235876
sce_ngs = 2
/