52 VPUBLIC
void VbuildG(
int *nxf,
int *nyf,
int *nzf,
55 double *pcFF,
double *acFF,
double *ac) {
58 MAT2(acFF, *nxf * *nyf * *nzf, 27);
68 RAT2(pcFF, 1, 1), RAT2(pcFF, 1, 2), RAT2(pcFF, 1, 3), RAT2(pcFF, 1, 4), RAT2(pcFF, 1, 5),
69 RAT2(pcFF, 1, 6), RAT2(pcFF, 1, 7), RAT2(pcFF, 1, 8), RAT2(pcFF, 1, 9),
70 RAT2(pcFF, 1, 10), RAT2(pcFF, 1, 11), RAT2(pcFF, 1, 12), RAT2(pcFF, 1, 13), RAT2(pcFF, 1, 14),
71 RAT2(pcFF, 1, 15), RAT2(pcFF, 1, 16), RAT2(pcFF, 1, 17), RAT2(pcFF, 1, 18),
72 RAT2(pcFF, 1, 19), RAT2(pcFF, 1, 20), RAT2(pcFF, 1, 21), RAT2(pcFF, 1, 22), RAT2(pcFF, 1, 23),
73 RAT2(pcFF, 1, 24), RAT2(pcFF, 1, 25), RAT2(pcFF, 1, 26), RAT2(pcFF, 1, 27),
77 RAT2(ac, 1, 1), RAT2(ac, 1, 2), RAT2(ac, 1, 3),
79 RAT2(ac, 1, 5), RAT2(ac, 1, 6),
80 RAT2(ac, 1, 7), RAT2(ac, 1, 8), RAT2(ac, 1, 9), RAT2(ac, 1, 10),
81 RAT2(ac, 1, 11), RAT2(ac, 1, 12), RAT2(ac, 1, 13), RAT2(ac, 1, 14)
85 }
else if (*numdia == 7) {
92 RAT2(pcFF, 1, 1), RAT2(pcFF, 1, 2), RAT2(pcFF, 1, 3), RAT2(pcFF, 1, 4), RAT2(pcFF, 1, 5),
93 RAT2(pcFF, 1, 6), RAT2(pcFF, 1, 7), RAT2(pcFF, 1, 8), RAT2(pcFF, 1, 9),
94 RAT2(pcFF, 1, 10), RAT2(pcFF, 1, 11), RAT2(pcFF, 1, 12), RAT2(pcFF, 1, 13), RAT2(pcFF, 1, 14),
95 RAT2(pcFF, 1, 15), RAT2(pcFF, 1, 16), RAT2(pcFF, 1, 17), RAT2(pcFF, 1, 18),
96 RAT2(pcFF, 1, 19), RAT2(pcFF, 1, 20), RAT2(pcFF, 1, 21), RAT2(pcFF, 1, 22), RAT2(pcFF, 1, 23),
97 RAT2(pcFF, 1, 24), RAT2(pcFF, 1, 25), RAT2(pcFF, 1, 26), RAT2(pcFF, 1, 27),
99 RAT2(acFF, 1, 1), RAT2(acFF, 1, 2), RAT2(acFF, 1, 3), RAT2(acFF, 1, 4),
101 RAT2(ac, 1, 1), RAT2(ac, 1, 2), RAT2(ac, 1, 3),
103 RAT2(ac, 1, 5), RAT2(ac, 1, 6),
104 RAT2(ac, 1, 7), RAT2(ac, 1, 8), RAT2(ac, 1, 9), RAT2(ac, 1, 10),
105 RAT2(ac, 1, 11), RAT2(ac, 1, 12), RAT2(ac, 1, 13), RAT2(ac, 1, 14)
109 }
else if (*numdia == 27) {
116 RAT2(pcFF, 1, 1), RAT2(pcFF, 1, 2), RAT2(pcFF, 1, 3), RAT2(pcFF, 1, 4), RAT2(pcFF, 1, 5),
117 RAT2(pcFF, 1, 6), RAT2(pcFF, 1, 7), RAT2(pcFF, 1, 8), RAT2(pcFF, 1, 9),
118 RAT2(pcFF, 1, 10), RAT2(pcFF, 1, 11), RAT2(pcFF, 1, 12), RAT2(pcFF, 1, 13), RAT2(pcFF, 1, 14),
119 RAT2(pcFF, 1, 15), RAT2(pcFF, 1, 16), RAT2(pcFF, 1, 17), RAT2(pcFF, 1, 18),
120 RAT2(pcFF, 1, 19), RAT2(pcFF, 1, 20), RAT2(pcFF, 1, 21), RAT2(pcFF, 1, 22), RAT2(pcFF, 1, 23),
121 RAT2(pcFF, 1, 24), RAT2(pcFF, 1, 25), RAT2(pcFF, 1, 26), RAT2(pcFF, 1, 27),
123 RAT2(acFF, 1, 1), RAT2(acFF, 1, 2), RAT2(acFF, 1, 3), RAT2(acFF, 1, 4),
124 RAT2(acFF, 1, 5), RAT2(acFF, 1, 6), RAT2(acFF, 1, 7), RAT2(acFF, 1, 8), RAT2(acFF, 1, 9),
125 RAT2(acFF, 1, 10), RAT2(acFF, 1, 11), RAT2(acFF, 1, 12), RAT2(acFF, 1, 13), RAT2(acFF, 1, 14),
127 RAT2(ac, 1, 1), RAT2(ac, 1, 2), RAT2(ac, 1, 3),
129 RAT2(ac, 1, 5), RAT2(ac, 1, 6),
130 RAT2(ac, 1, 7), RAT2(ac, 1, 8), RAT2(ac, 1, 9), RAT2(ac, 1, 10),
131 RAT2(ac, 1, 11), RAT2(ac, 1, 12), RAT2(ac, 1, 13), RAT2(ac, 1, 14)
136 Vnm_print(2,
"BUILDG: invalid stencil type given...\n");
141 int *
nx,
int *
ny,
int *
nz,
142 double *oPC,
double *oPN,
double *oPS,
double *oPE,
double *oPW,
143 double *oPNE,
double *oPNW,
double *oPSE,
double *oPSW,
144 double *uPC,
double *uPN,
double *uPS,
double *uPE,
double *uPW,
145 double *uPNE,
double *uPNW,
double *uPSE,
double *uPSW,
146 double *dPC,
double *dPN,
double *dPS,
double *dPE,
double *dPW,
147 double *dPNE,
double *dPNW,
double *dPSE,
double *dPSW,
149 double *XoC,
double *XoE,
double *XoN,
151 double *XoNE,
double *XoNW,
152 double *XuE,
double *XuW,
double *XuN,
double *XuS,
153 double *XuNE,
double *XuNW,
double *XuSE,
double *XuSW
156 int i, j, k, ii, jj, kk;
157 int im1, ip1, jm1, jp1, km1, kp1;
158 int iim1, iip1, jjm1, jjp1, kkm1, kkp1;
159 int nxm1, nym1, nzm1;
161 double TMP1_XOC, TMP2_XOC, TMP3_XOC;
163 MAT3( oC, *nxf, *nyf, *nzf);
165 MAT3( XoC, *
nx, *
ny, *
nz);
166 MAT3( XoE, *
nx, *
ny, *
nz);
167 MAT3( XoN, *
nx, *
ny, *
nz);
169 MAT3( XuC, *
nx, *
ny, *
nz);
171 MAT3(XoNE, *
nx, *
ny, *
nz);
172 MAT3(XoNW, *
nx, *
ny, *
nz);
174 MAT3( XuE, *
nx, *
ny, *
nz);
175 MAT3( XuW, *
nx, *
ny, *
nz);
176 MAT3( XuN, *
nx, *
ny, *
nz);
177 MAT3( XuS, *
nx, *
ny, *
nz);
178 MAT3(XuNE, *
nx, *
ny, *
nz);
179 MAT3(XuNW, *
nx, *
ny, *
nz);
180 MAT3(XuSE, *
nx, *
ny, *
nz);
181 MAT3(XuSW, *
nx, *
ny, *
nz);
183 MAT3( oPC, *
nx, *
ny, *
nz);
184 MAT3( oPN, *
nx, *
ny, *
nz);
185 MAT3( oPS, *
nx, *
ny, *
nz);
186 MAT3( oPE, *
nx, *
ny, *
nz);
187 MAT3( oPW, *
nx, *
ny, *
nz);
188 MAT3(oPNE, *
nx, *
ny, *
nz);
189 MAT3(oPNW, *
nx, *
ny, *
nz);
190 MAT3(oPSE, *
nx, *
ny, *
nz);
191 MAT3(oPSW, *
nx, *
ny, *
nz);
193 MAT3( uPC, *
nx, *
ny, *
nz);
194 MAT3( uPN, *
nx, *
ny, *
nz);
195 MAT3( uPS, *
nx, *
ny, *
nz);
196 MAT3( uPE, *
nx, *
ny, *
nz);
197 MAT3( uPW, *
nx, *
ny, *
nz);
198 MAT3(uPNE, *
nx, *
ny, *
nz);
199 MAT3(uPNW, *
nx, *
ny, *
nz);
200 MAT3(uPSE, *
nx, *
ny, *
nz);
201 MAT3(uPSW, *
nx, *
ny, *
nz);
203 MAT3( dPC, *
nx, *
ny, *
nz);
204 MAT3( dPN, *
nx, *
ny, *
nz);
205 MAT3( dPS, *
nx, *
ny, *
nz);
206 MAT3( dPE, *
nx, *
ny, *
nz);
207 MAT3( dPW, *
nx, *
ny, *
nz);
208 MAT3(dPNE, *
nx, *
ny, *
nz);
209 MAT3(dPNW, *
nx, *
ny, *
nz);
210 MAT3(dPSE, *
nx, *
ny, *
nz);
211 MAT3(dPSW, *
nx, *
ny, *
nz);
221 for(kk=2; kk<=*
nz-1; kk++) {
224 for(jj=2; jj<=*
ny-1; jj++) {
227 for(ii=2; ii<=*
nx-1; i++) {
250 VAT3( uPS, ii,jj,kk) * VAT3( uPS, ii,jj,kk) * VAT3( oC, i,jm1,kp1)
251 + VAT3(dPSW, ii,jj,kk) * VAT3(dPSW, ii,jj,kk) * VAT3( oC, im1,jm1,km1)
252 + VAT3(oPSW, ii,jj,kk) * VAT3(oPSW, ii,jj,kk) * VAT3( oC, im1,jm1,k)
253 + VAT3(uPSW, ii,jj,kk) * VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1)
254 + VAT3( dPW, ii,jj,kk) * VAT3( dPW, ii,jj,kk) * VAT3( oC, im1,j,km1)
255 + VAT3( oPW, ii,jj,kk) * VAT3( oPW, ii,jj,kk) * VAT3( oC, im1,j,k)
256 + VAT3(uPNW, ii,jj,kk) * VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1)
257 + VAT3( dPS, ii,jj,kk) * VAT3( dPS, ii,jj,kk) * VAT3( oC, i,jm1,km1)
258 + VAT3( oPS, ii,jj,kk) * VAT3( oPS, ii,jj,kk) * VAT3( oC, i,jm1,k);
261 VAT3( dPC, ii,jj,kk) * VAT3( dPC, ii,jj,kk) * VAT3( oC, i,j,km1)
262 + VAT3( oPC, ii,jj,kk) * VAT3( oPC, ii,jj,kk) * VAT3( oC, i,j,k)
263 + VAT3( uPC, ii,jj,kk) * VAT3( uPC, ii,jj,kk) * VAT3( oC, i,j,kp1)
264 + VAT3( dPN, ii,jj,kk) * VAT3( dPN, ii,jj,kk) * VAT3( oC, i,jp1,km1)
265 + VAT3( oPN, ii,jj,kk) * VAT3( oPN, ii,jj,kk) * VAT3( oC, i,jp1,k)
266 + VAT3( uPW, ii,jj,kk) * VAT3( uPW, ii,jj,kk) * VAT3( oC, im1,j,kp1)
267 + VAT3(dPNW, ii,jj,kk) * VAT3(dPNW, ii,jj,kk) * VAT3( oC, im1,jp1,km1)
268 + VAT3(oPNW, ii,jj,kk) * VAT3(oPNW, ii,jj,kk) * VAT3( oC, im1,jp1,k)
269 + VAT3( oPE, ii,jj,kk) * VAT3( oPE, ii,jj,kk) * VAT3( oC, ip1,j,k);
272 VAT3( uPE, ii,jj,kk) * VAT3( uPE, ii,jj,kk) * VAT3( oC, ip1,j,kp1)
273 + VAT3(dPNE, ii,jj,kk) * VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1)
274 + VAT3(oPNE, ii,jj,kk) * VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,k)
275 + VAT3(uPNE, ii,jj,kk) * VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1)
276 + VAT3( uPN, ii,jj,kk) * VAT3( uPN, ii,jj,kk) * VAT3( oC, i,jp1,kp1)
277 + VAT3(dPSE, ii,jj,kk) * VAT3(dPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,km1)
278 + VAT3(oPSE, ii,jj,kk) * VAT3(oPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,k)
279 + VAT3(uPSE, ii,jj,kk) * VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1)
280 + VAT3( dPE, ii,jj,kk) * VAT3( dPE, ii,jj,kk) * VAT3( oC, ip1,j,km1);
282 VAT3( XoC, ii,jj,kk) = TMP1_XOC + TMP2_XOC + TMP3_XOC;
290 VAT3( XoE, ii,jj,kk) =
291 - VAT3(dPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
292 - VAT3(oPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
293 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
294 - VAT3( dPE, ii,jj,kk) * VAT3( oC, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
295 - VAT3( oPE, ii,jj,kk) * VAT3( oC, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
296 - VAT3( uPE, ii,jj,kk) * VAT3( oC, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
297 - VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
298 - VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
299 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPNW, iip1,jj,kk);
307 VAT3( XoN, ii,jj,kk) =
308 - VAT3(dPNW, ii,jj,kk) * VAT3( oC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
309 - VAT3(oPNW, ii,jj,kk) * VAT3( oC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
310 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
311 - VAT3( dPN, ii,jj,kk) * VAT3( oC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
312 - VAT3( oPN, ii,jj,kk) * VAT3( oC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
313 - VAT3( uPN, ii,jj,kk) * VAT3( oC, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
314 - VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
315 - VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1, k) * VAT3(oPSE, ii,jjp1,kk)
316 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk);
324 VAT3( XuC, ii,jj,kk) =
325 - VAT3(dPSW, ii,jj,kkp1) * VAT3( oC, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
326 - VAT3( dPW, ii,jj,kkp1) * VAT3( oC, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
327 - VAT3(dPNW, ii,jj,kkp1) * VAT3( oC, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
328 - VAT3( dPS, ii,jj,kkp1) * VAT3( oC, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
329 - VAT3( dPC, ii,jj,kkp1) * VAT3( oC, i,j,kp1) * VAT3( uPC, ii,jj,kk)
330 - VAT3( dPN, ii,jj,kkp1) * VAT3( oC, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
331 - VAT3(dPSE, ii,jj,kkp1) * VAT3( oC, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
332 - VAT3( dPE, ii,jj,kkp1) * VAT3( oC, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
333 - VAT3(dPNE, ii,jj,kkp1) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPNE, ii,jj,kk);
341 VAT3(XoNE, ii,jj,kk) =
342 - VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
343 - VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
344 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk);
353 VAT3(XoNW, ii,jj,kk) =
354 - VAT3(dPNW, ii,jj,kk) * VAT3( oC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
355 - VAT3(oPNW, ii,jj,kk) * VAT3( oC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
356 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk);
364 VAT3( XuE, ii,jj,kk) =
365 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
366 - VAT3( uPE, ii,jj,kk) * VAT3( oC, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
367 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1);
375 VAT3( XuW, ii,jj,kk) =
376 - VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
377 - VAT3( uPW, ii,jj,kk) * VAT3( oC, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
378 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1);
386 VAT3( XuN, ii,jj,kk) =
387 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
388 - VAT3( uPN, ii,jj,kk) * VAT3( oC, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
389 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1);
397 VAT3( XuS, ii,jj,kk) =
398 - VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
399 - VAT3( uPS, ii,jj,kk) * VAT3( oC, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
400 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1);
408 VAT3(XuNE, ii,jj,kk) =
409 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1);
417 VAT3(XuNW, ii,jj,kk) =
418 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1);
426 VAT3(XuSE, ii,jj,kk) =
427 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1);
435 VAT3(XuSW, ii,jj,kk) =
436 - VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1);
446 int *
nx,
int *
ny,
int *
nz,
447 double *oPC,
double *oPN,
double *oPS,
double *oPE,
double *oPW,
448 double *oPNE,
double *oPNW,
double *oPSE,
double *oPSW,
449 double *uPC,
double *uPN,
double *uPS,
double *uPE,
double *uPW,
450 double *uPNE,
double *uPNW,
double *uPSE,
double *uPSW,
451 double *dPC,
double *dPN,
double *dPS,
double *dPE,
double *dPW,
452 double *dPNE,
double *dPNW,
double *dPSE,
double *dPSW,
453 double *oC,
double *oE,
double *oN,
double *uC,
454 double *XoC,
double *XoE,
double *XoN,
456 double *XoNE,
double *XoNW,
457 double *XuE,
double *XuW,
double *XuN,
double *XuS,
458 double *XuNE,
double *XuNW,
double *XuSE,
double *XuSW) {
462 int im1, ip1, im2, ip2;
463 int jm1, jp1, jm2, jp2;
464 int km1, kp1, km2, kp2;
468 int nxm1, nym1, nzm1;
470 double TMP1_XOC, TMP2_XOC, TMP3_XOC, TMP4_XOC;
471 double TMP5_XOC, TMP6_XOC, TMP7_XOC, TMP8_XOC;
473 double TMP1_XOE, TMP2_XOE, TMP3_XOE, TMP4_XOE;
474 double TMP1_XON, TMP2_XON, TMP3_XON, TMP4_XON;
475 double TMP1_XUC, TMP2_XUC, TMP3_XUC, TMP4_XUC;
477 MAT3( oC, *nxf, *nyf, *nzf);
478 MAT3( oE, *nxf, *nyf, *nzf);
479 MAT3( oN, *nxf, *nyf, *nzf);
481 MAT3( uC, *nxf, *nyf, *nzf);
483 MAT3( XoC, *
nx, *
ny, *
nz);
484 MAT3( XoE, *
nx, *
ny, *
nz);
485 MAT3( XoN, *
nx, *
ny, *
nz);
486 MAT3(XoNE, *
nx, *
ny, *
nz);
487 MAT3(XoNW, *
nx, *
ny, *
nz);
489 MAT3( XuC, *
nx, *
ny, *
nz);
490 MAT3( XuE, *
nx, *
ny, *
nz);
491 MAT3( XuW, *
nx, *
ny, *
nz);
492 MAT3( XuN, *
nx, *
ny, *
nz);
493 MAT3( XuS, *
nx, *
ny, *
nz);
494 MAT3(XuNE, *
nx, *
ny, *
nz);
495 MAT3(XuNW, *
nx, *
ny, *
nz);
496 MAT3(XuSE, *
nx, *
ny, *
nz);
497 MAT3(XuSW, *
nx, *
ny, *
nz);
499 MAT3( oPC, *
nx, *
ny, *
nz);
500 MAT3( oPN, *
nx, *
ny, *
nz);
501 MAT3( oPS, *
nx, *
ny, *
nz);
502 MAT3( oPE, *
nx, *
ny, *
nz);
503 MAT3( oPW, *
nx, *
ny, *
nz);
504 MAT3(oPNE, *
nx, *
ny, *
nz);
505 MAT3(oPNW, *
nx, *
ny, *
nz);
506 MAT3(oPSE, *
nx, *
ny, *
nz);
507 MAT3(oPSW, *
nx, *
ny, *
nz);
509 MAT3( uPC, *
nx, *
ny, *
nz);
510 MAT3( uPN, *
nx, *
ny, *
nz);
511 MAT3( uPS, *
nx, *
ny, *
nz);
512 MAT3( uPE, *
nx, *
ny, *
nz);
513 MAT3( uPW, *
nx, *
ny, *
nz);
514 MAT3(uPNE, *
nx, *
ny, *
nz);
515 MAT3(uPNW, *
nx, *
ny, *
nz);
516 MAT3(uPSE, *
nx, *
ny, *
nz);
517 MAT3(uPSW, *
nx, *
ny, *
nz);
519 MAT3( dPC, *
nx, *
ny, *
nz);
520 MAT3( dPN, *
nx, *
ny, *
nz);
521 MAT3( dPS, *
nx, *
ny, *
nz);
522 MAT3( dPE, *
nx, *
ny, *
nz);
523 MAT3( dPW, *
nx, *
ny, *
nz);
524 MAT3(dPNE, *
nx, *
ny, *
nz);
525 MAT3(dPNW, *
nx, *
ny, *
nz);
526 MAT3(dPSE, *
nx, *
ny, *
nz);
527 MAT3(dPSW, *
nx, *
ny, *
nz);
537 for(kk=2; kk<=*
nz-1; kk++) {
540 for(jj=2; jj<=*
ny-1; jj++) {
543 for(ii=2; ii<=*
nx-1; ii++) {
572 VAT3(dPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
573 - VAT3( uC, im1,jm1,km1) * VAT3(oPSW, ii,jj,kk)
574 - VAT3( oN, im1,jm1,km1) * VAT3( dPW, ii,jj,kk)
575 - VAT3( oE, im1,jm1,km1) * VAT3( dPS, ii,jj,kk))
577 + VAT3(oPNE, ii,jj,kk) * (- VAT3( oE, i,jp1, k) * VAT3( oPN, ii,jj,kk)
578 - VAT3( oN, ip1, j, k) * VAT3( oPE, ii,jj,kk)
579 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
580 + VAT3( oC, ip1,jp1, k) * VAT3(oPNE, ii,jj,kk)
581 - VAT3( uC, ip1,jp1, k) * VAT3(uPNE, ii,jj,kk))
583 + VAT3(dPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
584 - VAT3( oN, ip1, j,km1) * VAT3( dPE, ii,jj,kk)
585 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
586 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNE, ii,jj,kk));
589 VAT3(dPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
590 + VAT3( oC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
591 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSE, ii,jj,kk)
592 - VAT3( oN, ip1,jm1,km1) * VAT3( dPE, ii,jj,kk))
594 + VAT3( uPE, ii,jj,kk) * (- VAT3( oE, i, j,kp1) * VAT3( uPC, ii,jj,kk)
595 - VAT3( oN, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
596 - VAT3( uC, ip1, j, k) * VAT3( oPE, ii,jj,kk)
597 + VAT3( oC, ip1, j,kp1) * VAT3( uPE, ii,jj,kk)
598 - VAT3( oN, ip1, j,kp1) * VAT3(uPNE, ii,jj,kk))
600 + VAT3( oPE, ii,jj,kk) * (- VAT3( oE, i, j, k) * VAT3( oPC, ii,jj,kk)
601 - VAT3( oN, ip1,jm1, k) * VAT3(oPSE, ii,jj,kk)
602 - VAT3( uC, ip1, j,km1) * VAT3( dPE, ii,jj,kk)
603 + VAT3( oC, ip1, j, k) * VAT3( oPE, ii,jj,kk)
604 - VAT3( uC, ip1, j, k) * VAT3( uPE, ii,jj,kk)
605 - VAT3( oN, ip1, j, k) * VAT3(oPNE, ii,jj,kk));
608 + VAT3( dPE, ii,jj,kk) * (- VAT3( oE, i, j,km1) * VAT3( dPC, ii,jj,kk)
609 - VAT3( oN, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
610 + VAT3( oC, ip1, j,km1) * VAT3( dPE, ii,jj,kk)
611 - VAT3( uC, ip1, j,km1) * VAT3( oPE, ii,jj,kk)
612 - VAT3( oN, ip1, j,km1) * VAT3(dPNE, ii,jj,kk))
614 + VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
615 - VAT3( uC, ip1,jm1, k) * VAT3(oPSE, ii,jj,kk)
616 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
617 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPE, ii,jj,kk))
619 + VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
620 - VAT3( oN, ip1, j,kp1) * VAT3( uPE, ii,jj,kk)
621 - VAT3( uC, ip1,jp1, k) * VAT3(oPNE, ii,jj,kk)
622 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNE, ii,jj,kk));
625 + VAT3( oPS, ii,jj,kk) * (- VAT3( oE, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
626 - VAT3( uC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
627 + VAT3( oC, i,jm1, k) * VAT3( oPS, ii,jj,kk)
628 - VAT3( uC, i,jm1, k) * VAT3( uPS, ii,jj,kk)
629 - VAT3( oN, i,jm1, k) * VAT3( oPC, ii,jj,kk)
630 - VAT3( oE, i,jm1, k) * VAT3(oPSE, ii,jj,kk))
632 + VAT3( dPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
633 + VAT3( oC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
634 - VAT3( uC, i,jm1,km1) * VAT3( oPS, ii,jj,kk)
635 - VAT3( oN, i,jm1,km1) * VAT3( dPC, ii,jj,kk)
636 - VAT3( oE, i,jm1,km1) * VAT3(dPSE, ii,jj,kk))
638 + VAT3(oPSE, ii,jj,kk) * (- VAT3( oE, i,jm1, k) * VAT3( oPS, ii,jj,kk)
639 - VAT3( uC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
640 + VAT3( oC, ip1,jm1, k) * VAT3(oPSE, ii,jj,kk)
641 - VAT3( uC, ip1,jm1, k) * VAT3(uPSE, ii,jj,kk)
642 - VAT3( oN, ip1,jm1, k) * VAT3( oPE, ii,jj,kk));
645 + VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
646 - VAT3( oN, i, j,km1) * VAT3( dPC, ii,jj,kk)
647 + VAT3( oC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
648 - VAT3( uC, i,jp1,km1) * VAT3( oPN, ii,jj,kk)
649 - VAT3( oE, i,jp1,km1) * VAT3(dPNE, ii,jj,kk))
651 + VAT3( uPC, ii,jj,kk) * (- VAT3( oE, im1, j,kp1) * VAT3( uPW, ii,jj,kk)
652 - VAT3( oN, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
653 - VAT3( uC, i, j, k) * VAT3( oPC, ii,jj,kk)
654 + VAT3( oC, i, j,kp1) * VAT3( uPC, ii,jj,kk)
655 - VAT3( oN, i, j,kp1) * VAT3( uPN, ii,jj,kk)
656 - VAT3( oE, i, j,kp1) * VAT3( uPE, ii,jj,kk))
658 + VAT3( oPC, ii,jj,kk) * (- VAT3( oE, im1, j, k) * VAT3( oPW, ii,jj,kk)
659 - VAT3( oN, i,jm1, k) * VAT3( oPS, ii,jj,kk)
660 - VAT3( uC, i, j,km1) * VAT3( dPC, ii,jj,kk)
661 + VAT3( oC, i, j, k) * VAT3( oPC, ii,jj,kk)
662 - VAT3( uC, i, j, k) * VAT3( uPC, ii,jj,kk)
663 - VAT3( oN, i, j, k) * VAT3( oPN, ii,jj,kk)
664 - VAT3( oE, i, j, k) * VAT3( oPE, ii,jj,kk));
667 + VAT3( dPC, ii,jj,kk) * (- VAT3( oE, im1, j,km1) * VAT3( dPW, ii,jj,kk)
668 - VAT3( oN, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
669 + VAT3( oC, i, j,km1) * VAT3( dPC, ii,jj,kk)
670 - VAT3( uC, i, j,km1) * VAT3( oPC, ii,jj,kk)
671 - VAT3( oN, i, j,km1) * VAT3( dPN, ii,jj,kk)
672 - VAT3( oE, i, j,km1) * VAT3( dPE, ii,jj,kk))
674 + VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
675 - VAT3( uC, i,jm1, k) * VAT3( oPS, ii,jj,kk)
676 + VAT3( oC, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
677 - VAT3( oN, i,jm1,kp1) * VAT3( uPC, ii,jj,kk)
678 - VAT3( oE, i,jm1,kp1) * VAT3(uPSE, ii,jj,kk))
680 + VAT3(uPNW, ii,jj,kk) * (- VAT3( oN, im1, j,kp1) * VAT3( uPW, ii,jj,kk)
681 - VAT3( uC, im1,jp1, k) * VAT3(oPNW, ii,jj,kk)
682 + VAT3( oC, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
683 - VAT3( oE, im1,jp1,kp1) * VAT3( uPN, ii,jj,kk));
686 + VAT3( dPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
687 + VAT3( oC, im1, j,km1) * VAT3( dPW, ii,jj,kk)
688 - VAT3( uC, im1, j,km1) * VAT3( oPW, ii,jj,kk)
689 - VAT3( oN, im1, j,km1) * VAT3(dPNW, ii,jj,kk)
690 - VAT3( oE, im1, j,km1) * VAT3( dPC, ii,jj,kk))
692 + VAT3(uPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
693 + VAT3( oC, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
694 - VAT3( oN, im1,jm1,kp1) * VAT3( uPW, ii,jj,kk)
695 - VAT3( oE, im1,jm1,kp1) * VAT3( uPS, ii,jj,kk))
697 + VAT3(oPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
698 + VAT3( oC, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
699 - VAT3( uC, im1,jm1, k) * VAT3(uPSW, ii,jj,kk)
700 - VAT3( oN, im1,jm1, k) * VAT3( oPW, ii,jj,kk)
701 - VAT3( oE, im1,jm1, k) * VAT3( oPS, ii,jj,kk));
704 + VAT3(oPNW, ii,jj,kk) * (- VAT3( oN, im1, j, k) * VAT3( oPW, ii,jj,kk)
705 - VAT3( uC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
706 + VAT3( oC, im1,jp1, k) * VAT3(oPNW, ii,jj,kk)
707 - VAT3( uC, im1,jp1, k) * VAT3(uPNW, ii,jj,kk)
708 - VAT3( oE, im1,jp1, k) * VAT3( oPN, ii,jj,kk))
710 + VAT3(dPNW, ii,jj,kk) * (- VAT3( oN, im1, j,km1) * VAT3( dPW, ii,jj,kk)
711 + VAT3( oC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
712 - VAT3( uC, im1,jp1,km1) * VAT3(oPNW, ii,jj,kk)
713 - VAT3( oE, im1,jp1,km1) * VAT3( dPN, ii,jj,kk))
715 + VAT3( oPW, ii,jj,kk) * (- VAT3( oN, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
716 - VAT3( uC, im1, j,km1) * VAT3( dPW, ii,jj,kk)
717 + VAT3( oC, im1, j, k) * VAT3( oPW, ii,jj,kk)
718 - VAT3( uC, im1, j, k) * VAT3( uPW, ii,jj,kk)
719 - VAT3( oN, im1, j, k) * VAT3(oPNW, ii,jj,kk)
720 - VAT3( oE, im1, j, k) * VAT3( oPC, ii,jj,kk));
723 + VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
724 - VAT3( uC, im1, j, k) * VAT3( oPW, ii,jj,kk)
725 + VAT3( oC, im1, j,kp1) * VAT3( uPW, ii,jj,kk)
726 - VAT3( oN, im1, j,kp1) * VAT3(uPNW, ii,jj,kk)
727 - VAT3( oE, im1, j,kp1) * VAT3( uPC, ii,jj,kk))
729 + VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
730 - VAT3( oN, i, j,kp1) * VAT3( uPC, ii,jj,kk)
731 - VAT3( uC, i,jp1, k) * VAT3( oPN, ii,jj,kk)
732 + VAT3( oC, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
733 - VAT3( oE, i,jp1,kp1) * VAT3(uPNE, ii,jj,kk))
735 + VAT3( oPN, ii,jj,kk) * (- VAT3( oE, im1,jp1, k) * VAT3(oPNW, ii,jj,kk)
736 - VAT3( oN, i, j, k) * VAT3( oPC, ii,jj,kk)
737 - VAT3( uC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
738 + VAT3( oC, i,jp1, k) * VAT3( oPN, ii,jj,kk)
739 - VAT3( uC, i,jp1, k) * VAT3( uPN, ii,jj,kk)
740 - VAT3( oE, i,jp1, k) * VAT3(oPNE, ii,jj,kk));
742 VAT3( XoC, ii,jj,kk) = TMP1_XOC + TMP2_XOC + TMP3_XOC
743 + TMP4_XOC + TMP5_XOC + TMP6_XOC
744 + TMP7_XOC + TMP8_XOC + TMP9_XOC;
754 VAT3( dPS, ii,jj,kk) * VAT3( oE, i,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
755 + VAT3( oPS, ii,jj,kk) * VAT3( oE, i,jm1, k) * VAT3(oPSW, iip1,jj,kk)
756 + VAT3( uPS, ii,jj,kk) * VAT3( oE, i,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
757 + VAT3( dPC, ii,jj,kk) * VAT3( oE, i, j,km1) * VAT3( dPW, iip1,jj,kk)
758 + VAT3( oPC, ii,jj,kk) * VAT3( oE, i, j, k) * VAT3( oPW, iip1,jj,kk)
759 + VAT3( uPC, ii,jj,kk) * VAT3( oE, i, j,kp1) * VAT3( uPW, iip1,jj,kk)
760 + VAT3( dPN, ii,jj,kk) * VAT3( oE, i,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
761 + VAT3( oPN, ii,jj,kk) * VAT3( oE, i,jp1, k) * VAT3(oPNW, iip1,jj,kk)
762 + VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(uPNW, iip1,jj,kk)
764 - VAT3(dPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
765 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSW, iip1,jj,kk)
766 - VAT3( oN, ip1,jm1,km1) * VAT3( dPW, iip1,jj,kk)
767 - VAT3( oE, ip1,jm1,km1) * VAT3( dPS, iip1,jj,kk));
770 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
771 + VAT3( oC, ip1,jm1, k) * VAT3(oPSW, iip1,jj,kk)
772 - VAT3( uC, ip1,jm1, k) * VAT3(uPSW, iip1,jj,kk)
773 - VAT3( oN, ip1,jm1, k) * VAT3( oPW, iip1,jj,kk)
774 - VAT3( oE, ip1,jm1, k) * VAT3( oPS, iip1,jj,kk))
776 - VAT3(uPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1, k) * VAT3(oPSW, iip1,jj,kk)
777 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
778 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPW, iip1,jj,kk)
779 - VAT3( oE, ip1,jm1,kp1) * VAT3( uPS, iip1,jj,kk))
781 - VAT3( dPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
782 + VAT3( oC, ip1, j,km1) * VAT3( dPW, iip1,jj,kk)
783 - VAT3( uC, ip1, j,km1) * VAT3( oPW, iip1,jj,kk)
784 - VAT3( oN, ip1, j,km1) * VAT3(dPNW, iip1,jj,kk)
785 - VAT3( oE, ip1, j,km1) * VAT3( dPC, iip1,jj,kk));
788 - VAT3( oPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1, k) * VAT3(oPSW, iip1,jj,kk)
789 - VAT3( uC, ip1, j,km1) * VAT3( dPW, iip1,jj,kk)
790 + VAT3( oC, ip1, j, k) * VAT3( oPW, iip1,jj,kk)
791 - VAT3( uC, ip1, j, k) * VAT3( uPW, iip1,jj,kk)
792 - VAT3( oN, ip1, j, k) * VAT3(oPNW, iip1,jj,kk)
793 - VAT3( oE, ip1, j, k) * VAT3( oPC, iip1,jj,kk))
795 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
796 - VAT3( uC, ip1, j, k) * VAT3( oPW, iip1,jj,kk)
797 + VAT3( oC, ip1, j,kp1) * VAT3( uPW, iip1,jj,kk)
798 - VAT3( oN, ip1, j,kp1) * VAT3(uPNW, iip1,jj,kk)
799 - VAT3( oE, ip1, j,kp1) * VAT3( uPC, iip1,jj,kk))
801 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oN, ip1, j,km1) * VAT3( dPW, iip1,jj,kk)
802 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
803 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNW, iip1,jj,kk)
804 - VAT3( oE, ip1,jp1,km1) * VAT3( dPN, iip1,jj,kk));
807 - VAT3(oPNE, ii,jj,kk) * (- VAT3( oN, ip1, j, k) * VAT3( oPW, iip1,jj,kk)
808 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
809 + VAT3( oC, ip1,jp1, k) * VAT3(oPNW, iip1,jj,kk)
810 - VAT3( uC, ip1,jp1, k) * VAT3(uPNW, iip1,jj,kk)
811 - VAT3( oE, ip1,jp1, k) * VAT3( oPN, iip1,jj,kk))
813 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oN, ip1, j,kp1) * VAT3( uPW, iip1,jj,kk)
814 - VAT3( uC, ip1,jp1, k) * VAT3(oPNW, iip1,jj,kk)
815 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNW, iip1,jj,kk)
816 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPN, iip1,jj,kk));
818 VAT3( XoE, ii,jj,kk) = TMP1_XOE + TMP2_XOE + TMP3_XOE + TMP4_XOE;
830 VAT3( dPW, ii,jj,kk) * VAT3( oN, im1,j,km1) * VAT3(dPSW, ii,jjp1,kk)
831 + VAT3( oPW, ii,jj,kk) * VAT3( oN, im1,j,k) * VAT3(oPSW, ii,jjp1,kk)
832 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(uPSW, ii,jjp1,kk)
834 - VAT3(dPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
835 - VAT3( uC, im1,jp1,km1) * VAT3(oPSW, ii,jjp1,kk)
836 - VAT3( oN, im1,jp1,km1) * VAT3( dPW, ii,jjp1,kk)
837 - VAT3( oE, im1,jp1,km1) * VAT3( dPS, ii,jjp1,kk))
839 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
840 + VAT3( oC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
841 - VAT3( uC, im1,jp1,k) * VAT3(uPSW, ii,jjp1,kk)
842 - VAT3( oN, im1,jp1,k) * VAT3( oPW, ii,jjp1,kk)
843 - VAT3( oE, im1,jp1,k) * VAT3( oPS, ii,jjp1,kk));
846 - VAT3(uPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
847 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
848 - VAT3( oN, im1,jp1,kp1) * VAT3( uPW, ii,jjp1,kk)
849 - VAT3( oE, im1,jp1,kp1) * VAT3( uPS, ii,jjp1,kk))
851 + VAT3( dPC, ii,jj,kk) * VAT3( oN, i,j,km1) * VAT3( dPS, ii,jjp1,kk)
852 + VAT3( oPC, ii,jj,kk) * VAT3( oN, i,j,k) * VAT3( oPS, ii,jjp1,kk)
853 + VAT3( uPC, ii,jj,kk) * VAT3( oN, i,j,kp1) * VAT3( uPS, ii,jjp1,kk)
855 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
856 + VAT3( oC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
857 - VAT3( uC, i,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
858 - VAT3( oN, i,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
859 - VAT3( oE, i,jp1,km1) * VAT3(dPSE, ii,jjp1,kk));
862 - VAT3( oPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
863 - VAT3( uC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
864 + VAT3( oC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
865 - VAT3( uC, i,jp1,k) * VAT3( uPS, ii,jjp1,kk)
866 - VAT3( oN, i,jp1,k) * VAT3( oPC, ii,jjp1,kk)
867 - VAT3( oE, i,jp1,k) * VAT3(oPSE, ii,jjp1,kk))
869 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
870 - VAT3( uC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
871 + VAT3( oC, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
872 - VAT3( oN, i,jp1,kp1) * VAT3( uPC, ii,jjp1,kk)
873 - VAT3( oE, i,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk))
875 + VAT3( dPE, ii,jj,kk) * VAT3( oN, ip1,j,km1) * VAT3(dPSE, ii,jjp1,kk)
876 + VAT3( oPE, ii,jj,kk) * VAT3( oN, ip1,j,k) * VAT3(oPSE, ii,jjp1,kk)
877 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(uPSE, ii,jjp1,kk);
880 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
881 + VAT3( oC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
882 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSE, ii,jjp1,kk)
883 - VAT3( oN, ip1,jp1,km1) * VAT3( dPE, ii,jjp1,kk))
885 - VAT3(oPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
886 - VAT3( uC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
887 + VAT3( oC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
888 - VAT3( uC, ip1,jp1,k) * VAT3(uPSE, ii,jjp1,kk)
889 - VAT3( oN, ip1,jp1,k) * VAT3( oPE, ii,jjp1,kk))
891 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
892 - VAT3( uC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
893 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk)
894 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPE, ii,jjp1,kk));
896 VAT3( XoN, ii,jj,kk) = TMP1_XON + TMP2_XON + TMP3_XON + TMP4_XON;
906 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPSW, ii,jj,kkp1)
908 - VAT3(uPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
909 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSW, ii,jj,kkp1)
910 - VAT3( oN, im1,jm1,kp1) * VAT3( dPW, ii,jj,kkp1)
911 - VAT3( oE, im1,jm1,kp1) * VAT3( dPS, ii,jj,kkp1))
913 + VAT3( oPW, ii,jj,kk) * VAT3( uC, im1,j,k) * VAT3( dPW, ii,jj,kkp1)
915 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
916 + VAT3( oC, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
917 - VAT3( uC, im1,j,kp1) * VAT3( oPW, ii,jj,kkp1)
918 - VAT3( oN, im1,j,kp1) * VAT3(dPNW, ii,jj,kkp1)
919 - VAT3( oE, im1,j,kp1) * VAT3( dPC, ii,jj,kkp1))
921 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPNW, ii,jj,kkp1);
924 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oN, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
925 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
926 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNW, ii,jj,kkp1)
927 - VAT3( oE, im1,jp1,kp1) * VAT3( dPN, ii,jj,kkp1))
929 + VAT3( oPS, ii,jj,kk) * VAT3( uC, i,jm1,k) * VAT3( dPS, ii,jj,kkp1)
931 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
932 + VAT3( oC, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
933 - VAT3( uC, i,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
934 - VAT3( oN, i,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
935 - VAT3( oE, i,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1))
937 + VAT3( oPC, ii,jj,kk) * VAT3( uC, i,j,k) * VAT3( dPC, ii,jj,kkp1)
939 - VAT3( uPC, ii,jj,kk) * (- VAT3( oE, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
940 - VAT3( oN, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
941 + VAT3( oC, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
942 - VAT3( uC, i,j,kp1) * VAT3( oPC, ii,jj,kkp1)
943 - VAT3( oN, i,j,kp1) * VAT3( dPN, ii,jj,kkp1)
944 - VAT3( oE, i,j,kp1) * VAT3( dPE, ii,jj,kkp1));
947 + VAT3( oPN, ii,jj,kk) * VAT3( uC, i,jp1,k) * VAT3( dPN, ii,jj,kkp1)
949 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
950 - VAT3( oN, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
951 + VAT3( oC, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
952 - VAT3( uC, i,jp1,kp1) * VAT3( oPN, ii,jj,kkp1)
953 - VAT3( oE, i,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1))
955 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPSE, ii,jj,kkp1)
957 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
958 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
959 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSE, ii,jj,kkp1)
960 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPE, ii,jj,kkp1));
963 + VAT3( oPE, ii,jj,kk) * VAT3( uC, ip1,j,k) * VAT3( dPE, ii,jj,kkp1)
965 - VAT3( uPE, ii,jj,kk) * (- VAT3( oE, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
966 - VAT3( oN, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
967 + VAT3( oC, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
968 - VAT3( uC, ip1,j,kp1) * VAT3( oPE, ii,jj,kkp1)
969 - VAT3( oN, ip1,j,kp1) * VAT3(dPNE, ii,jj,kkp1))
971 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPNE, ii,jj,kkp1)
973 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
974 - VAT3( oN, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
975 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1)
976 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNE, ii,jj,kkp1));
978 VAT3( XuC, ii,jj,kk) = TMP1_XUC + TMP2_XUC + TMP3_XUC + TMP4_XUC;
986 VAT3(XoNE, ii,jj,kk) =
987 VAT3( dPN, ii,jj,kk) * VAT3( oE, i,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
988 + VAT3( oPN, ii,jj,kk) * VAT3( oE, i,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
989 + VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
990 + VAT3( dPE, ii,jj,kk) * VAT3( oN, ip1,j,km1) * VAT3(dPSW, iip1,jjp1,kk)
991 + VAT3( oPE, ii,jj,kk) * VAT3( oN, ip1,j,k) * VAT3(oPSW, iip1,jjp1,kk)
992 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(uPSW, iip1,jjp1,kk)
994 - VAT3(dPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
995 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSW, iip1,jjp1,kk)
996 - VAT3( oN, ip1,jp1,km1) * VAT3( dPW, iip1,jjp1,kk)
997 - VAT3( oE, ip1,jp1,km1) * VAT3( dPS, iip1,jjp1,kk))
999 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
1000 + VAT3( oC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
1001 - VAT3( uC, ip1,jp1,k) * VAT3(uPSW, iip1,jjp1,kk)
1002 - VAT3( oN, ip1,jp1,k) * VAT3( oPW, iip1,jjp1,kk)
1003 - VAT3( oE, ip1,jp1,k) * VAT3( oPS, iip1,jjp1,kk))
1005 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
1006 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
1007 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPW, iip1,jjp1,kk)
1008 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPS, iip1,jjp1,kk));
1016 VAT3(XoNW, ii,jj,kk) =
1017 VAT3( dPW, ii,jj,kk) * VAT3( oN, im1,j,km1) * VAT3(dPSE, iim1,jjp1,kk)
1018 + VAT3( oPW, ii,jj,kk) * VAT3( oN, im1,j,k) * VAT3(oPSE, iim1,jjp1,kk)
1019 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(uPSE, iim1,jjp1,kk)
1021 - VAT3(dPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
1022 + VAT3( oC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
1023 - VAT3( uC, im1,jp1,km1) * VAT3(oPSE, iim1,jjp1,kk)
1024 - VAT3( oN, im1,jp1,km1) * VAT3( dPE, iim1,jjp1,kk))
1026 - VAT3(oPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
1027 - VAT3( uC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
1028 + VAT3( oC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
1029 - VAT3( uC, im1,jp1,k) * VAT3(uPSE, iim1,jjp1,kk)
1030 - VAT3( oN, im1,jp1,k) * VAT3( oPE, iim1,jjp1,kk))
1032 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( uPS, iim1,jjp1,kk)
1033 - VAT3( uC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
1034 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk)
1035 - VAT3( oN, im1,jp1,kp1) * VAT3( uPE, iim1,jjp1,kk))
1037 + VAT3( dPN, ii,jj,kk) * VAT3( oE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
1038 + VAT3( oPN, ii,jj,kk) * VAT3( oE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
1039 + VAT3( uPN, ii,jj,kk) * VAT3( oE, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk);
1047 VAT3( XuE, ii,jj,kk) =
1048 VAT3( uPS, ii,jj,kk) * VAT3( oE, i,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
1049 + VAT3( uPC, ii,jj,kk) * VAT3( oE, i,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
1050 + VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
1051 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPSW, iip1,jj,kkp1)
1053 - VAT3(uPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
1054 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSW, iip1,jj,kkp1)
1055 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPW, iip1,jj,kkp1)
1056 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPS, iip1,jj,kkp1))
1058 + VAT3( oPE, ii,jj,kk) * VAT3( uC, ip1,j,k) * VAT3( dPW, iip1,jj,kkp1)
1060 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
1061 + VAT3( oC, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
1062 - VAT3( uC, ip1,j,kp1) * VAT3( oPW, iip1,jj,kkp1)
1063 - VAT3( oN, ip1,j,kp1) * VAT3(dPNW, iip1,jj,kkp1)
1064 - VAT3( oE, ip1,j,kp1) * VAT3( dPC, iip1,jj,kkp1))
1066 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPNW, iip1,jj,kkp1)
1068 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oN, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
1069 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
1070 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNW, iip1,jj,kkp1)
1071 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPN, iip1,jj,kkp1));
1079 VAT3( XuW, ii,jj,kk) =
1080 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPSE, iim1,jj,kkp1)
1082 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oE, im2,jm1,kp1) * VAT3( dPS, iim1,jj,kkp1)
1083 + VAT3( oC, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
1084 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSE, iim1,jj,kkp1)
1085 - VAT3( oN, im1,jm1,kp1) * VAT3( dPE, iim1,jj,kkp1))
1087 + VAT3( oPW, ii,jj,kk) * VAT3( uC, im1,j,k) * VAT3( dPE, iim1,jj,kkp1)
1089 - VAT3( uPW, ii,jj,kk) * (- VAT3( oE, im2,j,kp1) * VAT3( dPC, iim1,jj,kkp1)
1090 - VAT3( oN, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
1091 + VAT3( oC, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
1092 - VAT3( uC, im1,j,kp1) * VAT3( oPE, iim1,jj,kkp1)
1093 - VAT3( oN, im1,j,kp1) * VAT3(dPNE, iim1,jj,kkp1))
1095 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPNE, iim1,jj,kkp1)
1097 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( dPN, iim1,jj,kkp1)
1098 - VAT3( oN, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
1099 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1)
1100 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNE, iim1,jj,kkp1))
1102 + VAT3( uPS, ii,jj,kk) * VAT3( oE, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
1103 + VAT3( uPC, ii,jj,kk) * VAT3( oE, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
1104 + VAT3( uPN, ii,jj,kk) * VAT3( oE, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1);
1112 VAT3( XuN, ii,jj,kk) =
1113 VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
1114 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPSW, ii,jjp1,kkp1)
1116 - VAT3(uPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
1117 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
1118 - VAT3( oN, im1,jp1,kp1) * VAT3( dPW, ii,jjp1,kkp1)
1119 - VAT3( oE, im1,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1))
1121 + VAT3( uPC, ii,jj,kk) * VAT3( oN, i,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
1122 + VAT3( oPN, ii,jj,kk) * VAT3( uC, i,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
1124 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
1125 + VAT3( oC, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
1126 - VAT3( uC, i,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
1127 - VAT3( oN, i,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
1128 - VAT3( oE, i,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1))
1130 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
1131 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPSE, ii,jjp1,kkp1)
1133 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
1134 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
1135 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSE, ii,jjp1,kkp1)
1136 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPE, ii,jjp1,kkp1));
1144 VAT3( XuS, ii,jj,kk) =
1145 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPNW, ii,jjm1,kkp1)
1147 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oN, im1,jm2,kp1) * VAT3( dPW, ii,jjm1,kkp1)
1148 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
1149 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
1150 - VAT3( oE, im1,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1))
1152 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
1153 + VAT3( oPS, ii,jj,kk) * VAT3( uC, i,jm1,k) * VAT3( dPN, ii,jjm1,kkp1)
1155 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
1156 - VAT3( oN, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
1157 + VAT3( oC, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
1158 - VAT3( uC, i,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1)
1159 - VAT3( oE, i,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1))
1161 + VAT3( uPC, ii,jj,kk) * VAT3( oN, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
1162 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPNE, ii,jjm1,kkp1)
1164 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
1165 - VAT3( oN, ip1,jm2,kp1) * VAT3( dPE, ii,jjm1,kkp1)
1166 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
1167 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNE, ii,jjm1,kkp1))
1169 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1);
1177 VAT3(XuNE, ii,jj,kk) =
1178 VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
1179 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
1180 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPSW, iip1,jjp1,kkp1)
1182 - VAT3(uPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
1183 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
1184 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPW, iip1,jjp1,kkp1)
1185 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPS, iip1,jjp1,kkp1));
1193 VAT3(XuNW, ii,jj,kk) =
1194 VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
1195 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPSE, iim1,jjp1,kkp1)
1197 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( dPS, iim1,jjp1,kkp1)
1198 + VAT3( oC, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
1199 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSE, iim1,jjp1,kkp1)
1200 - VAT3( oN, im1,jp1,kp1) * VAT3( dPE, iim1,jjp1,kkp1))
1202 + VAT3( uPN, ii,jj,kk) * VAT3( oE, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1);
1210 VAT3(XuSE, ii,jj,kk) =
1211 VAT3( uPS, ii,jj,kk) * VAT3( oE, i,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
1212 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPNW, iip1,jjm1,kkp1)
1214 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oN, ip1,jm2,kp1) * VAT3( dPW, iip1,jjm1,kkp1)
1215 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
1216 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNW, iip1,jjm1,kkp1)
1217 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPN, iip1,jjm1,kkp1))
1219 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1);
1227 VAT3(XuSW, ii,jj,kk) =
1228 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPNE, iim1,jjm1,kkp1)
1230 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oE, im2,jm1,kp1) * VAT3( dPN, iim1,jjm1,kkp1)
1231 - VAT3( oN, im1,jm2,kp1) * VAT3( dPE, iim1,jjm1,kkp1)
1232 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
1233 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
1235 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
1236 + VAT3( uPS, ii,jj,kk) * VAT3( oE, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1);
1248 int *
nx,
int *
ny,
int *
nz,
1249 double *oPC,
double *oPN,
double *oPS,
double *oPE,
double *oPW,
1250 double *oPNE,
double *oPNW,
double *oPSE,
double *oPSW,
1251 double *uPC,
double *uPN,
double *uPS,
double *uPE,
double *uPW,
1252 double *uPNE,
double *uPNW,
double *uPSE,
double *uPSW,
1253 double *dPC,
double *dPN,
double *dPS,
double *dPE,
double *dPW,
1254 double *dPNE,
double *dPNW,
double *dPSE,
double *dPSW,
1255 double *oC,
double *oE,
double *oN,
double *uC,
1256 double *oNE,
double *oNW,
double *uE,
double *uW,
double *uN,
1257 double *uS,
double *uNE,
double *uNW,
double *uSE,
double *uSW,
1258 double *XoC,
double *XoE,
double *XoN,
1260 double *XoNE,
double *XoNW,
1261 double *XuE,
double *XuW,
double *XuN,
double *XuS,
1262 double *XuNE,
double *XuNW,
double *XuSE,
double *XuSW) {
1266 int im1, ip1, im2, ip2;
1267 int jm1, jp1, jm2, jp2;
1268 int km1, kp1, km2, kp2;
1272 int nxm1, nym1, nzm1;
1274 double TMP1_XOC, TMP2_XOC, TMP3_XOC, TMP4_XOC;
1275 double TMP5_XOC, TMP6_XOC, TMP7_XOC, TMP8_XOC;
1276 double TMP9_XOC, TMP10_XOC, TMP11_XOC, TMP12_XOC;
1277 double TMP13_XOC, TMP14_XOC, TMP15_XOC, TMP16_XOC;
1278 double TMP17_XOC, TMP18_XOC, TMP19_XOC, TMP20_XOC;
1279 double TMP21_XOC, TMP22_XOC, TMP23_XOC, TMP24_XOC;
1280 double TMP25_XOC, TMP26_XOC, TMP27_XOC;
1282 double TMP1_XOE, TMP2_XOE, TMP3_XOE, TMP4_XOE;
1283 double TMP5_XOE, TMP6_XOE, TMP7_XOE, TMP8_XOE;
1284 double TMP9_XOE, TMP10_XOE, TMP11_XOE, TMP12_XOE;
1286 double TMP1_XON, TMP2_XON, TMP3_XON, TMP4_XON;
1287 double TMP5_XON, TMP6_XON, TMP7_XON, TMP8_XON;
1288 double TMP9_XON, TMP10_XON, TMP11_XON, TMP12_XON;
1290 double TMP1_XUC, TMP2_XUC, TMP3_XUC, TMP4_XUC;
1291 double TMP5_XUC, TMP6_XUC, TMP7_XUC, TMP8_XUC;
1292 double TMP9_XUC, TMP10_XUC, TMP11_XUC, TMP12_XUC;
1294 double TMP1_XONE, TMP2_XONE, TMP3_XONE, TMP4_XONE;
1295 double TMP5_XONE, TMP6_XONE;
1297 double TMP1_XONW, TMP2_XONW, TMP3_XONW, TMP4_XONW;
1298 double TMP5_XONW, TMP6_XONW;
1300 double TMP1_XUE, TMP2_XUE, TMP3_XUE, TMP4_XUE;
1301 double TMP5_XUE, TMP6_XUE;
1303 double TMP1_XUW, TMP2_XUW, TMP3_XUW, TMP4_XUW;
1304 double TMP5_XUW, TMP6_XUW;
1306 double TMP1_XUN, TMP2_XUN, TMP3_XUN, TMP4_XUN;
1307 double TMP5_XUN, TMP6_XUN;
1309 double TMP1_XUS, TMP2_XUS, TMP3_XUS, TMP4_XUS;
1310 double TMP5_XUS, TMP6_XUS;
1312 double TMP1_XUNE, TMP2_XUNE, TMP1_XUNW, TMP2_XUNW;
1313 double TMP1_XUSE, TMP2_XUSE, TMP1_XUSW, TMP2_XUSW;
1315 MAT3( oC, *nxf, *nyf, *nzf);
1316 MAT3( oE, *nxf, *nyf, *nzf);
1317 MAT3( oN, *nxf, *nyf, *nzf);
1318 MAT3( oNE, *nxf, *nyf, *nzf);
1319 MAT3( oNW, *nxf, *nyf, *nzf);
1321 MAT3( uC, *nxf, *nyf, *nzf);
1322 MAT3( uE, *nxf, *nyf, *nzf);
1323 MAT3( uW, *nxf, *nyf, *nzf);
1324 MAT3( uN, *nxf, *nyf, *nzf);
1325 MAT3( uS, *nxf, *nyf, *nzf);
1326 MAT3( uNE, *nxf, *nyf, *nzf);
1327 MAT3( uNW, *nxf, *nyf, *nzf);
1328 MAT3( uSE, *nxf, *nyf, *nzf);
1329 MAT3( uSW, *nxf, *nyf, *nzf);
1331 MAT3( XoC, *
nx, *
ny, *
nz);
1332 MAT3( XoE, *
nx, *
ny, *
nz);
1333 MAT3( XoN, *
nx, *
ny, *
nz);
1334 MAT3(XoNE, *
nx, *
ny, *
nz);
1335 MAT3(XoNW, *
nx, *
ny, *
nz);
1337 MAT3( XuC, *
nx, *
ny, *
nz);
1338 MAT3( XuE, *
nx, *
ny, *
nz);
1339 MAT3( XuW, *
nx, *
ny, *
nz);
1340 MAT3( XuN, *
nx, *
ny, *
nz);
1341 MAT3( XuS, *
nx, *
ny, *
nz);
1342 MAT3(XuNE, *
nx, *
ny, *
nz);
1343 MAT3(XuNW, *
nx, *
ny, *
nz);
1344 MAT3(XuSE, *
nx, *
ny, *
nz);
1345 MAT3(XuSW, *
nx, *
ny, *
nz);
1347 MAT3( oPC, *
nx, *
ny, *
nz);
1348 MAT3( oPN, *
nx, *
ny, *
nz);
1349 MAT3( oPS, *
nx, *
ny, *
nz);
1350 MAT3( oPE, *
nx, *
ny, *
nz);
1351 MAT3( oPW, *
nx, *
ny, *
nz);
1352 MAT3(oPNE, *
nx, *
ny, *
nz);
1353 MAT3(oPNW, *
nx, *
ny, *
nz);
1354 MAT3(oPSE, *
nx, *
ny, *
nz);
1355 MAT3(oPSW, *
nx, *
ny, *
nz);
1357 MAT3( uPC, *
nx, *
ny, *
nz);
1358 MAT3( uPN, *
nx, *
ny, *
nz);
1359 MAT3( uPS, *
nx, *
ny, *
nz);
1360 MAT3( uPE, *
nx, *
ny, *
nz);
1361 MAT3( uPW, *
nx, *
ny, *
nz);
1362 MAT3(uPNE, *
nx, *
ny, *
nz);
1363 MAT3(uPNW, *
nx, *
ny, *
nz);
1364 MAT3(uPSE, *
nx, *
ny, *
nz);
1365 MAT3(uPSW, *
nx, *
ny, *
nz);
1367 MAT3( dPC, *
nx, *
ny, *
nz);
1368 MAT3( dPN, *
nx, *
ny, *
nz);
1369 MAT3( dPS, *
nx, *
ny, *
nz);
1370 MAT3( dPE, *
nx, *
ny, *
nz);
1371 MAT3( dPW, *
nx, *
ny, *
nz);
1372 MAT3(dPNE, *
nx, *
ny, *
nz);
1373 MAT3(dPNW, *
nx, *
ny, *
nz);
1374 MAT3(dPSE, *
nx, *
ny, *
nz);
1375 MAT3(dPSW, *
nx, *
ny, *
nz);
1385 for(kk=2; kk<=*
nz-1; kk++) {
1388 for(jj=2; jj<=*
ny-1; jj++) {
1391 for(ii=2; ii<=*
nx-1; ii++) {
1422 VAT3( oPN, ii,jj,kk) * (- VAT3( uNE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1423 - VAT3( oNE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1424 - VAT3( uSW, i,jp1,k) * VAT3( uPW, ii,jj,kk)
1425 - VAT3( uE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1426 - VAT3( oE, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1427 - VAT3( uW, i,jp1,k) * VAT3(uPNW, ii,jj,kk)
1428 - VAT3( uN, i,j,km1) * VAT3( dPC, ii,jj,kk)
1429 - VAT3( oN, i,j,k) * VAT3( oPC, ii,jj,kk)
1430 - VAT3( uS, i,jp1,k) * VAT3( uPC, ii,jj,kk)
1431 - VAT3( uC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1432 + VAT3( oC, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1433 - VAT3( uC, i,jp1,k) * VAT3( uPN, ii,jj,kk)
1434 - VAT3( uNW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1435 - VAT3( oNW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1436 - VAT3( uSE, i,jp1,k) * VAT3( uPE, ii,jj,kk)
1437 - VAT3( uW, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1438 - VAT3( oE, i,jp1,k) * VAT3(oPNE, ii,jj,kk)
1439 - VAT3( uE, i,jp1,k) * VAT3(uPNE, ii,jj,kk));
1442 + VAT3( dPN, ii,jj,kk) * (- VAT3( oNE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1443 - VAT3( uSW, i,jp1,km1) * VAT3( oPW, ii,jj,kk)
1444 - VAT3( oE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1445 - VAT3( uW, i,jp1,km1) * VAT3(oPNW, ii,jj,kk)
1446 - VAT3( oN, i,j,km1) * VAT3( dPC, ii,jj,kk)
1447 - VAT3( uS, i,jp1,km1) * VAT3( oPC, ii,jj,kk)
1448 + VAT3( oC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1449 - VAT3( uC, i,jp1,km1) * VAT3( oPN, ii,jj,kk)
1450 - VAT3( oNW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1451 - VAT3( uSE, i,jp1,km1) * VAT3( oPE, ii,jj,kk)
1452 - VAT3( oE, i,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1453 - VAT3( uE, i,jp1,km1) * VAT3(oPNE, ii,jj,kk));
1456 + VAT3( dPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1457 - VAT3( uSW, i,j,km1) * VAT3(oPSW, ii,jj,kk)
1458 - VAT3( oE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1459 - VAT3( uW, i,j,km1) * VAT3( oPW, ii,jj,kk)
1460 - VAT3( oNW, i,j,km1) * VAT3(dPNW, ii,jj,kk)
1461 - VAT3( uNW, i,j,km1) * VAT3(oPNW, ii,jj,kk)
1462 - VAT3( oN, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1463 - VAT3( uS, i,j,km1) * VAT3( oPS, ii,jj,kk)
1464 + VAT3( oC, i,j,km1) * VAT3( dPC, ii,jj,kk)
1465 - VAT3( uC, i,j,km1) * VAT3( oPC, ii,jj,kk)
1466 - VAT3( oN, i,j,km1) * VAT3( dPN, ii,jj,kk)
1467 - VAT3( uN, i,j,km1) * VAT3( oPN, ii,jj,kk)
1468 - VAT3( oNW, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1469 - VAT3( uSE, i,j,km1) * VAT3(oPSE, ii,jj,kk)
1470 - VAT3( oE, i,j,km1) * VAT3( dPE, ii,jj,kk)
1471 - VAT3( uE, i,j,km1) * VAT3( oPE, ii,jj,kk)
1472 - VAT3( oNE, i,j,km1) * VAT3(dPNE, ii,jj,kk)
1473 - VAT3( uNE, i,j,km1) * VAT3(oPNE, ii,jj,kk));
1476 + VAT3( uPC, ii,jj,kk) * (- VAT3( uNE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1477 - VAT3( oNE, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1478 - VAT3( uE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1479 - VAT3( oE, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1480 - VAT3( uSE, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1481 - VAT3( oNW, i,j,kp1) * VAT3(uPNW, ii,jj,kk)
1482 - VAT3( uN, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1483 - VAT3( oN, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1484 - VAT3( uC, i,j,k) * VAT3( oPC, ii,jj,kk)
1485 + VAT3( oC, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1486 - VAT3( uS, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1487 - VAT3( oN, i,j,kp1) * VAT3( uPN, ii,jj,kk)
1488 - VAT3( uNW, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1489 - VAT3( oNW, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1490 - VAT3( uW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1491 - VAT3( oE, i,j,kp1) * VAT3( uPE, ii,jj,kk)
1492 - VAT3( uSW, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1493 - VAT3( oNE, i,j,kp1) * VAT3(uPNE, ii,jj,kk));
1496 + VAT3( oPC, ii,jj,kk) * (- VAT3( uW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1497 - VAT3( oE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1498 - VAT3( uSE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1499 - VAT3( uNE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1500 - VAT3( uN, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1501 - VAT3( oNE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1502 - VAT3( uE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1503 - VAT3( oNW, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1504 - VAT3( uC, i,j,km1) * VAT3( dPC, ii,jj,kk)
1505 - VAT3( uNW, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1506 - VAT3( uSW, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1507 - VAT3( uS, i,jp1,km1) * VAT3( dPN, ii,jj,kk) - VAT3( oN, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1508 - VAT3( uNE, i,j,k) * VAT3(uPNE, ii,jj,kk) - VAT3( oNE, i,j,k) * VAT3(oPNE, ii,jj,kk)
1509 - VAT3( uE, i,j,k) * VAT3( uPE, ii,jj,kk) - VAT3( uSE, i,j,k) * VAT3(uPSE, ii,jj,kk)
1510 - VAT3( oN, i,j,k) * VAT3( oPN, ii,jj,kk) - VAT3( oE, i,j,k) * VAT3( oPE, ii,jj,kk)
1511 - VAT3( uS, i,j,k) * VAT3( uPS, ii,jj,kk) + VAT3( oC, i,j,k) * VAT3( oPC, ii,jj,kk)
1512 - VAT3( uSW, i,j,k) * VAT3(uPSW, ii,jj,kk) - VAT3( uN, i,j,k) * VAT3( uPN, ii,jj,kk)
1513 - VAT3( uC, i,j,k) * VAT3( uPC, ii,jj,kk) - VAT3( uW, i,j,k) * VAT3( uPW, ii,jj,kk)
1514 - VAT3( oNW, i,j,k) * VAT3(oPNW, ii,jj,kk) - VAT3( uNW, i,j,k) * VAT3(uPNW, ii,jj,kk));
1517 + VAT3( uPS, ii,jj,kk) * (- VAT3( uE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1518 - VAT3( oE, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1519 - VAT3( uSE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1520 - VAT3( oNW, i,jm1,kp1) * VAT3( uPW, ii,jj,kk)
1521 - VAT3( uC, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1522 + VAT3( oC, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1523 - VAT3( uS, i,j,k) * VAT3( oPC, ii,jj,kk)
1524 - VAT3( oN, i,jm1,kp1) * VAT3( uPC, ii,jj,kk)
1525 - VAT3( uW, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1526 - VAT3( oE, i,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1527 - VAT3( uSW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1528 - VAT3( oNE, i,jm1,kp1) * VAT3( uPE, ii,jj,kk));
1531 + VAT3( oPS, ii,jj,kk) * (- VAT3( uE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1532 - VAT3( oE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1533 - VAT3( uW, i,jm1,k) * VAT3(uPSW, ii,jj,kk)
1534 - VAT3( uSE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1535 - VAT3( oNW, i,jm1,k) * VAT3( oPW, ii,jj,kk)
1536 - VAT3( uNW, i,jm1,k) * VAT3( uPW, ii,jj,kk)
1537 - VAT3( uC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1538 + VAT3( oC, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1539 - VAT3( uC, i,jm1,k) * VAT3( uPS, ii,jj,kk)
1540 - VAT3( uS, i,j,km1) * VAT3( dPC, ii,jj,kk)
1541 - VAT3( oN, i,jm1,k) * VAT3( oPC, ii,jj,kk)
1542 - VAT3( uN, i,jm1,k) * VAT3( uPC, ii,jj,kk)
1543 - VAT3( uW, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1544 - VAT3( oE, i,jm1,k) * VAT3(oPSE, ii,jj,kk)
1545 - VAT3( uE, i,jm1,k) * VAT3(uPSE, ii,jj,kk)
1546 - VAT3( uSW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1547 - VAT3( oNE, i,jm1,k) * VAT3( oPE, ii,jj,kk)
1548 - VAT3( uNE, i,jm1,k) * VAT3( uPE, ii,jj,kk));
1551 + VAT3( dPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1552 - VAT3( uW, i,jm1,km1) * VAT3(oPSW, ii,jj,kk)
1553 - VAT3( oNW, i,jm1,km1) * VAT3( dPW, ii,jj,kk)
1554 - VAT3( uNW, i,jm1,km1) * VAT3( oPW, ii,jj,kk)
1555 + VAT3( oC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1556 - VAT3( uC, i,jm1,km1) * VAT3( oPS, ii,jj,kk)
1557 - VAT3( oN, i,jm1,km1) * VAT3( dPC, ii,jj,kk)
1558 - VAT3( uN, i,jm1,km1) * VAT3( oPC, ii,jj,kk)
1559 - VAT3( oE, i,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1560 - VAT3( uE, i,jm1,km1) * VAT3(oPSE, ii,jj,kk)
1561 - VAT3( oNE, i,jm1,km1) * VAT3( dPE, ii,jj,kk)
1562 - VAT3( uNE, i,jm1,km1) * VAT3( oPE, ii,jj,kk));
1565 + VAT3(uPNW, ii,jj,kk) * (- VAT3( uN, im1,j,k) * VAT3( oPW, ii,jj,kk)
1566 - VAT3( oN, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1567 - VAT3( uC, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1568 + VAT3( oC, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
1569 - VAT3( uNW, i,j,k) * VAT3( oPC, ii,jj,kk)
1570 - VAT3( oNW, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1571 - VAT3( uW, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1572 - VAT3( oE, im1,jp1,kp1) * VAT3( uPN, ii,jj,kk));
1575 + VAT3(oPNW, ii,jj,kk) * (- VAT3( uN, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1576 - VAT3( oN, im1,j,k) * VAT3( oPW, ii,jj,kk)
1577 - VAT3( uS, im1,jp1,k) * VAT3( uPW, ii,jj,kk)
1578 - VAT3( uC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1579 + VAT3( oC, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1580 - VAT3( uC, im1,jp1,k) * VAT3(uPNW, ii,jj,kk)
1581 - VAT3( uNW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1582 - VAT3( oNW, i,j,k) * VAT3( oPC, ii,jj,kk)
1583 - VAT3( uSE, im1,jp1,k) * VAT3( uPC, ii,jj,kk)
1584 - VAT3( uW, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1585 - VAT3( oE, im1,jp1,k) * VAT3( oPN, ii,jj,kk)
1586 - VAT3( uE, im1,jp1,k) * VAT3( uPN, ii,jj,kk));
1589 + VAT3( uPW, ii,jj,kk) * (- VAT3( uN, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1590 - VAT3( oN, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1591 - VAT3( uC, im1,j,k) * VAT3( oPW, ii,jj,kk)
1592 + VAT3( oC, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1593 - VAT3( uS, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1594 - VAT3( oN, im1,j,kp1) * VAT3(uPNW, ii,jj,kk)
1595 - VAT3( uNW, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1596 - VAT3( oNW, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1597 - VAT3( uW, i,j,k) * VAT3( oPC, ii,jj,kk)
1598 - VAT3( oE, im1,j,kp1) * VAT3( uPC, ii,jj,kk)
1599 - VAT3( uSW, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1600 - VAT3( oNE, im1,j,kp1) * VAT3( uPN, ii,jj,kk));
1603 + VAT3(dPNW, ii,jj,kk) * (- VAT3( oN, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1604 - VAT3( uS, im1,jp1,km1) * VAT3( oPW, ii,jj,kk)
1605 + VAT3( oC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1606 - VAT3( uC, im1,jp1,km1) * VAT3(oPNW, ii,jj,kk)
1607 - VAT3( oNW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1608 - VAT3( uSE, im1,jp1,km1) * VAT3( oPC, ii,jj,kk)
1609 - VAT3( oE, im1,jp1,km1) * VAT3( dPN, ii,jj,kk)
1610 - VAT3( uE, im1,jp1,km1) * VAT3( oPN, ii,jj,kk));
1613 + VAT3( oPW, ii,jj,kk) * (- VAT3( uN, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1614 - VAT3( oN, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1615 - VAT3( uS, im1,j,k) * VAT3(uPSW, ii,jj,kk)
1616 - VAT3( uC, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1617 + VAT3( oC, im1,j,k) * VAT3( oPW, ii,jj,kk)
1618 - VAT3( uC, im1,j,k) * VAT3( uPW, ii,jj,kk)
1619 - VAT3( uS, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1620 - VAT3( oN, im1,j,k) * VAT3(oPNW, ii,jj,kk)
1621 - VAT3( uN, im1,j,k) * VAT3(uPNW, ii,jj,kk)
1622 - VAT3( uNW, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1623 - VAT3( oNW, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1624 - VAT3( uSE, im1,j,k) * VAT3( uPS, ii,jj,kk)
1625 - VAT3( uW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1626 - VAT3( oE, im1,j,k) * VAT3( oPC, ii,jj,kk)
1627 - VAT3( uE, im1,j,k) * VAT3( uPC, ii,jj,kk)
1628 - VAT3( uSW, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1629 - VAT3( oNE, im1,j,k) * VAT3( oPN, ii,jj,kk)
1630 - VAT3( uNE, im1,j,k) * VAT3( uPN, ii,jj,kk));
1633 + VAT3(uPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1634 + VAT3( oC, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1635 - VAT3( uS, im1,j,k) * VAT3( oPW, ii,jj,kk)
1636 - VAT3( oN, im1,jm1,kp1) * VAT3( uPW, ii,jj,kk)
1637 - VAT3( uW, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1638 - VAT3( oE, im1,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1639 - VAT3( uSW, i,j,k) * VAT3( oPC, ii,jj,kk)
1640 - VAT3( oNE, im1,jm1,kp1) * VAT3( uPC, ii,jj,kk));
1643 + VAT3(oPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1644 + VAT3( oC, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1645 - VAT3( uC, im1,jm1,k) * VAT3(uPSW, ii,jj,kk)
1646 - VAT3( uS, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1647 - VAT3( oN, im1,jm1,k) * VAT3( oPW, ii,jj,kk)
1648 - VAT3( uN, im1,jm1,k) * VAT3( uPW, ii,jj,kk)
1649 - VAT3( uW, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1650 - VAT3( oE, im1,jm1,k) * VAT3( oPS, ii,jj,kk)
1651 - VAT3( uE, im1,jm1,k) * VAT3( uPS, ii,jj,kk)
1652 - VAT3( uSW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1653 - VAT3( oNE, im1,jm1,k) * VAT3( oPC, ii,jj,kk)
1654 - VAT3( uNE, im1,jm1,k) * VAT3( uPC, ii,jj,kk));
1657 + VAT3( dPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1658 - VAT3( uS, im1,j,km1) * VAT3(oPSW, ii,jj,kk)
1659 + VAT3( oC, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1660 - VAT3( uC, im1,j,km1) * VAT3( oPW, ii,jj,kk)
1661 - VAT3( oN, im1,j,km1) * VAT3(dPNW, ii,jj,kk)
1662 - VAT3( uN, im1,j,km1) * VAT3(oPNW, ii,jj,kk)
1663 - VAT3( oNW, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1664 - VAT3( uSE, im1,j,km1) * VAT3( oPS, ii,jj,kk)
1665 - VAT3( oE, im1,j,km1) * VAT3( dPC, ii,jj,kk)
1666 - VAT3( uE, im1,j,km1) * VAT3( oPC, ii,jj,kk)
1667 - VAT3( oNE, im1,j,km1) * VAT3( dPN, ii,jj,kk)
1668 - VAT3( uNE, im1,j,km1) * VAT3( oPN, ii,jj,kk));
1671 + VAT3(uPNE, ii,jj,kk) * (- VAT3( uNE, i,j,k) * VAT3( oPC, ii,jj,kk)
1672 - VAT3( oNE, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1673 - VAT3( uE, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1674 - VAT3( oE, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
1675 - VAT3( uN, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1676 - VAT3( oN, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
1677 - VAT3( uC, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1678 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNE, ii,jj,kk));
1681 + VAT3( uPE, ii,jj,kk) * (- VAT3( uNE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1682 - VAT3( oNE, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1683 - VAT3( uE, i,j,k) * VAT3( oPC, ii,jj,kk)
1684 - VAT3( oE, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1685 - VAT3( uSE, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1686 - VAT3( oNW, ip1,j,kp1) * VAT3( uPN, ii,jj,kk)
1687 - VAT3( uN, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1688 - VAT3( oN, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1689 - VAT3( uC, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1690 + VAT3( oC, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
1691 - VAT3( uS, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1692 - VAT3( oN, ip1,j,kp1) * VAT3(uPNE, ii,jj,kk));
1695 + VAT3(dPNE, ii,jj,kk) * (- VAT3( oNE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1696 - VAT3( uSW, ip1,jp1,km1) * VAT3( oPC, ii,jj,kk)
1697 - VAT3( oE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1698 - VAT3( uW, ip1,jp1,km1) * VAT3( oPN, ii,jj,kk)
1699 - VAT3( oN, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1700 - VAT3( uS, ip1,jp1,km1) * VAT3( oPE, ii,jj,kk)
1701 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1702 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNE, ii,jj,kk));
1705 + VAT3(oPNE, ii,jj,kk) * (- VAT3( uNE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1706 - VAT3( oNE, i,j,k) * VAT3( oPC, ii,jj,kk)
1707 - VAT3( uSW, ip1,jp1,k) * VAT3( uPC, ii,jj,kk)
1708 - VAT3( uE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1709 - VAT3( oE, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1710 - VAT3( uW, ip1,jp1,k) * VAT3( uPN, ii,jj,kk)
1711 - VAT3( uN, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1712 - VAT3( oN, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1713 - VAT3( uS, ip1,jp1,k) * VAT3( uPE, ii,jj,kk)
1714 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1715 + VAT3( oC, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1716 - VAT3( uC, ip1,jp1,k) * VAT3(uPNE, ii,jj,kk));
1719 + VAT3(oPSE, ii,jj,kk) * (- VAT3( uE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1720 - VAT3( oE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1721 - VAT3( uW, ip1,jm1,k) * VAT3( uPS, ii,jj,kk)
1722 - VAT3( uSE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1723 - VAT3( oNW, ip1,jm1,k) * VAT3( oPC, ii,jj,kk)
1724 - VAT3( uNW, ip1,jm1,k) * VAT3( uPC, ii,jj,kk)
1725 - VAT3( uC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1726 + VAT3( oC, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1727 - VAT3( uC, ip1,jm1,k) * VAT3(uPSE, ii,jj,kk)
1728 - VAT3( uS, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1729 - VAT3( oN, ip1,jm1,k) * VAT3( oPE, ii,jj,kk)
1730 - VAT3( uN, ip1,jm1,k) * VAT3( uPE, ii,jj,kk));
1733 + VAT3(dPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1734 - VAT3( uW, ip1,jm1,km1) * VAT3( oPS, ii,jj,kk)
1735 - VAT3( oNW, ip1,jm1,km1) * VAT3( dPC, ii,jj,kk)
1736 - VAT3( uNW, ip1,jm1,km1) * VAT3( oPC, ii,jj,kk)
1737 + VAT3( oC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1738 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSE, ii,jj,kk)
1739 - VAT3( oN, ip1,jm1,km1) * VAT3( dPE, ii,jj,kk)
1740 - VAT3( uN, ip1,jm1,km1) * VAT3( oPE, ii,jj,kk));
1743 + VAT3(uPSE, ii,jj,kk) * (- VAT3( uE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1744 - VAT3( oE, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1745 - VAT3( uSE, i,j,k) * VAT3( oPC, ii,jj,kk)
1746 - VAT3( oNW, ip1,jm1,kp1) * VAT3( uPC, ii,jj,kk)
1747 - VAT3( uC, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1748 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1749 - VAT3( uS, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1750 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPE, ii,jj,kk));
1753 + VAT3( oPE, ii,jj,kk) * (- VAT3( uNE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1754 - VAT3( oNE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1755 - VAT3( uSW, ip1,j,k) * VAT3( uPS, ii,jj,kk)
1756 - VAT3( uE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1757 - VAT3( oE, i,j,k) * VAT3( oPC, ii,jj,kk)
1758 - VAT3( uW, ip1,j,k) * VAT3( uPC, ii,jj,kk)
1759 - VAT3( uSE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1760 - VAT3( oNW, ip1,j,k) * VAT3( oPN, ii,jj,kk)
1761 - VAT3( uNW, ip1,j,k) * VAT3( uPN, ii,jj,kk)
1762 - VAT3( uN, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1763 - VAT3( oN, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1764 - VAT3( uS, ip1,j,k) * VAT3(uPSE, ii,jj,kk)
1765 - VAT3( uC, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1766 + VAT3( oC, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1767 - VAT3( uC, ip1,j,k) * VAT3( uPE, ii,jj,kk)
1768 - VAT3( uS, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1769 - VAT3( oN, ip1,j,k) * VAT3(oPNE, ii,jj,kk)
1770 - VAT3( uN, ip1,j,k) * VAT3(uPNE, ii,jj,kk));
1773 + VAT3( dPE, ii,jj,kk) * (- VAT3( oNE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1774 - VAT3( uSW, ip1,j,km1) * VAT3( oPS, ii,jj,kk)
1775 - VAT3( oE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1776 - VAT3( uW, ip1,j,km1) * VAT3( oPC, ii,jj,kk)
1777 - VAT3( oNW, ip1,j,km1) * VAT3( dPN, ii,jj,kk)
1778 - VAT3( uNW, ip1,j,km1) * VAT3( oPN, ii,jj,kk)
1779 - VAT3( oN, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1780 - VAT3( uS, ip1,j,km1) * VAT3(oPSE, ii,jj,kk)
1781 + VAT3( oC, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1782 - VAT3( uC, ip1,j,km1) * VAT3( oPE, ii,jj,kk)
1783 - VAT3( oN, ip1,j,km1) * VAT3(dPNE, ii,jj,kk)
1784 - VAT3( uN, ip1,j,km1) * VAT3(oPNE, ii,jj,kk));
1787 + VAT3( uPN, ii,jj,kk) * (- VAT3( uNE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1788 - VAT3( oNE, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1789 - VAT3( uE, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1790 - VAT3( oE, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
1791 - VAT3( uN, i,j,k) * VAT3( oPC, ii,jj,kk)
1792 - VAT3( oN, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1793 - VAT3( uC, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1794 + VAT3( oC, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
1795 - VAT3( uNW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1796 - VAT3( oNW, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
1797 - VAT3( uW, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1798 - VAT3( oE, i,jp1,kp1) * VAT3(uPNE, ii,jj,kk));
1801 + VAT3(dPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1802 - VAT3( uC, im1,jm1,km1) * VAT3(oPSW, ii,jj,kk)
1803 - VAT3( oN, im1,jm1,km1) * VAT3( dPW, ii,jj,kk)
1804 - VAT3( uN, im1,jm1,km1) * VAT3( oPW, ii,jj,kk)
1805 - VAT3( oE, im1,jm1,km1) * VAT3( dPS, ii,jj,kk)
1806 - VAT3( uE, im1,jm1,km1) * VAT3( oPS, ii,jj,kk)
1807 - VAT3( oNE, im1,jm1,km1) * VAT3( dPC, ii,jj,kk)
1808 - VAT3( uNE, im1,jm1,km1) * VAT3( oPC, ii,jj,kk));
1810 VAT3( XoC, ii,jj,kk) = TMP1_XOC + TMP2_XOC + TMP3_XOC + TMP4_XOC
1811 + TMP5_XOC + TMP6_XOC + TMP7_XOC + TMP8_XOC + TMP9_XOC
1812 + TMP10_XOC + TMP11_XOC + TMP12_XOC + TMP13_XOC + TMP14_XOC
1813 + TMP15_XOC + TMP16_XOC + TMP17_XOC + TMP18_XOC + TMP19_XOC
1814 + TMP20_XOC + TMP21_XOC + TMP22_XOC + TMP23_XOC + TMP24_XOC
1815 + TMP25_XOC + TMP26_XOC + TMP27_XOC;
1825 - VAT3( dPS, ii,jj,kk) * (- VAT3( oE, i,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1826 - VAT3( uE, i,jm1,km1) * VAT3(oPSW, iip1,jj,kk)
1827 - VAT3( oNE, i,jm1,km1) * VAT3( dPW, iip1,jj,kk)
1828 - VAT3( uNE, i,jm1,km1) * VAT3( oPW, iip1,jj,kk))
1830 - VAT3( oPS, ii,jj,kk) * (- VAT3( uW, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1831 - VAT3( oE, i,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1832 - VAT3( uE, i,jm1,k) * VAT3(uPSW, iip1,jj,kk)
1833 - VAT3( uSW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1834 - VAT3( oNE, i,jm1,k) * VAT3( oPW, iip1,jj,kk)
1835 - VAT3( uNE, i,jm1,k) * VAT3( uPW, iip1,jj,kk))
1837 - VAT3( uPS, ii,jj,kk) * (- VAT3( uW, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1838 - VAT3( oE, i,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1839 - VAT3( uSW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1840 - VAT3( oNE, i,jm1,kp1) * VAT3( uPW, iip1,jj,kk));
1843 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNW, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1844 - VAT3( uSE, i,j,km1) * VAT3(oPSW, iip1,jj,kk)
1845 - VAT3( oE, i,j,km1) * VAT3( dPW, iip1,jj,kk)
1846 - VAT3( uE, i,j,km1) * VAT3( oPW, iip1,jj,kk)
1847 - VAT3( oNE, i,j,km1) * VAT3(dPNW, iip1,jj,kk)
1848 - VAT3( uNE, i,j,km1) * VAT3(oPNW, iip1,jj,kk))
1850 - VAT3( oPC, ii,jj,kk) * (- VAT3( uNW, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1851 - VAT3( oNW, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1852 - VAT3( uSE, i,j,k) * VAT3(uPSW, iip1,jj,kk)
1853 - VAT3( uW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1854 - VAT3( oE, i,j,k) * VAT3( oPW, iip1,jj,kk)
1855 - VAT3( uE, i,j,k) * VAT3( uPW, iip1,jj,kk)
1856 - VAT3( uSW, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1857 - VAT3( oNE, i,j,k) * VAT3(oPNW, iip1,jj,kk)
1858 - VAT3( uNE, i,j,k) * VAT3(uPNW, iip1,jj,kk));
1861 - VAT3( uPC, ii,jj,kk) * (- VAT3( uNW, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1862 - VAT3( oNW, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1863 - VAT3( uW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1864 - VAT3( oE, i,j,kp1) * VAT3( uPW, iip1,jj,kk)
1865 - VAT3( uSW, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1866 - VAT3( oNE, i,j,kp1) * VAT3(uPNW, iip1,jj,kk))
1868 - VAT3( dPN, ii,jj,kk) * (- VAT3( oNW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1869 - VAT3( uSE, i,jp1,km1) * VAT3( oPW, iip1,jj,kk)
1870 - VAT3( oE, i,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1871 - VAT3( uE, i,jp1,km1) * VAT3(oPNW, iip1,jj,kk))
1873 - VAT3( oPN, ii,jj,kk) * (- VAT3( uNW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1874 - VAT3( oNW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1875 - VAT3( uSE, i,jp1,k) * VAT3( uPW, iip1,jj,kk)
1876 - VAT3( uW, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1877 - VAT3( oE, i,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1878 - VAT3( uE, i,jp1,k) * VAT3(uPNW, iip1,jj,kk));
1881 - VAT3( uPN, ii,jj,kk) * (- VAT3( uNW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1882 - VAT3( oNW, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
1883 - VAT3( uW, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1884 - VAT3( oE, i,jp1,kp1) * VAT3(uPNW, iip1,jj,kk))
1886 - VAT3(dPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1887 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSW, iip1,jj,kk)
1888 - VAT3( oN, ip1,jm1,km1) * VAT3( dPW, iip1,jj,kk)
1889 - VAT3( uN, ip1,jm1,km1) * VAT3( oPW, iip1,jj,kk)
1890 - VAT3( oE, ip1,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1891 - VAT3( uE, ip1,jm1,km1) * VAT3( oPS, iip1,jj,kk)
1892 - VAT3( oNE, ip1,jm1,km1) * VAT3( dPC, iip1,jj,kk)
1893 - VAT3( uNE, ip1,jm1,km1) * VAT3( oPC, iip1,jj,kk));
1896 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1897 + VAT3( oC, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1898 - VAT3( uC, ip1,jm1,k) * VAT3(uPSW, iip1,jj,kk)
1899 - VAT3( uS, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1900 - VAT3( oN, ip1,jm1,k) * VAT3( oPW, iip1,jj,kk)
1901 - VAT3( uN, ip1,jm1,k) * VAT3( uPW, iip1,jj,kk)
1902 - VAT3( uW, ip2,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1903 - VAT3( oE, ip1,jm1,k) * VAT3( oPS, iip1,jj,kk)
1904 - VAT3( uE, ip1,jm1,k) * VAT3( uPS, iip1,jj,kk)
1905 - VAT3( uSW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1906 - VAT3( oNE, ip1,jm1,k) * VAT3( oPC, iip1,jj,kk)
1907 - VAT3( uNE, ip1,jm1,k) * VAT3( uPC, iip1,jj,kk));
1910 - VAT3(uPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1911 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1912 - VAT3( uS, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1913 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPW, iip1,jj,kk)
1914 - VAT3( uW, ip2,jm1,k) * VAT3( oPS, iip1,jj,kk)
1915 - VAT3( oE, ip1,jm1,kp1) * VAT3( uPS, iip1,jj,kk)
1916 - VAT3( uSW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
1917 - VAT3( oNE, ip1,jm1,kp1) * VAT3( uPC, iip1,jj,kk));
1920 - VAT3( dPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1921 - VAT3( uS, ip1,j,km1) * VAT3(oPSW, iip1,jj,kk)
1922 + VAT3( oC, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1923 - VAT3( uC, ip1,j,km1) * VAT3( oPW, iip1,jj,kk)
1924 - VAT3( oN, ip1,j,km1) * VAT3(dPNW, iip1,jj,kk)
1925 - VAT3( uN, ip1,j,km1) * VAT3(oPNW, iip1,jj,kk)
1926 - VAT3( oNW, ip2,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1927 - VAT3( uSE, ip1,j,km1) * VAT3( oPS, iip1,jj,kk)
1928 - VAT3( oE, ip1,j,km1) * VAT3( dPC, iip1,jj,kk)
1929 - VAT3( uE, ip1,j,km1) * VAT3( oPC, iip1,jj,kk)
1930 - VAT3( oNE, ip1,j,km1) * VAT3( dPN, iip1,jj,kk)
1931 - VAT3( uNE, ip1,j,km1) * VAT3( oPN, iip1,jj,kk));
1934 - VAT3( oPE, ii,jj,kk) * (- VAT3( uN, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1935 - VAT3( oN, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1936 - VAT3( uS, ip1,j,k) * VAT3(uPSW, iip1,jj,kk)
1937 - VAT3( uC, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1938 + VAT3( oC, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1939 - VAT3( uC, ip1,j,k) * VAT3( uPW, iip1,jj,kk)
1940 - VAT3( uS, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1941 - VAT3( oN, ip1,j,k) * VAT3(oPNW, iip1,jj,kk)
1942 - VAT3( uN, ip1,j,k) * VAT3(uPNW, iip1,jj,kk)
1943 - VAT3( uNW, ip2,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1944 - VAT3( oNW, ip2,jm1,k) * VAT3( oPS, iip1,jj,kk)
1945 - VAT3( uSE, ip1,j,k) * VAT3( uPS, iip1,jj,kk)
1946 - VAT3( uW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1947 - VAT3( oE, ip1,j,k) * VAT3( oPC, iip1,jj,kk)
1948 - VAT3( uE, ip1,j,k) * VAT3( uPC, iip1,jj,kk)
1949 - VAT3( uSW, ip2,jp1,km1) * VAT3( dPN, iip1,jj,kk)
1950 - VAT3( oNE, ip1,j,k) * VAT3( oPN, iip1,jj,kk)
1951 - VAT3( uNE, ip1,j,k) * VAT3( uPN, iip1,jj,kk));
1954 - VAT3( uPE, ii,jj,kk) * (- VAT3( uN, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1955 - VAT3( oN, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1956 - VAT3( uC, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1957 + VAT3( oC, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
1958 - VAT3( uS, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1959 - VAT3( oN, ip1,j,kp1) * VAT3(uPNW, iip1,jj,kk)
1960 - VAT3( uNW, ip2,jm1,k) * VAT3( oPS, iip1,jj,kk)
1961 - VAT3( oNW, ip2,jm1,kp1) * VAT3( uPS, iip1,jj,kk)
1962 - VAT3( uW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
1963 - VAT3( oE, ip1,j,kp1) * VAT3( uPC, iip1,jj,kk)
1964 - VAT3( uSW, ip2,jp1,k) * VAT3( oPN, iip1,jj,kk)
1965 - VAT3( oNE, ip1,j,kp1) * VAT3( uPN, iip1,jj,kk));
1968 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oN, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1969 - VAT3( uS, ip1,jp1,km1) * VAT3( oPW, iip1,jj,kk)
1970 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1971 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNW, iip1,jj,kk)
1972 - VAT3( oNW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1973 - VAT3( uSE, ip1,jp1,km1) * VAT3( oPC, iip1,jj,kk)
1974 - VAT3( oE, ip1,jp1,km1) * VAT3( dPN, iip1,jj,kk)
1975 - VAT3( uE, ip1,jp1,km1) * VAT3( oPN, iip1,jj,kk));
1978 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uN, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1979 - VAT3( oN, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1980 - VAT3( uS, ip1,jp1,k) * VAT3( uPW, iip1,jj,kk)
1981 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1982 + VAT3( oC, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1983 - VAT3( uC, ip1,jp1,k) * VAT3(uPNW, iip1,jj,kk)
1984 - VAT3( uNW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1985 - VAT3( oNW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
1986 - VAT3( uSE, ip1,jp1,k) * VAT3( uPC, iip1,jj,kk)
1987 - VAT3( uW, ip2,jp1,km1) * VAT3( dPN, iip1,jj,kk)
1988 - VAT3( oE, ip1,jp1,k) * VAT3( oPN, iip1,jj,kk)
1989 - VAT3( uE, ip1,jp1,k) * VAT3( uPN, iip1,jj,kk));
1992 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uN, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1993 - VAT3( oN, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
1994 - VAT3( uC, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1995 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNW, iip1,jj,kk)
1996 - VAT3( uNW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
1997 - VAT3( oNW, ip2,j,kp1) * VAT3( uPC, iip1,jj,kk)
1998 - VAT3( uW, ip2,jp1,k) * VAT3( oPN, iip1,jj,kk)
1999 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPN, iip1,jj,kk));
2001 VAT3( XoE, ii,jj,kk) = TMP1_XOE + TMP2_XOE + TMP3_XOE + TMP4_XOE
2002 + TMP5_XOE + TMP6_XOE + TMP7_XOE + TMP8_XOE + TMP9_XOE
2003 + TMP10_XOE + TMP11_XOE + TMP12_XOE;
2013 - VAT3( dPW, ii,jj,kk) * (- VAT3( oN, im1,j,km1) * VAT3(dPSW, ii,jjp1,kk)
2014 - VAT3( uN, im1,j,km1) * VAT3(oPSW, ii,jjp1,kk)
2015 - VAT3( oNE, im1,j,km1) * VAT3( dPS, ii,jjp1,kk)
2016 - VAT3( uNE, im1,j,km1) * VAT3( oPS, ii,jjp1,kk))
2018 - VAT3( oPW, ii,jj,kk) * (- VAT3( uS, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2019 - VAT3( oN, im1,j,k) * VAT3(oPSW, ii,jjp1,kk)
2020 - VAT3( uN, im1,j,k) * VAT3(uPSW, ii,jjp1,kk)
2021 - VAT3( uSW, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2022 - VAT3( oNE, im1,j,k) * VAT3( oPS, ii,jjp1,kk)
2023 - VAT3( uNE, im1,j,k) * VAT3( uPS, ii,jjp1,kk))
2025 - VAT3( uPW, ii,jj,kk) * (- VAT3( uS, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2026 - VAT3( oN, im1,j,kp1) * VAT3(uPSW, ii,jjp1,kk)
2027 - VAT3( uSW, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2028 - VAT3( oNE, im1,j,kp1) * VAT3( uPS, ii,jjp1,kk));
2031 - VAT3(dPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2032 - VAT3( uC, im1,jp1,km1) * VAT3(oPSW, ii,jjp1,kk)
2033 - VAT3( oN, im1,jp1,km1) * VAT3( dPW, ii,jjp1,kk)
2034 - VAT3( uN, im1,jp1,km1) * VAT3( oPW, ii,jjp1,kk)
2035 - VAT3( oE, im1,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2036 - VAT3( uE, im1,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
2037 - VAT3( oNE, im1,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
2038 - VAT3( uNE, im1,jp1,km1) * VAT3( oPC, ii,jjp1,kk));
2041 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2042 + VAT3( oC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2043 - VAT3( uC, im1,jp1,k) * VAT3(uPSW, ii,jjp1,kk)
2044 - VAT3( uS, im1,jp2,km1) * VAT3( dPW, ii,jjp1,kk)
2045 - VAT3( oN, im1,jp1,k) * VAT3( oPW, ii,jjp1,kk)
2046 - VAT3( uN, im1,jp1,k) * VAT3( uPW, ii,jjp1,kk)
2047 - VAT3( uW, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2048 - VAT3( oE, im1,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2049 - VAT3( uE, im1,jp1,k) * VAT3( uPS, ii,jjp1,kk)
2050 - VAT3( uSW, i,jp2,km1) * VAT3( dPC, ii,jjp1,kk)
2051 - VAT3( oNE, im1,jp1,k) * VAT3( oPC, ii,jjp1,kk)
2052 - VAT3( uNE, im1,jp1,k) * VAT3( uPC, ii,jjp1,kk));
2055 - VAT3(uPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2056 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
2057 - VAT3( uS, im1,jp2,k) * VAT3( oPW, ii,jjp1,kk)
2058 - VAT3( oN, im1,jp1,kp1) * VAT3( uPW, ii,jjp1,kk)
2059 - VAT3( uW, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2060 - VAT3( oE, im1,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
2061 - VAT3( uSW, i,jp2,k) * VAT3( oPC, ii,jjp1,kk)
2062 - VAT3( oNE, im1,jp1,kp1) * VAT3( uPC, ii,jjp1,kk))
2064 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNW, i,j,km1) * VAT3(dPSW, ii,jjp1,kk)
2065 - VAT3( uNW, i,j,km1) * VAT3(oPSW, ii,jjp1,kk)
2066 - VAT3( oN, i,j,km1) * VAT3( dPS, ii,jjp1,kk)
2067 - VAT3( uN, i,j,km1) * VAT3( oPS, ii,jjp1,kk)
2068 - VAT3( oNE, i,j,km1) * VAT3(dPSE, ii,jjp1,kk)
2069 - VAT3( uNE, i,j,km1) * VAT3(oPSE, ii,jjp1,kk));
2072 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2073 - VAT3( oNW, i,j,k) * VAT3(oPSW, ii,jjp1,kk)
2074 - VAT3( uNW, i,j,k) * VAT3(uPSW, ii,jjp1,kk)
2075 - VAT3( uS, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2076 - VAT3( oN, i,j,k) * VAT3( oPS, ii,jjp1,kk)
2077 - VAT3( uN, i,j,k) * VAT3( uPS, ii,jjp1,kk)
2078 - VAT3( uSW, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2079 - VAT3( oNE, i,j,k) * VAT3(oPSE, ii,jjp1,kk)
2080 - VAT3( uNE, i,j,k) * VAT3(uPSE, ii,jjp1,kk))
2082 - VAT3( uPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2083 - VAT3( oNW, i,j,kp1) * VAT3(uPSW, ii,jjp1,kk)
2084 - VAT3( uS, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2085 - VAT3( oN, i,j,kp1) * VAT3( uPS, ii,jjp1,kk)
2086 - VAT3( uSW, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2087 - VAT3( oNE, i,j,kp1) * VAT3(uPSE, ii,jjp1,kk));
2090 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2091 - VAT3( uW, i,jp1,km1) * VAT3(oPSW, ii,jjp1,kk)
2092 - VAT3( oNW, i,jp1,km1) * VAT3( dPW, ii,jjp1,kk)
2093 - VAT3( uNW, i,jp1,km1) * VAT3( oPW, ii,jjp1,kk)
2094 + VAT3( oC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2095 - VAT3( uC, i,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
2096 - VAT3( oN, i,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
2097 - VAT3( uN, i,jp1,km1) * VAT3( oPC, ii,jjp1,kk)
2098 - VAT3( oE, i,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2099 - VAT3( uE, i,jp1,km1) * VAT3(oPSE, ii,jjp1,kk)
2100 - VAT3( oNE, i,jp1,km1) * VAT3( dPE, ii,jjp1,kk)
2101 - VAT3( uNE, i,jp1,km1) * VAT3( oPE, ii,jjp1,kk));
2104 - VAT3( oPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2105 - VAT3( oE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2106 - VAT3( uW, i,jp1,k) * VAT3(uPSW, ii,jjp1,kk)
2107 - VAT3( uSE, im1,jp2,km1) * VAT3( dPW, ii,jjp1,kk)
2108 - VAT3( oNW, i,jp1,k) * VAT3( oPW, ii,jjp1,kk)
2109 - VAT3( uNW, i,jp1,k) * VAT3( uPW, ii,jjp1,kk)
2110 - VAT3( uC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2111 + VAT3( oC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2112 - VAT3( uC, i,jp1,k) * VAT3( uPS, ii,jjp1,kk)
2113 - VAT3( uS, i,jp2,km1) * VAT3( dPC, ii,jjp1,kk)
2114 - VAT3( oN, i,jp1,k) * VAT3( oPC, ii,jjp1,kk)
2115 - VAT3( uN, i,jp1,k) * VAT3( uPC, ii,jjp1,kk)
2116 - VAT3( uW, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2117 - VAT3( oE, i,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2118 - VAT3( uE, i,jp1,k) * VAT3(uPSE, ii,jjp1,kk)
2119 - VAT3( uSW, ip1,jp2,km1) * VAT3( dPE, ii,jjp1,kk)
2120 - VAT3( oNE, i,jp1,k) * VAT3( oPE, ii,jjp1,kk)
2121 - VAT3( uNE, i,jp1,k) * VAT3( uPE, ii,jjp1,kk));
2124 - VAT3( uPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2125 - VAT3( oE, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
2126 - VAT3( uSE, im1,jp2,k) * VAT3( oPW, ii,jjp1,kk)
2127 - VAT3( oNW, i,jp1,kp1) * VAT3( uPW, ii,jjp1,kk)
2128 - VAT3( uC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2129 + VAT3( oC, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
2130 - VAT3( uS, i,jp2,k) * VAT3( oPC, ii,jjp1,kk)
2131 - VAT3( oN, i,jp1,kp1) * VAT3( uPC, ii,jjp1,kk)
2132 - VAT3( uW, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2133 - VAT3( oE, i,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk)
2134 - VAT3( uSW, ip1,jp2,k) * VAT3( oPE, ii,jjp1,kk)
2135 - VAT3( oNE, i,jp1,kp1) * VAT3( uPE, ii,jjp1,kk))
2137 - VAT3( dPE, ii,jj,kk) * (- VAT3( oNW, ip1,j,km1) * VAT3( dPS, ii,jjp1,kk)
2138 - VAT3( uNW, ip1,j,km1) * VAT3( oPS, ii,jjp1,kk)
2139 - VAT3( oN, ip1,j,km1) * VAT3(dPSE, ii,jjp1,kk)
2140 - VAT3( uN, ip1,j,km1) * VAT3(oPSE, ii,jjp1,kk));
2143 - VAT3( oPE, ii,jj,kk) * (- VAT3( uSE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2144 - VAT3( oNW, ip1,j,k) * VAT3( oPS, ii,jjp1,kk)
2145 - VAT3( uNW, ip1,j,k) * VAT3( uPS, ii,jjp1,kk)
2146 - VAT3( uS, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2147 - VAT3( oN, ip1,j,k) * VAT3(oPSE, ii,jjp1,kk)
2148 - VAT3( uN, ip1,j,k) * VAT3(uPSE, ii,jjp1,kk))
2150 - VAT3( uPE, ii,jj,kk) * (- VAT3( uSE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2151 - VAT3( oNW, ip1,j,kp1) * VAT3( uPS, ii,jjp1,kk)
2152 - VAT3( uS, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2153 - VAT3( oN, ip1,j,kp1) * VAT3(uPSE, ii,jjp1,kk));
2156 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2157 - VAT3( uW, ip1,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
2158 - VAT3( oNW, ip1,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
2159 - VAT3( uNW, ip1,jp1,km1) * VAT3( oPC, ii,jjp1,kk)
2160 + VAT3( oC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2161 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSE, ii,jjp1,kk)
2162 - VAT3( oN, ip1,jp1,km1) * VAT3( dPE, ii,jjp1,kk)
2163 - VAT3( uN, ip1,jp1,km1) * VAT3( oPE, ii,jjp1,kk));
2166 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2167 - VAT3( oE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2168 - VAT3( uW, ip1,jp1,k) * VAT3( uPS, ii,jjp1,kk)
2169 - VAT3( uSE, i,jp2,km1) * VAT3( dPC, ii,jjp1,kk)
2170 - VAT3( oNW, ip1,jp1,k) * VAT3( oPC, ii,jjp1,kk)
2171 - VAT3( uNW, ip1,jp1,k) * VAT3( uPC, ii,jjp1,kk)
2172 - VAT3( uC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2173 + VAT3( oC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2174 - VAT3( uC, ip1,jp1,k) * VAT3(uPSE, ii,jjp1,kk)
2175 - VAT3( uS, ip1,jp2,km1) * VAT3( dPE, ii,jjp1,kk)
2176 - VAT3( oN, ip1,jp1,k) * VAT3( oPE, ii,jjp1,kk)
2177 - VAT3( uN, ip1,jp1,k) * VAT3( uPE, ii,jjp1,kk));
2180 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2181 - VAT3( oE, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
2182 - VAT3( uSE, i,jp2,k) * VAT3( oPC, ii,jjp1,kk)
2183 - VAT3( oNW, ip1,jp1,kp1) * VAT3( uPC, ii,jjp1,kk)
2184 - VAT3( uC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2185 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk)
2186 - VAT3( uS, ip1,jp2,k) * VAT3( oPE, ii,jjp1,kk)
2187 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPE, ii,jjp1,kk));
2189 VAT3( XoN, ii,jj,kk) = TMP1_XON + TMP2_XON + TMP3_XON + TMP4_XON
2190 + TMP5_XON + TMP6_XON + TMP7_XON + TMP8_XON + TMP9_XON
2191 + TMP10_XON + TMP11_XON + TMP12_XON;
2202 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,k) * VAT3(dPSW, ii,jj,kkp1)
2203 - VAT3( uN, im1,jm1,k) * VAT3( dPW, ii,jj,kkp1)
2204 - VAT3( uE, im1,jm1,k) * VAT3( dPS, ii,jj,kkp1)
2205 - VAT3( uNE, im1,jm1,k) * VAT3( dPC, ii,jj,kkp1))
2207 - VAT3(uPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2208 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSW, ii,jj,kkp1)
2209 - VAT3( oN, im1,jm1,kp1) * VAT3( dPW, ii,jj,kkp1)
2210 - VAT3( uN, im1,jm1,kp1) * VAT3( oPW, ii,jj,kkp1)
2211 - VAT3( oE, im1,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2212 - VAT3( uE, im1,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
2213 - VAT3( oNE, im1,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
2214 - VAT3( uNE, im1,jm1,kp1) * VAT3( oPC, ii,jj,kkp1));
2217 - VAT3( oPW, ii,jj,kk) * (- VAT3( uS, im1,j,k) * VAT3(dPSW, ii,jj,kkp1)
2218 - VAT3( uC, im1,j,k) * VAT3( dPW, ii,jj,kkp1)
2219 - VAT3( uN, im1,j,k) * VAT3(dPNW, ii,jj,kkp1)
2220 - VAT3( uSE, im1,j,k) * VAT3( dPS, ii,jj,kkp1)
2221 - VAT3( uE, im1,j,k) * VAT3( dPC, ii,jj,kkp1)
2222 - VAT3( uNE, im1,j,k) * VAT3( dPN, ii,jj,kkp1));
2225 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2226 - VAT3( uS, im1,j,kp1) * VAT3(oPSW, ii,jj,kkp1)
2227 + VAT3( oC, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2228 - VAT3( uC, im1,j,kp1) * VAT3( oPW, ii,jj,kkp1)
2229 - VAT3( oN, im1,j,kp1) * VAT3(dPNW, ii,jj,kkp1)
2230 - VAT3( uN, im1,j,kp1) * VAT3(oPNW, ii,jj,kkp1)
2231 - VAT3( oNW, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2232 - VAT3( uSE, im1,j,kp1) * VAT3( oPS, ii,jj,kkp1)
2233 - VAT3( oE, im1,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2234 - VAT3( uE, im1,j,kp1) * VAT3( oPC, ii,jj,kkp1)
2235 - VAT3( oNE, im1,j,kp1) * VAT3( dPN, ii,jj,kkp1)
2236 - VAT3( uNE, im1,j,kp1) * VAT3( oPN, ii,jj,kkp1))
2238 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uS, im1,jp1,k) * VAT3( dPW, ii,jj,kkp1)
2239 - VAT3( uC, im1,jp1,k) * VAT3(dPNW, ii,jj,kkp1)
2240 - VAT3( uSE, im1,jp1,k) * VAT3( dPC, ii,jj,kkp1)
2241 - VAT3( uE, im1,jp1,k) * VAT3( dPN, ii,jj,kkp1));
2244 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oN, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2245 - VAT3( uS, im1,jp1,kp1) * VAT3( oPW, ii,jj,kkp1)
2246 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
2247 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNW, ii,jj,kkp1)
2248 - VAT3( oNW, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2249 - VAT3( uSE, im1,jp1,kp1) * VAT3( oPC, ii,jj,kkp1)
2250 - VAT3( oE, im1,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
2251 - VAT3( uE, im1,jp1,kp1) * VAT3( oPN, ii,jj,kkp1))
2253 - VAT3( oPS, ii,jj,kk) * (- VAT3( uW, i,jm1,k) * VAT3(dPSW, ii,jj,kkp1)
2254 - VAT3( uNW, i,jm1,k) * VAT3( dPW, ii,jj,kkp1)
2255 - VAT3( uC, i,jm1,k) * VAT3( dPS, ii,jj,kkp1)
2256 - VAT3( uN, i,jm1,k) * VAT3( dPC, ii,jj,kkp1)
2257 - VAT3( uE, i,jm1,k) * VAT3(dPSE, ii,jj,kkp1)
2258 - VAT3( uNE, i,jm1,k) * VAT3( dPE, ii,jj,kkp1));
2261 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2262 - VAT3( uW, i,jm1,kp1) * VAT3(oPSW, ii,jj,kkp1)
2263 - VAT3( oNW, i,jm1,kp1) * VAT3( dPW, ii,jj,kkp1)
2264 - VAT3( uNW, i,jm1,kp1) * VAT3( oPW, ii,jj,kkp1)
2265 + VAT3( oC, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2266 - VAT3( uC, i,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
2267 - VAT3( oN, i,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
2268 - VAT3( uN, i,jm1,kp1) * VAT3( oPC, ii,jj,kkp1)
2269 - VAT3( oE, i,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2270 - VAT3( uE, i,jm1,kp1) * VAT3(oPSE, ii,jj,kkp1)
2271 - VAT3( oNE, i,jm1,kp1) * VAT3( dPE, ii,jj,kkp1)
2272 - VAT3( uNE, i,jm1,kp1) * VAT3( oPE, ii,jj,kkp1));
2275 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, i,j,k) * VAT3(dPSW, ii,jj,kkp1)
2276 - VAT3( uW, i,j,k) * VAT3( dPW, ii,jj,kkp1)
2277 - VAT3( uNW, i,j,k) * VAT3(dPNW, ii,jj,kkp1)
2278 - VAT3( uS, i,j,k) * VAT3( dPS, ii,jj,kkp1)
2279 - VAT3( uC, i,j,k) * VAT3( dPC, ii,jj,kkp1)
2280 - VAT3( uN, i,j,k) * VAT3( dPN, ii,jj,kkp1)
2281 - VAT3( uSE, i,j,k) * VAT3(dPSE, ii,jj,kkp1)
2282 - VAT3( uE, i,j,k) * VAT3( dPE, ii,jj,kkp1)
2283 - VAT3( uNE, i,j,k) * VAT3(dPNE, ii,jj,kkp1));
2286 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2287 - VAT3( uSW, i,j,kp1) * VAT3(oPSW, ii,jj,kkp1)
2288 - VAT3( oE, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2289 - VAT3( uW, i,j,kp1) * VAT3( oPW, ii,jj,kkp1)
2290 - VAT3( oNW, i,j,kp1) * VAT3(dPNW, ii,jj,kkp1)
2291 - VAT3( uNW, i,j,kp1) * VAT3(oPNW, ii,jj,kkp1)
2292 - VAT3( oN, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2293 - VAT3( uS, i,j,kp1) * VAT3( oPS, ii,jj,kkp1)
2294 + VAT3( oC, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2295 - VAT3( uC, i,j,kp1) * VAT3( oPC, ii,jj,kkp1)
2296 - VAT3( oN, i,j,kp1) * VAT3( dPN, ii,jj,kkp1)
2297 - VAT3( uN, i,j,kp1) * VAT3( oPN, ii,jj,kkp1)
2298 - VAT3( oNW, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2299 - VAT3( uSE, i,j,kp1) * VAT3(oPSE, ii,jj,kkp1)
2300 - VAT3( oE, i,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2301 - VAT3( uE, i,j,kp1) * VAT3( oPE, ii,jj,kkp1)
2302 - VAT3( oNE, i,j,kp1) * VAT3(dPNE, ii,jj,kkp1)
2303 - VAT3( uNE, i,j,kp1) * VAT3(oPNE, ii,jj,kkp1));
2306 - VAT3( oPN, ii,jj,kk) * (- VAT3( uSW, i,jp1,k) * VAT3( dPW, ii,jj,kkp1)
2307 - VAT3( uW, i,jp1,k) * VAT3(dPNW, ii,jj,kkp1)
2308 - VAT3( uS, i,jp1,k) * VAT3( dPC, ii,jj,kkp1)
2309 - VAT3( uC, i,jp1,k) * VAT3( dPN, ii,jj,kkp1)
2310 - VAT3( uSE, i,jp1,k) * VAT3( dPE, ii,jj,kkp1)
2311 - VAT3( uE, i,jp1,k) * VAT3(dPNE, ii,jj,kkp1));
2314 - VAT3( uPN, ii,jj,kk) * (- VAT3( oNE, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2315 - VAT3( uSW, i,jp1,kp1) * VAT3( oPW, ii,jj,kkp1)
2316 - VAT3( oE, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
2317 - VAT3( uW, i,jp1,kp1) * VAT3(oPNW, ii,jj,kkp1)
2318 - VAT3( oN, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2319 - VAT3( uS, i,jp1,kp1) * VAT3( oPC, ii,jj,kkp1)
2320 + VAT3( oC, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
2321 - VAT3( uC, i,jp1,kp1) * VAT3( oPN, ii,jj,kkp1)
2322 - VAT3( oNW, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2323 - VAT3( uSE, i,jp1,kp1) * VAT3( oPE, ii,jj,kkp1)
2324 - VAT3( oE, i,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1)
2325 - VAT3( uE, i,jp1,kp1) * VAT3(oPNE, ii,jj,kkp1))
2327 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uW, ip1,jm1,k) * VAT3( dPS, ii,jj,kkp1)
2328 - VAT3( uNW, ip1,jm1,k) * VAT3( dPC, ii,jj,kkp1)
2329 - VAT3( uC, ip1,jm1,k) * VAT3(dPSE, ii,jj,kkp1)
2330 - VAT3( uN, ip1,jm1,k) * VAT3( dPE, ii,jj,kkp1));
2333 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2334 - VAT3( uW, ip1,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
2335 - VAT3( oNW, ip1,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
2336 - VAT3( uNW, ip1,jm1,kp1) * VAT3( oPC, ii,jj,kkp1)
2337 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2338 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSE, ii,jj,kkp1)
2339 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPE, ii,jj,kkp1)
2340 - VAT3( uN, ip1,jm1,kp1) * VAT3( oPE, ii,jj,kkp1))
2342 - VAT3( oPE, ii,jj,kk) * (- VAT3( uSW, ip1,j,k) * VAT3( dPS, ii,jj,kkp1)
2343 - VAT3( uW, ip1,j,k) * VAT3( dPC, ii,jj,kkp1)
2344 - VAT3( uNW, ip1,j,k) * VAT3( dPN, ii,jj,kkp1)
2345 - VAT3( uS, ip1,j,k) * VAT3(dPSE, ii,jj,kkp1)
2346 - VAT3( uC, ip1,j,k) * VAT3( dPE, ii,jj,kkp1)
2347 - VAT3( uN, ip1,j,k) * VAT3(dPNE, ii,jj,kkp1));
2350 - VAT3( uPE, ii,jj,kk) * (- VAT3( oNE, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2351 - VAT3( uSW, ip1,j,kp1) * VAT3( oPS, ii,jj,kkp1)
2352 - VAT3( oE, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2353 - VAT3( uW, ip1,j,kp1) * VAT3( oPC, ii,jj,kkp1)
2354 - VAT3( oNW, ip1,j,kp1) * VAT3( dPN, ii,jj,kkp1)
2355 - VAT3( uNW, ip1,j,kp1) * VAT3( oPN, ii,jj,kkp1)
2356 - VAT3( oN, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2357 - VAT3( uS, ip1,j,kp1) * VAT3(oPSE, ii,jj,kkp1)
2358 + VAT3( oC, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2359 - VAT3( uC, ip1,j,kp1) * VAT3( oPE, ii,jj,kkp1)
2360 - VAT3( oN, ip1,j,kp1) * VAT3(dPNE, ii,jj,kkp1)
2361 - VAT3( uN, ip1,j,kp1) * VAT3(oPNE, ii,jj,kkp1));
2364 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uSW, ip1,jp1,k) * VAT3( dPC, ii,jj,kkp1)
2365 - VAT3( uW, ip1,jp1,k) * VAT3( dPN, ii,jj,kkp1)
2366 - VAT3( uS, ip1,jp1,k) * VAT3( dPE, ii,jj,kkp1)
2367 - VAT3( uC, ip1,jp1,k) * VAT3(dPNE, ii,jj,kkp1))
2369 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oNE, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2370 - VAT3( uSW, ip1,jp1,kp1) * VAT3( oPC, ii,jj,kkp1)
2371 - VAT3( oE, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
2372 - VAT3( uW, ip1,jp1,kp1) * VAT3( oPN, ii,jj,kkp1)
2373 - VAT3( oN, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2374 - VAT3( uS, ip1,jp1,kp1) * VAT3( oPE, ii,jj,kkp1)
2375 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1)
2376 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNE, ii,jj,kkp1));
2378 VAT3( XuC, ii,jj,kk) = TMP1_XUC + TMP2_XUC + TMP3_XUC + TMP4_XUC
2379 + TMP5_XUC + TMP6_XUC + TMP7_XUC + TMP8_XUC + TMP9_XUC
2380 + TMP10_XUC + TMP11_XUC + TMP12_XUC;
2390 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNE, i,j,km1) * VAT3(dPSW, iip1,jjp1,kk)
2391 - VAT3( uNE, i,j,km1) * VAT3(oPSW, iip1,jjp1,kk))
2393 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2394 - VAT3( oNE, i,j,k) * VAT3(oPSW, iip1,jjp1,kk)
2395 - VAT3( uNE, i,j,k) * VAT3(uPSW, iip1,jjp1,kk))
2397 - VAT3( uPC, ii,jj,kk) * (- VAT3( uSW, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2398 - VAT3( oNE, i,j,kp1) * VAT3(uPSW, iip1,jjp1,kk))
2400 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2401 - VAT3( uE, i,jp1,km1) * VAT3(oPSW, iip1,jjp1,kk)
2402 - VAT3( oNE, i,jp1,km1) * VAT3( dPW, iip1,jjp1,kk)
2403 - VAT3( uNE, i,jp1,km1) * VAT3( oPW, iip1,jjp1,kk));
2406 - VAT3( oPN, ii,jj,kk) * (- VAT3( uW, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2407 - VAT3( oE, i,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2408 - VAT3( uE, i,jp1,k) * VAT3(uPSW, iip1,jjp1,kk)
2409 - VAT3( uSW, ip1,jp2,km1) * VAT3( dPW, iip1,jjp1,kk)
2410 - VAT3( oNE, i,jp1,k) * VAT3( oPW, iip1,jjp1,kk)
2411 - VAT3( uNE, i,jp1,k) * VAT3( uPW, iip1,jjp1,kk))
2413 - VAT3( uPN, ii,jj,kk) * (- VAT3( uW, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2414 - VAT3( oE, i,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
2415 - VAT3( uSW, ip1,jp2,k) * VAT3( oPW, iip1,jjp1,kk)
2416 - VAT3( oNE, i,jp1,kp1) * VAT3( uPW, iip1,jjp1,kk))
2418 - VAT3( dPE, ii,jj,kk) * (- VAT3( oN, ip1,j,km1) * VAT3(dPSW, iip1,jjp1,kk)
2419 - VAT3( uN, ip1,j,km1) * VAT3(oPSW, iip1,jjp1,kk)
2420 - VAT3( oNE, ip1,j,km1) * VAT3( dPS, iip1,jjp1,kk)
2421 - VAT3( uNE, ip1,j,km1) * VAT3( oPS, iip1,jjp1,kk));
2424 - VAT3( oPE, ii,jj,kk) * (- VAT3( uS, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2425 - VAT3( oN, ip1,j,k) * VAT3(oPSW, iip1,jjp1,kk)
2426 - VAT3( uN, ip1,j,k) * VAT3(uPSW, iip1,jjp1,kk)
2427 - VAT3( uSW, ip2,jp1,km1) * VAT3( dPS, iip1,jjp1,kk)
2428 - VAT3( oNE, ip1,j,k) * VAT3( oPS, iip1,jjp1,kk)
2429 - VAT3( uNE, ip1,j,k) * VAT3( uPS, iip1,jjp1,kk))
2431 - VAT3( uPE, ii,jj,kk) * (- VAT3( uS, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2432 - VAT3( oN, ip1,j,kp1) * VAT3(uPSW, iip1,jjp1,kk)
2433 - VAT3( uSW, ip2,jp1,k) * VAT3( oPS, iip1,jjp1,kk)
2434 - VAT3( oNE, ip1,j,kp1) * VAT3( uPS, iip1,jjp1,kk));
2437 - VAT3(dPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2438 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSW, iip1,jjp1,kk)
2439 - VAT3( oN, ip1,jp1,km1) * VAT3( dPW, iip1,jjp1,kk)
2440 - VAT3( uN, ip1,jp1,km1) * VAT3( oPW, iip1,jjp1,kk)
2441 - VAT3( oE, ip1,jp1,km1) * VAT3( dPS, iip1,jjp1,kk)
2442 - VAT3( uE, ip1,jp1,km1) * VAT3( oPS, iip1,jjp1,kk)
2443 - VAT3( oNE, ip1,jp1,km1) * VAT3( dPC, iip1,jjp1,kk)
2444 - VAT3( uNE, ip1,jp1,km1) * VAT3( oPC, iip1,jjp1,kk));
2447 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2448 + VAT3( oC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2449 - VAT3( uC, ip1,jp1,k) * VAT3(uPSW, iip1,jjp1,kk)
2450 - VAT3( uS, ip1,jp2,km1) * VAT3( dPW, iip1,jjp1,kk)
2451 - VAT3( oN, ip1,jp1,k) * VAT3( oPW, iip1,jjp1,kk)
2452 - VAT3( uN, ip1,jp1,k) * VAT3( uPW, iip1,jjp1,kk)
2453 - VAT3( uW, ip2,jp1,km1) * VAT3( dPS, iip1,jjp1,kk)
2454 - VAT3( oE, ip1,jp1,k) * VAT3( oPS, iip1,jjp1,kk)
2455 - VAT3( uE, ip1,jp1,k) * VAT3( uPS, iip1,jjp1,kk)
2456 - VAT3( uSW, ip2,jp2,km1) * VAT3( dPC, iip1,jjp1,kk)
2457 - VAT3( oNE, ip1,jp1,k) * VAT3( oPC, iip1,jjp1,kk)
2458 - VAT3( uNE, ip1,jp1,k) * VAT3( uPC, iip1,jjp1,kk));
2461 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2462 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
2463 - VAT3( uS, ip1,jp2,k) * VAT3( oPW, iip1,jjp1,kk)
2464 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPW, iip1,jjp1,kk)
2465 - VAT3( uW, ip2,jp1,k) * VAT3( oPS, iip1,jjp1,kk)
2466 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPS, iip1,jjp1,kk)
2467 - VAT3( uSW, ip2,jp2,k) * VAT3( oPC, iip1,jjp1,kk)
2468 - VAT3( oNE, ip1,jp1,kp1) * VAT3( uPC, iip1,jjp1,kk));
2470 VAT3(XoNE, ii,jj,kk) = TMP1_XONE + TMP2_XONE + TMP3_XONE + TMP4_XONE
2471 + TMP5_XONE + TMP6_XONE;
2481 - VAT3( dPW, ii,jj,kk) * (- VAT3( oNW, im1,j,km1) * VAT3( dPS, iim1,jjp1,kk)
2482 - VAT3( uNW, im1,j,km1) * VAT3( oPS, iim1,jjp1,kk)
2483 - VAT3( oN, im1,j,km1) * VAT3(dPSE, iim1,jjp1,kk)
2484 - VAT3( uN, im1,j,km1) * VAT3(oPSE, iim1,jjp1,kk))
2486 - VAT3( oPW, ii,jj,kk) * (- VAT3( uSE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
2487 - VAT3( oNW, im1,j,k) * VAT3( oPS, iim1,jjp1,kk)
2488 - VAT3( uNW, im1,j,k) * VAT3( uPS, iim1,jjp1,kk)
2489 - VAT3( uS, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2490 - VAT3( oN, im1,j,k) * VAT3(oPSE, iim1,jjp1,kk)
2491 - VAT3( uN, im1,j,k) * VAT3(uPSE, iim1,jjp1,kk));
2494 - VAT3( uPW, ii,jj,kk) * (- VAT3( uSE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
2495 - VAT3( oNW, im1,j,kp1) * VAT3( uPS, iim1,jjp1,kk)
2496 - VAT3( uS, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2497 - VAT3( oN, im1,j,kp1) * VAT3(uPSE, iim1,jjp1,kk))
2499 - VAT3(dPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
2500 - VAT3( uW, im1,jp1,km1) * VAT3( oPS, iim1,jjp1,kk)
2501 - VAT3( oNW, im1,jp1,km1) * VAT3( dPC, iim1,jjp1,kk)
2502 - VAT3( uNW, im1,jp1,km1) * VAT3( oPC, iim1,jjp1,kk)
2503 + VAT3( oC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2504 - VAT3( uC, im1,jp1,km1) * VAT3(oPSE, iim1,jjp1,kk)
2505 - VAT3( oN, im1,jp1,km1) * VAT3( dPE, iim1,jjp1,kk)
2506 - VAT3( uN, im1,jp1,km1) * VAT3( oPE, iim1,jjp1,kk));
2509 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
2510 - VAT3( oE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
2511 - VAT3( uW, im1,jp1,k) * VAT3( uPS, iim1,jjp1,kk)
2512 - VAT3( uSE, im2,jp2,km1) * VAT3( dPC, iim1,jjp1,kk)
2513 - VAT3( oNW, im1,jp1,k) * VAT3( oPC, iim1,jjp1,kk)
2514 - VAT3( uNW, im1,jp1,k) * VAT3( uPC, iim1,jjp1,kk)
2515 - VAT3( uC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2516 + VAT3( oC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2517 - VAT3( uC, im1,jp1,k) * VAT3(uPSE, iim1,jjp1,kk)
2518 - VAT3( uS, im1,jp2,km1) * VAT3( dPE, iim1,jjp1,kk)
2519 - VAT3( oN, im1,jp1,k) * VAT3( oPE, iim1,jjp1,kk)
2520 - VAT3( uN, im1,jp1,k) * VAT3( uPE, iim1,jjp1,kk));
2523 - VAT3(uPNW, ii,jj,kk) * (- VAT3( uE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
2524 - VAT3( oE, im2,jp1,kp1) * VAT3( uPS, iim1,jjp1,kk)
2525 - VAT3( uSE, im2,jp2,k) * VAT3( oPC, iim1,jjp1,kk)
2526 - VAT3( oNW, im1,jp1,kp1) * VAT3( uPC, iim1,jjp1,kk)
2527 - VAT3( uC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2528 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk)
2529 - VAT3( uS, im1,jp2,k) * VAT3( oPE, iim1,jjp1,kk)
2530 - VAT3( oN, im1,jp1,kp1) * VAT3( uPE, iim1,jjp1,kk))
2532 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNW, i,j,km1) * VAT3(dPSE, iim1,jjp1,kk)
2533 - VAT3( uNW, i,j,km1) * VAT3(oPSE, iim1,jjp1,kk));
2536 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2537 - VAT3( oNW, i,j,k) * VAT3(oPSE, iim1,jjp1,kk)
2538 - VAT3( uNW, i,j,k) * VAT3(uPSE, iim1,jjp1,kk))
2540 - VAT3( uPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2541 - VAT3( oNW, i,j,kp1) * VAT3(uPSE, iim1,jjp1,kk))
2543 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2544 - VAT3( uW, i,jp1,km1) * VAT3(oPSE, iim1,jjp1,kk)
2545 - VAT3( oNW, i,jp1,km1) * VAT3( dPE, iim1,jjp1,kk)
2546 - VAT3( uNW, i,jp1,km1) * VAT3( oPE, iim1,jjp1,kk));
2549 - VAT3( oPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2550 - VAT3( oE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2551 - VAT3( uW, i,jp1,k) * VAT3(uPSE, iim1,jjp1,kk)
2552 - VAT3( uSE, im1,jp2,km1) * VAT3( dPE, iim1,jjp1,kk)
2553 - VAT3( oNW, i,jp1,k) * VAT3( oPE, iim1,jjp1,kk)
2554 - VAT3( uNW, i,jp1,k) * VAT3( uPE, iim1,jjp1,kk))
2556 - VAT3( uPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2557 - VAT3( oE, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk)
2558 - VAT3( uSE, im1,jp2,k) * VAT3( oPE, iim1,jjp1,kk)
2559 - VAT3( oNW, i,jp1,kp1) * VAT3( uPE, iim1,jjp1,kk));
2561 VAT3(XoNW, ii,jj,kk) = TMP1_XONW + TMP2_XONW + TMP3_XONW + TMP4_XONW
2562 + TMP5_XONW + TMP6_XONW;
2571 - VAT3( oPS, ii,jj,kk) * (- VAT3( uE, i,jm1,k) * VAT3(dPSW, iip1,jj,kkp1)
2572 - VAT3( uNE, i,jm1,k) * VAT3( dPW, iip1,jj,kkp1))
2574 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2575 - VAT3( uE, i,jm1,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2576 - VAT3( oNE, i,jm1,kp1) * VAT3( dPW, iip1,jj,kkp1)
2577 - VAT3( uNE, i,jm1,kp1) * VAT3( oPW, iip1,jj,kkp1))
2579 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSE, i,j,k) * VAT3(dPSW, iip1,jj,kkp1)
2580 - VAT3( uE, i,j,k) * VAT3( dPW, iip1,jj,kkp1)
2581 - VAT3( uNE, i,j,k) * VAT3(dPNW, iip1,jj,kkp1));
2584 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2585 - VAT3( uSE, i,j,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2586 - VAT3( oE, i,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2587 - VAT3( uE, i,j,kp1) * VAT3( oPW, iip1,jj,kkp1)
2588 - VAT3( oNE, i,j,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2589 - VAT3( uNE, i,j,kp1) * VAT3(oPNW, iip1,jj,kkp1))
2591 - VAT3( oPN, ii,jj,kk) * (- VAT3( uSE, i,jp1,k) * VAT3( dPW, iip1,jj,kkp1)
2592 - VAT3( uE, i,jp1,k) * VAT3(dPNW, iip1,jj,kkp1))
2594 - VAT3( uPN, ii,jj,kk) * (- VAT3( oNW, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2595 - VAT3( uSE, i,jp1,kp1) * VAT3( oPW, iip1,jj,kkp1)
2596 - VAT3( oE, i,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2597 - VAT3( uE, i,jp1,kp1) * VAT3(oPNW, iip1,jj,kkp1));
2600 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,k) * VAT3(dPSW, iip1,jj,kkp1)
2601 - VAT3( uN, ip1,jm1,k) * VAT3( dPW, iip1,jj,kkp1)
2602 - VAT3( uE, ip1,jm1,k) * VAT3( dPS, iip1,jj,kkp1)
2603 - VAT3( uNE, ip1,jm1,k) * VAT3( dPC, iip1,jj,kkp1))
2605 - VAT3(uPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2606 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2607 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPW, iip1,jj,kkp1)
2608 - VAT3( uN, ip1,jm1,kp1) * VAT3( oPW, iip1,jj,kkp1)
2609 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPS, iip1,jj,kkp1)
2610 - VAT3( uE, ip1,jm1,kp1) * VAT3( oPS, iip1,jj,kkp1)
2611 - VAT3( oNE, ip1,jm1,kp1) * VAT3( dPC, iip1,jj,kkp1)
2612 - VAT3( uNE, ip1,jm1,kp1) * VAT3( oPC, iip1,jj,kkp1));
2615 - VAT3( oPE, ii,jj,kk) * (- VAT3( uS, ip1,j,k) * VAT3(dPSW, iip1,jj,kkp1)
2616 - VAT3( uC, ip1,j,k) * VAT3( dPW, iip1,jj,kkp1)
2617 - VAT3( uN, ip1,j,k) * VAT3(dPNW, iip1,jj,kkp1)
2618 - VAT3( uSE, ip1,j,k) * VAT3( dPS, iip1,jj,kkp1)
2619 - VAT3( uE, ip1,j,k) * VAT3( dPC, iip1,jj,kkp1)
2620 - VAT3( uNE, ip1,j,k) * VAT3( dPN, iip1,jj,kkp1));
2623 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2624 - VAT3( uS, ip1,j,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2625 + VAT3( oC, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2626 - VAT3( uC, ip1,j,kp1) * VAT3( oPW, iip1,jj,kkp1)
2627 - VAT3( oN, ip1,j,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2628 - VAT3( uN, ip1,j,kp1) * VAT3(oPNW, iip1,jj,kkp1)
2629 - VAT3( oNW, ip2,jm1,kp1) * VAT3( dPS, iip1,jj,kkp1)
2630 - VAT3( uSE, ip1,j,kp1) * VAT3( oPS, iip1,jj,kkp1)
2631 - VAT3( oE, ip1,j,kp1) * VAT3( dPC, iip1,jj,kkp1)
2632 - VAT3( uE, ip1,j,kp1) * VAT3( oPC, iip1,jj,kkp1)
2633 - VAT3( oNE, ip1,j,kp1) * VAT3( dPN, iip1,jj,kkp1)
2634 - VAT3( uNE, ip1,j,kp1) * VAT3( oPN, iip1,jj,kkp1));
2637 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uS, ip1,jp1,k) * VAT3( dPW, iip1,jj,kkp1)
2638 - VAT3( uC, ip1,jp1,k) * VAT3(dPNW, iip1,jj,kkp1)
2639 - VAT3( uSE, ip1,jp1,k) * VAT3( dPC, iip1,jj,kkp1)
2640 - VAT3( uE, ip1,jp1,k) * VAT3( dPN, iip1,jj,kkp1))
2642 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oN, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2643 - VAT3( uS, ip1,jp1,kp1) * VAT3( oPW, iip1,jj,kkp1)
2644 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2645 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNW, iip1,jj,kkp1)
2646 - VAT3( oNW, ip2,j,kp1) * VAT3( dPC, iip1,jj,kkp1)
2647 - VAT3( uSE, ip1,jp1,kp1) * VAT3( oPC, iip1,jj,kkp1)
2648 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPN, iip1,jj,kkp1)
2649 - VAT3( uE, ip1,jp1,kp1) * VAT3( oPN, iip1,jj,kkp1));
2651 VAT3( XuE, ii,jj,kk) = TMP1_XUE + TMP2_XUE + TMP3_XUE + TMP4_XUE
2652 + TMP5_XUE + TMP6_XUE;
2661 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uW, im1,jm1,k) * VAT3( dPS, iim1,jj,kkp1)
2662 - VAT3( uNW, im1,jm1,k) * VAT3( dPC, iim1,jj,kkp1)
2663 - VAT3( uC, im1,jm1,k) * VAT3(dPSE, iim1,jj,kkp1)
2664 - VAT3( uN, im1,jm1,k) * VAT3( dPE, iim1,jj,kkp1))
2666 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oE, im2,jm1,kp1) * VAT3( dPS, iim1,jj,kkp1)
2667 - VAT3( uW, im1,jm1,kp1) * VAT3( oPS, iim1,jj,kkp1)
2668 - VAT3( oNW, im1,jm1,kp1) * VAT3( dPC, iim1,jj,kkp1)
2669 - VAT3( uNW, im1,jm1,kp1) * VAT3( oPC, iim1,jj,kkp1)
2670 + VAT3( oC, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2671 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2672 - VAT3( oN, im1,jm1,kp1) * VAT3( dPE, iim1,jj,kkp1)
2673 - VAT3( uN, im1,jm1,kp1) * VAT3( oPE, iim1,jj,kkp1));
2676 - VAT3( oPW, ii,jj,kk) * (- VAT3( uSW, im1,j,k) * VAT3( dPS, iim1,jj,kkp1)
2677 - VAT3( uW, im1,j,k) * VAT3( dPC, iim1,jj,kkp1)
2678 - VAT3( uNW, im1,j,k) * VAT3( dPN, iim1,jj,kkp1)
2679 - VAT3( uS, im1,j,k) * VAT3(dPSE, iim1,jj,kkp1)
2680 - VAT3( uC, im1,j,k) * VAT3( dPE, iim1,jj,kkp1)
2681 - VAT3( uN, im1,j,k) * VAT3(dPNE, iim1,jj,kkp1));
2684 - VAT3( uPW, ii,jj,kk) * (- VAT3( oNE, im2,jm1,kp1) * VAT3( dPS, iim1,jj,kkp1)
2685 - VAT3( uSW, im1,j,kp1) * VAT3( oPS, iim1,jj,kkp1)
2686 - VAT3( oE, im2,j,kp1) * VAT3( dPC, iim1,jj,kkp1)
2687 - VAT3( uW, im1,j,kp1) * VAT3( oPC, iim1,jj,kkp1)
2688 - VAT3( oNW, im1,j,kp1) * VAT3( dPN, iim1,jj,kkp1)
2689 - VAT3( uNW, im1,j,kp1) * VAT3( oPN, iim1,jj,kkp1)
2690 - VAT3( oN, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2691 - VAT3( uS, im1,j,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2692 + VAT3( oC, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2693 - VAT3( uC, im1,j,kp1) * VAT3( oPE, iim1,jj,kkp1)
2694 - VAT3( oN, im1,j,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2695 - VAT3( uN, im1,j,kp1) * VAT3(oPNE, iim1,jj,kkp1));
2698 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uSW, im1,jp1,k) * VAT3( dPC, iim1,jj,kkp1)
2699 - VAT3( uW, im1,jp1,k) * VAT3( dPN, iim1,jj,kkp1)
2700 - VAT3( uS, im1,jp1,k) * VAT3( dPE, iim1,jj,kkp1)
2701 - VAT3( uC, im1,jp1,k) * VAT3(dPNE, iim1,jj,kkp1))
2703 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oNE, im2,j,kp1) * VAT3( dPC, iim1,jj,kkp1)
2704 - VAT3( uSW, im1,jp1,kp1) * VAT3( oPC, iim1,jj,kkp1)
2705 - VAT3( oE, im2,jp1,kp1) * VAT3( dPN, iim1,jj,kkp1)
2706 - VAT3( uW, im1,jp1,kp1) * VAT3( oPN, iim1,jj,kkp1)
2707 - VAT3( oN, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2708 - VAT3( uS, im1,jp1,kp1) * VAT3( oPE, iim1,jj,kkp1)
2709 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2710 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNE, iim1,jj,kkp1));
2713 - VAT3( oPS, ii,jj,kk) * (- VAT3( uW, i,jm1,k) * VAT3(dPSE, iim1,jj,kkp1)
2714 - VAT3( uNW, i,jm1,k) * VAT3( dPE, iim1,jj,kkp1))
2716 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2717 - VAT3( uW, i,jm1,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2718 - VAT3( oNW, i,jm1,kp1) * VAT3( dPE, iim1,jj,kkp1)
2719 - VAT3( uNW, i,jm1,kp1) * VAT3( oPE, iim1,jj,kkp1))
2721 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, i,j,k) * VAT3(dPSE, iim1,jj,kkp1)
2722 - VAT3( uW, i,j,k) * VAT3( dPE, iim1,jj,kkp1)
2723 - VAT3( uNW, i,j,k) * VAT3(dPNE, iim1,jj,kkp1));
2726 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2727 - VAT3( uSW, i,j,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2728 - VAT3( oE, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2729 - VAT3( uW, i,j,kp1) * VAT3( oPE, iim1,jj,kkp1)
2730 - VAT3( oNW, i,j,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2731 - VAT3( uNW, i,j,kp1) * VAT3(oPNE, iim1,jj,kkp1))
2733 - VAT3( oPN, ii,jj,kk) * (- VAT3( uSW, i,jp1,k) * VAT3( dPE, iim1,jj,kkp1)
2734 - VAT3( uW, i,jp1,k) * VAT3(dPNE, iim1,jj,kkp1))
2736 - VAT3( uPN, ii,jj,kk) * (- VAT3( oNE, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2737 - VAT3( uSW, i,jp1,kp1) * VAT3( oPE, iim1,jj,kkp1)
2738 - VAT3( oE, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2739 - VAT3( uW, i,jp1,kp1) * VAT3(oPNE, iim1,jj,kkp1));
2741 VAT3( XuW, ii,jj,kk) = TMP1_XUW + TMP2_XUW + TMP3_XUW + TMP4_XUW
2742 + TMP5_XUW + TMP6_XUW;
2751 - VAT3( oPW, ii,jj,kk) * (- VAT3( uN, im1,j,k) * VAT3(dPSW, ii,jjp1,kkp1)
2752 - VAT3( uNE, im1,j,k) * VAT3( dPS, ii,jjp1,kkp1))
2754 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,j,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2755 - VAT3( uN, im1,j,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2756 - VAT3( oNE, im1,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2757 - VAT3( uNE, im1,j,kp1) * VAT3( oPS, ii,jjp1,kkp1))
2759 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,k) * VAT3(dPSW, ii,jjp1,kkp1)
2760 - VAT3( uN, im1,jp1,k) * VAT3( dPW, ii,jjp1,kkp1)
2761 - VAT3( uE, im1,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
2762 - VAT3( uNE, im1,jp1,k) * VAT3( dPC, ii,jjp1,kkp1));
2765 - VAT3(uPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2766 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2767 - VAT3( oN, im1,jp1,kp1) * VAT3( dPW, ii,jjp1,kkp1)
2768 - VAT3( uN, im1,jp1,kp1) * VAT3( oPW, ii,jjp1,kkp1)
2769 - VAT3( oE, im1,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2770 - VAT3( uE, im1,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2771 - VAT3( oNE, im1,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
2772 - VAT3( uNE, im1,jp1,kp1) * VAT3( oPC, ii,jjp1,kkp1))
2774 - VAT3( oPC, ii,jj,kk) * (- VAT3( uNW, i,j,k) * VAT3(dPSW, ii,jjp1,kkp1)
2775 - VAT3( uN, i,j,k) * VAT3( dPS, ii,jjp1,kkp1)
2776 - VAT3( uNE, i,j,k) * VAT3(dPSE, ii,jjp1,kkp1));
2779 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, i,j,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2780 - VAT3( uNW, i,j,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2781 - VAT3( oN, i,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2782 - VAT3( uN, i,j,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2783 - VAT3( oNE, i,j,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2784 - VAT3( uNE, i,j,kp1) * VAT3(oPSE, ii,jjp1,kkp1))
2786 - VAT3( oPN, ii,jj,kk) * (- VAT3( uW, i,jp1,k) * VAT3(dPSW, ii,jjp1,kkp1)
2787 - VAT3( uNW, i,jp1,k) * VAT3( dPW, ii,jjp1,kkp1)
2788 - VAT3( uC, i,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
2789 - VAT3( uN, i,jp1,k) * VAT3( dPC, ii,jjp1,kkp1)
2790 - VAT3( uE, i,jp1,k) * VAT3(dPSE, ii,jjp1,kkp1)
2791 - VAT3( uNE, i,jp1,k) * VAT3( dPE, ii,jjp1,kkp1));
2794 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2795 - VAT3( uW, i,jp1,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2796 - VAT3( oNW, i,jp1,kp1) * VAT3( dPW, ii,jjp1,kkp1)
2797 - VAT3( uNW, i,jp1,kp1) * VAT3( oPW, ii,jjp1,kkp1)
2798 + VAT3( oC, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2799 - VAT3( uC, i,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2800 - VAT3( oN, i,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
2801 - VAT3( uN, i,jp1,kp1) * VAT3( oPC, ii,jjp1,kkp1)
2802 - VAT3( oE, i,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2803 - VAT3( uE, i,jp1,kp1) * VAT3(oPSE, ii,jjp1,kkp1)
2804 - VAT3( oNE, i,jp1,kp1) * VAT3( dPE, ii,jjp1,kkp1)
2805 - VAT3( uNE, i,jp1,kp1) * VAT3( oPE, ii,jjp1,kkp1));
2808 - VAT3( oPE, ii,jj,kk) * (- VAT3( uNW, ip1,j,k) * VAT3( dPS, ii,jjp1,kkp1)
2809 - VAT3( uN, ip1,j,k) * VAT3(dPSE, ii,jjp1,kkp1))
2811 - VAT3( uPE, ii,jj,kk) * (- VAT3( oNW, ip1,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2812 - VAT3( uNW, ip1,j,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2813 - VAT3( oN, ip1,j,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2814 - VAT3( uN, ip1,j,kp1) * VAT3(oPSE, ii,jjp1,kkp1))
2816 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uW, ip1,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
2817 - VAT3( uNW, ip1,jp1,k) * VAT3( dPC, ii,jjp1,kkp1)
2818 - VAT3( uC, ip1,jp1,k) * VAT3(dPSE, ii,jjp1,kkp1)
2819 - VAT3( uN, ip1,jp1,k) * VAT3( dPE, ii,jjp1,kkp1));
2822 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2823 - VAT3( uW, ip1,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2824 - VAT3( oNW, ip1,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
2825 - VAT3( uNW, ip1,jp1,kp1) * VAT3( oPC, ii,jjp1,kkp1)
2826 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2827 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSE, ii,jjp1,kkp1)
2828 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPE, ii,jjp1,kkp1)
2829 - VAT3( uN, ip1,jp1,kp1) * VAT3( oPE, ii,jjp1,kkp1));
2831 VAT3( XuN, ii,jj,kk) = TMP1_XUN + TMP2_XUN + TMP3_XUN + TMP4_XUN
2832 + TMP5_XUN + TMP6_XUN;
2842 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uS, im1,jm1,k) * VAT3( dPW, ii,jjm1,kkp1)
2843 - VAT3( uC, im1,jm1,k) * VAT3(dPNW, ii,jjm1,kkp1)
2844 - VAT3( uSE, im1,jm1,k) * VAT3( dPC, ii,jjm1,kkp1)
2845 - VAT3( uE, im1,jm1,k) * VAT3( dPN, ii,jjm1,kkp1))
2847 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oN, im1,jm2,kp1) * VAT3( dPW, ii,jjm1,kkp1)
2848 - VAT3( uS, im1,jm1,kp1) * VAT3( oPW, ii,jjm1,kkp1)
2849 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2850 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2851 - VAT3( oNW, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
2852 - VAT3( uSE, im1,jm1,kp1) * VAT3( oPC, ii,jjm1,kkp1)
2853 - VAT3( oE, im1,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2854 - VAT3( uE, im1,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1));
2857 - VAT3( oPW, ii,jj,kk) * (- VAT3( uS, im1,j,k) * VAT3(dPNW, ii,jjm1,kkp1)
2858 - VAT3( uSE, im1,j,k) * VAT3( dPN, ii,jjm1,kkp1))
2860 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2861 - VAT3( uS, im1,j,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2862 - VAT3( oNW, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2863 - VAT3( uSE, im1,j,kp1) * VAT3( oPN, ii,jjm1,kkp1))
2865 - VAT3( oPS, ii,jj,kk) * (- VAT3( uSW, i,jm1,k) * VAT3( dPW, ii,jjm1,kkp1)
2866 - VAT3( uW, i,jm1,k) * VAT3(dPNW, ii,jjm1,kkp1)
2867 - VAT3( uS, i,jm1,k) * VAT3( dPC, ii,jjm1,kkp1)
2868 - VAT3( uC, i,jm1,k) * VAT3( dPN, ii,jjm1,kkp1)
2869 - VAT3( uSE, i,jm1,k) * VAT3( dPE, ii,jjm1,kkp1)
2870 - VAT3( uE, i,jm1,k) * VAT3(dPNE, ii,jjm1,kkp1));
2873 - VAT3( uPS, ii,jj,kk) * (- VAT3( oNE, im1,jm2,kp1) * VAT3( dPW, ii,jjm1,kkp1)
2874 - VAT3( uSW, i,jm1,kp1) * VAT3( oPW, ii,jjm1,kkp1)
2875 - VAT3( oE, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2876 - VAT3( uW, i,jm1,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2877 - VAT3( oN, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
2878 - VAT3( uS, i,jm1,kp1) * VAT3( oPC, ii,jjm1,kkp1)
2879 + VAT3( oC, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2880 - VAT3( uC, i,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2881 - VAT3( oNW, ip1,jm2,kp1) * VAT3( dPE, ii,jjm1,kkp1)
2882 - VAT3( uSE, i,jm1,kp1) * VAT3( oPE, ii,jjm1,kkp1)
2883 - VAT3( oE, i,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2884 - VAT3( uE, i,jm1,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2887 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, i,j,k) * VAT3(dPNW, ii,jjm1,kkp1)
2888 - VAT3( uS, i,j,k) * VAT3( dPN, ii,jjm1,kkp1)
2889 - VAT3( uSE, i,j,k) * VAT3(dPNE, ii,jjm1,kkp1))
2891 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2892 - VAT3( uSW, i,j,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2893 - VAT3( oN, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2894 - VAT3( uS, i,j,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2895 - VAT3( oNW, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2896 - VAT3( uSE, i,j,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2899 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uSW, ip1,jm1,k) * VAT3( dPC, ii,jjm1,kkp1)
2900 - VAT3( uW, ip1,jm1,k) * VAT3( dPN, ii,jjm1,kkp1)
2901 - VAT3( uS, ip1,jm1,k) * VAT3( dPE, ii,jjm1,kkp1)
2902 - VAT3( uC, ip1,jm1,k) * VAT3(dPNE, ii,jjm1,kkp1))
2904 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oNE, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
2905 - VAT3( uSW, ip1,jm1,kp1) * VAT3( oPC, ii,jjm1,kkp1)
2906 - VAT3( oE, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2907 - VAT3( uW, ip1,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2908 - VAT3( oN, ip1,jm2,kp1) * VAT3( dPE, ii,jjm1,kkp1)
2909 - VAT3( uS, ip1,jm1,kp1) * VAT3( oPE, ii,jjm1,kkp1)
2910 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2911 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2914 - VAT3( oPE, ii,jj,kk) * (- VAT3( uSW, ip1,j,k) * VAT3( dPN, ii,jjm1,kkp1)
2915 - VAT3( uS, ip1,j,k) * VAT3(dPNE, ii,jjm1,kkp1))
2917 - VAT3( uPE, ii,jj,kk) * (- VAT3( oNE, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2918 - VAT3( uSW, ip1,j,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2919 - VAT3( oN, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2920 - VAT3( uS, ip1,j,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2922 VAT3( XuS, ii,jj,kk) = TMP1_XUS + TMP2_XUS + TMP3_XUS + TMP4_XUS
2923 + TMP5_XUS + TMP6_XUS;
2932 VAT3( oPC, ii,jj,kk) * VAT3( uNE, i,j,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2934 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, i,j,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2935 - VAT3( uNE, i,j,kp1) * VAT3(oPSW, iip1,jjp1,kkp1))
2937 - VAT3( oPN, ii,jj,kk) * (- VAT3( uE, i,jp1,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2938 - VAT3( uNE, i,jp1,k) * VAT3( dPW, iip1,jjp1,kkp1))
2940 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2941 - VAT3( uE, i,jp1,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
2942 - VAT3( oNE, i,jp1,kp1) * VAT3( dPW, iip1,jjp1,kkp1)
2943 - VAT3( uNE, i,jp1,kp1) * VAT3( oPW, iip1,jjp1,kkp1))
2945 - VAT3( oPE, ii,jj,kk) * (- VAT3( uN, ip1,j,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2946 - VAT3( uNE, ip1,j,k) * VAT3( dPS, iip1,jjp1,kkp1));
2949 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,j,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2950 - VAT3( uN, ip1,j,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
2951 - VAT3( oNE, ip1,j,kp1) * VAT3( dPS, iip1,jjp1,kkp1)
2952 - VAT3( uNE, ip1,j,kp1) * VAT3( oPS, iip1,jjp1,kkp1))
2954 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2955 - VAT3( uN, ip1,jp1,k) * VAT3( dPW, iip1,jjp1,kkp1)
2956 - VAT3( uE, ip1,jp1,k) * VAT3( dPS, iip1,jjp1,kkp1)
2957 - VAT3( uNE, ip1,jp1,k) * VAT3( dPC, iip1,jjp1,kkp1))
2959 - VAT3(uPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2960 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
2961 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPW, iip1,jjp1,kkp1)
2962 - VAT3( uN, ip1,jp1,kp1) * VAT3( oPW, iip1,jjp1,kkp1)
2963 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPS, iip1,jjp1,kkp1)
2964 - VAT3( uE, ip1,jp1,kp1) * VAT3( oPS, iip1,jjp1,kkp1)
2965 - VAT3( oNE, ip1,jp1,kp1) * VAT3( dPC, iip1,jjp1,kkp1)
2966 - VAT3( uNE, ip1,jp1,kp1) * VAT3( oPC, iip1,jjp1,kkp1));
2967 VAT3(XuNE, ii,jj,kk) = TMP1_XUNE + TMP2_XUNE;
2977 - VAT3( oPW, ii,jj,kk) * (- VAT3( uNW, im1,j,k) * VAT3( dPS, iim1,jjp1,kkp1)
2978 - VAT3( uN, im1,j,k) * VAT3(dPSE, iim1,jjp1,kkp1))
2980 - VAT3( uPW, ii,jj,kk) * (- VAT3( oNW, im1,j,kp1) * VAT3( dPS, iim1,jjp1,kkp1)
2981 - VAT3( uNW, im1,j,kp1) * VAT3( oPS, iim1,jjp1,kkp1)
2982 - VAT3( oN, im1,j,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
2983 - VAT3( uN, im1,j,kp1) * VAT3(oPSE, iim1,jjp1,kkp1))
2985 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uW, im1,jp1,k) * VAT3( dPS, iim1,jjp1,kkp1)
2986 - VAT3( uNW, im1,jp1,k) * VAT3( dPC, iim1,jjp1,kkp1)
2987 - VAT3( uC, im1,jp1,k) * VAT3(dPSE, iim1,jjp1,kkp1)
2988 - VAT3( uN, im1,jp1,k) * VAT3( dPE, iim1,jjp1,kkp1));
2991 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( dPS, iim1,jjp1,kkp1)
2992 - VAT3( uW, im1,jp1,kp1) * VAT3( oPS, iim1,jjp1,kkp1)
2993 - VAT3( oNW, im1,jp1,kp1) * VAT3( dPC, iim1,jjp1,kkp1)
2994 - VAT3( uNW, im1,jp1,kp1) * VAT3( oPC, iim1,jjp1,kkp1)
2995 + VAT3( oC, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
2996 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSE, iim1,jjp1,kkp1)
2997 - VAT3( oN, im1,jp1,kp1) * VAT3( dPE, iim1,jjp1,kkp1)
2998 - VAT3( uN, im1,jp1,kp1) * VAT3( oPE, iim1,jjp1,kkp1))
3000 + VAT3( oPC, ii,jj,kk) * VAT3( uNW, i,j,k) * VAT3(dPSE, iim1,jjp1,kkp1)
3002 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, i,j,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
3003 - VAT3( uNW, i,j,kp1) * VAT3(oPSE, iim1,jjp1,kkp1))
3005 - VAT3( oPN, ii,jj,kk) * (- VAT3( uW, i,jp1,k) * VAT3(dPSE, iim1,jjp1,kkp1)
3006 - VAT3( uNW, i,jp1,k) * VAT3( dPE, iim1,jjp1,kkp1))
3008 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
3009 - VAT3( uW, i,jp1,kp1) * VAT3(oPSE, iim1,jjp1,kkp1)
3010 - VAT3( oNW, i,jp1,kp1) * VAT3( dPE, iim1,jjp1,kkp1)
3011 - VAT3( uNW, i,jp1,kp1) * VAT3( oPE, iim1,jjp1,kkp1));
3013 VAT3(XuNW, ii,jj,kk) = TMP1_XUNW + TMP2_XUNW;
3023 - VAT3( oPS, ii,jj,kk) * (- VAT3( uSE, i,jm1,k) * VAT3( dPW, iip1,jjm1,kkp1)
3024 - VAT3( uE, i,jm1,k) * VAT3(dPNW, iip1,jjm1,kkp1))
3026 - VAT3( uPS, ii,jj,kk) * (- VAT3( oNW, ip1,jm2,kp1) * VAT3( dPW, iip1,jjm1,kkp1)
3027 - VAT3( uSE, i,jm1,kp1) * VAT3( oPW, iip1,jjm1,kkp1)
3028 - VAT3( oE, i,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3029 - VAT3( uE, i,jm1,kp1) * VAT3(oPNW, iip1,jjm1,kkp1))
3031 + VAT3( oPC, ii,jj,kk) * VAT3( uSE, i,j,k) * VAT3(dPNW, iip1,jjm1,kkp1)
3033 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3034 - VAT3( uSE, i,j,kp1) * VAT3(oPNW, iip1,jjm1,kkp1));
3037 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uS, ip1,jm1,k) * VAT3( dPW, iip1,jjm1,kkp1)
3038 - VAT3( uC, ip1,jm1,k) * VAT3(dPNW, iip1,jjm1,kkp1)
3039 - VAT3( uSE, ip1,jm1,k) * VAT3( dPC, iip1,jjm1,kkp1)
3040 - VAT3( uE, ip1,jm1,k) * VAT3( dPN, iip1,jjm1,kkp1))
3042 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oN, ip1,jm2,kp1) * VAT3( dPW, iip1,jjm1,kkp1)
3043 - VAT3( uS, ip1,jm1,kp1) * VAT3( oPW, iip1,jjm1,kkp1)
3044 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3045 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNW, iip1,jjm1,kkp1)
3046 - VAT3( oNW, ip2,jm2,kp1) * VAT3( dPC, iip1,jjm1,kkp1)
3047 - VAT3( uSE, ip1,jm1,kp1) * VAT3( oPC, iip1,jjm1,kkp1)
3048 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPN, iip1,jjm1,kkp1)
3049 - VAT3( uE, ip1,jm1,kp1) * VAT3( oPN, iip1,jjm1,kkp1))
3051 - VAT3( oPE, ii,jj,kk) * (- VAT3( uS, ip1,j,k) * VAT3(dPNW, iip1,jjm1,kkp1)
3052 - VAT3( uSE, ip1,j,k) * VAT3( dPN, iip1,jjm1,kkp1))
3054 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3055 - VAT3( uS, ip1,j,kp1) * VAT3(oPNW, iip1,jjm1,kkp1)
3056 - VAT3( oNW, ip2,jm1,kp1) * VAT3( dPN, iip1,jjm1,kkp1)
3057 - VAT3( uSE, ip1,j,kp1) * VAT3( oPN, iip1,jjm1,kkp1));
3058 VAT3(XuSE, ii,jj,kk) = TMP1_XUSE + TMP2_XUSE;
3068 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uSW, im1,jm1,k) * VAT3( dPC, iim1,jjm1,kkp1)
3069 - VAT3( uW, im1,jm1,k) * VAT3( dPN, iim1,jjm1,kkp1)
3070 - VAT3( uS, im1,jm1,k) * VAT3( dPE, iim1,jjm1,kkp1)
3071 - VAT3( uC, im1,jm1,k) * VAT3(dPNE, iim1,jjm1,kkp1))
3073 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oNE, im2,jm2,kp1) * VAT3( dPC, iim1,jjm1,kkp1)
3074 - VAT3( uSW, im1,jm1,kp1) * VAT3( oPC, iim1,jjm1,kkp1)
3075 - VAT3( oE, im2,jm1,kp1) * VAT3( dPN, iim1,jjm1,kkp1)
3076 - VAT3( uW, im1,jm1,kp1) * VAT3( oPN, iim1,jjm1,kkp1)
3077 - VAT3( oN, im1,jm2,kp1) * VAT3( dPE, iim1,jjm1,kkp1)
3078 - VAT3( uS, im1,jm1,kp1) * VAT3( oPE, iim1,jjm1,kkp1)
3079 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3080 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
3082 - VAT3( oPW, ii,jj,kk) * (- VAT3( uSW, im1,j,k) * VAT3( dPN, iim1,jjm1,kkp1)
3083 - VAT3( uS, im1,j,k) * VAT3(dPNE, iim1,jjm1,kkp1));
3086 - VAT3( uPW, ii,jj,kk) * (- VAT3( oNE, im2,jm1,kp1) * VAT3( dPN, iim1,jjm1,kkp1)
3087 - VAT3( uSW, im1,j,kp1) * VAT3( oPN, iim1,jjm1,kkp1)
3088 - VAT3( oN, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3089 - VAT3( uS, im1,j,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
3091 - VAT3( oPS, ii,jj,kk) * (- VAT3( uSW, i,jm1,k) * VAT3( dPE, iim1,jjm1,kkp1)
3092 - VAT3( uW, i,jm1,k) * VAT3(dPNE, iim1,jjm1,kkp1))
3094 - VAT3( uPS, ii,jj,kk) * (- VAT3( oNE, im1,jm2,kp1) * VAT3( dPE, iim1,jjm1,kkp1)
3095 - VAT3( uSW, i,jm1,kp1) * VAT3( oPE, iim1,jjm1,kkp1)
3096 - VAT3( oE, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3097 - VAT3( uW, i,jm1,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
3099 + VAT3( oPC, ii,jj,kk) * VAT3( uSW, i,j,k) * VAT3(dPNE, iim1,jjm1,kkp1)
3101 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3102 - VAT3( uSW, i,j,kp1) * VAT3(oPNE, iim1,jjm1,kkp1));
3103 VAT3(XuSW, ii,jj,kk) = TMP1_XUSW + TMP2_XUSW;
VPUBLIC void VbuildG_7(int *nxf, int *nyf, int *nzf, int *nx, int *ny, int *nz, double *oPC, double *oPN, double *oPS, double *oPE, double *oPW, double *oPNE, double *oPNW, double *oPSE, double *oPSW, double *uPC, double *uPN, double *uPS, double *uPE, double *uPW, double *uPNE, double *uPNW, double *uPSE, double *uPSW, double *dPC, double *dPN, double *dPS, double *dPE, double *dPW, double *dPNE, double *dPNW, double *dPSE, double *dPSW, double *oC, double *oE, double *oN, double *uC, double *XoC, double *XoE, double *XoN, double *XuC, double *XoNE, double *XoNW, double *XuE, double *XuW, double *XuN, double *XuS, double *XuNE, double *XuNW, double *XuSE, double *XuSW)
Computes a 27-point galerkin coarse grid matrix from a 7-point fine grid matrix.
VPUBLIC void VbuildG_27(int *nxf, int *nyf, int *nzf, int *nx, int *ny, int *nz, double *oPC, double *oPN, double *oPS, double *oPE, double *oPW, double *oPNE, double *oPNW, double *oPSE, double *oPSW, double *uPC, double *uPN, double *uPS, double *uPE, double *uPW, double *uPNE, double *uPNW, double *uPSE, double *uPSW, double *dPC, double *dPN, double *dPS, double *dPE, double *dPW, double *dPNE, double *dPNW, double *dPSE, double *dPSW, double *oC, double *oE, double *oN, double *uC, double *oNE, double *oNW, double *uE, double *uW, double *uN, double *uS, double *uNE, double *uNW, double *uSE, double *uSW, double *XoC, double *XoE, double *XoN, double *XuC, double *XoNE, double *XoNW, double *XuE, double *XuW, double *XuN, double *XuS, double *XuNE, double *XuNW, double *XuSE, double *XuSW)
Compute a 27-point galerkin coarse grid matrix from a 27-point fine grid matrix.
VPUBLIC void VbuildG(int *nxf, int *nyf, int *nzf, int *nxc, int *nyc, int *nzc, int *numdia, double *pcFF, double *acFF, double *ac)
Build Galerkin matrix structures.
VPUBLIC void VbuildG_1(int *nxf, int *nyf, int *nzf, int *nx, int *ny, int *nz, double *oPC, double *oPN, double *oPS, double *oPE, double *oPW, double *oPNE, double *oPNW, double *oPSE, double *oPSW, double *uPC, double *uPN, double *uPS, double *uPE, double *uPW, double *uPNE, double *uPNW, double *uPSE, double *uPSW, double *dPC, double *dPN, double *dPS, double *dPE, double *dPW, double *dPNE, double *dPNW, double *dPSE, double *dPSW, double *oC, double *XoC, double *XoE, double *XoN, double *XuC, double *XoNE, double *XoNW, double *XuE, double *XuW, double *XuN, double *XuS, double *XuNE, double *XuNW, double *XuSE, double *XuSW)
Computes a 27-point galerkin coarse grid matrix from a 1-point (i.e., diagonal) fine grid matrix.