/[MITgcm]/MITgcm/pkg/smooth/smooth_rhs.F
ViewVC logotype

Diff of /MITgcm/pkg/smooth/smooth_rhs.F

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

revision 1.1 by gforget, Mon Feb 15 23:46:04 2010 UTC revision 1.2 by gforget, Tue Aug 21 13:46:30 2012 UTC
# Line 98  c 2nd k loop: flux computation Line 98  c 2nd k loop: flux computation
98         DO i=1-OLx,sNx+OLx         DO i=1-OLx,sNx+OLx
99          df(i,j,bi,bj)        = 0. _d 0          df(i,j,bi,bj)        = 0. _d 0
100          xA(i,j,bi,bj) = _dyG(i,j,bi,bj)          xA(i,j,bi,bj) = _dyG(i,j,bi,bj)
101       &   *drF(k)*_hFacW(i,j,k,bi,bj)       &   *drF(k)*smooth_hFacW(i,j,k,bi,bj)
102          yA(i,j,bi,bj) = _dxG(i,j,bi,bj)          yA(i,j,bi,bj) = _dxG(i,j,bi,bj)
103       &   *drF(k)*_hFacS(i,j,k,bi,bj)       &   *drF(k)*smooth_hFacS(i,j,k,bi,bj)
104        IF (K .EQ. 1) THEN        IF (K .EQ. 1) THEN
105            maskUp(i,j,bi,bj) = 0.            maskUp(i,j,bi,bj) = 0.
106        ELSE        ELSE
# Line 344  c 3rd k loop: Divergence of fluxes Line 344  c 3rd k loop: Divergence of fluxes
344           DO j=jMin,jMax           DO j=jMin,jMax
345            DO i=iMin,iMax            DO i=iMin,iMax
346          gt_in(i,j,k,bi,bj)=gt_in(i,j,k,bi,bj)          gt_in(i,j,k,bi,bj)=gt_in(i,j,k,bi,bj)
347       &   -_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)*recip_rA(i,j,bi,bj)       &   -smooth_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
348         &   *recip_rA(i,j,bi,bj)
349       &   *( (fZon(i+1,j,k,bi,bj)-fZon(i,j,k,bi,bj))       &   *( (fZon(i+1,j,k,bi,bj)-fZon(i,j,k,bi,bj))
350       &     +(fMer(i,j+1,k,bi,bj)-fMer(i,j,k,bi,bj))       &     +(fMer(i,j+1,k,bi,bj)-fMer(i,j,k,bi,bj))
351       &     +(fVerT(i,j,k,bi,bj)-fVerT(i,j,k-1,bi,bj))*rkSign       &     +(fVerT(i,j,k,bi,bj)-fVerT(i,j,k-1,bi,bj))*rkSign
# Line 355  c 3rd k loop: Divergence of fluxes Line 356  c 3rd k loop: Divergence of fluxes
356           DO j=jMin,jMax           DO j=jMin,jMax
357            DO i=iMin,iMax            DO i=iMin,iMax
358          gt_in(i,j,k,bi,bj)=gt_in(i,j,k,bi,bj)          gt_in(i,j,k,bi,bj)=gt_in(i,j,k,bi,bj)
359       &   -_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)*recip_rA(i,j,bi,bj)       &   -smooth_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
360         &   *recip_rA(i,j,bi,bj)
361       &   *( (fZon(i+1,j,k,bi,bj)-fZon(i,j,k,bi,bj))       &   *( (fZon(i+1,j,k,bi,bj)-fZon(i,j,k,bi,bj))
362       &     +(fMer(i,j+1,k,bi,bj)-fMer(i,j,k,bi,bj))       &     +(fMer(i,j+1,k,bi,bj)-fMer(i,j,k,bi,bj))
363       &     +(fVerT(i,j,k,bi,bj))*rkSign       &     +(fVerT(i,j,k,bi,bj))*rkSign
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22