IndexRundeck: P1SoM40Created: Wed May 8 03:30:18 EDT 2024

LAKES

File: LAKES.f
Summary: LAKES subroutines for Lakes and Rivers
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 |

Simplex Website Curator: Igor Aleinov — NASA Official: Gavin A. Schmidt

Contact GISS NASA Privacy PolicyAccessibility