#---------------------------------------------------------------------------------------
# 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 2.E+0
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
PBLpar 3
waveenh 0
momflxpart 0
c_d -999.
kwe 100.
relwind 0
eos 5
lindens 0
nmeltpoint 1
Turbpar 2
stabfunc 1
kepsbc 1
soiltype 5
soil_depth 10.
soilswitch 1
saltice 1
nsoilcols 5
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
#.
#----------------------------------------------------------------------------------------
# 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