/[MITgcm]/MITgcm/model/src/calc_gw.F
ViewVC logotype

Diff of /MITgcm/model/src/calc_gw.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by adcroft, Mon Mar 22 15:54:03 1999 UTC revision 1.2 by adcroft, Tue Mar 14 17:47:25 2000 UTC
# Line 177  caja           and an hFacUS (above V po Line 177  caja           and an hFacUS (above V po
177         ENDDO         ENDDO
178        ENDDO        ENDDO
179    
 C     Add div(W*) RHS of elliptic equation  
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
         DO K=1,1  
          DO j=1,sNy  
           DO i=1,sNx  
            IF ( _hFacC(i,j,k,bi,bj) .GT. 0. ) THEN  
             cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) +  
      &       (  
      &        -wVel(i,j,k+1,bi,bj)  
      &         )*_rA(i,j,bi,bj)/deltatMom  
            ENDIF  
           ENDDO  
          ENDDO  
         ENDDO  
         DO K=2,Nr-1  
          DO j=1,sNy  
           DO i=1,sNx  
            IF ( _hFacC(i,j,k,bi,bj) .GT. 0. ) THEN  
             cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) +  
      &       ( wVel(i,j,k  ,bi,bj)  
      &        -wVel(i,j,k+1,bi,bj)  
      &       )*_rA(i,j,bi,bj)/deltatMom  
            ENDIF  
           ENDDO  
          ENDDO  
         ENDDO  
         DO K=Nr,Nr  
          DO j=1,sNy  
           DO i=1,sNx  
            IF ( _hFacC(i,j,k,bi,bj) .GT. 0. ) THEN  
             cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) +  
      &       wVel(i,j,k  ,bi,bj)  
      &       *_rA(i,j,bi,bj)/deltatMom  
            ENDIF  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
         
180  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
181    
182        RETURN        RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22