/[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.4 by cnh, Mon May 25 20:05:55 1998 UTC revision 1.5 by cnh, Wed May 27 21:01:47 1998 UTC
# Line 104  C     Advective component of zonal flux Line 104  C     Advective component of zonal flux
104  C     Zonal tracer gradient  C     Zonal tracer gradient
105        DO j=jMin,jMax        DO j=jMin,jMax
106         DO i=iMin,iMax         DO i=iMin,iMax
107          dTdx(i,j) = rdxC(i,j,bi,bj)*          dTdx(i,j) = _rdxC(i,j,bi,bj)*
108       &  (theta(i,j,k,bi,bj)-theta(i-1,j,k,bi,bj))       &  (theta(i,j,k,bi,bj)-theta(i-1,j,k,bi,bj))
109         ENDDO         ENDDO
110        ENDDO        ENDDO
# Line 156  C--   Interpolate terms for Redi/GM sche Line 156  C--   Interpolate terms for Redi/GM sche
156        DO j=jMin,jMax        DO j=jMin,jMax
157         DO i=iMin,iMax         DO i=iMin,iMax
158          dTdx(i,j) = 0.5*(          dTdx(i,j) = 0.5*(
159       &   +0.5*(maskW(i+1,j,k,bi,bj)*rdxC(i+1,j,bi,bj)*       &   +0.5*(maskW(i+1,j,k,bi,bj)*_rdxC(i+1,j,bi,bj)*
160       &           (theta(i+1,j,k,bi,bj)-theta(i,j,k,bi,bj))       &           (theta(i+1,j,k,bi,bj)-theta(i,j,k,bi,bj))
161       &        +maskW(i,j,k,bi,bj)*rdxC(i,j,bi,bj)*       &        +maskW(i,j,k,bi,bj)*_rdxC(i,j,bi,bj)*
162       &           (theta(i,j,k,bi,bj)-theta(i-1,j,k,bi,bj)))       &           (theta(i,j,k,bi,bj)-theta(i-1,j,k,bi,bj)))
163       &   +0.5*(maskW(i+1,j,km1,bi,bj)*rdxC(i+1,j,bi,bj)*       &   +0.5*(maskW(i+1,j,km1,bi,bj)*_rdxC(i+1,j,bi,bj)*
164       &           (theta(i+1,j,km1,bi,bj)-theta(i,j,km1,bi,bj))       &           (theta(i+1,j,km1,bi,bj)-theta(i,j,km1,bi,bj))
165       &        +maskW(i,j,km1,bi,bj)*rdxC(i,j,bi,bj)*       &        +maskW(i,j,km1,bi,bj)*_rdxC(i,j,bi,bj)*
166       &           (theta(i,j,km1,bi,bj)-theta(i-1,j,km1,bi,bj)))       &           (theta(i,j,km1,bi,bj)-theta(i-1,j,km1,bi,bj)))
167       &       )       &       )
168         ENDDO         ENDDO
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22