--- MITgcm/pkg/generic_advdiff/gad_dst3_adv_r.F	2001/09/10 13:09:04	1.1
+++ MITgcm/pkg/generic_advdiff/gad_dst3_adv_r.F	2002/01/17 16:48:59	1.2
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/generic_advdiff/gad_dst3_adv_r.F,v 1.1 2001/09/10 13:09:04 adcroft Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/generic_advdiff/gad_dst3_adv_r.F,v 1.2 2002/01/17 16:48:59 adcroft Exp $
 C $Name:  $
 
 #include "GAD_OPTIONS.h"
@@ -57,13 +57,13 @@
       DO j=1-Oly,sNy+Oly
        DO i=1-Olx,sNx+Olx
         Rjp=(tracer(i,j,k,bi,bj)-tracer(i,j,kp1,bi,bj))
-     &         *maskC(i,j,kp1,bi,bj)
+     &         *maskC(i,j,kp1,bi_arg,bj_arg)
         Rj =(tracer(i,j,km1,bi,bj)-tracer(i,j,k,bi,bj))
-     &         *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
+     &         *maskC(i,j,k,bi_arg,bj_arg)*maskC(i,j,km1,bi_arg,bj_arg)
         Rjm=(tracer(i,j,km2,bi,bj)-tracer(i,j,km1,bi,bj))
-     &         *maskC(i,j,km1,bi,bj)
+     &         *maskC(i,j,km1,bi_arg,bj_arg)
 
-        cfl=abs(wVel(i,j,k,bi,bj)*dTarg*recip_drc(k))
+        cfl=abs(wVel(i,j,k,bi_arg,bj_arg)*dTarg*recip_drc(k))
         d0=(2.-cfl)*(1.-cfl)*oneSixth
         d1=(1.-cfl*cfl)*oneSixth
 c       thetaP=0.