Compute interior fluxes

*INTERFACE:*

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

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

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)

Original author(s): Lars Umlauf

Karsten Bolding 2012-12-28