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 |
SPEED_SQ = U1**2 + V1**2 |
SPEED_SQ = U1**2 + V1**2 |
81 |
IF ( SPEED_SQ .LT. SEAICE_EPS_SQ ) THEN |
IF ( SPEED_SQ .LT. SEAICE_EPS_SQ ) THEN |
82 |
UG(I,J)=SEAICE_EPS |
UG(I,J)=SEAICE_EPS |