ip_z_interpol

INTERFACE:

    subroutine ip_z_interpol()
DESCRIPTION:

Here, the horizontal gradients of buoyancy, $(\partial_x^* b)_k$ and $(\partial_y^* b)_k$, are directly calculated in $z$-coordinates by linearly interpolating the buoyancies in the vertical to the evaluation point (which is the interface vertically located between the velocity points). In the case that extrapolations become necessary near the sloping surface (or more likely) near the sloping bottom, then the last regular buoyancy value (surface value or bottom value) is used. USES:

    use internal_pressure
  $ use omp_lib
    IMPLICIT NONE
REVISION HISTORY:
    Original author(s): Hans Burchard & Karsten Bolding
LOCAL VARIABLES:
    integer                   :: i,j,k, rc
    REALTYPE                  :: dxm1,dym1
    REALTYPE                  :: grdl,grdu,buoyl,prgr,dxz,dyz
    integer                   :: kplus,kminus
    REALTYPE, POINTER         :: zx(:)
    REALTYPE                  :: buoyplus,buoyminus