--- MITgcm/pkg/seaice/groatb.F 2003/02/28 09:47:31 1.4.4.4 +++ MITgcm/pkg/seaice/groatb.F 2003/02/18 05:33:55 1.5 @@ -18,7 +18,6 @@ #include "DYNVARS.h" #include "FFIELDS.h" #include "SEAICE.h" -#include "SEAICE_PARAMS.h" #include "SEAICE_FFIELDS.h" C === Routine arguments === @@ -34,7 +33,7 @@ INTEGER i, j, bi, bj INTEGER KOPEN - _RL U1, V1, SPEED_SQ + _RL U1, V1 _RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) _RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) @@ -73,12 +72,7 @@ 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 - ELSE - UG(I,J)=SQRT(SPEED_SQ) - ENDIF + UG(I,J)=SQRT(U1**2+V1**2) ENDDO ENDDO