--- MITgcm_contrib/dgoldberg/streamice/streamice_taub.F	2012/09/20 02:04:45	1.2
+++ MITgcm_contrib/dgoldberg/streamice/streamice_taub.F	2012/09/25 03:03:21	1.4
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/dgoldberg/streamice/streamice_taub.F,v 1.2 2012/09/20 02:04:45 heimbach Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/dgoldberg/streamice/streamice_taub.F,v 1.4 2012/09/25 03:03:21 dgoldberg Exp $
 C $Name:  $
 
 #include "STREAMICE_OPTIONS.h"
@@ -33,8 +33,9 @@
 
 
 C     LOCAL VARIABLES
-      INTEGER i, j, bi, bj, k, l, umid, vmid
+      INTEGER i, j, bi, bj, k, l
       _RL ux, uy, vx, vy, exx, eyy, exy, second_inv, unorm
+      _RL umid, vmid
 
       DO bj=myByLo(myThid),myByHi(myThid)
        DO bi=myBxLo(myThid),myBxHi(myThid)
@@ -46,10 +47,31 @@
 cph)
           IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN
   
+!            streamice_taubx (i,j,bi,bj) = 
+!      &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
+!            streamice_tauby (i,j,bi,bj) = 
+!      &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
+
+           umid = 0.0
+           vmid = 0.0
+
+           DO k=0,1
+            DO l=0,1
+             umid = umid + 0.25 * 
+     &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) * 
+     &        recip_rA(i,j,bi,bj) *
+     &        U_streamice(i+k,j+l,bi,bj) 
+             vmid = vmid + 0.25 * 
+     &        dxG(i,j+l,bi,bj)*dyG(i+k,j,bi,bj) * 
+     &        recip_rA(i,j,bi,bj) *
+     &        V_streamice(i+k,j+l,bi,bj)
+            ENDDO
+           ENDDO
+  
            streamice_taubx (i,j,bi,bj) = 
-     &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
+     &      tau_beta_eff_streamice(i,j,bi,bj) * umid
            streamice_tauby (i,j,bi,bj) = 
-     &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
+     &      tau_beta_eff_streamice(i,j,bi,bj) * vmid
            
           ENDIF
          ENDDO