/[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.7 by cnh, Thu May 28 15:09:29 1998 UTC revision 1.8 by cnh, Sat May 30 02:10:16 1998 UTC
# 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
# 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 197  C     Note: For K=1 then KM1=1 this give Line 197  C     Note: For K=1 then KM1=1 this give
197  C           boundary condition.  C           boundary condition.
198        DO j=jMin,jMax        DO j=jMin,jMax
199         DO i=iMin,iMax         DO i=iMin,iMax
200          df(i,j) = zA(i,j,bi,bj)*(          df(i,j) = _zA(i,j,bi,bj)*(
201       &   -(diffKzT+KapGM(i,j)*K33(i,j,k))*rdzC(k)       &   -(diffKzT+KapGM(i,j)*K33(i,j,k))*rdzC(k)
202       &   *(theta(i,j,kM1,bi,bj)-theta(i,j,k,bi,bj))       &   *(theta(i,j,kM1,bi,bj)-theta(i,j,k,bi,bj))
203       &   -KapGM(i,j)*K13(i,j,k)*dTdx(i,j)       &   -KapGM(i,j)*K13(i,j,k)*dTdx(i,j)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22