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/ |
23 |
|
& MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE, |
24 |
|
& SEAICE_EPS, SEAICE_EPS_SQ |
25 |
|
_RL SEAICE_EPS, SEAICE_EPS_SQ |
26 |
|
_RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE |
27 |
|
|
28 |
C === Routine arguments === |
C === Routine arguments === |
29 |
C myThid - Thread no. that called this routine. |
C myThid - Thread no. that called this routine. |
30 |
_RL A22 |
_RL A22 |
38 |
|
|
39 |
INTEGER i, j, bi, bj |
INTEGER i, j, bi, bj |
40 |
INTEGER KOPEN |
INTEGER KOPEN |
41 |
_RL U1, V1 |
_RL U1, V1, SPEED_SQ |
42 |
|
|
43 |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
44 |
_RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
73 |
|
|
74 |
DO J=1,sNy |
DO J=1,sNy |
75 |
DO I=1,sNx |
DO I=1,sNx |
76 |
U1=0.25 _d 0*(UWIND(I-1,J-1,bi,bj)+UWIND(I-1,J,bi,bj) |
U1=0.25 _d 0*(UWIND(I+1,J+1,bi,bj)+UWIND(I+1,J,bi,bj) |
77 |
1 +UWIND(I,J-1,bi,bj)+UWIND(I,J,bi,bj)) |
1 +UWIND(I,J+1,bi,bj)+UWIND(I,J,bi,bj)) |
78 |
V1=0.25 _d 0*(VWIND(I-1,J-1,bi,bj)+VWIND(I-1,J,bi,bj) |
V1=0.25 _d 0*(VWIND(I+1,J+1,bi,bj)+VWIND(I+1,J,bi,bj) |
79 |
1 +VWIND(I,J-1,bi,bj)+VWIND(I,J,bi,bj)) |
1 +VWIND(I,J+1,bi,bj)+VWIND(I,J,bi,bj)) |
80 |
UG(I,J)=SQRT(U1**2+V1**2) |
SPEED_SQ = U1**2 + V1**2 |
81 |
|
IF ( SPEED_SQ .LT. SEAICE_EPS_SQ ) THEN |
82 |
|
UG(I,J)=SEAICE_EPS |
83 |
|
ELSE |
84 |
|
UG(I,J)=SQRT(SPEED_SQ) |
85 |
|
ENDIF |
86 |
ENDDO |
ENDDO |
87 |
ENDDO |
ENDDO |
88 |
|
|