--- MITgcm/pkg/seaice/groatb.F 2003/03/03 06:58:38 1.1.2.5 +++ MITgcm/pkg/seaice/groatb.F 2003/02/13 19:28:39 1.4.4.1 @@ -19,12 +19,7 @@ #include "FFIELDS.h" #include "SEAICE.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 @@ -38,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) @@ -67,7 +62,7 @@ TMIX(I,J,bi,bj)=theta(I,J,1,bi,bj)+273.16 _d +00 #ifdef SEAICE_DEBUG TMIX(I,J,bi,bj)=MAX(TMIX(I,J,bi,bj),271.2 _d +00) -#endif +#endif SEAICE_DEBUG ENDDO ENDDO @@ -77,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 @@ -97,7 +87,7 @@ ENDDO ENDDO -#endif /* ALLOW_SEAICE */ +#endif ALLOW_SEAICE RETURN END