18 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
19 |
#include "FFIELDS.h" |
#include "FFIELDS.h" |
20 |
#include "SEAICE.h" |
#include "SEAICE.h" |
21 |
|
#include "SEAICE_PARAMS.h" |
22 |
#include "SEAICE_FFIELDS.h" |
#include "SEAICE_FFIELDS.h" |
23 |
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 |
|
|
|
|
24 |
C === Routine arguments === |
C === Routine arguments === |
25 |
C myThid - Thread no. that called this routine. |
C myThid - Thread no. that called this routine. |
26 |
_RL A22 |
_RL A22 |
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*(UWIND(I-1,J-1,bi,bj)+UWIND(I-1,J,bi,bj) |
73 |
1 +UWIND(I,J+1,bi,bj)+UWIND(I,J,bi,bj)) |
1 +UWIND(I,J-1,bi,bj)+UWIND(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*(VWIND(I-1,J-1,bi,bj)+VWIND(I-1,J,bi,bj) |
75 |
1 +VWIND(I,J+1,bi,bj)+VWIND(I,J,bi,bj)) |
1 +VWIND(I,J-1,bi,bj)+VWIND(I,J,bi,bj)) |
76 |
SPEED_SQ = U1**2 + V1**2 |
SPEED_SQ = U1**2 + V1**2 |
77 |
IF ( SPEED_SQ .LT. SEAICE_EPS_SQ ) THEN |
IF ( SPEED_SQ .LT. SEAICE_EPS_SQ ) THEN |
78 |
UG(I,J)=SEAICE_EPS |
UG(I,J)=SEAICE_EPS |