/[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.7 by jmc, Sat Oct 22 20:17:44 2005 UTC revision 1.8 by heimbach, Wed Jun 7 01:55:14 2006 UTC
# Line 113  c           rUpwind= (2. _d 0*rHigh - wC Line 113  c           rUpwind= (2. _d 0*rHigh - wC
113             a5d(i,j,k)   = a5d(i,j,k)             a5d(i,j,k)   = a5d(i,j,k)
114       &                  + rC4km       &                  + rC4km
115       &                   *deltaTarg(k)       &                   *deltaTarg(k)
116       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
117             b5d(i,j,k)   = b5d(i,j,k)             b5d(i,j,k)   = b5d(i,j,k)
118       &                  - ( (rCenter+rUpwind) + rC4km )       &                  - ( (rCenter+rUpwind) + rC4km )
119       &                   *deltaTarg(k)       &                   *deltaTarg(k)
120       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
121             c5d(i,j,k)   = c5d(i,j,k)             c5d(i,j,k)   = c5d(i,j,k)
122       &                  - ( (rCenter-rUpwind) + rC4kp )       &                  - ( (rCenter-rUpwind) + rC4kp )
123       &                   *deltaTarg(k)       &                   *deltaTarg(k)
124       &                    *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                    *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
125             d5d(i,j,k)   = d5d(i,j,k)             d5d(i,j,k)   = d5d(i,j,k)
126       &                  + rC4kp       &                  + rC4kp
127       &                   *deltaTarg(k)       &                   *deltaTarg(k)
128       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
129             b5d(i,j,k-1) = b5d(i,j,k-1)             b5d(i,j,k-1) = b5d(i,j,k-1)
130       &                  - rC4km       &                  - rC4km
131       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
132       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
133             c5d(i,j,k-1) = c5d(i,j,k-1)             c5d(i,j,k-1) = c5d(i,j,k-1)
134       &                  + ( (rCenter+rUpwind) + rC4km )       &                  + ( (rCenter+rUpwind) + rC4km )
135       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
136       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
137             d5d(i,j,k-1) = d5d(i,j,k-1)             d5d(i,j,k-1) = d5d(i,j,k-1)
138       &                  + ( (rCenter-rUpwind) + rC4kp )       &                  + ( (rCenter-rUpwind) + rC4kp )
139       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
140       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
141             e5d(i,j,k-1) = e5d(i,j,k-1)             e5d(i,j,k-1) = e5d(i,j,k-1)
142       &                  - rC4kp       &                  - rC4kp
143       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
144       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
145           ENDDO           ENDDO
146         ENDDO         ENDDO
147    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22