--- MITgcm/pkg/generic_advdiff/gad_dst3fl_adv_x.F	2001/11/08 23:39:34	1.4
+++ MITgcm/pkg/generic_advdiff/gad_dst3fl_adv_x.F	2002/03/06 01:29:36	1.5
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/generic_advdiff/gad_dst3fl_adv_x.F,v 1.4 2001/11/08 23:39:34 heimbach Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/generic_advdiff/gad_dst3fl_adv_x.F,v 1.5 2002/03/06 01:29:36 jmc Exp $
 C $Name:  $
 
 #include "GAD_OPTIONS.h"
@@ -31,8 +31,10 @@
       INTEGER myThid
 
 C     == Local variables ==
+C     uFld   :: velocity [m/s], zonal component
       INTEGER i,j
       _RL Rjm,Rj,Rjp,cfl,d0,d1,psiP,psiM,thetaP,thetaM
+      _RL uFld
 
       DO j=1-Oly,sNy+Oly
        uT(1-Olx,j)=0.D0
@@ -43,7 +45,10 @@
         Rj =(tracer(i,j)-tracer(i-1,j))*maskW(i,j,k,bi,bj)
         Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskW(i-1,j,k,bi,bj)
 
-        cfl=abs(uVel(i,j,k,bi,bj)*deltaT*recip_dxc(i,j,bi,bj))
+c       uFld = uVel(i,j,k,bi,bj)
+        uFld = uTrans(i,j)*recip_dyG(i,j,bi,bj)
+     &       *recip_drF(k)*recip_hFacW(i,j,k,bi,bj)
+        cfl=abs(uFld*deltaT*recip_dxC(i,j,bi,bj))
         d0=(2.D0-cfl)*(1.D0-cfl)*oneSixth
         d1=(1.D0-cfl*cfl)*oneSixth
 c       thetaP=0.D0