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

Diff of /MITgcm/pkg/generic_advdiff/gad_c4_adv_x.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 59  CEOP Line 59  CEOP
59         uT(sNx+Olx,j)=0.         uT(sNx+Olx,j)=0.
60         DO i=1-Olx+2,sNx+Olx-1         DO i=1-Olx+2,sNx+Olx-1
61          Rjp=(tracer(i+1,j)-tracer(i,j))          Rjp=(tracer(i+1,j)-tracer(i,j))
62  c    &     *maskW(i+1,j,k,bi,bj)       &     *maskW(i+1,j,k,bi,bj)
63          Rj =(tracer(i,j)-tracer(i-1,j))          Rj =(tracer(i,j)-tracer(i-1,j))
64  c    &     *maskW(i,j,k,bi,bj)       &     *maskW(i,j,k,bi,bj)
65          Rjm=(tracer(i-1,j)-tracer(i-2,j))          Rjm=(tracer(i-1,j)-tracer(i-2,j))
66  c    &     *maskW(i-1,j,k,bi,bj)       &     *maskW(i-1,j,k,bi,bj)
67  C- jmc: mask not needed here above if Rjjp & Rjjm & uTrans are masked          Rjjp=(Rjp-Rj)
68          Rjjp=(Rjp-Rj)*maskW(i+1,j,k,bi,bj)          Rjjm=(Rj-Rjm)
         Rjjm=(Rj-Rjm)*maskW(i-1,j,k,bi,bj)  
69          uT(i,j) =          uT(i,j) =
70       &   uTrans(i,j)*(       &   uTrans(i,j)*(
71       &     Tracer(i,j)+Tracer(i-1,j)-oneSixth*( Rjjp+Rjjm )       &     Tracer(i,j)+Tracer(i-1,j)-oneSixth*( Rjjp+Rjjm )
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22