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

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

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

revision 1.5 by jmc, Sat Dec 4 18:50:34 2004 UTC revision 1.6 by jmc, Wed Jun 22 00:27:47 2005 UTC
# Line 94  C--   Compute the upwind fraction: Line 94  C--   Compute the upwind fraction:
94              Cr=Rjp/Rj              Cr=Rjp/Rj
95            ENDIF            ENDIF
96            upwindFac(i,j) = 1. _d 0            upwindFac(i,j) = 1. _d 0
97       &                   - Limiter(Cr) * ( 1. _d 0 + abs(w_CFL) )       &                   - Limiter(Cr) * ( 1. _d 0 + ABS(w_CFL) )
98            upwindFac(i,j) = max( -1. _d 0, upwindFac(i,j) )            upwindFac(i,j) = MAX( -1. _d 0, upwindFac(i,j) )
99           ELSE           ELSE
100            upwindFac(i,j) = 0. _d 0            upwindFac(i,j) = 0. _d 0
101           ENDIF           ENDIF
# Line 105  C--   Compute the upwind fraction: Line 105  C--   Compute the upwind fraction:
105  C--    Add centered & upwind contributions  C--    Add centered & upwind contributions
106         DO j=jMin,jMax         DO j=jMin,jMax
107           DO i=iMin,iMax           DO i=iMin,iMax
108             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
109             rUpwind = abs(rCenter)*upwindFac(i,j)             rUpwind = ABS(rCenter)*upwindFac(i,j)
110             a3d(i,j,k)   = a3d(i,j,k)             a3d(i,j,k)   = a3d(i,j,k)
111       &                  + (rCenter-rUpwind)*deltaTarg(k)       &                  - (rCenter+rUpwind)*deltaTarg(k)
112       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                   *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
113             b3d(i,j,k)   = b3d(i,j,k)             b3d(i,j,k)   = b3d(i,j,k)
114       &                  + (rCenter+rUpwind)*deltaTarg(k)       &                  - (rCenter-rUpwind)*deltaTarg(k)
115       &                    *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)       &                    *recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
116             b3d(i,j,k-1) = b3d(i,j,k-1)             b3d(i,j,k-1) = b3d(i,j,k-1)
117       &                  - (rCenter-rUpwind)*deltaTarg(k-1)       &                  + (rCenter+rUpwind)*deltaTarg(k-1)
118       &                    *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                    *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
119             c3d(i,j,k-1) = c3d(i,j,k-1)             c3d(i,j,k-1) = c3d(i,j,k-1)
120       &                  - (rCenter+rUpwind)*deltaTarg(k-1)       &                  + (rCenter-rUpwind)*deltaTarg(k-1)
121       &                    *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)       &                    *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1)
122           ENDDO           ENDDO
123         ENDDO         ENDDO
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22