#---------------------------------------------------------------------------------------
# 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.
backdiff 2
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