28 |
#ifdef ALLOW_SEAICE |
#ifdef ALLOW_SEAICE |
29 |
|
|
30 |
C === Local variables === |
C === Local variables === |
31 |
C i,j,k,bi,bj - Loop counters |
C i,j,bi,bj - Loop counters |
32 |
|
|
33 |
INTEGER i, j, k, bi, bj |
INTEGER i, j, bi, bj |
34 |
_RL DELTXX1, DELTYY1, DELTXX, DELTYY |
_RL DELTXX1, DELTYY1, DELTXX, DELTYY |
35 |
_RL HEFF1 (1-OLx:sNx+OLx, 1-OLy:sNy+OLy) |
_RL HEFF1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
36 |
|
|
37 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
38 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
39 |
|
|
40 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
41 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
42 |
HEFF1(I,J)=0.0 |
HEFF1(I,J,bi,bj)=0.0 _d 0 |
43 |
ENDDO |
ENDDO |
44 |
ENDDO |
ENDDO |
45 |
|
|
52 |
& *CSTICE(I,J,bi,bj)*CSTICE(I,J,bi,bj)) |
& *CSTICE(I,J,bi,bj)*CSTICE(I,J,bi,bj)) |
53 |
DELTYY=DELTYY1/(DYTICE(I,J,bi,bj)*DYTICE(I,J,bi,bj) |
DELTYY=DELTYY1/(DYTICE(I,J,bi,bj)*DYTICE(I,J,bi,bj) |
54 |
& *CSTICE(I,J,bi,bj)) |
& *CSTICE(I,J,bi,bj)) |
55 |
HEFF1(I,J)=DELTXX*((HEFF(I+1,J,3,bi,bj)-HEFF(I,J,3,bi,bj)) |
HEFF1(I,J,bi,bj)=DELTXX*((HEFF(I+1,J,3,bi,bj) |
56 |
& *HEFFM(I+1,J,bi,bj) |
& -HEFF(I,J,3,bi,bj))*HEFFM(I+1,J,bi,bj) |
57 |
& -(HEFF(I,J,3,bi,bj)-HEFF(I-1,J,3,bi,bj))*HEFFM(I-1,J,bi,bj)) |
& -(HEFF(I,J,3,bi,bj)-HEFF(I-1,J,3,bi,bj))*HEFFM(I-1,J,bi,bj)) |
58 |
& +DELTYY*((HEFF(I,J+1,3,bi,bj)-HEFF(I,J,3,bi,bj)) |
& +DELTYY*((HEFF(I,J+1,3,bi,bj)-HEFF(I,J,3,bi,bj)) |
59 |
& *CSUICE(I,J,bi,bj)*HEFFM(I,J+1,bi,bj) |
& *CSUICE(I+1,J+1,bi,bj)*HEFFM(I,J+1,bi,bj) |
60 |
& -(HEFF(I,J,3,bi,bj)-HEFF(I,J-1,3,bi,bj)) |
& -(HEFF(I,J,3,bi,bj)-HEFF(I,J-1,3,bi,bj)) |
61 |
& *CSUICE(I,J-1,bi,bj)*HEFFM(I,J-1,bi,bj)) |
& *CSUICE(I+1,J,bi,bj)*HEFFM(I,J-1,bi,bj)) |
62 |
ENDDO |
ENDDO |
63 |
ENDDO |
ENDDO |
64 |
|
|
71 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
72 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
73 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
74 |
HEFF(I,J,3,bi,bj)=HEFF1(I,J) |
HEFF(I,J,3,bi,bj)=HEFF1(I,J,bi,bj) |
75 |
ENDDO |
ENDDO |
76 |
ENDDO |
ENDDO |
77 |
ENDDO |
ENDDO |
78 |
ENDDO |
ENDDO |
79 |
|
|
80 |
#endif ALLOW_SEAICE |
#endif /* ALLOW_SEAICE */ |
81 |
|
|
82 |
RETURN |
RETURN |
83 |
END |
END |