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

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

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

revision 1.4 by cnh, Tue Feb 6 04:22:36 2001 UTC revision 1.5 by adcroft, Wed Feb 7 03:00:00 2001 UTC
# Line 63  C     tracer cell face *above* this leve Line 63  C     tracer cell face *above* this leve
63            ENDDO            ENDDO
64          ENDDO          ENDDO
65        ELSEIF (k.eq.Nr) THEN        ELSEIF (k.eq.Nr) THEN
66  C       DO j=1-Oly,sNy+Oly-1          DO j=1-Oly,sNy+Oly-1
67  C         DO i=1-Olx,sNx+Olx-1            DO i=1-Olx,sNx+Olx-1
         DO j=1-1,sNy+1  
           DO i=1-1,sNx+1  
68              wFld(i,j,k,bi,bj) =              wFld(i,j,k,bi,bj) =
69       &         -(       &         -(
70       &            uTrans(i+1,j)-uTrans(i,j)       &            uTrans(i+1,j)-uTrans(i,j)
71       &           +vTrans(i,j+1)-vTrans(i,j)       &           +vTrans(i,j+1)-vTrans(i,j)
72       &          )/rA(i,j,bi,bj)       &          )*recip_rA(i,j,bi,bj)
73            ENDDO            ENDDO
74          ENDDO          ENDDO
75        ELSE        ELSE
76  C      DO j=1-Oly,sNy+Oly-1         DO j=1-Oly,sNy+Oly-1
77  C        DO i=1-Olx,sNx+Olx-1           DO i=1-Olx,sNx+Olx-1
        DO j=1-1,sNy+1  
          DO i=1-1,sNx+1  
78             wFld(i,j,k,bi,bj) = wFld(i,j,k+1,bi,bj)             wFld(i,j,k,bi,bj) = wFld(i,j,k+1,bi,bj)
79       &         -(       &         -(
80       &            uTrans(i+1,j)-uTrans(i,j)       &            uTrans(i+1,j)-uTrans(i,j)
81       &           +vTrans(i,j+1)-vTrans(i,j)       &           +vTrans(i,j+1)-vTrans(i,j)
82       &          )/rA(i,j,bi,bj)       &          )*recip_rA(i,j,bi,bj)
83            ENDDO            ENDDO
84          ENDDO          ENDDO
85        ENDIF        ENDIF
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22