--- MITgcm_contrib/dgoldberg/streamice/streamice_taub.F	2012/09/18 17:04:40	1.1
+++ 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.1 2012/09/18 17:04:40 dgoldberg 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"
@@ -40,12 +40,37 @@
        DO bi=myBxLo(myThid),myBxHi(myThid)
         DO j=1,sNy
          DO i=1,sNx
+cph(
+          streamice_taubx (i,j,bi,bj) = 0.
+          streamice_tauby (i,j,bi,bj) = 0.
+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