--- MITgcm/pkg/seaice/groatb.F 2002/12/05 08:06:17 1.1.2.2 +++ MITgcm/pkg/seaice/groatb.F 2002/11/12 20:47:27 1.2 @@ -43,16 +43,14 @@ _EXCH_XY_R8(GAIRX, myThid) _EXCH_XY_R8(GAIRY, myThid) -c$taf loop = parallel DO bj=myByLo(myThid),myByHi(myThid) -c$taf loop = parallel DO bi=myBxLo(myThid),myBxHi(myThid) C DETERMINE AMOUNT OF OPEN WATER AND ICE THICKNESS DO J=1,sNy DO I=1,sNx AREA(I,J,2,bi,bj)=MAX(A22,AREA(I,J,2,bi,bj)) - FHEFF(I,J,bi,bj)=0.0 _d 0 + FHEFF(I,J,bi,bj)=0.0 HICE(I,J)=HEFF(I,J,2,bi,bj)/AREA(I,J,2,bi,bj) ENDDO ENDDO @@ -60,18 +58,18 @@ C NOW DETERMINE MIXED LAYER TEMPERATURE DO J=1,sNy DO I=1,sNx - TMIX(I,J,bi,bj)=theta(I,J,1,bi,bj)+273.16 _d +00 + TMIX(I,J,bi,bj)=theta(I,J,1,bi,bj)+273.16E+00 #ifdef SEAICE_DEBUG - TMIX(I,J,bi,bj)=MAX(TMIX(I,J,bi,bj),271.2 _d +00) + TMIX(I,J,bi,bj)=MAX(TMIX(I,J,bi,bj),271.2E+00) #endif SEAICE_DEBUG ENDDO ENDDO DO J=1,sNy DO I=1,sNx - U1=0.25 _d 0*(GAIRX(I-1,J-1,bi,bj)+GAIRX(I-1,J,bi,bj) + U1=0.25*(GAIRX(I-1,J-1,bi,bj)+GAIRX(I-1,J,bi,bj) 1 +GAIRX(I,J-1,bi,bj)+GAIRX(I,J,bi,bj)) - V1=0.25 _d 0*(GAIRY(I-1,J-1,bi,bj)+GAIRY(I-1,J,bi,bj) + V1=0.25*(GAIRY(I-1,J-1,bi,bj)+GAIRY(I-1,J,bi,bj) 1 +GAIRY(I,J-1,bi,bj)+GAIRY(I,J,bi,bj)) UG(I,J)=SQRT(U1**2+V1**2) ENDDO