C **********************************************************************
C *                                                                    *
C *                                SOFTWARE LICENSING                  *
C *                                                                    *
C *                This program is free software; you can redistribute *
C *                it and/or modify it under the terms of the GNU      *
C *                General Public License as published by the Free     *
C *                Software Foundation, either Version 2 of the        *
C *                license, or (at your option) any later version.     *
C *                                                                    *
C *                This program is distributed in the hope that it     *
C *                will be useful, but without any warranty; without   *
C *                even the implied warranty of merchantability or     *
C *                fitness for a particular purpose. See the GNU       *
C *                General Public License for more details.            *
C *                                                                    *
C *                A copy of the GNU General Public License is         *
C *                available at http://www.gnu.org/copyleft/gpl.html   *
C *                or by writing to the Free Software Foundation, Inc.,*
C *                59 Temple Place - Suite 330, Boston, MA 02111, USA. *
C *                                                                    *
C **********************************************************************

      SUBROUTINE BAROPG(DRHOX,DRHOY,TRNU,TRNV,SGW)
C     VERSION(09/27/90)
C
      INCLUDE 'comdeck'
      SAVE
      DIMENSION DRHOX(IM,JM,KB),DRHOY(IM,JM,KB),TRNU(IM,JM),TRNV(IM,JM)
C
CC      DO 50 K=1,KB
CC      DO 50 J=1,JM
CC      DO 50 I=1,IM
CC50    RHO(I,J,K)=RHO(I,J,K)-RMEAN(I,J,K)
C
c Introduce SGW's scheme to get density (NKIM 12/17/02)
c
c      DATA SGW/0.5/
c
      DO K=1,KB
         DO J=1,JM
            DO I=1,IM
               RHO(I,J,K)=SGW*RHO(I,J,K)+(1-SGW)*RMEAN(I,J,K)
            ENDDO
         ENDDO
      ENDDO
cend
C-------- X COMPONENT OF BAROCLINIC PRESSURE GRADIENT ------------------
      DO 300 J=2,JMM1
      DO 300 I=2,IMM1
 300  DRHOX(I,J,1)=.5*GRAV*(-ZZ(1))*(DT(I,J)+DT(I-1,J))
     .    *(RHO(I,J,1)-RHO(I-1,J,1))
      DO 310 J=2,JMM1
      DO 310 K=2,KBM1
      DO 310 I=2,IMM1
 310  DRHOX(I,J,K)=DRHOX(I,J,K-1)
     .      +GRAV*.25*(ZZ(K-1)-ZZ(K))*(DT(I,J)+DT(I-1,J))
     .      *(RHO(I,J,K)-RHO(I-1,J,K)+RHO(I,J,K-1)-RHO(I-1,J,K-1))
     .      +GRAV*.25*(ZZ(K-1)+ZZ(K))*(DT(I,J)-DT(I-1,J))
     .      *(RHO(I,J,K)+RHO(I-1,J,K)-RHO(I,J,K-1)-RHO(I-1,J,K-1))
C
      DO 360 J=2,JMM1
      DO 360 K=1,KBM1
      DO 360 I=2,IMM1
 360  DRHOX(I,J,K)=.25*(DT(I,J)+DT(I-1,J))*DRHOX(I,J,K)*DUM(I,J)
     .     *DUMWAD(I,J)*(H2(I,J)+H2(I-1,J)) ! Not in original NYHOPS W&D
!     .     *(H2(I,J)+H2(I-1,J))
C
C-------- Y COMPONENT OF BAROCLINIC PRESSURE GRADIENT ------------------
      DO 500 J=3,JMM1
      DO 500 I=2,IMM1
 500  DRHOY(I,J,1)=.5*GRAV*(-ZZ(1))*(DT(I,J)+DT(I,J-1))
     .                *(RHO(I,J,1)-RHO(I,J-1,1))
      DO 510 J=3,JMM1
      DO 510 K=2,KBM1
      DO 510 I=2,IMM1
 510  DRHOY(I,J,K)=DRHOY(I,J,K-1)
     .      +GRAV*.25*(ZZ(K-1)-ZZ(K))*(DT(I,J)+DT(I,J-1))
     .        *(RHO(I,J,K)-RHO(I,J-1,K)+RHO(I,J,K-1)-RHO(I,J-1,K-1))
     .      +GRAV*.25*(ZZ(K-1)+ZZ(K))*(DT(I,J)-DT(I,J-1))
     .        *(RHO(I,J,K)+RHO(I,J-1,K)-RHO(I,J,K-1)-RHO(I,J-1,K-1))
C
      DO 560 J=3,JMM1
      DO 560 K=1,KBM1
      DO 560 I=2,IMM1
 560  DRHOY(I,J,K)=.25*(DT(I,J)+DT(I,J-1))*DRHOY(I,J,K)*DVM(I,J)
     .       *DVMWAD(I,J)*(H1(I,J)+H1(I,J-1)) ! Not in original NYHOPS W&D
!     .       *(H1(I,J)+H1(I,J-1))
      DO 580 J=1,JM
      DO 580 K=1,KB
      DO 580 I=1,IM
      DRHOX(I,J,K)=RAMP*DRHOX(I,J,K)*DUM(I,J)*DUMWAD(I,J) ! Not in original NYHOPS W&D
 580  DRHOY(I,J,K)=RAMP*DRHOY(I,J,K)*DVM(I,J)*DVMWAD(I,J) ! Not in original NYHOPS W&D
C
C-------- VERTICALLY INTEGRATE PRESSURE GRADIENT -----------------------
      DO 590 J=1,JM
      DO 590 I=1,IM
      TRNU(I,J)=0.0
 590  TRNV(I,J)=0.0
      DO 570 J=1,JM
      DO 570 K=1,KBM1
      DO 570 I=1,IM
      TRNU(I,J)=TRNU(I,J)+DRHOX(I,J,K)*DZ(K)
 570  TRNV(I,J)=TRNV(I,J)+DRHOY(I,J,K)*DZ(K)
C
CC      DO 600 K=1,KB
CC      DO 600 J=1,JM
CC      DO 600 I=1,IM
CC 600  RHO(I,J,K)=RHO(I,J,K)+RMEAN(I,J,K)
CCC
      RETURN
      END