Here, the interior diffusivities (defined as the diffusivities outside the surface and bottom boundary layers) are computed. The algorithms are identical to those suggested by Large et al. (1994). For numerical efficiency, the algorithms for different physical processes are active only if certain pre-processor macros are defined in cppdefs.h.
IMPLICIT NONEINPUT PARAMETERS:
number of grid cells integer :: nlev square of buoyancy frequency (1/s^2) REALTYPE :: NN(0:nlev) square of buoyancy frequencies caused by temperature and salinity stratification REALTYPE :: NNT(0:nlev),NNS(0:nlev) square of shear frequency (1/s^2) REALTYPE :: SS(0:nlev)REVISION HISTORY:
Original author(s): Lars Umlauf