*ID WMCOIHF
*/
*/ First, ICE_CTL
*/
*DECLARE ICE_CTL
*I ICE_CTL.219

c ITT doesn't seem to get passed in at vn4.5...
c     IF (itt .eq. 1 .and. o_mype .eq. 0) THEN
c       print *,'wmc: AICEMIZFRY, L_HADCM4O2I: ',AICEMIZFRY,L_HADCM4O2I
c     ENDIF

*/ Define, make "Magnitude of U Dashed" (mud) and pass it to ICEFLOE
*I ICE_CTL.83
      real mud(imt,jmt)       ! Mag of o-i velocity shear (on a-grid)
      real mud_uv(imt,jmtm1)  ! Mag of o-i velocity shear (on u-grid)

*I ORH1F305.5731 
      DO i=1,imt
        DO j=1,jmtm1
          mud_uv(i,j)=sqrt((ucurrent(i,j)-uice(i,j))**2
     &                    +(vcurrent(i,j)-vice(i,j))**2)
        ENDDO
      ENDDO
      CALL uv_to_h( 
*CALL ARGOINDX                       
     & mud_uv,mud,imt,jmt,jmtm1)  
*I ODC1F405.341
     &,mud                 ! Magnitude of o-i vel shear (on a-grid)
*/ Add diags
*I ODC1F405.326
C         IF (SF(150,32)) CALL COPYODIAGL(IMT,JMTM1,1,.TRUE.,0.
C     &           ,uice,ocean_uv,STASHWORK(SI(150,32,im_index)))
C
C         IF (SF(151,32)) CALL COPYODIAGL(IMT,JMTM1,1,.TRUE.,0.
C     &           ,vice,ocean_uv,STASHWORK(SI(151,32,im_index)))
*/
*/ Now onto ICEFLOE
*/
*DECLARE ICEFLOE
*I ODC1F405.439
     &,mud                 ! Magnitude of o-i vel shear (on a-grid)
*I ODC1F405.440
      real mud(imt,jmt)
*I ICEFLOE.272
C
C Make OI heat flux proportional to shear between ocean and ice.
C Constant of 0.1 (m/s) means no change at 10cm/s, increases above
C this and decreases below...
C
      OCEANFLX(I,J) = OCEANFLX(I,J) * mud(i,j) / 0.1