... | ... | @@ -109,8 +109,7 @@ by at least one space. Some keywords are followed by a group of numbers. |
|
|
|*area_lake* | Lake area, m^2 |
|
|
|
|*cellipt* | The ratio of x dimension to y dimension of a lake, n/d |
|
|
|
|*lakeform* | The form of the lake's horizontal cross-section: 1 -- ellipse, 2 -- rectangle |
|
|
|
|*trib_inflow* | Total tributaries' inflow, m^3/s; if -9999., then
|
|
|
it is dynamically adjusted to keep the lake depth nearly constant |
|
|
|
|*trib_inflow* | Total tributaries' inflow, m^3/s; if -9999., then it is dynamically adjusted to keep the lake depth nearly constant |
|
|
|
|*morphometry* | Depth - lake cross-section area table, m^2; if not given, the area is constant with depth |
|
|
|
|*effl_outflow* | effluent discharge parameters group, specifying polynomial dependence of discharge on water level; the last value is relative altitude of effluent bottom over lake bottom (at deepest points, respectively); the first N values are polynomial coefficients, where N stands after *effl_outflow* keyword |
|
|
|
|
... | ... | @@ -160,6 +159,7 @@ Setup file has the same format, as the driver file. The list of parameters is gi |
|
|
|
|
|
| Keyword | Description and acceptable values set |
|
|
|
|:---------------|:-----------------------------------------------------------|
|
|
|
|*nstep_keps* | Number of tim esteps of k-epsilon parameterization per on model timestep, recommended value is 1 |
|
|
|
|*M* | Number of layers in water |
|
|
|
|*ns* | Number of layers in soil |
|
|
|
|*nsoilcols* | Number of soil columns (1 - single soil column contacting with the lowest numerical layer (of water or deep ice) only, \> 1 - additional soil columns at the bottom slope, contacting with all numerical layers of water, ice and deep ice; morphometry information is necessary in driver file, if *nsoilcols* \> 1) |
|
... | ... | @@ -181,7 +181,6 @@ Setup file has the same format, as the driver file. The list of parameters is gi |
|
|
|*relwind* | In momentum flux calculation the currents speed is substracted from wind speed: 1 - on, 0 - off |
|
|
|
|*eos* | Equation of state: 1 - from Hostetler model (freshwater), 2 - TEOS-2010 (currently not operational), 3 - for Kivu lake (Schmid et al., 2003), 4 - UNESCO formula (freshwater), 5 - McCatcheon et al., 1993 |
|
|
|
|*lindens* | Switch for using linearized equation of state respective to temperature and salinity: 1 - on, 0 - off |
|
|
|
|*lindens* | Linear dependency of density on temperature and salinity: 0 - off, 1 - on |
|
|
|
|*nmeltpoint* | Melting point: 0 - constant, 0 °C, 1 - linearly dependent on salinity, 2 - TEOS-2010 formula |
|
|
|
|*Turbpar* | Turbulence model: 2 - $` K-\epsilon `$ closure; 8 - Hendersson-Sellers diffusivity (1985); 9 - Henderson-Sellers diffusivity adopted for shallow lakes |
|
|
|
|*stabfunc* | Stability function in $` K-\epsilon `$ : 1 - constant coefficients (standard K-\epsilon model); 2 - stability functions according to (Canuto et al., 2001); 3 - stability functions according to (Galperin et al., 1988) |
|
... | ... | @@ -199,6 +198,12 @@ Setup file has the same format, as the driver file. The list of parameters is gi |
|
|
|*sedim* | Gravitational sedimentation of tracer in water: 0 - off, 1 - on |
|
|
|
|*salsoil* | Switch for salinity transport in soil: 0 - off (zero flux at the bottom), 1 - on |
|
|
|
|*dyn_pgrad* | Dynamic pressure gradient: 0 - off, 1 - barotropic pressure gradient, 2 - two-layer approximation for pressure gradient, 3 - multilayer model for pressure gradient |
|
|
|
|*botfric* | The law for friction at the sloping bottom: 0 - no bottom friction, 1 - linear friction, 2 - quadratic friction |
|
|
|
|*horvisc* | Horizontal viscosity, m^2/s, used to damp the 1-st horizontal ve-
|
|
|
locity mode |
|
|
|
|*backdiff* | Switch for background diffusivity by Hondzo and Stefan (1993); 1
|
|
|
– off, 1 – on |
|
|
|
|*nManning* | Manning’s coefficient, relevant for dyn_pgrad =5 |
|
|
|
|*zero_model* | Switch for zero-dimensional model: 0 -off, 1 - on |
|
|
|
|
|
|
|
*thermokarst_meth_prod* | Switch for old organics methane production under thermokarst lakes: 1. - on, 0. - off |
|
... | ... | @@ -214,7 +219,7 @@ Setup file has the same format, as the driver file. The list of parameters is gi |
|
|
|
|
|
| Keyword | Description and acceptable values set |
|
|
|
|:---------------|:-----------------------------------------------------------|
|
|
|
|*T_profile* | Initial temperature (°C, col.2), salinity (kg/kg, col.3), CH_4 (mol/m^3, col.4), CO_2 (mol/m^3, col.5), O_2 (mol/m^3, col.6) profiles, at depth levels (m, col.1) |
|
|
|
|*T_profile* | Initial temperature (°C, col.2), salinity (kg/kg, col.3), CH_4 (mol/m^3, col.4), CO_2 (mol/m^3, col.5), O_2 (mol/m^3, col.6), DIP (mol/m^3, col.7) profiles, at depth levels (m, col.1) |
|
|
|
|
|
|
**Tributaries and effluents**
|
|
|
|
... | ... | @@ -224,6 +229,10 @@ Setup file has the same format, as the driver file. The list of parameters is gi |
|
|
|*N_tribin* | The number of tributaries, the next line should be a sequence of *N_tribin* numbers, from 1 to 4, describing the location of each tributary |
|
|
|
|*N_triblev* | The number of levels in tributaries data, relevant if tribheat = 1 or 2 |
|
|
|
|*iefflloc* | from 1 to 4, the location of an effluent |
|
|
|
|*fileinflow* | The name of file in data folder with time-dependent inflow profiles. Relevant if tribheat = 2,3 |
|
|
|
|*fileoutflow* | The name of file in data folder with time-dependent outflow profiles. Relevant if tribheat = 2,3 |
|
|
|
|*dttribupdate* | Time step for profiles update in fileinflow and fileoutflow, days. Relevant if tribheat = 2,3 |
|
|
|
|
|
|
|
|
|
## Running the model
|
|
|
The model is launched by typing in the terminal
|
... | ... | |