Newer
Older
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
#---------------------------------------------------------------------------------------
# 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
M 20
ns 10
Mice 10
d_surf 1.E-2
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 pure 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 1
PBLpar 1
waveenh 0
momflxpart 1
c_d -999
kwe 100.
relwind 0
eos 5
lindens 0
nmeltpoint 1
Turbpar 2
stabfunc 2
kepsbc 1
soiltype 5
soil_depth 10.
soilswitch 1
saltice 0
nsoilcols 5
tricemethhydr 0.
carbon_model 1
skin 0
massflux 0
ifrad 1
ifbubble 1
sedim 0
salsoil 0
dyn_pgrad 0
pgrad 0.
nManning 5.E-2
horvisc 0.
Victor Stepanenko
committed
backdiff 2
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
backdiff0 -999.
botfric 1
zero_model 0
thermokarst_meth_prod 0.
soil_meth_prod 1.
outflpar 0
#
sensflux0 100.
momflux0 1.e-15
soilbotflx 0.
cuette 0
#
VmaxCH4aeroboxid -999.
khsCH4 -999.
khsO2 -999.
r0methprod -999.
#
deadvol 0.
#
#----------------------------------------------------------------------------------------
# INITIAL CONDITIONS FOR TEMPERATURE
#----------------------------------------------------------------------------------------
#
#T_profile 15
#0.1528798623 15.6593492318 1.0967626609E-3 0. 0. 0.
#0.4586201021 15.3698204159 1.1244681301E-3 0. 0. 0.
#0.7643487883 15.2880392733 1.1389363045E-3 0. 0. 0.
#1.0700665454 15.3784681078 1.2157976331E-3 0. 0. 0.
#1.3757244583 15.4705468887 1.6825468273E-3 0. 0. 0.
#1.6811083613 15.3023162147 3.5392949352E-3 0. 0. 0.
#1.9859062169 14.8535425863 6.5826901114E-3 0. 0. 0.
#2.2900303837 14.5124893839 9.1725234989E-3 0. 0. 0.
#2.5934888392 13.6183860236 12.0747625176E-3 0. 0. 0.
#2.8960785319 13.2161902403 16.431655897E-3 0. 0. 0.
#3.1976916336 12.5317665948 20.3271090163E-3 0. 0. 0.
#3.4987146359 11.6504498125 21.2290960955E-3 0. 0. 0.
#3.7994820663 11.5276435987 22.3359795778E-3 0. 0. 0.
#4.0946210218 11.6325191677 23.7094310675E-3 0. 0. 0.
#7. 5. 30.E-3 0. 0. 0.
T_profile 8
0.5 10.59 1.0E-3 0. 0. 0. 0.
1.0 10.56 1.0E-3 0. 0. 0. 0.
1.5 10.56 1.0E-3 0. 0. 0. 0.
2.0 10.57 1.0E-3 0. 0. 0. 0.
2.5 10.56 1.0E-3 0. 0. 0. 0.
3.0 10.54 1.0E-3 0. 0. 0. 0.
3.5 10.49 1.0E-3 0. 0. 0. 0.
4.0 10.47 1.0E-3 0. 0. 0. 0.
#
T_soilprofile 2
0. 10.
7. 10.
#.
#----------------------------------------------------------------------------------------
# BOUNDARY CONDITIONS: TRIBUTARIES AND EFFLUENTS
#----------------------------------------------------------------------------------------
# DESCRIPTION
# tribheat --- the switch for thermal effect of tributaries and effluents, 0 - OFF
# 1 - ON
# 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 0
iefflloc 1
fileinflow 'BolshoiVilui20153_inflows.dat'
fileoutflow 'BolshoiVilui20153_outflow.dat'
dttribupdate 0.25
#
#----------------------------------------------------------------------------------------
# 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 0
daily 1
hourly 1
everystep 0
time_series 1
dt_out 1
nscreen 1000
scale_output 0
accum_begin 2003060100
accum_end 2004060100
zserout -999.
rtemp 1
-999. -999. -999.
#
ngrid_out 8
0.5
1.0
1.5
2.0
2.5
3.0
3.5
3.95
#
ngridice_out 2
0.
0.2
#
ngridsoil_out 8
0.
1.
2.
3.
4.
5.
6.
7.
#
#----------------------------------------------------------------------------------------
# NOTE: VARIABLES, DENOTED BY ASTERISK *, ARE USED ONLY IN STANDALONE RUNS OF THE MODEL
#----------------------------------------------------------------------------------------
end