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

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

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

revision 1.4 by adcroft, Fri Sep 21 13:11:43 2001 UTC revision 1.5 by jmc, Thu May 1 16:12:36 2003 UTC
# Line 61  CEOP Line 61  CEOP
61        DO j=1-Oly+2,sNy+Oly-1        DO j=1-Oly+2,sNy+Oly-1
62         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
63          Rjp=(tracer(i,j+1)-tracer(i,j))          Rjp=(tracer(i,j+1)-tracer(i,j))
64  c    &     *maskS(i,j+1,k,bi,bj)       &     *maskS(i,j+1,k,bi,bj)
65          Rj =(tracer(i,j)-tracer(i,j-1))          Rj =(tracer(i,j)-tracer(i,j-1))
66  c    &     *maskS(i,j,k,bi,bj)       &     *maskS(i,j,k,bi,bj)
67          Rjm=(tracer(i,j-1)-tracer(i,j-2))          Rjm=(tracer(i,j-1)-tracer(i,j-2))
68  c    &     *maskS(i,j-1,k,bi,bj)       &     *maskS(i,j-1,k,bi,bj)
69  C- jmc: mask not needed here above if Rjjp & Rjjm & vTrans are masked          Rjjp=(Rjp-Rj)
70          Rjjp=(Rjp-Rj)*maskS(i,j+1,k,bi,bj)          Rjjm=(Rj-Rjm)
         Rjjm=(Rj-Rjm)*maskS(i,j-1,k,bi,bj)  
71          vT(i,j) =          vT(i,j) =
72       &   vTrans(i,j)*(       &   vTrans(i,j)*(
73       &     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.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22