44 |
# ifdef ALLOW_NONHYDROSTATIC |
# ifdef ALLOW_NONHYDROSTATIC |
45 |
# include "CG3D.h" |
# include "CG3D.h" |
46 |
# endif |
# endif |
47 |
|
# ifdef EXACT_CONSERV |
48 |
|
# include "SURFACE.h" |
49 |
|
# endif |
50 |
#endif |
#endif |
51 |
|
|
52 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
154 |
CADJ STORE vvel = tapelev3, key = ilev_3 |
CADJ STORE vvel = tapelev3, key = ilev_3 |
155 |
CADJ STORE wvel = tapelev3, key = ilev_3 |
CADJ STORE wvel = tapelev3, key = ilev_3 |
156 |
CADJ STORE etan = tapelev3, key = ilev_3 |
CADJ STORE etan = tapelev3, key = ilev_3 |
|
CADJ STORE etanm1 = tapelev3, key = ilev_3 |
|
157 |
CADJ STORE gu = tapelev3, key = ilev_3 |
CADJ STORE gu = tapelev3, key = ilev_3 |
158 |
CADJ STORE gv = tapelev3, key = ilev_3 |
CADJ STORE gv = tapelev3, key = ilev_3 |
159 |
#ifdef INCLUDE_CD_CODE |
#ifdef INCLUDE_CD_CODE |
160 |
|
CADJ STORE etanm1 = tapelev3, key = ilev_3 |
161 |
CADJ STORE uveld = tapelev3, key = ilev_3 |
CADJ STORE uveld = tapelev3, key = ilev_3 |
162 |
CADJ STORE vveld = tapelev3, key = ilev_3 |
CADJ STORE vveld = tapelev3, key = ilev_3 |
163 |
CADJ STORE unm1 = tapelev3, key = ilev_3 |
CADJ STORE unm1 = tapelev3, key = ilev_3 |
169 |
CADJ STORE tr1 = tapelev3, key = ilev_3 |
CADJ STORE tr1 = tapelev3, key = ilev_3 |
170 |
CADJ STORE gtr1nm1 = tapelev3, key = ilev_3 |
CADJ STORE gtr1nm1 = tapelev3, key = ilev_3 |
171 |
#endif |
#endif |
172 |
|
#ifdef EXACT_CONSERV |
173 |
|
CADJ STORE etaH = tapelev3, key = ilev_3 |
174 |
|
CADJ STORE hDivFlow = tapelev3, key = ilev_3 |
175 |
|
#endif |
176 |
|
|
177 |
c-- Initialise storage for the middle loop. |
c-- Initialise storage for the middle loop. |
178 |
CADJ INIT tapelev2 = USER |
CADJ INIT tapelev2 = USER |
189 |
CADJ STORE vvel = tapelev2, key = ilev_2 |
CADJ STORE vvel = tapelev2, key = ilev_2 |
190 |
CADJ STORE wvel = tapelev2, key = ilev_2 |
CADJ STORE wvel = tapelev2, key = ilev_2 |
191 |
CADJ STORE etan = tapelev2, key = ilev_2 |
CADJ STORE etan = tapelev2, key = ilev_2 |
|
CADJ STORE etanm1 = tapelev2, key = ilev_2 |
|
192 |
CADJ STORE gu = tapelev2, key = ilev_2 |
CADJ STORE gu = tapelev2, key = ilev_2 |
193 |
CADJ STORE gv = tapelev2, key = ilev_2 |
CADJ STORE gv = tapelev2, key = ilev_2 |
194 |
#ifdef INCLUDE_CD_CODE |
#ifdef INCLUDE_CD_CODE |
195 |
|
CADJ STORE etanm1 = tapelev2, key = ilev_2 |
196 |
CADJ STORE uveld = tapelev2, key = ilev_2 |
CADJ STORE uveld = tapelev2, key = ilev_2 |
197 |
CADJ STORE vveld = tapelev2, key = ilev_2 |
CADJ STORE vveld = tapelev2, key = ilev_2 |
198 |
CADJ STORE unm1 = tapelev2, key = ilev_2 |
CADJ STORE unm1 = tapelev2, key = ilev_2 |
204 |
CADJ STORE tr1 = tapelev2, key = ilev_2 |
CADJ STORE tr1 = tapelev2, key = ilev_2 |
205 |
CADJ STORE gtr1nm1 = tapelev2, key = ilev_2 |
CADJ STORE gtr1nm1 = tapelev2, key = ilev_2 |
206 |
#endif |
#endif |
207 |
|
#ifdef EXACT_CONSERV |
208 |
|
CADJ STORE etaH = tapelev2, key = ilev_2 |
209 |
|
CADJ STORE hDivFlow = tapelev2, key = ilev_2 |
210 |
|
#endif |
211 |
|
|
212 |
c-- Initialize storage for the innermost loop. |
c-- Initialize storage for the innermost loop. |
213 |
c-- Always check common block sizes for the checkpointing! |
c-- Always check common block sizes for the checkpointing! |