INTERFACE:
subroutine ip_z_interpol()DESCRIPTION:
Here, the horizontal gradients of buoyancy, and , are directly calculated in -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 NONEREVISION HISTORY:
Original author(s): Hans Burchard & Karsten BoldingLOCAL 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