43 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
44 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
45 |
|
|
46 |
uVel(I,J,K,bi,bj)=0. |
uVel(I,J,K,bi,bj)=0. _d 0 |
47 |
vVel(I,J,K,bi,bj)=0. |
vVel(I,J,K,bi,bj)=0. _d 0 |
48 |
wVel(I,J,K,bi,bj)=0. |
wVel(I,J,K,bi,bj)=0. _d 0 |
49 |
theta(I,J,K,bi,bj)=0. |
theta(I,J,K,bi,bj)=0. _d 0 |
50 |
salt(I,J,K,bi,bj)=0. |
salt(I,J,K,bi,bj)=0. _d 0 |
51 |
|
|
52 |
gU(I,J,K,bi,bj)=0. |
gU(I,J,K,bi,bj)=0. _d 0 |
53 |
gV(I,J,K,bi,bj)=0. |
gV(I,J,K,bi,bj)=0. _d 0 |
54 |
gT(I,J,K,bi,bj)=0. |
gT(I,J,K,bi,bj)=0. _d 0 |
55 |
gS(I,J,K,bi,bj)=0. |
gS(I,J,K,bi,bj)=0. _d 0 |
56 |
gUnm1(I,J,K,bi,bj)=0. |
gUnm1(I,J,K,bi,bj)=0. _d 0 |
57 |
gVnm1(I,J,K,bi,bj)=0. |
gVnm1(I,J,K,bi,bj)=0. _d 0 |
58 |
gTnm1(I,J,K,bi,bj)=0. |
gTnm1(I,J,K,bi,bj)=0. _d 0 |
59 |
gSnm1(I,J,K,bi,bj)=0. |
gSnm1(I,J,K,bi,bj)=0. _d 0 |
60 |
|
|
61 |
#ifdef INCLUDE_CD_CODE |
#ifdef INCLUDE_CD_CODE |
62 |
guCD(I,J,K,bi,bj)=0. |
guCD(I,J,K,bi,bj)=0. _d 0 |
63 |
gvCD(I,J,K,bi,bj)=0. |
gvCD(I,J,K,bi,bj)=0. _d 0 |
64 |
uNM1(I,J,K,bi,bj)=0. |
uNM1(I,J,K,bi,bj)=0. _d 0 |
65 |
vNM1(I,J,K,bi,bj)=0. |
vNM1(I,J,K,bi,bj)=0. _d 0 |
66 |
uVeld(I,J,K,bi,bj)=0. |
uVeld(I,J,K,bi,bj)=0. _d 0 |
67 |
vVeld(I,J,K,bi,bj)=0. |
vVeld(I,J,K,bi,bj)=0. _d 0 |
68 |
#endif |
#endif |
69 |
|
|
70 |
#ifdef ALLOW_NONHYDROSTATIC |
#ifdef ALLOW_NONHYDROSTATIC |
71 |
gW(I,J,K,bi,bj)=0. |
gW(I,J,K,bi,bj)=0. _d 0 |
72 |
gWnm1(I,J,K,bi,bj)=0. |
gWnm1(I,J,K,bi,bj)=0. _d 0 |
73 |
#endif |
#endif |
74 |
|
|
75 |
ENDDO |
ENDDO |
79 |
C- 2D arrays |
C- 2D arrays |
80 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
81 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
82 |
etaN(I,J,bi,bj)=0. |
etaN(I,J,bi,bj)=0. _d 0 |
83 |
etaNm1(I,J,bi,bj)=0. |
etaNm1(I,J,bi,bj)=0. _d 0 |
84 |
ENDDO |
ENDDO |
85 |
ENDDO |
ENDDO |
86 |
|
|