ip_blumberg_mellor_lin

INTERFACE:

    subroutine ip_blumberg_mellor_lin()
DESCRIPTION:

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.

\begin{displaymath}\begin{array}{l} \displaystyle \frac12(h_{i,j,k}+h_{i,j,k+1})...
..._{i,j,k}}{z^c_{i,j,k+1}-z^c_{i,j,k}} \right)\Bigg], \end{array}\end{displaymath} (117)

where $ z^c_{i,j,k}$ is the $ z$-coordinate of the centre of the grid box with the index $ (i,j,k)$.

The discretisation of $ (\partial_y^* b)_k$ for the $ v$-equation is done accordingly. USES:

    use internal_pressure
    use variables_3d, only: kumin_pmz,kvmin_pmz
  $ use omp_lib
    IMPLICIT NONE
REVISION HISTORY:
    Original author(s): Hans Burchard & Karsten Bolding
LOCAL VARIABLES:
    integer                   :: i,j,k
    REALTYPE                  :: dxm1,dym1
    REALTYPE                  :: prgr,dxzu,dxzl,dyzu,dyzl
    REALTYPE                  :: dzr2,dzr1,dxru,dxrl,dyru,dyrl,aa,bb,cc



kklingbe 2017-10-02