#include "dims.h"
      SUBROUTINE GLPL(ZL,ZP)
      use cons_module,only: len1,imax,imaxp2,grav,re_inv
      implicit none
#include "params.h"
#include "fcom.h"
#include "index.h"
#include "buff.h"
#include "phys.h"
#include "vscr.h"
!
! Args:
      real,intent(out) :: ZL(1),ZP(1)
!
! Local:
      integer :: i
!
      CALL DLDP(NZ,ZL,ZP,LEN1-4)
      DO 1 I = 3,LEN1-2
        ZL(I) = ZL(I)*grav*RACS
        ZP(I) = ZP(I)*grav*re_inv
    1 CONTINUE
C     ****
C     ****     PERIODIC POINTS
C     ****
      DO 2 I = 1,2
        ZL(I) = ZL(I+IMAX)
        ZL(I+IMAXP2) = ZL(I+2)
        ZP(I) = ZP(I+IMAX)
        ZP(I+IMAXP2) = ZP(I+2)
    2 CONTINUE
      RETURN
      END
C