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

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

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

revision 1.1 by adcroft, Wed May 30 19:34:48 2001 UTC revision 1.2 by jmc, Wed Jul 11 22:46:32 2001 UTC
# Line 9  C $Name$ Line 9  C $Name$
9       I           tracer,       I           tracer,
10       O           uT,       O           uT,
11       I           myThid )       I           myThid )
12    C     /==========================================================\
13    C     | SUBROUTINE GAD_FLUXLIMIT_ADV_X                           |
14    C     | o Compute Zonal advective Flux of Tracer using           |
15    C     |   Flux Limiter Scheme                                    |
16    C     |==========================================================|
17        IMPLICIT NONE        IMPLICIT NONE
18    
19  C     == GLobal variables ==  C     == GLobal variables ==
# Line 31  C     == Local variables == Line 36  C     == Local variables ==
36    
37        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
38         uT(1-Olx,j)=0.         uT(1-Olx,j)=0.
39         DO i=1-Olx+1,sNx+Olx         uT(2-Olx,j)=0.
40           uT(sNx+Olx,j)=0.
41           DO i=1-Olx+2,sNx+Olx-1
42          Rjp=(tracer(i+1,j)-tracer(i,j))*maskW(i+1,j,k,bi,bj)          Rjp=(tracer(i+1,j)-tracer(i,j))*maskW(i+1,j,k,bi,bj)
43          Rj=(tracer(i,j)-tracer(i-1,j))*maskW(i,j,k,bi,bj)          Rj=(tracer(i,j)-tracer(i-1,j))*maskW(i,j,k,bi,bj)
44          Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskW(i-1,j,k,bi,bj)          Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskW(i-1,j,k,bi,bj)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22