25 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
26 |
#include "PARAMS.h" |
#include "PARAMS.h" |
27 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
|
#ifdef ALLOW_NONHYDROSTATIC |
|
|
#include "GW.h" |
|
|
#endif |
|
28 |
|
|
29 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
30 |
C == Routine arguments == |
C == Routine arguments == |
58 |
gV(I,J,K,bi,bj)=0. _d 0 |
gV(I,J,K,bi,bj)=0. _d 0 |
59 |
gT(I,J,K,bi,bj)=0. _d 0 |
gT(I,J,K,bi,bj)=0. _d 0 |
60 |
gS(I,J,K,bi,bj)=0. _d 0 |
gS(I,J,K,bi,bj)=0. _d 0 |
61 |
gUnm1(I,J,K,bi,bj)=0. _d 0 |
#ifdef ALLOW_ADAMSBASHFORTH_3 |
62 |
gVnm1(I,J,K,bi,bj)=0. _d 0 |
guNm(I,J,K,bi,bj,1)=0. _d 0 |
63 |
gTnm1(I,J,K,bi,bj)=0. _d 0 |
gvNm(I,J,K,bi,bj,1)=0. _d 0 |
64 |
gSnm1(I,J,K,bi,bj)=0. _d 0 |
gtNm(I,J,K,bi,bj,1)=0. _d 0 |
65 |
|
gsNm(I,J,K,bi,bj,1)=0. _d 0 |
66 |
#ifdef ALLOW_NONHYDROSTATIC |
guNm(I,J,K,bi,bj,2)=0. _d 0 |
67 |
ceh3 needs an IF ( useNONHYDROSTATIC ) THEN |
gvNm(I,J,K,bi,bj,2)=0. _d 0 |
68 |
gW(I,J,K,bi,bj)=0. _d 0 |
gtNm(I,J,K,bi,bj,2)=0. _d 0 |
69 |
gWnm1(I,J,K,bi,bj)=0. _d 0 |
gsNm(I,J,K,bi,bj,2)=0. _d 0 |
70 |
|
#else |
71 |
|
guNm1(I,J,K,bi,bj) =0. _d 0 |
72 |
|
gvNm1(I,J,K,bi,bj) =0. _d 0 |
73 |
|
gtNm1(I,J,K,bi,bj) =0. _d 0 |
74 |
|
gsNm1(I,J,K,bi,bj) =0. _d 0 |
75 |
#endif |
#endif |
76 |
|
|
77 |
totPhiHyd(I,J,K,bi,bj)=0. _d 0 |
totPhiHyd(I,J,K,bi,bj)=0. _d 0 |
86 |
etaN(I,J,bi,bj)=0. _d 0 |
etaN(I,J,bi,bj)=0. _d 0 |
87 |
etaH(I,J,bi,bj)=0. _d 0 |
etaH(I,J,bi,bj)=0. _d 0 |
88 |
phiHydLow(I,J,bi,bj)=0. _d 0 |
phiHydLow(I,J,bi,bj)=0. _d 0 |
89 |
|
hMixLayer(I,J,bi,bj)=0. _d 0 |
90 |
ENDDO |
ENDDO |
91 |
ENDDO |
ENDDO |
92 |
|
|