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

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

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

revision 1.1 by jmc, Thu Jul 12 00:31:59 2001 UTC revision 1.2 by jmc, Mon Aug 20 20:44:01 2001 UTC
# Line 13  C     /================================= Line 13  C     /=================================
13  C     | SUBROUTINE GAD_U3_ADV_Y                                  |  C     | SUBROUTINE GAD_U3_ADV_Y                                  |
14  C     | o Compute Meridional advective Flux of Tracer using      |  C     | o Compute Meridional advective Flux of Tracer using      |
15  C     |   3rd Order Upwind Scheme                                |  C     |   3rd Order Upwind Scheme                                |
16    C     | o reduced to 2nd or 3rd Order near a boundary            |
17  C     |==========================================================|  C     |==========================================================|
18        IMPLICIT NONE        IMPLICIT NONE
19    
# Line 40  C     == Local variables == Line 41  C     == Local variables ==
41        ENDDO        ENDDO
42        DO j=1-Oly+2,sNy+Oly-1        DO j=1-Oly+2,sNy+Oly-1
43         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
44          Rjp=(tracer(i,j+1)-tracer(i,j))*maskS(i,j+1,k,bi,bj)          Rjp=(tracer(i,j+1)-tracer(i,j))
45          Rj =(tracer(i,j)-tracer(i,j-1))*maskS(i,j,k,bi,bj)  c    &     *maskS(i,j+1,k,bi,bj)
46          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskS(i,j-1,k,bi,bj)          Rj =(tracer(i,j)-tracer(i,j-1))
47          Rjjp=Rjp-Rj  c    &     *maskS(i,j,k,bi,bj)
48          Rjjm=Rj-Rjm          Rjm=(tracer(i,j-1)-tracer(i,j-2))
49    c    &     *maskS(i,j-1,k,bi,bj)
50    C- jmc: mask not needed here above if Rjjp & Rjjm & vTrans are masked
51            Rjjp=(Rjp-Rj)*maskS(i,j+1,k,bi,bj)
52            Rjjm=(Rj-Rjm)*maskS(i,j-1,k,bi,bj)
53          vT(i,j) =          vT(i,j) =
54       &   vTrans(i,j)*(       &   vTrans(i,j)*(
55       &     Tracer(i,j)+Tracer(i,j-1)-oneSixth*( Rjjp+Rjjm )       &     Tracer(i,j)+Tracer(i,j-1)-oneSixth*( Rjjp+Rjjm )
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22