/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_u3c4_impl_r.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_u3c4_impl_r.F

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

revision 1.4 by jmc, Sat Dec 4 18:50:34 2004 UTC revision 1.5 by jmc, Wed Jun 22 00:27:47 2005 UTC
# Line 88  C--   process interior interface only: Line 88  C--   process interior interface only:
88  C--    Add centered & upwind contributions  C--    Add centered & upwind contributions
89         DO j=jMin,jMax         DO j=jMin,jMax
90           DO i=iMin,iMax           DO i=iMin,iMax
91             rCenter= 0.5 _d 0 *rTrans(i,j)*recip_rA(i,j,bi,bj)*rkFac             rCenter= 0.5 _d 0 *rTrans(i,j)*recip_rA(i,j,bi,bj)*rkSign
92             mskM   = maskC(i,j,km2,bi,bj)*maskM2             mskM   = maskC(i,j,km2,bi,bj)*maskM2
93             mskP   = maskC(i,j,kp1,bi,bj)*maskP1             mskP   = maskC(i,j,kp1,bi,bj)*maskP1
94             rC4km  = oneSixth*rCenter*mskM             rC4km  = oneSixth*rCenter*mskM
# Line 98  C--    Add centered & upwind contributio Line 98  C--    Add centered & upwind contributio
98              rU3km  = 0. _d 0              rU3km  = 0. _d 0
99              rU3kp  = 0. _d 0              rU3kp  = 0. _d 0
100             ELSE             ELSE
101              rU1k   = oneSixth*abs(rCenter)              rU1k   = oneSixth*ABS(rCenter)
102              rUpwind= rU1k+rU1k              rUpwind= rU1k+rU1k
103              rU3km  = rU1k*mskM              rU3km  = rU1k*mskM
104              rU3kp  = rU1k*mskP              rU3kp  = rU1k*mskP
105             ENDIF             ENDIF
106             a5d(i,j,k)   = a5d(i,j,k)             a5d(i,j,k)   = a5d(i,j,k)
107       &                  - (rC4km - rU3km)       &                  + (rC4km + rU3km)
108       &                   *deltaTarg(k)       &                   *deltaTarg(k)
109       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
110             b5d(i,j,k)   = b5d(i,j,k)             b5d(i,j,k)   = b5d(i,j,k)
111       &                  + (rCenter + rC4km - rUpwind - rU3km)       &                  - (rCenter + rC4km + rUpwind + rU3km)
112       &                   *deltaTarg(k)       &                   *deltaTarg(k)
113       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
114             c5d(i,j,k)   = c5d(i,j,k)             c5d(i,j,k)   = c5d(i,j,k)
115       &                  + (rCenter + rC4kp + rUpwind + rU3kp)       &                  - (rCenter + rC4kp - rUpwind - rU3kp)
116       &                   *deltaTarg(k)       &                   *deltaTarg(k)
117       &                    *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                    *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
118             d5d(i,j,k)   = d5d(i,j,k)             d5d(i,j,k)   = d5d(i,j,k)
119       &                  - (rC4kp + rU3kp)       &                  + (rC4kp - rU3kp)
120       &                   *deltaTarg(k)       &                   *deltaTarg(k)
121       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
122             b5d(i,j,k-1) = b5d(i,j,k-1)             b5d(i,j,k-1) = b5d(i,j,k-1)
123       &                  + (rC4km - rU3km)       &                  - (rC4km + rU3km)
124       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
125       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
126             c5d(i,j,k-1) = c5d(i,j,k-1)             c5d(i,j,k-1) = c5d(i,j,k-1)
127       &                  - (rCenter + rC4km - rUpwind - rU3km)       &                  + (rCenter + rC4km + rUpwind + rU3km)
128       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
129       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
130             d5d(i,j,k-1) = d5d(i,j,k-1)             d5d(i,j,k-1) = d5d(i,j,k-1)
131       &                  - (rCenter + rC4kp + rUpwind + rU3kp)       &                  + (rCenter + rC4kp - rUpwind - rU3kp)
132       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
133       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
134             e5d(i,j,k-1) = e5d(i,j,k-1)             e5d(i,j,k-1) = e5d(i,j,k-1)
135       &                  + (rC4kp + rU3kp)       &                  - (rC4kp - rU3kp)
136       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
137       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
138           ENDDO           ENDDO
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22