Here, the internal part of the pressure gradient is discretised according to Mellor et al. (1994). The crucial part of this term, which is (in the case of the -equation), is discretised between two vertically adjacent velocity points:
where is the -coordinate of the interface in the T-point above the grid box with the index .
The discretisation of for the -equation is done accordingly.
In this routine, as a first step, the interface heights are calculated in the T-points, in order to allow for the calculation of the coordinate slopes in the U- and V-points. In a second step, the expression (116) equivalent formulation for the -direction are integrated up downwards, beginning from the surface. USES:
use internal_pressure $ use omp_lib IMPLICIT NONEREVISION HISTORY:
Original author(s): Hans Burchard, Adolf Stips, Karsten BoldingLOCAL VARIABLES:
integer :: i,j,k REALTYPE :: dxm1,dym1 REALTYPE :: grdl,grdu,buoyl,buoyu,prgr,dxz,dyz