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

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

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

revision 1.14 by cnh, Sat Aug 22 17:51:07 1998 UTC revision 1.15 by cnh, Wed Oct 28 03:11:36 1998 UTC
# Line 161  C--   Interpolate terms for Redi/GM sche Line 161  C--   Interpolate terms for Redi/GM sche
161        DO j=jMin,jMax        DO j=jMin,jMax
162         DO i=iMin,iMax         DO i=iMin,iMax
163          dSdx(i,j) = 0.5*(          dSdx(i,j) = 0.5*(
164       &   +0.5*(_maskW(i+1,j,k,bi,bj)*_recip_dxC(i+1,j,bi,bj)*       &   +0.5*(_maskW(i+1,j,k,bi,bj)
165         &         *_recip_dxC(i+1,j,bi,bj)*
166       &           (salt(i+1,j,k,bi,bj)-salt(i,j,k,bi,bj))       &           (salt(i+1,j,k,bi,bj)-salt(i,j,k,bi,bj))
167       &        +_maskW(i,j,k,bi,bj)*_recip_dxC(i,j,bi,bj)*       &        +_maskW(i,j,k,bi,bj)
168         &         *_recip_dxC(i,j,bi,bj)*
169       &           (salt(i,j,k,bi,bj)-salt(i-1,j,k,bi,bj)))       &           (salt(i,j,k,bi,bj)-salt(i-1,j,k,bi,bj)))
170       &   +0.5*(_maskW(i+1,j,km1,bi,bj)*_recip_dxC(i+1,j,bi,bj)*       &   +0.5*(_maskW(i+1,j,km1,bi,bj)
171         &         *_recip_dxC(i+1,j,bi,bj)*
172       &           (salt(i+1,j,km1,bi,bj)-salt(i,j,km1,bi,bj))       &           (salt(i+1,j,km1,bi,bj)-salt(i,j,km1,bi,bj))
173       &        +_maskW(i,j,km1,bi,bj)*_recip_dxC(i,j,bi,bj)*       &        +_maskW(i,j,km1,bi,bj)
174         &         *_recip_dxC(i,j,bi,bj)*
175       &           (salt(i,j,km1,bi,bj)-salt(i-1,j,km1,bi,bj)))       &           (salt(i,j,km1,bi,bj)-salt(i-1,j,km1,bi,bj)))
176       &       )       &       )
177         ENDDO         ENDDO
# Line 175  C--   Interpolate terms for Redi/GM sche Line 179  C--   Interpolate terms for Redi/GM sche
179        DO j=jMin,jMax        DO j=jMin,jMax
180         DO i=iMin,iMax         DO i=iMin,iMax
181          dSdy(i,j) = 0.5*(          dSdy(i,j) = 0.5*(
182       &   +0.5*(_maskS(i,j,k,bi,bj)*_recip_dyC(i,j,bi,bj)*       &   +0.5*(_maskS(i,j,k,bi,bj)
183         &         *_recip_dyC(i,j,bi,bj)*
184       &           (salt(i,j,k,bi,bj)-salt(i,j-1,k,bi,bj))       &           (salt(i,j,k,bi,bj)-salt(i,j-1,k,bi,bj))
185       &        +_maskS(i,j+1,k,bi,bj)*_recip_dyC(i,j+1,bi,bj)*       &        +_maskS(i,j+1,k,bi,bj)
186         &         *_recip_dyC(i,j+1,bi,bj)*
187       &           (salt(i,j+1,k,bi,bj)-salt(i,j,k,bi,bj)))       &           (salt(i,j+1,k,bi,bj)-salt(i,j,k,bi,bj)))
188       &   +0.5*(_maskS(i,j,km1,bi,bj)*_recip_dyC(i,j,bi,bj)*       &   +0.5*(_maskS(i,j,km1,bi,bj)
189         &         *_recip_dyC(i,j,bi,bj)*
190       &           (salt(i,j,km1,bi,bj)-salt(i,j-1,km1,bi,bj))       &           (salt(i,j,km1,bi,bj)-salt(i,j-1,km1,bi,bj))
191       &        +_maskS(i,j+1,km1,bi,bj)*_recip_dyC(i,j+1,bi,bj)*       &        +_maskS(i,j+1,km1,bi,bj)
192         &         *_recip_dyC(i,j+1,bi,bj)*
193       &           (salt(i,j+1,km1,bi,bj)-salt(i,j,km1,bi,bj)))       &           (salt(i,j+1,km1,bi,bj)-salt(i,j,km1,bi,bj)))
194       &       )       &       )
195         ENDDO         ENDDO
# Line 240  C           they are not algorithmically Line 248  C           they are not algorithmically
248  C           are not used.  C           are not used.
249        DO j=jMin,jMax        DO j=jMin,jMax
250         DO i=iMin,iMax         DO i=iMin,iMax
251  #define _recip_VolS(i,j,k,bi,bj) _recip_hFacC(i,j,k,bi,bj)*recip_drF(k)/_rA(i,j,bi,bj)  #define _recip_VolS1(i,j,k,bi,bj) _recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
252    #define _recip_VolS2(i,j,k,bi,bj) /_rA(i,j,bi,bj)
253          gS(i,j,k,bi,bj)=          gS(i,j,k,bi,bj)=
254       &   -_recip_VolS(i,j,k,bi,bj)       &   -_recip_VolS1(i,j,k,bi,bj)
255         &    _recip_VolS2(i,j,k,bi,bj)
256       &   *(       &   *(
257       &    +( fZon(i+1,j)-fZon(i,j) )       &    +( fZon(i+1,j)-fZon(i,j) )
258       &    +( fMer(i,j+1)-fMer(i,j) )       &    +( fMer(i,j+1)-fMer(i,j) )
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22