This subroutine updates the explicit solution of (72) and (73) with shape indicated by (76). In addition to the simplifications discussed in section 4.25, is assumed in (73) to eliminate the dependency on according to (74). As discussed in section 4.4, this implies that the last of (73) is replaced by (75). Thus, the -term in (76) drops out, and the solution is characterized by and only.
As a consequence, the numerators and the denominator appearing in (79) are of somewhat different form compared to the result in section 4.24. They can be written as
The coefficients of are given by
These polynomials correspond to a slightly generalized form of the solution suggested by Canuto et al. (2001) and Cheng et al. (2002). For cases with unstable stratification, the same clipping conditions on is applied as described in section 4.27. For the cases of extreme shear, the limiter described in the context of (85) is active.
use turbulence, only: an,as,at use turbulence, only: cmue1,cmue2 use turbulence, only: cm0 use turbulence, only: cc1 use turbulence, only: ct1,ctt use turbulence, only: a1,a2,a3,a4,a5 use turbulence, only: at1,at2,at3,at4,at5 IMPLICIT NONEINPUT PARAMETERS:
number of vertical layers integer, intent(in) :: nlevDEFINED PARAMETERS:
REALTYPE, parameter :: asLimitFact=1.0d0 REALTYPE, parameter :: anLimitFact=0.5d0REVISION HISTORY:
Original author(s): Lars Umlauf
Karsten Bolding 2012-12-28