Here, the internal pressure gradient calculation is carried out on the basis of the same buoyancy stencil than in the method according to Mellor et al. (1994) (see routine ip_blumberg_mellor), but in such a way that the pressure gradient numerically vanishes for linear stratification without horizontal gradients.
where is the -coordinate of the centre of the grid box with the index .
The discretisation of for the -equation is done accordingly. USES:
use internal_pressure use variables_3d, only: kumin_pmz,kvmin_pmz $ use omp_lib IMPLICIT NONEREVISION HISTORY:
Original author(s): Hans Burchard & Karsten BoldingLOCAL VARIABLES:
integer :: i,j,k REALTYPE :: dxm1,dym1 REALTYPE :: prgr,dxzu,dxzl,dyzu,dyzl REALTYPE :: dzr2,dzr1,dxru,dxrl,dyru,dyrl,aa,bb,cc