!&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
    SUBROUTINE PDTE
!     ******************************************************************
!$$$  SUBPROGRAM DOCUMENTATION BLOCK
!                .      .    .
! SUBPROGRAM:    PDTE        SURFACE PRESSURE TENDENCY CALC
!   PRGRMMR: JANJIC          ORG: W/NP22     DATE: 94-03-08

! ABSTRACT:
!     PDTE VERTICALLY INTEGRATES THE MASS FLUX DIVERGENCE TO
!     OBTAIN THE SURFACE PRESSURE TENDENCY AND ETADOT ON THE
!     LAYER INTERFACES.

! PROGRAM HISTORY LOG:
!   87-06-??  JANJIC     - ORIGINATOR
!   95-03-25  BLACK      - CONVERSION FROM 1-D TO 2-D IN HORIZONTAL
!   95-11-20  ABELES     - PARALLEL OPTIMIZATION
!   96-03-29  BLACK      - REMOVED SCRCH
!   98-10-30  BLACK      - MODIFIED FOR DISTRIBUTED MEMORY

! USAGE: CALL PDTE FROM MAIN PROGRAM EBU
!   INPUT ARGUMENT LIST:
!       NONE

!   OUTPUT ARGUMENT LIST:
!     NONE

!   OUTPUT FILES:
!     NONE

!   SUBPROGRAMS CALLED:

!     UNIQUE: NONE

!     LIBRARY: NONE

!   COMMON BLOCKS: CTLBLK
!                  MASKS
!                  DYNAM
!                  VRBLS

! ATTRIBUTES:
!   LANGUAGE: FORTRAN 90
!   MACHINE : IBM SP
!$$$
!***********************************************************************
    PARAMETER &
    (CP=1004.6)
!-----------------------------------------------------------------------
    INCLUDE "parmeta.f90"
    INCLUDE "mpp.h"
#include "sp.h"
!-----------------------------------------------------------------------
    PARAMETER &
    (IMJM=IM*JM-JM/2,JAM=6+2*(JM-10),LM1=LM-1,LP1=LM+1)
!-----------------------------------------------------------------------
    LOGICAL :: &
    RUN,FIRST,RESTRT,SIGMA
!-----------------------------------------------------------------------
    INCLUDE "COMM_CTLBLK.f90"
!-----------------------------------------------------------------------
    INCLUDE "COMM_MASKS.f90"
!-----------------------------------------------------------------------
    INCLUDE "COMM_DYNAM.f90"
!-----------------------------------------------------------------------
    INCLUDE "COMM_VRBLS.f90"
    INCLUDE "COMM_PVRBLS.f90"
    INCLUDE "COMM_CLDWTR.f90"
!-----------------------------------------------------------------------
    INCLUDE "COMM_CONTIN.f90"
!-----------------------------------------------------------------------
    DIMENSION &
    PRET  (IDIM1:IDIM2,JDIM1:JDIM2),RPSL  (IDIM1:IDIM2,JDIM1:JDIM2)
!    &,TT    (IDIM1:IDIM2,JDIM1:JDIM2),CPEV  (IDIM1:IDIM2,JDIM1:JDIM2)
!    &,DPDE  (IDIM1:IDIM2,JDIM1:JDIM2)
!-----------------------------------------------------------------------
!-----------------------------------------------------------------------
!--------------COMPUTATION OF PRESSURE TENDENCY & PREPARATIONS----------
    DO 100 L=2,LM
    ! omp parallel do
        DO 100 J=MYJS,MYJE2
            DO 100 I=MYIS,MYIE
                DIV(I,J,L)=DIV(I,J,L-1)+DIV(I,J,L)
    100 END DO

! omp parallel do
    DO 110 J=MYJS2,MYJE2
        DO 110 I=MYIS,MYIE
            PSDT(I,J)=-DIV(I,J,LM)
            PRET(I,J)=PSDT(I,J)*RES(I,J)
            RPSL(I,J)=1./PDSL(I,J)
    110 END DO
!--------------COMPUTATION OF ETADT-------------------------------------
! omp parallel do
    DO 120 L=1,LM1
        DO 120 J=MYJS2,MYJE2
            DO 120 I=MYIS,MYIE
                ETADT(I,J,L)=-(PRET(I,J)*ETA(L+1)+DIV(I,J,L)) &
                *HTM(I,J,L+1)*RPSL(I,J)
    120 END DO
!--------------KINETIC ENERGY GENERATION TERMS IN T EQUATION------------
! omp parallel do
    DO 130 J=MYJS2,MYJE2
        DO 130 I=MYIS,MYIE
        !     DPDE(I,J)=DETA(1)*PDSL(I,J)
        !     TT(I,J)=T(I,J,1)
            OMGALF(I,J,1)=OMGALF(I,J,1)-DIV(I,J,1)*RTOP(I,J,1)*EF4T
            T(I,J,1)=T(I,J,1)-DIV(I,J,1)*RTOP(I,J,1)*EF4T
        !     CPEV(I,J)=CP*(T(I,J,1)-TT(I,J))*DPDE(I,J)
    130 END DO
!-----------------------------------------------------------------------
! omp parallel do
    DO 140 J=MYJS2,MYJE2
        DO 140 I=MYIS,MYIE
        !     DPDE(I,J)=DETA(L)*PDSL(I,J)
        !     TT(I,J)=T(I,J,L)
            OMGALF(I,J,L)=OMGALF(I,J,L)-(DIV(I,J,L-1)+DIV(I,J,L)) &
            *RTOP(I,J,L)*EF4T &
            *HTM(I,J,L)
            T(I,J,L)=T(I,J,L)-(DIV(I,J,L-1)+DIV(I,J,L))*RTOP(I,J,L)*EF4T
        !     CPEV(I,J)=CP*(T(I,J,L)-TT(I,J))*DPDE(I,J)
    140 END DO
    145 CONTINUE
!-----------------------------------------------------------------------
! omp parallel do
    DO 150 J=MYJS2,MYJE2
        DO 150 I=MYIS,MYIE
        !     DPDE(I,J)=DETA(LM)*PDSL(I,J)
        !     TT(I,J)=T(I,J,LM)
            OMGALF(I,J,LM)=OMGALF(I,J,LM)+(PRET(I,J)-DIV(I,J,LM1)) &
            *RTOP(I,J,LM)*EF4T &
            *HTM(I,J,LM)
            T(I,J,LM)=T(I,J,LM)+(PRET(I,J)-DIV(I,J,LM1))*RTOP(I,J,LM)*EF4T
        !     CPEV(I,J)=CP*(T(I,J,LM)-TT(I,J))*DPDE(I,J)
    150 END DO

    DO 160 L=LM,2,-1
    ! omp parallel do
        DO 160 J=MYJS,MYJE2
            DO 160 I=MYIS,MYIE
                DIV(I,J,L)=DIV(I,J,L)-DIV(I,J,L-1)
    160 END DO
!-----------------------------------------------------------------------
    RETURN
    END SUBROUTINE PDTE