Index | Rundeck: E6F40 | Created: Thu May 16 03:30:27 EDT 2024 |
ATURB.f
Summary: Author :
Version:
Modules:
Global Subroutines:
- apply_fluxes_to_atm
- a dummy subroutine that replaces the real one needed by DRYCNV.
- atm_diffus
- atm_diffus updates u,v,t,q due to turbulent transport throughout
all GCM layers using a non-local turbulence model
- de_solver_edge
- differential equation solver using tridiagonal method.
The differential equation is expressed as
d/dt x = d/dz (P1 d/dz x) - P3 x + P4
where x is the unknown to be solved,
x, P3 and P4 are at the layer edge ze, while
P1 is at the layer middle z.
- de_solver_main
- differential equation solver using tridiagonal method.
The differential equation is expressed as
d/dt x = d/dz (P1 d/dz x) + P4
where x is the unknown to be solved,
x and P4 are at the layer middle z, while
P1 is at the layer edge ze.
- dout
- dout writes out diagnostics at (i,j)
- e_gcm
- finds the turbulent kinetic energy (e, in m^2/s^2).
Within the PBL, e is determined according to the
parameterization of the Large Eddy Simulation (LES) data
(Moeng and Sullivan, 1994), above the PBL, e is calculated
by the second order closure model of Cheng et al. (2002).
- find_pbl_top
- finds the PBL height (dbl, in meters)
and the main level index immediately above (ldbl),
using the bulk Richardson number criterion
(Holtslag and Boville, 1993).
- getdz
- getdz computes the 3-d finite difference dz and dze
as well as the 3-d density rho and rhoe
called at the primary grid (A-grid)
- k_gcm
- computes the turbulent stability functions Km (for momentum)
and Kh (for heat and moisture), as well as the fluxes
(local and non-local).
Within the convective PBL, it is according to
Holtslag and Boville (1993); within the stable PBL or above
the PBL, it is according to Cheng et al. (2002).
- l_gcm
- calculates the turbulent length scale
(lscale, in meters). Within the PBL, it is according
to Nakanishi(2001); above the PBL, we generalized and
employed a formula by Holtslag and Boville (1993).
- zze
- finds the layer middle and edge heights, z and ze
Note that z(L) is between ze(L) and ze(L+1).
Depends on the following files:
Used by the following files: