--- MITgcm_contrib/dgoldberg/streamice/streamice_taub.F	2012/09/18 17:04:40	1.1
+++ MITgcm_contrib/dgoldberg/streamice/streamice_taub.F	2012/12/10 02:34:45	1.5
@@ -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.5 2012/12/10 02:34:45 dgoldberg Exp $
 C $Name:  $
 
 #include "STREAMICE_OPTIONS.h"
@@ -30,22 +30,49 @@
 !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
 
 #ifdef ALLOW_STREAMICE
+#ifdef STREAMICE_HYBRID_STRESS
 
 
 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)
         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
@@ -55,5 +82,6 @@
 
 
 #endif
+#endif
       RETURN
       END