Index | Rundeck: P1SoM40 | Created: Wed May 8 03:30:18 EDT 2024 |
LAKES | File: LAKES.f |
Author : Gavin Schmidt/Gary Russell
Version: 2010/08/04 (based on LB265); enhanced in June 2021 if the preprocessor option RVR_ELEV is selected: River speed is variable, emergency direction outflows are replaced by a more stable scheme, the River Direction file is now optional, but if one is used, select for the 144x90 grid RVR=RD2HX2E.nc, not RD_Fd.nc.
Subroutines:
- lkmix
- LKMIX calculates mixing and entrainment in lakes
- lksourc
- LKSOURC applies fluxes to lake in ice-covered and ice-free areas
Global Variables:
- grav : used from constant
- bygrav : used from constant
- shw : used from constant
- rhow : used from constant
- lhm : used from constant
- shi : used from constant
- teeny : used from constant
- undef : used from constant
- kdirec :
integer, allocatable, dimension(:, :)
- directions for river flow
Used by: | :alloc_lakes | :riverf |
- dhorz :
real*8, allocatable, dimension(:, :)
- horizontal distance to downstream box (m)
Used by: | :alloc_lakes |
- xyzc :
real*8, allocatable, dimension(:, :, :)
- (x,y,z) unit vectors on sphere of primary cell centers
Used by: | :alloc_lakes |
- jflow :
integer, allocatable, dimension(:, :)
- grid box indexes for downstream direction
Used by: | :alloc_lakes | :riverf |
- iflow :
integer, allocatable, dimension(:, :)
- grid box indexes for downstream direction
Used by: | :alloc_lakes | :riverf |
- jfl911 :
integer, allocatable, dimension(:, :)
- grid box indexes for emergency downstream direction
Used by: | :alloc_lakes | :riverf |
- ifl911 :
integer, allocatable, dimension(:, :)
- grid box indexes for emergency downstream direction
Used by: | :alloc_lakes | :riverf |
- kd911 :
integer, allocatable, dimension(:, :)
- emergency directions for river flow
Used by: | :alloc_lakes | :riverf |
- rate :
real*8, allocatable, dimension(:, :)
- rate of river flow downslope (fraction)
Used by: | :alloc_lakes | :riverf |
- minmld :
real*8, parameter
- minimum mixed layer depth in lake (m)
Initial Value= 1.
Used by: | :daily_lake | :ground_lk | :surface |
- hlake_min :
real*8, parameter
- minimum sill depth for lake (m)
Initial Value= 1.
Used by: | :daily_lake |
- tmaxrho :
real*8, parameter
- temperature of maximum density (pure water) (C)
Initial Value= 4.
- kvlake :
real*8, parameter
- lake diffusion constant at mixed layer depth (m^2/s)
Initial Value= 1d-5
- tfl :
real*8, parameter
- freezing temperature for lakes (=0 C)
Initial Value= 0.
- ac2lmin :
real*8, parameter
- minimum ice thickness for lake ice (kg/m^2)
Initial Value= 0.1
- ac1lmin :
real*8, parameter
- minimum ice thickness for lake ice (kg/m^2)
Initial Value= 0.1
- fleadlk :
real*8, parameter
- lead fraction for lakes
Initial Value= 0.
- byzeta :
real*8, parameter
- reciprocal of solar rad. extinction depth for lake (1/m)
Initial Value= 1./0.35d0
Used by: | :ground_lk |
- river_fac :
real*8
- Factor to multiply runoff by to balance sea level
Initial Value= 1.
Used by: | :riverf |
- init_flake :
integer
- used to make sure FLAKE is properly initialised
Initial Value= 1
- variable_lk :
integer
- 1 if lakes are to be variable
Initial Value= 0
Used by: | :daily_lake |
- lake_rise_max :
real*8
- amount of lake rise (m) over sill level before
Initial Value= 1d2
Used by: | :riverf |