Newer
Older
#---------------------------------------------------------------------------------------
# GENERAL CONTROLS
#---------------------------------------------------------------------------------------
# DESCRIPTION
# path --- the directory, in which the model is launched
# (required to be set in some UNIX-systems)
# runmode 1 --- stand alone run
# runmode 2 --- running in atmospheric model as lake parametrization
# omp 0 --- OpenMP is not used
# omp 1 --- OpenMP is used
#---------------------------------------------------------------------------------------
#
path ''
runmode 1
omp 0
#
#----------------------------------------------------------------------------------------
# SPATIAL RESOLUTION OF THE MODEL
#----------------------------------------------------------------------------------------
# DESCRIPTION
# nstep_keps --- number of timesteps of k-epsilon parameterization per on model timestep
# M --- number of layers in water layer
# Mice --- number of layers in upper and deep layers
# ns --- number of levels in soil
# d_surf --- grid zooming parameter at the surface, n/d
# d_bot --- grid zooming parameter at the bottom, n/d
#----------------------------------------------------------------------------------------
#
nstep_keps 1
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
d_bot 1.E-2
#
#----------------------------------------------------------------------------------------
# CONTROLS FOR PHYSICS OF THE MODEL
#----------------------------------------------------------------------------------------
# DESCRIPTION
# PBL parameterization
# PBLpar -1 --- sensible, latent heat and momentum fluxes are given as input for the model
# PBLpar 0 --- the latent heat flux is set to zero, while sensible heat and momentum fluxes
# are constant in time, specified by sensflux0 and momflux0
# PBLpar 1 --- Businger-Dayer formulas (Monin-Obukhov theory) for exchange coefficients
# PBLpar 2 --- formulation from NH3d
# PBLpar 3 --- formulation from FLake
# PBLpar 4 --- formulation implemented by M.Chechin
# c_d --- the momentum exchange coefficient, n/d
# (if -999, momentum flux is calculated by surface flux scheme)
# waveenh 0 --- the shallow water correction of surface fluxes (Panin et al., 1996) is OFF
# waveenh 1 --- the shallow water correction of surface fluxes (Panin et al., 1996) is ON
# momflxpart 0 --- all momentum flux from the atmosphere is consumed by currents acceleration
# momflxpart 1 --- momentum flux from the atmosphere is partitioned between wave developemnt
# (controlled by fetch) and currents acceleration,
# following Lin et al. (2002, J. Phys. Ocean.)
# kwe --- the factor of turbulence enhancement by wave breaking (wave energy factor), n/d
# Relative to water currents wind
# relwind 0 --- relative wind is off
# relwind 1 --- relative wind is on
# Equation of state
# eos 1 --- from Hostetler model
# eos 2 --- from TEOS-2010
# eos 3 --- for Kivu lake including salinity
# nmeltpoint 1 --- melting point linearly dependent on salinity
# nmeltpoint 2 --- TEOS-2010 formula
# Turbulent mixing parameterization
# Turbpar 1 --- "Empirical" parametrization: Stepanenko, Lykosov (2005)
# Turbpar 2 --- "E-epsilon"("K-epsilon") parameterization: k=E**2/eps with
# prognostic equations for E and eps
# Turbpar 3 --- Nickuradze (NICK) formulation: Rodi (1993)
# Turbpar 4 --- Parabolic (PARAB) formulation: Engelund (1976)
# Turbpar 7 --- RNG (re-normalization group) formulation: Simoes (1998)
# stabfunc 1 --- constant stability functions (standard k-epsilon model)
# stabfunc 2 --- stability functions according to (Canuto et al., 2001)
# stabfunc 3 --- stability functions according to (Galperin et al., 1988)
# kepsbc 1 --- Neuman boundary conditions for unstratified sheared flow (Burchard, 2002)
# kepsbc 2 --- Neuman boundary conditions for unstratified non-sheared flow with wave breaking (Burchard, 2002)
# kepsbc 3 --- Neuman boundary conditions unstratified sheared flow with wave breaking (Burchard, 2002)
# kepsbc 4 --- Neuman boundary conditions for free convection
# Water surface albedo: variable or constant
# varalb 0 --- constant
# varalb 1 --- sun height dependent
# soiltype 1 --- the soil type is "sand"
# soiltype 2 --- the soil type is "loamy sand"
# soiltype 3 --- the soil type is "sandy loam"
# soiltype 4 --- the soil type is "loam"
# soiltype 5 --- the soil type is "silt loam"
# soiltype 6 --- the soil type is "sandy clay loam"
# soiltype 7 --- the soil type is "clay loam"
# soiltype 8 --- the soil type is "silty clay loam"
# soiltype 9 --- the soil type is "sandy clay"
# soiltype 10 --- the soil type is "silty clay"
# soiltype 11 --- the soil type is "clay"
# soil_depth --- depth of the soil layer, m
# thermokarst_meth_prod 0. --- switch for old organics methane production under thermokarst lakes is OFF
# thermokarst_meth_prod 1. --- switch for old organics methane production under thermokarst lakes is ON
# soil_meth_prod 0. --- switch for new organics methane production under lakes is OFF
# soil_meth_prod 1. --- switch for new organics methane production under lakes is ON
# tricemethhydr 0. --- ice in soil pores is treated as purPBLpare ice
# tricemethhydr 1. --- ice in soil pores is treated as methane hydrate
# skin 0 --- the skin temperature parameterization is off
# skin 1 --- the skin temperature parameterisation is on
# sedim 0 --- gravitational sedimentation of tracer is NOT taken into account
# sedim 1 --- gravitational sedimentation of tracer is taken into account
# massflux 0 --- the massflux parameterization of convection (Siebesma et al., 2007) if OFF
# massflux 1 --- the massflux parameterization of convection (Siebesma et al., 2007) if ON
# sensflux0 --- sensible heat flux upwards, constant in time (relevant if PBLpar = 0), W/m**2
# momflux0 --- momentum flux downwards (positive), constant in time (relevant if PBLpar = 0), N/m**2
# ifrad 1 --- all radiation fluxes at the water surface are taken into account
# ifrad 0 --- all radiation fluxes are set to zero
# dyn_pgrad 0 --- dynamic pressure gradient is OFF
# dyn_pgrad 1 --- dynamic pressure gradient is ON
# zero_model 0 --- zero-dimensional model is ON
# zero_model 1 --- zero-dimensional model is OFF
# outflpar 0 --- variables value at the outflow = cross-section mean
# outflpar 1 --- the cross-section mean = 0.5*(inflow value + outflow value)
# outflpar 2 --- variables at the outflow are calculated using Lagrangian approach
# Note: zero-dimensional model is now implemented only for open water season and one-point simulation
# deadvol --- the depth (m) corresponding to "dead volume" - the minimal allowed reservoir volume
#----------------------------------------------------------------------------------------
#
varalb 0
waveenh 0
momflxpart 0
c_d -999.
kwe 100.
relwind 0
nmeltpoint 1
Turbpar 2
stabfunc 1
kepsbc 1
soiltype 5
soil_depth 10.
soilswitch 1
saltice 1
tricemethhydr 0.
skin 0
massflux 0
ifrad 1
ifbubble 1
sedim 0
salsoil 0
dyn_pgrad 0
zero_model 0
thermokarst_meth_prod 0.
soil_meth_prod 1.
outflpar 0 #2 #0
nManning 0.
horvisc 0.
pgrad 0.
backdiff0 0.
backdiff 0
botfric 1
#
VmaxCH4aeroboxid 0.
khsCH4 0.
khsO2 0.
r0methprod 0.
#
sensflux0 100.
momflux0 1.e-15
soilbotflx 0.
cuette 0
carbon_model 2
#
deadvol 0.
#
#----------------------------------------------------------------------------------------
# INITIAL CONDITIONS FOR TEMPERATURE
#----------------------------------------------------------------------------------------
#
T_profile 8
0. 15.23 0. 0.0012 0.080 2.94E-1 0.
1. 15.33 0. 0.0012 0.080 2.93E-1 0.
3. 15.29 0. 0.0013 0.080 2.63E-1 0.
5. 14.92 0. 0.0016 0.082 2.63E-1 0.
7. 11.88 0. 0.0010 0.144 2.45E-1 0.
9. 9.92 0. 0.0023 0.171 1.0E-1 0.
11. 9.12 0. 0.1076 0.175 0.0E-0 0.
12. 8.96 0. 0.1513 0.173 0.0E-0 0.
T_soilprofile -1
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
#.
#----------------------------------------------------------------------------------------
# BOUNDARY CONDITIONS: TRIBUTARIES AND EFFLUENTS
#----------------------------------------------------------------------------------------
# DESCRIPTION
# tribheat --- the switch for thermal effect of tributaries and effluents, 0 - OFF
# 1 - ON
# U_tribin --- the velocity of flow in tributary, m/s
# U_tribout --- the velocity of flow in effluent, m/s
# T_tribin --- the temperature of tributary flow, deg. C
# width_tribin --- the transversal width of tributary, m
# width_tribout --- the transversal width of effluent, m
# inflowprof --- the profile of inflow:
# width [m] velocity [m] temperature [C]
# outflowprof --- the profile of outflow:
# width [m] velocity [m]
#-----------------------------------------------------------------------------------------
#
#
tribheat 0
N_tribin 0
N_triblev 41
iefflloc 0
#fileinflow 'Kuj2012-13_inflow.dat'
#fileoutflow 'Kuj2012-13_outflow.dat'
dttribupdate 15.2
#inflowprof
#1 100. 0. 8.
#2 100. 0. 8.
#3 100. 0. 8.
#3 100. 0. 8.
#5 100. 0. 8.
#6 100. 0. 8.
#7 100. 0. 8.
#8 100. 0. 8.
#9 100. 0. 8.
#10 100. 0. 8.
#11 100. 0. 8.
#12 100. 0. 8.
#13 100. 0. 8.
#14 100. 0. 8.
#15 100. 0. 8.
#16 100. 0. 8.
#17 100. 0. 8.
#18 100. 0. 8.
#19 100. 0. 8.
#20 100. 0. 8.
#21 100. 0. 8.
#22 100. 0. 8.
#23 100. 0. 8.
#24 100. 0. 8.
#25 100. 0. 8.
#26 100. 0. 8.
#27 100. 0. 8.
#28 100. 0. 8.
#29 100. 0. 8.
#30 100. 0. 8.
#31 100. 0. 8.
#32 100. 0. 8.
#33 100. 0. 8.
#34 100. 0. 8.
#35 100. 0. 8.
#36 100. 0. 8.
#37 100. 0. 8.
#38 100. 0. 8.
#39 100. 0. 8.
#40 100. 0. 8.
#41 100. 0. 8.
#42 100. 0. 8.
#43 100. 0. 8.
#44 100. 0. 8.
#45 100. 0. 8.
#46 100. 0. 8.
#47 100. 0. 8.
#48 100. 0. 8.
#49 100. 0. 8.
#50 100. 0. 8.
#51 100. 0. 8.
#52 100. 0. 8.
#53 100. 0. 8.
#54 100. 0. 8.
#55 100. 0. 8.
#56 100. 0. 8.
#57 100. 0. 8.
#58 100. 0. 8.
#59 100. 0. 8.
#60 100. 0. 8.
#61 100. 0. 8.
#62 100. 0. 8.
#63 100. 0. 8.
#64 100. 0. 8.
#65 100. 0. 8.
#66 100. 0. 8.
#67 100. 0. 8.
#68 100. 0. 8.
#69 100. 0. 8.
#70 100. 0. 8.
#71 100. 0. 8.
#72 100. 0. 8.
#73 100. 0. 8.
#74 100. 0. 8.
#75 100. 0. 8.
#76 100. 0. 8.
#77 100. 0. 8.
#78 100. 0. 8.
#79 100. 0. 8.
#80 100. 0. 8.
#81 100. 0. 8.
#
#outflowprof
#1 100. 0.
#2 100. 0.
#3 100. 0.
#4 100. 0.
#5 100. 0.
#6 100. 0.
#7 100. 0.
#8 100. 0.
#9 100. 0.
#10 100. 0.
#11 100. 0.
#12 100. 0.
#13 100. 0.
#14 100. 0.
#15 100. 0.
#16 100. 0.
#17 100. 0.
#18 100. 0.
#19 100. 0.
#20 100. 0.
#21 100. 0.
#22 100. 0.
#23 100. 0.
#24 100. 0.
#25 100. 0.
#26 100. 0.
#27 100. 0.
#28 100. 0.
#29 100. 0.
#30 100. 0.
#31 100. 0.
#32 100. 0.
#33 100. 0.
#34 100. 0.
#35 100. 0.
#36 100. 0.
#37 100. 0.
#38 100. 0.
#39 100. 0.
#40 100. 0.
#41 100. 0.
#42 100. 0.
#43 100. 0.
#44 100. 0.
#45 100. 0.
#46 100. 0.
#47 100. 0.
#48 100. 0.
#49 100. 0.
#50 100. 0.
#51 100. 0.
#52 100. 0.
#53 100. 0.
#54 100. 0.
#55 100. 0.
#56 100. 0.
#57 100. 0.
#58 100. 0.
#59 100. 0.
#60 100. 0.
#61 100. 0.
#62 100. 0.
#63 100. 0.
#64 100. 0.
#65 100. 0.
#66 100. 0.
#67 100. 0.
#68 100. 0.
#69 100. 0.
#70 100. 0.
#71 100. 0.
#72 100. 0.
#73 100. 0.
#74 100. 0.
#75 100. 0.
#76 100. 0.
#77 100. 0.
#78 100. 0.
#79 100. 0.
#80 100. 0.
#81 100. 0.
#
#----------------------------------------------------------------------------------------
# DATA ASSIMILATION CONTROLS (NOT OPERATIONAL: PUT EVERYTHING TO 0)
#----------------------------------------------------------------------------------------
# assim --- data assimilation technique: 0 - no data assimilation
# 1 -
# 2 - Raleigh damping towards observations
# 3 - Cressman weighting
# 4 - ?
# as_window --- assimilation window: 1 - as_window is spinup period
#
#----------------------------------------------------------------------------------------
#
error_cov 0
assim 0
#as_window 1
#
#----------------------------------------------------------------------------------------
# OUTPUT CONTROLS (for ASCII files)
#----------------------------------------------------------------------------------------
# DESCRIPTION
# turb_out --- output turbulence characteristics: 1 - on, 0 - off
# monthly* --- monthly mean profiles output: 1 - on, 0 - off
# daily* --- daily mean profiles output: 1 - on, 0 - off
# hourly* --- hourly mean profiles output: 1 - on, 0 - off
# everystep* --- every time step profiles output: 1 - on, 2 - on but without profiles, 0 - off
# time_series* --- output of time series of layer thickness and surface values:
# 1 - on, 0 - off
# dt_out* --- time interval for time series output, hours
# nscreen --- the period of screen output, timesteps
# scale_output --- the switch for scaling of output of turbulent characteristics: 1 - on, 0 - off
# (set 0 for simulations of ice-covered lakes)
# ngrid_out --- the number of output levels for vertical water temperature profiles:
# -1 - use numerical grid levels
# >0 - use ngrid_out levels (in meters) given below
# ngridsoil_out--- the number of output levels for vertical soil temperature profiles:
# -1 - use numerical grid levels
# >0 - use ngridsoil_out levels (in meters) given below
#----------------------------------------------------------------------------------------
#
turb_out 0
monthly 1
daily 0
hourly 0
everystep 0
time_series 1
dt_out 3.
nscreen 1000
scale_output 0
accum_begin 2003060100
accum_end 2004060100
rtemp 1
-999. -999. -999.
zserout -999.
#
ngrid_out 12
0.2
0.5
1.
2.
3.
5.
7.
10.
13.
16.
19.
22.
#
ngridsoil_out 11
0.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
ngridice_out -1
#
#----------------------------------------------------------------------------------------
# NOTE: VARIABLES, DENOTED BY ASTERISK *, ARE USED ONLY IN STANDALONE RUNS OF THE MODEL
#----------------------------------------------------------------------------------------
end