--- MITgcm/pkg/seaice/groatb.F	2003/04/30 02:24:55	1.1.2.6
+++ MITgcm/pkg/seaice/groatb.F	2003/02/28 09:47:31	1.4.4.4
@@ -18,13 +18,9 @@
 #include "DYNVARS.h"
 #include "FFIELDS.h"
 #include "SEAICE.h"
+#include "SEAICE_PARAMS.h"
 #include "SEAICE_FFIELDS.h"
-       COMMON /SEAICE_BOUND_RL/
-     &    MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE,
-     &    SEAICE_EPS, SEAICE_EPS_SQ
-      _RL SEAICE_EPS, SEAICE_EPS_SQ
-      _RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE
-
+ 
 C     === Routine arguments ===
 C     myThid - Thread no. that called this routine.
       _RL A22
@@ -73,10 +69,10 @@
 
         DO J=1,sNy
          DO I=1,sNx
-          U1=0.25 _d 0*(UWIND(I+1,J+1,bi,bj)+UWIND(I+1,J,bi,bj)
-     1      +UWIND(I,J+1,bi,bj)+UWIND(I,J,bi,bj))
-          V1=0.25 _d 0*(VWIND(I+1,J+1,bi,bj)+VWIND(I+1,J,bi,bj)
-     1      +VWIND(I,J+1,bi,bj)+VWIND(I,J,bi,bj))
+          U1=0.25 _d 0*(UWIND(I-1,J-1,bi,bj)+UWIND(I-1,J,bi,bj)
+     1      +UWIND(I,J-1,bi,bj)+UWIND(I,J,bi,bj))
+          V1=0.25 _d 0*(VWIND(I-1,J-1,bi,bj)+VWIND(I-1,J,bi,bj)
+     1      +VWIND(I,J-1,bi,bj)+VWIND(I,J,bi,bj))
           SPEED_SQ = U1**2 + V1**2
           IF ( SPEED_SQ .LT. SEAICE_EPS_SQ ) THEN
              UG(I,J)=SEAICE_EPS