Compute interior fluxes

INTERFACE:

    subroutine interior(nlev,NN,NNT,NNS,SS)

DESCRIPTION:

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.

• The shear instability algorithm is active if the macro KPP_SHEAR is defined.
• The internal wave algorithm is active if the macro KPP_INTERNAL_WAVE is defined.
• The convective instability algorithm is active if the macro KPP_CONVEC is defined.
• The double-diffusion algorithm is active if the macro KPP_DDMIX is defined. Note that in this case, the macro SALINITY has to be defined as well.

USES:

    IMPLICIT NONE

INPUT 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