56 |
_RL aC, aCw, aCs |
_RL aC, aCw, aCs |
57 |
CEOP |
CEOP |
58 |
|
|
59 |
|
C-- Initialize arrays in common blocs (CG2D.h) ; not really necessary |
60 |
|
C but safer when EXCH do not fill all the overlap regions. |
61 |
|
DO bj=myByLo(myThid),myByHi(myThid) |
62 |
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
63 |
|
DO J=1-OLy,sNy+OLy |
64 |
|
DO I=1-OLx,sNx+OLx |
65 |
|
aW2d(I,J,bi,bj) = 0. _d 0 |
66 |
|
aS2d(I,J,bi,bj) = 0. _d 0 |
67 |
|
pW(I,J,bi,bj) = 0. _d 0 |
68 |
|
pS(I,J,bi,bj) = 0. _d 0 |
69 |
|
pC(I,J,bi,bj) = 0. _d 0 |
70 |
|
cg2d_q(I,J,bi,bj) = 0. _d 0 |
71 |
|
ENDDO |
72 |
|
ENDDO |
73 |
|
DO J=1-1,sNy+1 |
74 |
|
DO I=1-1,sNx+1 |
75 |
|
cg2d_r(I,J,bi,bj) = 0. _d 0 |
76 |
|
cg2d_s(I,J,bi,bj) = 0. _d 0 |
77 |
|
ENDDO |
78 |
|
ENDDO |
79 |
|
ENDDO |
80 |
|
ENDDO |
81 |
|
|
82 |
C-- Initialise laplace operator |
C-- Initialise laplace operator |
83 |
C aW2d: integral in Z Ax/dX |
C aW2d: integral in Z Ax/dX |
84 |
C aS2d: integral in Z Ay/dY |
C aS2d: integral in Z Ay/dY |