19 |
#include "FFIELDS.h" |
#include "FFIELDS.h" |
20 |
#include "SEAICE.h" |
#include "SEAICE.h" |
21 |
#include "SEAICE_FFIELDS.h" |
#include "SEAICE_FFIELDS.h" |
22 |
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 |
|
|
|
|
23 |
C === Routine arguments === |
C === Routine arguments === |
24 |
C myThid - Thread no. that called this routine. |
C myThid - Thread no. that called this routine. |
25 |
_RL A22 |
_RL A22 |
29 |
#ifdef ALLOW_SEAICE |
#ifdef ALLOW_SEAICE |
30 |
|
|
31 |
C === Local variables === |
C === Local variables === |
32 |
C i,j,bi,bj - Loop counters |
C i,j,k,bi,bj - Loop counters |
33 |
|
|
34 |
INTEGER i, j, bi, bj |
INTEGER i, j, k, bi, bj |
35 |
INTEGER KOPEN |
INTEGER KOPEN |
36 |
_RL U1, V1, SPEED_SQ |
_RL U1, V1 |
37 |
|
|
38 |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
39 |
|
_RL AR (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
40 |
_RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
41 |
|
|
42 |
C if this is not done before then do it here |
C if this is not done before then do it here |
43 |
_EXCH_XY_R8(UWIND, myThid) |
_EXCH_XY_R8(GAIRX, myThid) |
44 |
_EXCH_XY_R8(VWIND, myThid) |
_EXCH_XY_R8(GAIRY, myThid) |
45 |
|
|
46 |
c$taf loop = parallel |
c$taf loop = parallel |
47 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
63 |
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.16 _d +00 |
64 |
#ifdef SEAICE_DEBUG |
#ifdef SEAICE_DEBUG |
65 |
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.2 _d +00) |
66 |
#endif |
#endif SEAICE_DEBUG |
67 |
ENDDO |
ENDDO |
68 |
ENDDO |
ENDDO |
69 |
|
|
70 |
DO J=1,sNy |
DO J=1,sNy |
71 |
DO I=1,sNx |
DO I=1,sNx |
72 |
U1=0.25 _d 0*(UWIND(I+1,J+1,bi,bj)+UWIND(I+1,J,bi,bj) |
U1=0.25 _d 0*(GAIRX(I-1,J-1,bi,bj)+GAIRX(I-1,J,bi,bj) |
73 |
1 +UWIND(I,J+1,bi,bj)+UWIND(I,J,bi,bj)) |
1 +GAIRX(I,J-1,bi,bj)+GAIRX(I,J,bi,bj)) |
74 |
V1=0.25 _d 0*(VWIND(I+1,J+1,bi,bj)+VWIND(I+1,J,bi,bj) |
V1=0.25 _d 0*(GAIRY(I-1,J-1,bi,bj)+GAIRY(I-1,J,bi,bj) |
75 |
1 +VWIND(I,J+1,bi,bj)+VWIND(I,J,bi,bj)) |
1 +GAIRY(I,J-1,bi,bj)+GAIRY(I,J,bi,bj)) |
76 |
SPEED_SQ = U1**2 + V1**2 |
UG(I,J)=SQRT(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 |
|
77 |
ENDDO |
ENDDO |
78 |
ENDDO |
ENDDO |
79 |
|
|
88 |
ENDDO |
ENDDO |
89 |
ENDDO |
ENDDO |
90 |
|
|
91 |
#endif /* ALLOW_SEAICE */ |
#endif ALLOW_SEAICE |
92 |
|
|
93 |
RETURN |
RETURN |
94 |
END |
END |