1 |
C |
C |
2 |
|
|
3 |
#include "SEAICE_OPTIONS.h" |
#include "SEAICE_OPTIONS.h" |
4 |
|
|
5 |
CStartOfInterface |
CStartOfInterface |
6 |
SUBROUTINE groatb( A22, myThid ) |
SUBROUTINE groatb( A22, myThid ) |
7 |
C /==========================================================\ |
C /==========================================================\ |
43 |
|
|
44 |
INTEGER i, j, bi, bj |
INTEGER i, j, bi, bj |
45 |
INTEGER KOPEN |
INTEGER KOPEN |
46 |
_RL U1, V1, SPEED_SQ |
_RL SPEED_SQ |
47 |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL HICE (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
48 |
_RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL UG (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
49 |
|
|
54 |
_RL FICEP(1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL FICEP(1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
55 |
#endif |
#endif |
56 |
|
|
|
C if this is not done before then do it here |
|
|
_EXCH_XY_R8(UWIND, myThid) |
|
|
_EXCH_XY_R8(VWIND, myThid) |
|
|
|
|
57 |
CADJ loop = parallel |
CADJ loop = parallel |
58 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
59 |
CADJ loop = parallel |
CADJ loop = parallel |
110 |
|
|
111 |
DO J=1,sNy |
DO J=1,sNy |
112 |
DO I=1,sNx |
DO I=1,sNx |
113 |
U1=0.25 _d 0*(UWIND(I+1,J+1,bi,bj)+UWIND(I+1,J,bi,bj) |
SPEED_SQ = UWIND(I,J,bi,bj)**2 + VWIND(I,J,bi,bj)**2 |
|
1 +UWIND(I,J+1,bi,bj)+UWIND(I,J,bi,bj)) |
|
|
V1=0.25 _d 0*(VWIND(I+1,J+1,bi,bj)+VWIND(I+1,J,bi,bj) |
|
|
1 +VWIND(I,J+1,bi,bj)+VWIND(I,J,bi,bj)) |
|
|
SPEED_SQ = U1**2 + V1**2 |
|
114 |
IF ( SPEED_SQ .LE. SEAICE_EPS_SQ ) THEN |
IF ( SPEED_SQ .LE. SEAICE_EPS_SQ ) THEN |
115 |
UG(I,J)=SEAICE_EPS |
UG(I,J)=SEAICE_EPS |
116 |
ELSE |
ELSE |