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

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

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

revision 1.2 by jmc, Fri Nov 4 01:31:36 2005 UTC revision 1.3 by heimbach, Wed Jun 7 01:55:14 2006 UTC
# Line 149  C--    Add centered & upwind contributio Line 149  C--    Add centered & upwind contributio
149             a5d(i,j,k)   = a5d(i,j,k)             a5d(i,j,k)   = a5d(i,j,k)
150       &                  + rC4km       &                  + rC4km
151       &                   *deltaTarg(k)       &                   *deltaTarg(k)
152       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
153             b5d(i,j,k)   = b5d(i,j,k)             b5d(i,j,k)   = b5d(i,j,k)
154       &                  - ( (rCenter+rUpwind) + rC4km )       &                  - ( (rCenter+rUpwind) + rC4km )
155       &                   *deltaTarg(k)       &                   *deltaTarg(k)
156       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
157             c5d(i,j,k)   = c5d(i,j,k)             c5d(i,j,k)   = c5d(i,j,k)
158       &                  - ( (rCenter-rUpwind) + rC4kp )       &                  - ( (rCenter-rUpwind) + rC4kp )
159       &                   *deltaTarg(k)       &                   *deltaTarg(k)
160       &                    *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                    *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
161             d5d(i,j,k)   = d5d(i,j,k)             d5d(i,j,k)   = d5d(i,j,k)
162       &                  + rC4kp       &                  + rC4kp
163       &                   *deltaTarg(k)       &                   *deltaTarg(k)
164       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
165             b5d(i,j,k-1) = b5d(i,j,k-1)             b5d(i,j,k-1) = b5d(i,j,k-1)
166       &                  - rC4km       &                  - rC4km
167       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
168       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
169             c5d(i,j,k-1) = c5d(i,j,k-1)             c5d(i,j,k-1) = c5d(i,j,k-1)
170       &                  + ( (rCenter+rUpwind) + rC4km )       &                  + ( (rCenter+rUpwind) + rC4km )
171       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
172       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
173             d5d(i,j,k-1) = d5d(i,j,k-1)             d5d(i,j,k-1) = d5d(i,j,k-1)
174       &                  + ( (rCenter-rUpwind) + rC4kp )       &                  + ( (rCenter-rUpwind) + rC4kp )
175       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
176       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
177             e5d(i,j,k-1) = e5d(i,j,k-1)             e5d(i,j,k-1) = e5d(i,j,k-1)
178       &                  - rC4kp       &                  - rC4kp
179       &                   *deltaTarg(k-1)       &                   *deltaTarg(k-1)
180       &                   *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                   *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
181           ENDDO           ENDDO
182         ENDDO         ENDDO
183    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22