--- MITgcm_contrib/dgoldberg/streamice/streamice_taub.F	2012/09/20 02:04:45	1.2
+++ MITgcm_contrib/dgoldberg/streamice/streamice_taub.F	2012/09/20 15:45:13	1.3
@@ -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.3 2012/09/20 15:45:13 dgoldberg Exp $
 C $Name:  $
 
 #include "STREAMICE_OPTIONS.h"
@@ -46,10 +46,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