43 |
_EXCH_XY_R8(GAIRX, myThid) |
_EXCH_XY_R8(GAIRX, myThid) |
44 |
_EXCH_XY_R8(GAIRY, myThid) |
_EXCH_XY_R8(GAIRY, myThid) |
45 |
|
|
|
c$taf loop = parallel |
|
46 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
|
c$taf loop = parallel |
|
47 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
48 |
|
|
49 |
C DETERMINE AMOUNT OF OPEN WATER AND ICE THICKNESS |
C DETERMINE AMOUNT OF OPEN WATER AND ICE THICKNESS |
50 |
DO J=1,sNy |
DO J=1,sNy |
51 |
DO I=1,sNx |
DO I=1,sNx |
52 |
AREA(I,J,2,bi,bj)=MAX(A22,AREA(I,J,2,bi,bj)) |
AREA(I,J,2,bi,bj)=MAX(A22,AREA(I,J,2,bi,bj)) |
53 |
FHEFF(I,J,bi,bj)=0.0 _d 0 |
FHEFF(I,J,bi,bj)=0.0 |
54 |
HICE(I,J)=HEFF(I,J,2,bi,bj)/AREA(I,J,2,bi,bj) |
HICE(I,J)=HEFF(I,J,2,bi,bj)/AREA(I,J,2,bi,bj) |
55 |
ENDDO |
ENDDO |
56 |
ENDDO |
ENDDO |
58 |
C NOW DETERMINE MIXED LAYER TEMPERATURE |
C NOW DETERMINE MIXED LAYER TEMPERATURE |
59 |
DO J=1,sNy |
DO J=1,sNy |
60 |
DO I=1,sNx |
DO I=1,sNx |
61 |
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 |
62 |
#ifdef SEAICE_DEBUG |
#ifdef SEAICE_DEBUG |
63 |
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) |
64 |
#endif SEAICE_DEBUG |
#endif SEAICE_DEBUG |
65 |
ENDDO |
ENDDO |
66 |
ENDDO |
ENDDO |
67 |
|
|
68 |
DO J=1,sNy |
DO J=1,sNy |
69 |
DO I=1,sNx |
DO I=1,sNx |
70 |
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) |
71 |
1 +GAIRX(I,J-1,bi,bj)+GAIRX(I,J,bi,bj)) |
1 +GAIRX(I,J-1,bi,bj)+GAIRX(I,J,bi,bj)) |
72 |
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) |
73 |
1 +GAIRY(I,J-1,bi,bj)+GAIRY(I,J,bi,bj)) |
1 +GAIRY(I,J-1,bi,bj)+GAIRY(I,J,bi,bj)) |
74 |
UG(I,J)=SQRT(U1**2+V1**2) |
UG(I,J)=SQRT(U1**2+V1**2) |
75 |
ENDDO |
ENDDO |