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

Diff of /MITgcm/pkg/generic_advdiff/gad_fluxlimit_adv_r.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, Thu Jul 12 00:30:27 2001 UTC
# Line 9  C $Name$ Line 9  C $Name$
9       I           tracer,       I           tracer,
10       O           wT,       O           wT,
11       I           myThid )       I           myThid )
12    C     /==========================================================\
13    C     | SUBROUTINE GAD_FLUXLIMIT_ADV_R                           |
14    C     | o Compute vertical 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 45  C     == Local variables == Line 50  C     == Local variables ==
50           Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj))           Rjp=(tracer(i,j,kp1,bi,bj)-tracer(i,j,k,bi,bj))
51       &        *maskC(i,j,kp1,bi,bj)       &        *maskC(i,j,kp1,bi,bj)
52           Rj=(tracer(i,j,k,bi,bj)-tracer(i,j,kM1,bi,bj))           Rj=(tracer(i,j,k,bi,bj)-tracer(i,j,kM1,bi,bj))
      &        *maskC(i,j,k,bi,bj)  
53           Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,kM2,bi,bj))           Rjm=(tracer(i,j,km1,bi,bj)-tracer(i,j,kM2,bi,bj))
54       &        *maskC(i,j,km1,bi,bj)       &        *maskC(i,j,km2,bi,bj)
55           IF (Rj.NE.0.) THEN           IF (Rj.NE.0.) THEN
56            IF (rTrans(i,j).LT.0) THEN            IF (rTrans(i,j).LT.0) THEN
57              Cr=Rjm/Rj              Cr=Rjm/Rj
# Line 62  C     == Local variables == Line 66  C     == Local variables ==
66            ENDIF            ENDIF
67           ENDIF           ENDIF
68           Cr=Limiter(Cr)           Cr=Limiter(Cr)
69           wT(i,j) =           wT(i,j) = maskC(i,j,kM1,bi,bj)*(
70       &     rTrans(i,j)       &     rTrans(i,j)*
      &      *(maskC(i,j,kM1,bi,bj)*  
71       &        (Tracer(i,j,k,bi,bj)+Tracer(i,j,kM1,bi,bj))*0.5 _d 0       &        (Tracer(i,j,k,bi,bj)+Tracer(i,j,kM1,bi,bj))*0.5 _d 0
72       &       +(maskC(i,j,k,bi,bj)-maskC(i,j,kM1,bi,bj))*       &    +(ABS(rTrans(i,j))*(1-Cr)
73       &        Tracer(i,j,k,bi,bj)       &      +rTrans(i,j)*wVel(i,j,k,bi,bj)*deltaT*recip_drC(k)
74       &       )       &                  *Cr
75       &     +maskC(i,j,km1,bi,bj)*0.5*(       &     )*Rj*0.5 _d 0                )
      &        (1-Cr)*ABS(rTrans(i,j))  
      &        +rTrans(i,j)*wVel(i,j,k,bi,bj)*deltaT  
      &         *recip_drC(k)*Cr  
      &        )*Rj  
76          ENDDO          ENDDO
77         ENDDO         ENDDO
78        ENDIF        ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22