--- MITgcm/pkg/generic_advdiff/gad_dst3_adv_y.F	2001/09/10 00:14:05	1.2
+++ MITgcm/pkg/generic_advdiff/gad_dst3_adv_y.F	2002/03/06 01:29:36	1.3
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/generic_advdiff/gad_dst3_adv_y.F,v 1.2 2001/09/10 00:14:05 adcroft Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/generic_advdiff/gad_dst3_adv_y.F,v 1.3 2002/03/06 01:29:36 jmc Exp $
 C $Name:  $
 
 #include "GAD_OPTIONS.h"
@@ -31,9 +31,11 @@
       INTEGER myThid
 
 C     == Local variables ==
+C     vFld    :: velocity [m/s], meridional component
       INTEGER i,j
       _RL Rjm,Rj,Rjp,cfl,d0,d1
       _RL psiP,psiM,thetaP,thetaM
+      _RL vFld
 
       DO i=1-Olx,sNx+Olx
        vT(i,1-Oly)=0.
@@ -46,7 +48,10 @@
         Rj =(tracer(i,j)-tracer(i,j-1))*maskS(i,j,k,bi,bj)
         Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskS(i,j-1,k,bi,bj)
 
-        cfl=abs(vVel(i,j,k,bi,bj)*deltaT*recip_dyc(i,j,bi,bj))
+c       vFld = vVel(i,j,k,bi,bj)
+        vFld = vTrans(i,j)*recip_dxG(i,j,bi,bj)
+     &       *recip_drF(k)*recip_hFacS(i,j,k,bi,bj)
+        cfl=abs(vFld*deltaT*recip_dyC(i,j,bi,bj))
         d0=(2.-cfl)*(1.-cfl)*oneSixth
         d1=(1.-cfl*cfl)*oneSixth
 c       thetaP=0.