In the ISPRAMIX ocean model (see Eifler and Schrimpf (1992)), another approach is used for considering stability effects on vertical mixing. The stability functions in this model are of the form:
The neutral Prandtl number used there is . The function is assumed to lay between the values 0.18 (corresponding to a supercritically stratified situation) and 2.0 (preventing it from growing too much under unstable conditions).
A formulation for can be derived from the definition of the flux Richardson number
and (202), see Beckers (1995):
where is the gradient Richardson number.
use turbulence, only: cm0_fix,Prandtl0_fix,xRF use turbulence, only: cmue1,cmue2,an,as IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: nlevREVISION HISTORY:
Original author(s): Manuel Ruiz Villarreal, Hans Burchard
Karsten Bolding 2012-12-28