ip_song_wright

INTERFACE:

    subroutine ip_song_wright()
DESCRIPTION:

Here, the pressure gradient is calculating according to an energy-conserving method suggested by Song (1998), which for the pressure gradient in $ x$-direction looks as:

\begin{displaymath}\begin{array}{l} \displaystyle \frac12(h_{i,j,k}+h_{i,j,k+1})...
...z^c_{i+1,j,k}-z^c_{i,j,k}}{\Delta x^u_{i,j}}\Bigg], \end{array}\end{displaymath} (118)

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 omp_lib
    IMPLICIT NONE
REVISION HISTORY:
    Original author(s): Hans Burchard & Karsten Bolding
LOCAL VARIABLES:
    integer                   :: i,j,k
    REALTYPE                  :: dxm1,dym1
    REALTYPE                  :: grdl,grdu,buoyl,buoyu,prgr,dxz,dyz



kklingbe 2017-10-02