C
C
      SUBROUTINE GETW(DX,DY,Z,DZ,H,U,V,W,IM,JM,KB)
C
C     Calc. vertical vel from U,V (interp. from w perpendicular to
C     sigma levels to real upward direction.)
C
      DIMENSION Z(KB),DZ(KB)
      DIMENSION DX(IM,JM),DY(IM,JM),H(IM,JM)
      DIMENSION U(IM,JM,KB),V(IM,JM,KB),W(IM,JM,KB)
      DIMENSION XFLUX(IM,JM,KB),YFLUX(IM,JM,KB)
C
      XFLUX=0.
      YFLUX=0.
C
C
C CALCULATE VERTICAL VELOCITY (RELATIVE TO SIGMA)
C
      DO 100 K=1,KB-1
      DO 100 J=1,JM
      DO 100 I=2,IM
 100  XFLUX(I,J,K)
     1 =.250*(DY(I,J)+DY(I-1,J))*(H(I,J)+H(I-1,J))*U(I,J,K)
      DO 120 K=1,KB-1
      DO 120 J=2,JM
      DO 120 I=1,IM
  120 YFLUX(I,J,K)
     1  =.250*(DX(I,J)+DX(I,J-1))*(H(I,J)+H(I,J-1))*V(I,J,K)
      W(:,:,1)=0.0
      DO 710 K=1,KB-1
      DO 710 J=2,JM-1
      DO 710 I=2,IM-1
 710  W(I,J,K+1)=W(I,J,K)
     1    +DZ(K)* (XFLUX(I+1,J,K)-XFLUX(I,J,K)
     2            +YFLUX(I,J+1,K)-YFLUX(I,J,K))/(DX(I,J)*DY(I,J))
C
C --------- CALC. REAL (UPWARD) VERTICAL VELOCITY ----------------
C (assuming EL<<H and D(EL)/Dt=0)
C
      DO 810 K=1,KB-1
      DO 810 J=2,JM-1
      DO 810 I=2,IM-1
 810  W(I,J,K+1)=W(I,J,K+1) + Z(K+1) * (
     1  (U(I,J,K)+U(I,J,K+1))*(H(I,J)-H(I-1,J))/(DX(I,J)+DX(I-1,J))
     2 +(V(I,J,K)+V(I,J,K+1))*(H(I,J)-H(I,J-1))/(DY(I,J)+DY(I,J-1)) )
C
      RETURN
      END
C