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

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

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

revision 1.5 by cnh, Wed May 27 21:01:47 1998 UTC revision 1.6 by cnh, Thu May 28 03:34:52 1998 UTC
# Line 134  C       Advective component of meridiona Line 134  C       Advective component of meridiona
134  C     Zonal tracer gradient  C     Zonal tracer gradient
135        DO j=jMin,jMax        DO j=jMin,jMax
136         DO i=iMin,iMax         DO i=iMin,iMax
137          dTdy(i,j) = rdyC(i,j,bi,bj)*          dTdy(i,j) = _rdyC(i,j,bi,bj)*
138       &  (theta(i,j,k,bi,bj)-theta(i,j-1,k,bi,bj))       &  (theta(i,j,k,bi,bj)-theta(i,j-1,k,bi,bj))
139         ENDDO         ENDDO
140        ENDDO        ENDDO
# Line 170  C--   Interpolate terms for Redi/GM sche Line 170  C--   Interpolate terms for Redi/GM sche
170        DO j=jMin,jMax        DO j=jMin,jMax
171         DO i=iMin,iMax         DO i=iMin,iMax
172          dTdy(i,j) = 0.5*(          dTdy(i,j) = 0.5*(
173       &   +0.5*(maskS(i,j,k,bi,bj)*rdyC(i,j,bi,bj)*       &   +0.5*(maskS(i,j,k,bi,bj)*_rdyC(i,j,bi,bj)*
174       &           (theta(i,j,k,bi,bj)-theta(i,j-1,k,bi,bj))       &           (theta(i,j,k,bi,bj)-theta(i,j-1,k,bi,bj))
175       &        +maskS(i,j+1,k,bi,bj)*rdyC(i,j+1,bi,bj)*       &        +maskS(i,j+1,k,bi,bj)*_rdyC(i,j+1,bi,bj)*
176       &           (theta(i,j+1,k,bi,bj)-theta(i,j,k,bi,bj)))       &           (theta(i,j+1,k,bi,bj)-theta(i,j,k,bi,bj)))
177       &   +0.5*(maskS(i,j,km1,bi,bj)*rdyC(i,j,bi,bj)*       &   +0.5*(maskS(i,j,km1,bi,bj)*_rdyC(i,j,bi,bj)*
178       &           (theta(i,j,km1,bi,bj)-theta(i,j-1,km1,bi,bj))       &           (theta(i,j,km1,bi,bj)-theta(i,j-1,km1,bi,bj))
179       &        +maskS(i,j+1,km1,bi,bj)*rdyC(i,j+1,bi,bj)*       &        +maskS(i,j+1,km1,bi,bj)*_rdyC(i,j+1,bi,bj)*
180       &           (theta(i,j+1,km1,bi,bj)-theta(i,j,km1,bi,bj)))       &           (theta(i,j+1,km1,bi,bj)-theta(i,j,km1,bi,bj)))
181       &       )       &       )
182         ENDDO         ENDDO
# Line 221  C           are not used. Line 221  C           are not used.
221        DO j=jMin,jMax        DO j=jMin,jMax
222         DO i=iMin,iMax         DO i=iMin,iMax
223          gT(i,j,k,bi,bj)=          gT(i,j,k,bi,bj)=
224       &   -rHFacC(i,j,k,bi,bj)*rdzF(k)*rDxF(i,j,bi,bj)*rDyF(i,j,bi,bj)       &   -rHFacC(i,j,k,bi,bj)*rdzF(k)*_rdxF(i,j,bi,bj)*_rdyF(i,j,bi,bj)
225       &   *(       &   *(
226       &    +( fZon(i+1,j)-fZon(i,j) )       &    +( fZon(i+1,j)-fZon(i,j) )
227       &    +( fMer(i,j+1)-fMer(i,j) )       &    +( fMer(i,j+1)-fMer(i,j) )
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22