Here, the pressure gradient is calculating according to an energy-conserving method suggested by Song (1998), which for the pressure gradient in -direction looks as:
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 omp_lib IMPLICIT NONEREVISION HISTORY:
Original author(s): Hans Burchard & Karsten BoldingLOCAL VARIABLES:
integer :: i,j,k REALTYPE :: dxm1,dym1 REALTYPE :: grdl,grdu,buoyl,buoyu,prgr,dxz,dyz