/[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.5 by adcroft, Fri Sep 21 13:11:43 2001 UTC revision 1.6 by adcroft, Thu Jan 17 16:48:59 2002 UTC
# Line 93  C      for maskC(...) and wVel(...) Line 93  C      for maskC(...) and wVel(...)
93         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
94          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
95           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))
96       &        *maskC(i,j,kp1,bi,bj)       &        *maskC(i,j,kp1,bi_arg,bj_arg)
97           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))
98           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))
99       &        *maskC(i,j,km2,bi,bj)       &        *maskC(i,j,km2,bi_arg,bj_arg)
100           IF (Rj.NE.0.) THEN           IF (Rj.NE.0.) THEN
101            IF (rTrans(i,j).LT.0.) THEN            IF (rTrans(i,j).LT.0.) THEN
102              Cr=Rjm/Rj              Cr=Rjm/Rj
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22