APBS 3.0.0
Loading...
Searching...
No Matches
buildGd.c
1
55#include "buildGd.h"
56
57VPUBLIC void VbuildG(int *nxf, int *nyf, int *nzf,
58 int *nxc, int *nyc, int *nzc,
59 int *numdia,
60 double *pcFF, double *acFF, double *ac) {
61
62 MAT2(pcFF, *nxc * *nyc * *nzc, 27);
63 MAT2(acFF, *nxf * *nyf * *nzf, 27);
64 MAT2( ac, *nxc * *nyc * *nzc, 27);
65
66 // Call the build routine ***
67 if (*numdia == 1) {
68
70
71 nxf, nyf, nzf, nxc, nyc, nzc,
72
73 RAT2(pcFF, 1, 1), RAT2(pcFF, 1, 2), RAT2(pcFF, 1, 3), RAT2(pcFF, 1, 4), RAT2(pcFF, 1, 5),
74 RAT2(pcFF, 1, 6), RAT2(pcFF, 1, 7), RAT2(pcFF, 1, 8), RAT2(pcFF, 1, 9),
75 RAT2(pcFF, 1, 10), RAT2(pcFF, 1, 11), RAT2(pcFF, 1, 12), RAT2(pcFF, 1, 13), RAT2(pcFF, 1, 14),
76 RAT2(pcFF, 1, 15), RAT2(pcFF, 1, 16), RAT2(pcFF, 1, 17), RAT2(pcFF, 1, 18),
77 RAT2(pcFF, 1, 19), RAT2(pcFF, 1, 20), RAT2(pcFF, 1, 21), RAT2(pcFF, 1, 22), RAT2(pcFF, 1, 23),
78 RAT2(pcFF, 1, 24), RAT2(pcFF, 1, 25), RAT2(pcFF, 1, 26), RAT2(pcFF, 1, 27),
79
80 RAT2(acFF, 1, 1),
81
82 RAT2(ac, 1, 1), RAT2(ac, 1, 2), RAT2(ac, 1, 3),
83 RAT2(ac, 1, 4),
84 RAT2(ac, 1, 5), RAT2(ac, 1, 6),
85 RAT2(ac, 1, 7), RAT2(ac, 1, 8), RAT2(ac, 1, 9), RAT2(ac, 1, 10),
86 RAT2(ac, 1, 11), RAT2(ac, 1, 12), RAT2(ac, 1, 13), RAT2(ac, 1, 14)
87
88 );
89
90 } else if (*numdia == 7) {
91
93
94 nxf, nyf, nzf,
95 nxc, nyc, nzc,
96
97 RAT2(pcFF, 1, 1), RAT2(pcFF, 1, 2), RAT2(pcFF, 1, 3), RAT2(pcFF, 1, 4), RAT2(pcFF, 1, 5),
98 RAT2(pcFF, 1, 6), RAT2(pcFF, 1, 7), RAT2(pcFF, 1, 8), RAT2(pcFF, 1, 9),
99 RAT2(pcFF, 1, 10), RAT2(pcFF, 1, 11), RAT2(pcFF, 1, 12), RAT2(pcFF, 1, 13), RAT2(pcFF, 1, 14),
100 RAT2(pcFF, 1, 15), RAT2(pcFF, 1, 16), RAT2(pcFF, 1, 17), RAT2(pcFF, 1, 18),
101 RAT2(pcFF, 1, 19), RAT2(pcFF, 1, 20), RAT2(pcFF, 1, 21), RAT2(pcFF, 1, 22), RAT2(pcFF, 1, 23),
102 RAT2(pcFF, 1, 24), RAT2(pcFF, 1, 25), RAT2(pcFF, 1, 26), RAT2(pcFF, 1, 27),
103
104 RAT2(acFF, 1, 1), RAT2(acFF, 1, 2), RAT2(acFF, 1, 3), RAT2(acFF, 1, 4),
105
106 RAT2(ac, 1, 1), RAT2(ac, 1, 2), RAT2(ac, 1, 3),
107 RAT2(ac, 1, 4),
108 RAT2(ac, 1, 5), RAT2(ac, 1, 6),
109 RAT2(ac, 1, 7), RAT2(ac, 1, 8), RAT2(ac, 1, 9), RAT2(ac, 1, 10),
110 RAT2(ac, 1, 11), RAT2(ac, 1, 12), RAT2(ac, 1, 13), RAT2(ac, 1, 14)
111
112 );
113
114 } else if (*numdia == 27) {
115
117
118 nxf, nyf, nzf,
119 nxc, nyc, nzc,
120
121 RAT2(pcFF, 1, 1), RAT2(pcFF, 1, 2), RAT2(pcFF, 1, 3), RAT2(pcFF, 1, 4), RAT2(pcFF, 1, 5),
122 RAT2(pcFF, 1, 6), RAT2(pcFF, 1, 7), RAT2(pcFF, 1, 8), RAT2(pcFF, 1, 9),
123 RAT2(pcFF, 1, 10), RAT2(pcFF, 1, 11), RAT2(pcFF, 1, 12), RAT2(pcFF, 1, 13), RAT2(pcFF, 1, 14),
124 RAT2(pcFF, 1, 15), RAT2(pcFF, 1, 16), RAT2(pcFF, 1, 17), RAT2(pcFF, 1, 18),
125 RAT2(pcFF, 1, 19), RAT2(pcFF, 1, 20), RAT2(pcFF, 1, 21), RAT2(pcFF, 1, 22), RAT2(pcFF, 1, 23),
126 RAT2(pcFF, 1, 24), RAT2(pcFF, 1, 25), RAT2(pcFF, 1, 26), RAT2(pcFF, 1, 27),
127
128 RAT2(acFF, 1, 1), RAT2(acFF, 1, 2), RAT2(acFF, 1, 3), RAT2(acFF, 1, 4),
129 RAT2(acFF, 1, 5), RAT2(acFF, 1, 6), RAT2(acFF, 1, 7), RAT2(acFF, 1, 8), RAT2(acFF, 1, 9),
130 RAT2(acFF, 1, 10), RAT2(acFF, 1, 11), RAT2(acFF, 1, 12), RAT2(acFF, 1, 13), RAT2(acFF, 1, 14),
131
132 RAT2(ac, 1, 1), RAT2(ac, 1, 2), RAT2(ac, 1, 3),
133 RAT2(ac, 1, 4),
134 RAT2(ac, 1, 5), RAT2(ac, 1, 6),
135 RAT2(ac, 1, 7), RAT2(ac, 1, 8), RAT2(ac, 1, 9), RAT2(ac, 1, 10),
136 RAT2(ac, 1, 11), RAT2(ac, 1, 12), RAT2(ac, 1, 13), RAT2(ac, 1, 14)
137
138 );
139
140 } else {
141 Vnm_print(2, "BUILDG: invalid stencil type given...\n");
142 }
143}
144
145VPUBLIC void VbuildG_1(int *nxf, int *nyf, int *nzf,
146 int *nx, int *ny, int *nz,
147 double *oPC, double *oPN, double *oPS, double *oPE, double *oPW,
148 double *oPNE, double *oPNW, double *oPSE, double *oPSW,
149 double *uPC, double *uPN, double *uPS, double *uPE, double *uPW,
150 double *uPNE, double *uPNW, double *uPSE, double *uPSW,
151 double *dPC, double *dPN, double *dPS, double *dPE, double *dPW,
152 double *dPNE, double *dPNW, double *dPSE, double *dPSW,
153 double *oC,
154 double *XoC, double *XoE, double *XoN,
155 double *XuC,
156 double *XoNE, double *XoNW,
157 double *XuE, double *XuW, double *XuN, double *XuS,
158 double *XuNE, double *XuNW, double *XuSE, double *XuSW
159 ) {
160
161 int i, j, k, ii, jj, kk;
162 int im1, ip1, jm1, jp1, km1, kp1;
163 int iim1, iip1, jjm1, jjp1, kkm1, kkp1;
164 int nxm1, nym1, nzm1;
165
166 double TMP1_XOC, TMP2_XOC, TMP3_XOC;
167
168 MAT3( oC, *nxf, *nyf, *nzf);
169
170 MAT3( XoC, *nx, *ny, *nz);
171 MAT3( XoE, *nx, *ny, *nz);
172 MAT3( XoN, *nx, *ny, *nz);
173
174 MAT3( XuC, *nx, *ny, *nz);
175
176 MAT3(XoNE, *nx, *ny, *nz);
177 MAT3(XoNW, *nx, *ny, *nz);
178
179 MAT3( XuE, *nx, *ny, *nz);
180 MAT3( XuW, *nx, *ny, *nz);
181 MAT3( XuN, *nx, *ny, *nz);
182 MAT3( XuS, *nx, *ny, *nz);
183 MAT3(XuNE, *nx, *ny, *nz);
184 MAT3(XuNW, *nx, *ny, *nz);
185 MAT3(XuSE, *nx, *ny, *nz);
186 MAT3(XuSW, *nx, *ny, *nz);
187
188 MAT3( oPC, *nx, *ny, *nz);
189 MAT3( oPN, *nx, *ny, *nz);
190 MAT3( oPS, *nx, *ny, *nz);
191 MAT3( oPE, *nx, *ny, *nz);
192 MAT3( oPW, *nx, *ny, *nz);
193 MAT3(oPNE, *nx, *ny, *nz);
194 MAT3(oPNW, *nx, *ny, *nz);
195 MAT3(oPSE, *nx, *ny, *nz);
196 MAT3(oPSW, *nx, *ny, *nz);
197
198 MAT3( uPC, *nx, *ny, *nz);
199 MAT3( uPN, *nx, *ny, *nz);
200 MAT3( uPS, *nx, *ny, *nz);
201 MAT3( uPE, *nx, *ny, *nz);
202 MAT3( uPW, *nx, *ny, *nz);
203 MAT3(uPNE, *nx, *ny, *nz);
204 MAT3(uPNW, *nx, *ny, *nz);
205 MAT3(uPSE, *nx, *ny, *nz);
206 MAT3(uPSW, *nx, *ny, *nz);
207
208 MAT3( dPC, *nx, *ny, *nz);
209 MAT3( dPN, *nx, *ny, *nz);
210 MAT3( dPS, *nx, *ny, *nz);
211 MAT3( dPE, *nx, *ny, *nz);
212 MAT3( dPW, *nx, *ny, *nz);
213 MAT3(dPNE, *nx, *ny, *nz);
214 MAT3(dPNW, *nx, *ny, *nz);
215 MAT3(dPSE, *nx, *ny, *nz);
216 MAT3(dPSW, *nx, *ny, *nz);
217
218 // Define n and determine number of mesh points
219 nxm1 = *nx - 1;
220 nym1 = *ny - 1;
221 nzm1 = *nz - 1;
222
223 //fprintf(data, "%s\n", PRINT_FUNC);
224
225 // Build the operator
226 for(kk=2; kk<=*nz-1; kk++) {
227 k = 2 * kk - 1;
228
229 for(jj=2; jj<=*ny-1; jj++) {
230 j = 2 * jj - 1;
231
232 for(ii=2; ii<=*nx-1; i++) {
233 i = 2 * ii - 1;
234
235 // Index computations
236 im1 = i - 1;
237 ip1 = i + 1;
238 jm1 = j - 1;
239 jp1 = j + 1;
240 km1 = k - 1;
241 kp1 = k + 1;
242 iim1 = ii - 1;
243 iip1 = ii + 1;
244 jjm1 = jj - 1;
245 jjp1 = jj + 1;
246 kkm1 = kk - 1;
247 kkp1 = kk + 1;
248
249 /* *************************************************************
250 * oC
251 * *************************************************************/
252
253 // XoC(ii,jj,kk) =
254 TMP1_XOC =
255 VAT3( uPS, ii,jj,kk) * VAT3( uPS, ii,jj,kk) * VAT3( oC, i,jm1,kp1)
256 + VAT3(dPSW, ii,jj,kk) * VAT3(dPSW, ii,jj,kk) * VAT3( oC, im1,jm1,km1)
257 + VAT3(oPSW, ii,jj,kk) * VAT3(oPSW, ii,jj,kk) * VAT3( oC, im1,jm1,k)
258 + VAT3(uPSW, ii,jj,kk) * VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1)
259 + VAT3( dPW, ii,jj,kk) * VAT3( dPW, ii,jj,kk) * VAT3( oC, im1,j,km1)
260 + VAT3( oPW, ii,jj,kk) * VAT3( oPW, ii,jj,kk) * VAT3( oC, im1,j,k)
261 + VAT3(uPNW, ii,jj,kk) * VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1)
262 + VAT3( dPS, ii,jj,kk) * VAT3( dPS, ii,jj,kk) * VAT3( oC, i,jm1,km1)
263 + VAT3( oPS, ii,jj,kk) * VAT3( oPS, ii,jj,kk) * VAT3( oC, i,jm1,k);
264
265 TMP2_XOC =
266 VAT3( dPC, ii,jj,kk) * VAT3( dPC, ii,jj,kk) * VAT3( oC, i,j,km1)
267 + VAT3( oPC, ii,jj,kk) * VAT3( oPC, ii,jj,kk) * VAT3( oC, i,j,k)
268 + VAT3( uPC, ii,jj,kk) * VAT3( uPC, ii,jj,kk) * VAT3( oC, i,j,kp1)
269 + VAT3( dPN, ii,jj,kk) * VAT3( dPN, ii,jj,kk) * VAT3( oC, i,jp1,km1)
270 + VAT3( oPN, ii,jj,kk) * VAT3( oPN, ii,jj,kk) * VAT3( oC, i,jp1,k)
271 + VAT3( uPW, ii,jj,kk) * VAT3( uPW, ii,jj,kk) * VAT3( oC, im1,j,kp1)
272 + VAT3(dPNW, ii,jj,kk) * VAT3(dPNW, ii,jj,kk) * VAT3( oC, im1,jp1,km1)
273 + VAT3(oPNW, ii,jj,kk) * VAT3(oPNW, ii,jj,kk) * VAT3( oC, im1,jp1,k)
274 + VAT3( oPE, ii,jj,kk) * VAT3( oPE, ii,jj,kk) * VAT3( oC, ip1,j,k);
275
276 TMP3_XOC =
277 VAT3( uPE, ii,jj,kk) * VAT3( uPE, ii,jj,kk) * VAT3( oC, ip1,j,kp1)
278 + VAT3(dPNE, ii,jj,kk) * VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1)
279 + VAT3(oPNE, ii,jj,kk) * VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,k)
280 + VAT3(uPNE, ii,jj,kk) * VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1)
281 + VAT3( uPN, ii,jj,kk) * VAT3( uPN, ii,jj,kk) * VAT3( oC, i,jp1,kp1)
282 + VAT3(dPSE, ii,jj,kk) * VAT3(dPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,km1)
283 + VAT3(oPSE, ii,jj,kk) * VAT3(oPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,k)
284 + VAT3(uPSE, ii,jj,kk) * VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1)
285 + VAT3( dPE, ii,jj,kk) * VAT3( dPE, ii,jj,kk) * VAT3( oC, ip1,j,km1);
286
287 VAT3( XoC, ii,jj,kk) = TMP1_XOC + TMP2_XOC + TMP3_XOC;
288
289 //fprintf(data, "%19.12E\n", VAT3(XoC, ii, jj, kk));
290
291 /* **************************************************************
292 * *** > oE;
293 * **************************************************************/
294
295 VAT3( XoE, ii,jj,kk) =
296 - VAT3(dPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
297 - VAT3(oPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
298 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
299 - VAT3( dPE, ii,jj,kk) * VAT3( oC, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
300 - VAT3( oPE, ii,jj,kk) * VAT3( oC, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
301 - VAT3( uPE, ii,jj,kk) * VAT3( oC, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
302 - VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
303 - VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
304 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPNW, iip1,jj,kk);
305
306 //fprintf(data, "%19.12E\n", VAT3(XoE, ii, jj, kk));
307
308 /* **************************************************************
309 * *** > oN;
310 * **************************************************************/
311
312 VAT3( XoN, ii,jj,kk) =
313 - VAT3(dPNW, ii,jj,kk) * VAT3( oC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
314 - VAT3(oPNW, ii,jj,kk) * VAT3( oC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
315 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
316 - VAT3( dPN, ii,jj,kk) * VAT3( oC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
317 - VAT3( oPN, ii,jj,kk) * VAT3( oC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
318 - VAT3( uPN, ii,jj,kk) * VAT3( oC, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
319 - VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
320 - VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1, k) * VAT3(oPSE, ii,jjp1,kk)
321 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk);
322
323 //fprintf(data, "%19.12E\n", VAT3(XoN, ii, jj, kk));
324
325 /* **************************************************************
326 * *** > uC;
327 * **************************************************************/
328
329 VAT3( XuC, ii,jj,kk) =
330 - VAT3(dPSW, ii,jj,kkp1) * VAT3( oC, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
331 - VAT3( dPW, ii,jj,kkp1) * VAT3( oC, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
332 - VAT3(dPNW, ii,jj,kkp1) * VAT3( oC, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
333 - VAT3( dPS, ii,jj,kkp1) * VAT3( oC, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
334 - VAT3( dPC, ii,jj,kkp1) * VAT3( oC, i,j,kp1) * VAT3( uPC, ii,jj,kk)
335 - VAT3( dPN, ii,jj,kkp1) * VAT3( oC, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
336 - VAT3(dPSE, ii,jj,kkp1) * VAT3( oC, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
337 - VAT3( dPE, ii,jj,kkp1) * VAT3( oC, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
338 - VAT3(dPNE, ii,jj,kkp1) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPNE, ii,jj,kk);
339
340 //fprintf(data, "%19.12E\n", VAT3(XuC, ii, jj, kk));
341
342 /* **************************************************************
343 * *** > oNE;
344 * **************************************************************/
345
346 VAT3(XoNE, ii,jj,kk) =
347 - VAT3(dPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
348 - VAT3(oPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
349 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk);
350
351 //fprintf(data, "%19.12E\n", VAT3(XoNE, ii, jj, kk));
352
353
354 /* **************************************************************
355 * *** > oNW;
356 * **************************************************************/
357
358 VAT3(XoNW, ii,jj,kk) =
359 - VAT3(dPNW, ii,jj,kk) * VAT3( oC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
360 - VAT3(oPNW, ii,jj,kk) * VAT3( oC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
361 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk);
362
363 //fprintf(data, "%19.12E\n", VAT3(XoNW, ii, jj, kk));
364
365 /* **************************************************************
366 * *** > uE;
367 * **************************************************************/
368
369 VAT3( XuE, ii,jj,kk) =
370 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
371 - VAT3( uPE, ii,jj,kk) * VAT3( oC, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
372 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1);
373
374 //fprintf(data, "%19.12E\n", VAT3(XuE, ii, jj, kk));
375
376 /* **************************************************************
377 * *** > uW;
378 * **************************************************************/
379
380 VAT3( XuW, ii,jj,kk) =
381 - VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
382 - VAT3( uPW, ii,jj,kk) * VAT3( oC, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
383 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1);
384
385 //fprintf(data, "%19.12E\n", VAT3(XuW, ii, jj, kk));
386
387 /* **************************************************************
388 * *** > uN;
389 * **************************************************************/
390
391 VAT3( XuN, ii,jj,kk) =
392 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
393 - VAT3( uPN, ii,jj,kk) * VAT3( oC, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
394 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1);
395
396 //fprintf(data, "%19.12E\n", VAT3(XuN, ii, jj, kk));
397
398 /* **************************************************************
399 * *** > uS;
400 * **************************************************************/
401
402 VAT3( XuS, ii,jj,kk) =
403 - VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
404 - VAT3( uPS, ii,jj,kk) * VAT3( oC, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
405 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1);
406
407 //fprintf(data, "%19.12E\n", VAT3(XuS, ii, jj, kk));
408
409 /* **************************************************************
410 * *** > uNE;
411 * **************************************************************/
412
413 VAT3(XuNE, ii,jj,kk) =
414 - VAT3(uPNE, ii,jj,kk) * VAT3( oC, ip1,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1);
415
416 //fprintf(data, "%19.12E\n", VAT3(XuNE, ii, jj, kk));
417
418 /* **************************************************************
419 * *** > uNW;
420 * **************************************************************/
421
422 VAT3(XuNW, ii,jj,kk) =
423 - VAT3(uPNW, ii,jj,kk) * VAT3( oC, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1);
424
425 //fprintf(data, "%19.12E\n", VAT3(XuNW, ii, jj, kk));
426
427 /* **************************************************************
428 * *** > uSE;
429 * **************************************************************/
430
431 VAT3(XuSE, ii,jj,kk) =
432 - VAT3(uPSE, ii,jj,kk) * VAT3( oC, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1);
433
434 //fprintf(data, "%19.12E\n", VAT3(XuSE, ii, jj, kk));
435
436 /* **************************************************************
437 * *** > uSW;
438 * **************************************************************/
439
440 VAT3(XuSW, ii,jj,kk) =
441 - VAT3(uPSW, ii,jj,kk) * VAT3( oC, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1);
442
443 //fprintf(data, "%19.12E\n", VAT3(XuSW, ii, jj, kk));
444
445 }
446 }
447 }
448}
449
450VPUBLIC void VbuildG_7(int *nxf, int *nyf, int *nzf,
451 int *nx, int *ny, int *nz,
452 double *oPC, double *oPN, double *oPS, double *oPE, double *oPW,
453 double *oPNE, double *oPNW, double *oPSE, double *oPSW,
454 double *uPC, double *uPN, double *uPS, double *uPE, double *uPW,
455 double *uPNE, double *uPNW, double *uPSE, double *uPSW,
456 double *dPC, double *dPN, double *dPS, double *dPE, double *dPW,
457 double *dPNE, double *dPNW, double *dPSE, double *dPSW,
458 double *oC, double *oE, double *oN, double *uC,
459 double *XoC, double *XoE, double *XoN,
460 double *XuC,
461 double *XoNE, double *XoNW,
462 double *XuE, double *XuW, double *XuN, double *XuS,
463 double *XuNE, double *XuNW, double *XuSE, double *XuSW) {
464
465 int i, j, k;
466 int ii, jj, kk;
467 int im1, ip1, im2, ip2;
468 int jm1, jp1, jm2, jp2;
469 int km1, kp1, km2, kp2;
470 int iim1, iip1;
471 int jjm1, jjp1;
472 int kkm1, kkp1;
473 int nxm1, nym1, nzm1;
474
475 double TMP1_XOC, TMP2_XOC, TMP3_XOC, TMP4_XOC;
476 double TMP5_XOC, TMP6_XOC, TMP7_XOC, TMP8_XOC;
477 double TMP9_XOC;
478 double TMP1_XOE, TMP2_XOE, TMP3_XOE, TMP4_XOE;
479 double TMP1_XON, TMP2_XON, TMP3_XON, TMP4_XON;
480 double TMP1_XUC, TMP2_XUC, TMP3_XUC, TMP4_XUC;
481
482 MAT3( oC, *nxf, *nyf, *nzf);
483 MAT3( oE, *nxf, *nyf, *nzf);
484 MAT3( oN, *nxf, *nyf, *nzf);
485
486 MAT3( uC, *nxf, *nyf, *nzf);
487
488 MAT3( XoC, *nx, *ny, *nz);
489 MAT3( XoE, *nx, *ny, *nz);
490 MAT3( XoN, *nx, *ny, *nz);
491 MAT3(XoNE, *nx, *ny, *nz);
492 MAT3(XoNW, *nx, *ny, *nz);
493
494 MAT3( XuC, *nx, *ny, *nz);
495 MAT3( XuE, *nx, *ny, *nz);
496 MAT3( XuW, *nx, *ny, *nz);
497 MAT3( XuN, *nx, *ny, *nz);
498 MAT3( XuS, *nx, *ny, *nz);
499 MAT3(XuNE, *nx, *ny, *nz);
500 MAT3(XuNW, *nx, *ny, *nz);
501 MAT3(XuSE, *nx, *ny, *nz);
502 MAT3(XuSW, *nx, *ny, *nz);
503
504 MAT3( oPC, *nx, *ny, *nz);
505 MAT3( oPN, *nx, *ny, *nz);
506 MAT3( oPS, *nx, *ny, *nz);
507 MAT3( oPE, *nx, *ny, *nz);
508 MAT3( oPW, *nx, *ny, *nz);
509 MAT3(oPNE, *nx, *ny, *nz);
510 MAT3(oPNW, *nx, *ny, *nz);
511 MAT3(oPSE, *nx, *ny, *nz);
512 MAT3(oPSW, *nx, *ny, *nz);
513
514 MAT3( uPC, *nx, *ny, *nz);
515 MAT3( uPN, *nx, *ny, *nz);
516 MAT3( uPS, *nx, *ny, *nz);
517 MAT3( uPE, *nx, *ny, *nz);
518 MAT3( uPW, *nx, *ny, *nz);
519 MAT3(uPNE, *nx, *ny, *nz);
520 MAT3(uPNW, *nx, *ny, *nz);
521 MAT3(uPSE, *nx, *ny, *nz);
522 MAT3(uPSW, *nx, *ny, *nz);
523
524 MAT3( dPC, *nx, *ny, *nz);
525 MAT3( dPN, *nx, *ny, *nz);
526 MAT3( dPS, *nx, *ny, *nz);
527 MAT3( dPE, *nx, *ny, *nz);
528 MAT3( dPW, *nx, *ny, *nz);
529 MAT3(dPNE, *nx, *ny, *nz);
530 MAT3(dPNW, *nx, *ny, *nz);
531 MAT3(dPSE, *nx, *ny, *nz);
532 MAT3(dPSW, *nx, *ny, *nz);
533
534 // Define n and determine number of mesh points
535 nxm1 = *nx - 1;
536 nym1 = *ny - 1;
537 nzm1 = *nz - 1;
538
539 //fprintf(data, "%s\n", PRINT_FUNC);
540
541 // Build the operator ***
542 for(kk=2; kk<=*nz-1; kk++) {
543 k = 2 * kk - 1;
544
545 for(jj=2; jj<=*ny-1; jj++) {
546 j = 2 * jj - 1;
547
548 for(ii=2; ii<=*nx-1; ii++) {
549 i = 2 * ii - 1;
550
551 // Index computations
552 im1 = i - 1;
553 ip1 = i + 1;
554 im2 = i - 2;
555 ip2 = i + 2;
556 jm1 = j - 1;
557 jp1 = j + 1;
558 jm2 = j - 2;
559 jp2 = j + 2;
560 km1 = k - 1;
561 kp1 = k + 1;
562 km2 = k - 2;
563 kp2 = k + 2;
564 iim1 = ii - 1;
565 iip1 = ii + 1;
566 jjm1 = jj - 1;
567 jjp1 = jj + 1;
568 kkm1 = kk - 1;
569 kkp1 = kk + 1;
570
571 /* *************************************************************
572 * *** > oC;
573 * *************************************************************/
574
575 // XoC(ii,jj,kk) =
576 TMP1_XOC =
577 VAT3(dPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
578 - VAT3( uC, im1,jm1,km1) * VAT3(oPSW, ii,jj,kk)
579 - VAT3( oN, im1,jm1,km1) * VAT3( dPW, ii,jj,kk)
580 - VAT3( oE, im1,jm1,km1) * VAT3( dPS, ii,jj,kk))
581
582 + VAT3(oPNE, ii,jj,kk) * (- VAT3( oE, i,jp1, k) * VAT3( oPN, ii,jj,kk)
583 - VAT3( oN, ip1, j, k) * VAT3( oPE, ii,jj,kk)
584 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
585 + VAT3( oC, ip1,jp1, k) * VAT3(oPNE, ii,jj,kk)
586 - VAT3( uC, ip1,jp1, k) * VAT3(uPNE, ii,jj,kk))
587
588 + VAT3(dPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
589 - VAT3( oN, ip1, j,km1) * VAT3( dPE, ii,jj,kk)
590 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
591 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNE, ii,jj,kk));
592
593 TMP2_XOC =
594 VAT3(dPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
595 + VAT3( oC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
596 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSE, ii,jj,kk)
597 - VAT3( oN, ip1,jm1,km1) * VAT3( dPE, ii,jj,kk))
598
599 + VAT3( uPE, ii,jj,kk) * (- VAT3( oE, i, j,kp1) * VAT3( uPC, ii,jj,kk)
600 - VAT3( oN, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
601 - VAT3( uC, ip1, j, k) * VAT3( oPE, ii,jj,kk)
602 + VAT3( oC, ip1, j,kp1) * VAT3( uPE, ii,jj,kk)
603 - VAT3( oN, ip1, j,kp1) * VAT3(uPNE, ii,jj,kk))
604
605 + VAT3( oPE, ii,jj,kk) * (- VAT3( oE, i, j, k) * VAT3( oPC, ii,jj,kk)
606 - VAT3( oN, ip1,jm1, k) * VAT3(oPSE, ii,jj,kk)
607 - VAT3( uC, ip1, j,km1) * VAT3( dPE, ii,jj,kk)
608 + VAT3( oC, ip1, j, k) * VAT3( oPE, ii,jj,kk)
609 - VAT3( uC, ip1, j, k) * VAT3( uPE, ii,jj,kk)
610 - VAT3( oN, ip1, j, k) * VAT3(oPNE, ii,jj,kk));
611
612 TMP3_XOC =
613 + VAT3( dPE, ii,jj,kk) * (- VAT3( oE, i, j,km1) * VAT3( dPC, ii,jj,kk)
614 - VAT3( oN, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
615 + VAT3( oC, ip1, j,km1) * VAT3( dPE, ii,jj,kk)
616 - VAT3( uC, ip1, j,km1) * VAT3( oPE, ii,jj,kk)
617 - VAT3( oN, ip1, j,km1) * VAT3(dPNE, ii,jj,kk))
618
619 + VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
620 - VAT3( uC, ip1,jm1, k) * VAT3(oPSE, ii,jj,kk)
621 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
622 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPE, ii,jj,kk))
623
624 + VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
625 - VAT3( oN, ip1, j,kp1) * VAT3( uPE, ii,jj,kk)
626 - VAT3( uC, ip1,jp1, k) * VAT3(oPNE, ii,jj,kk)
627 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNE, ii,jj,kk));
628
629 TMP4_XOC =
630 + VAT3( oPS, ii,jj,kk) * (- VAT3( oE, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
631 - VAT3( uC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
632 + VAT3( oC, i,jm1, k) * VAT3( oPS, ii,jj,kk)
633 - VAT3( uC, i,jm1, k) * VAT3( uPS, ii,jj,kk)
634 - VAT3( oN, i,jm1, k) * VAT3( oPC, ii,jj,kk)
635 - VAT3( oE, i,jm1, k) * VAT3(oPSE, ii,jj,kk))
636
637 + VAT3( dPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
638 + VAT3( oC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
639 - VAT3( uC, i,jm1,km1) * VAT3( oPS, ii,jj,kk)
640 - VAT3( oN, i,jm1,km1) * VAT3( dPC, ii,jj,kk)
641 - VAT3( oE, i,jm1,km1) * VAT3(dPSE, ii,jj,kk))
642
643 + VAT3(oPSE, ii,jj,kk) * (- VAT3( oE, i,jm1, k) * VAT3( oPS, ii,jj,kk)
644 - VAT3( uC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
645 + VAT3( oC, ip1,jm1, k) * VAT3(oPSE, ii,jj,kk)
646 - VAT3( uC, ip1,jm1, k) * VAT3(uPSE, ii,jj,kk)
647 - VAT3( oN, ip1,jm1, k) * VAT3( oPE, ii,jj,kk));
648
649 TMP5_XOC =
650 + VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
651 - VAT3( oN, i, j,km1) * VAT3( dPC, ii,jj,kk)
652 + VAT3( oC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
653 - VAT3( uC, i,jp1,km1) * VAT3( oPN, ii,jj,kk)
654 - VAT3( oE, i,jp1,km1) * VAT3(dPNE, ii,jj,kk))
655
656 + VAT3( uPC, ii,jj,kk) * (- VAT3( oE, im1, j,kp1) * VAT3( uPW, ii,jj,kk)
657 - VAT3( oN, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
658 - VAT3( uC, i, j, k) * VAT3( oPC, ii,jj,kk)
659 + VAT3( oC, i, j,kp1) * VAT3( uPC, ii,jj,kk)
660 - VAT3( oN, i, j,kp1) * VAT3( uPN, ii,jj,kk)
661 - VAT3( oE, i, j,kp1) * VAT3( uPE, ii,jj,kk))
662
663 + VAT3( oPC, ii,jj,kk) * (- VAT3( oE, im1, j, k) * VAT3( oPW, ii,jj,kk)
664 - VAT3( oN, i,jm1, k) * VAT3( oPS, ii,jj,kk)
665 - VAT3( uC, i, j,km1) * VAT3( dPC, ii,jj,kk)
666 + VAT3( oC, i, j, k) * VAT3( oPC, ii,jj,kk)
667 - VAT3( uC, i, j, k) * VAT3( uPC, ii,jj,kk)
668 - VAT3( oN, i, j, k) * VAT3( oPN, ii,jj,kk)
669 - VAT3( oE, i, j, k) * VAT3( oPE, ii,jj,kk));
670
671 TMP6_XOC =
672 + VAT3( dPC, ii,jj,kk) * (- VAT3( oE, im1, j,km1) * VAT3( dPW, ii,jj,kk)
673 - VAT3( oN, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
674 + VAT3( oC, i, j,km1) * VAT3( dPC, ii,jj,kk)
675 - VAT3( uC, i, j,km1) * VAT3( oPC, ii,jj,kk)
676 - VAT3( oN, i, j,km1) * VAT3( dPN, ii,jj,kk)
677 - VAT3( oE, i, j,km1) * VAT3( dPE, ii,jj,kk))
678
679 + VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
680 - VAT3( uC, i,jm1, k) * VAT3( oPS, ii,jj,kk)
681 + VAT3( oC, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
682 - VAT3( oN, i,jm1,kp1) * VAT3( uPC, ii,jj,kk)
683 - VAT3( oE, i,jm1,kp1) * VAT3(uPSE, ii,jj,kk))
684
685 + VAT3(uPNW, ii,jj,kk) * (- VAT3( oN, im1, j,kp1) * VAT3( uPW, ii,jj,kk)
686 - VAT3( uC, im1,jp1, k) * VAT3(oPNW, ii,jj,kk)
687 + VAT3( oC, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
688 - VAT3( oE, im1,jp1,kp1) * VAT3( uPN, ii,jj,kk));
689
690 TMP7_XOC =
691 + VAT3( dPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
692 + VAT3( oC, im1, j,km1) * VAT3( dPW, ii,jj,kk)
693 - VAT3( uC, im1, j,km1) * VAT3( oPW, ii,jj,kk)
694 - VAT3( oN, im1, j,km1) * VAT3(dPNW, ii,jj,kk)
695 - VAT3( oE, im1, j,km1) * VAT3( dPC, ii,jj,kk))
696
697 + VAT3(uPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
698 + VAT3( oC, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
699 - VAT3( oN, im1,jm1,kp1) * VAT3( uPW, ii,jj,kk)
700 - VAT3( oE, im1,jm1,kp1) * VAT3( uPS, ii,jj,kk))
701
702 + VAT3(oPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
703 + VAT3( oC, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
704 - VAT3( uC, im1,jm1, k) * VAT3(uPSW, ii,jj,kk)
705 - VAT3( oN, im1,jm1, k) * VAT3( oPW, ii,jj,kk)
706 - VAT3( oE, im1,jm1, k) * VAT3( oPS, ii,jj,kk));
707
708 TMP8_XOC =
709 + VAT3(oPNW, ii,jj,kk) * (- VAT3( oN, im1, j, k) * VAT3( oPW, ii,jj,kk)
710 - VAT3( uC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
711 + VAT3( oC, im1,jp1, k) * VAT3(oPNW, ii,jj,kk)
712 - VAT3( uC, im1,jp1, k) * VAT3(uPNW, ii,jj,kk)
713 - VAT3( oE, im1,jp1, k) * VAT3( oPN, ii,jj,kk))
714
715 + VAT3(dPNW, ii,jj,kk) * (- VAT3( oN, im1, j,km1) * VAT3( dPW, ii,jj,kk)
716 + VAT3( oC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
717 - VAT3( uC, im1,jp1,km1) * VAT3(oPNW, ii,jj,kk)
718 - VAT3( oE, im1,jp1,km1) * VAT3( dPN, ii,jj,kk))
719
720 + VAT3( oPW, ii,jj,kk) * (- VAT3( oN, im1,jm1, k) * VAT3(oPSW, ii,jj,kk)
721 - VAT3( uC, im1, j,km1) * VAT3( dPW, ii,jj,kk)
722 + VAT3( oC, im1, j, k) * VAT3( oPW, ii,jj,kk)
723 - VAT3( uC, im1, j, k) * VAT3( uPW, ii,jj,kk)
724 - VAT3( oN, im1, j, k) * VAT3(oPNW, ii,jj,kk)
725 - VAT3( oE, im1, j, k) * VAT3( oPC, ii,jj,kk));
726
727 TMP9_XOC =
728 + VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
729 - VAT3( uC, im1, j, k) * VAT3( oPW, ii,jj,kk)
730 + VAT3( oC, im1, j,kp1) * VAT3( uPW, ii,jj,kk)
731 - VAT3( oN, im1, j,kp1) * VAT3(uPNW, ii,jj,kk)
732 - VAT3( oE, im1, j,kp1) * VAT3( uPC, ii,jj,kk))
733
734 + VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
735 - VAT3( oN, i, j,kp1) * VAT3( uPC, ii,jj,kk)
736 - VAT3( uC, i,jp1, k) * VAT3( oPN, ii,jj,kk)
737 + VAT3( oC, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
738 - VAT3( oE, i,jp1,kp1) * VAT3(uPNE, ii,jj,kk))
739
740 + VAT3( oPN, ii,jj,kk) * (- VAT3( oE, im1,jp1, k) * VAT3(oPNW, ii,jj,kk)
741 - VAT3( oN, i, j, k) * VAT3( oPC, ii,jj,kk)
742 - VAT3( uC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
743 + VAT3( oC, i,jp1, k) * VAT3( oPN, ii,jj,kk)
744 - VAT3( uC, i,jp1, k) * VAT3( uPN, ii,jj,kk)
745 - VAT3( oE, i,jp1, k) * VAT3(oPNE, ii,jj,kk));
746
747 VAT3( XoC, ii,jj,kk) = TMP1_XOC + TMP2_XOC + TMP3_XOC
748 + TMP4_XOC + TMP5_XOC + TMP6_XOC
749 + TMP7_XOC + TMP8_XOC + TMP9_XOC;
750
751 //fprintf(data, "%19.12E\n", VAT3(XoC, ii, jj, kk));
752
753 /* *************************************************************
754 * *** > oE;
755 * *************************************************************/
756
757 // VAT3( XoE, ii,jj,kk) =
758 TMP1_XOE =
759 VAT3( dPS, ii,jj,kk) * VAT3( oE, i,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
760 + VAT3( oPS, ii,jj,kk) * VAT3( oE, i,jm1, k) * VAT3(oPSW, iip1,jj,kk)
761 + VAT3( uPS, ii,jj,kk) * VAT3( oE, i,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
762 + VAT3( dPC, ii,jj,kk) * VAT3( oE, i, j,km1) * VAT3( dPW, iip1,jj,kk)
763 + VAT3( oPC, ii,jj,kk) * VAT3( oE, i, j, k) * VAT3( oPW, iip1,jj,kk)
764 + VAT3( uPC, ii,jj,kk) * VAT3( oE, i, j,kp1) * VAT3( uPW, iip1,jj,kk)
765 + VAT3( dPN, ii,jj,kk) * VAT3( oE, i,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
766 + VAT3( oPN, ii,jj,kk) * VAT3( oE, i,jp1, k) * VAT3(oPNW, iip1,jj,kk)
767 + VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(uPNW, iip1,jj,kk)
768
769 - VAT3(dPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
770 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSW, iip1,jj,kk)
771 - VAT3( oN, ip1,jm1,km1) * VAT3( dPW, iip1,jj,kk)
772 - VAT3( oE, ip1,jm1,km1) * VAT3( dPS, iip1,jj,kk));
773
774 TMP2_XOE =
775 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
776 + VAT3( oC, ip1,jm1, k) * VAT3(oPSW, iip1,jj,kk)
777 - VAT3( uC, ip1,jm1, k) * VAT3(uPSW, iip1,jj,kk)
778 - VAT3( oN, ip1,jm1, k) * VAT3( oPW, iip1,jj,kk)
779 - VAT3( oE, ip1,jm1, k) * VAT3( oPS, iip1,jj,kk))
780
781 - VAT3(uPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1, k) * VAT3(oPSW, iip1,jj,kk)
782 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
783 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPW, iip1,jj,kk)
784 - VAT3( oE, ip1,jm1,kp1) * VAT3( uPS, iip1,jj,kk))
785
786 - VAT3( dPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
787 + VAT3( oC, ip1, j,km1) * VAT3( dPW, iip1,jj,kk)
788 - VAT3( uC, ip1, j,km1) * VAT3( oPW, iip1,jj,kk)
789 - VAT3( oN, ip1, j,km1) * VAT3(dPNW, iip1,jj,kk)
790 - VAT3( oE, ip1, j,km1) * VAT3( dPC, iip1,jj,kk));
791
792 TMP3_XOE =
793 - VAT3( oPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1, k) * VAT3(oPSW, iip1,jj,kk)
794 - VAT3( uC, ip1, j,km1) * VAT3( dPW, iip1,jj,kk)
795 + VAT3( oC, ip1, j, k) * VAT3( oPW, iip1,jj,kk)
796 - VAT3( uC, ip1, j, k) * VAT3( uPW, iip1,jj,kk)
797 - VAT3( oN, ip1, j, k) * VAT3(oPNW, iip1,jj,kk)
798 - VAT3( oE, ip1, j, k) * VAT3( oPC, iip1,jj,kk))
799
800 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
801 - VAT3( uC, ip1, j, k) * VAT3( oPW, iip1,jj,kk)
802 + VAT3( oC, ip1, j,kp1) * VAT3( uPW, iip1,jj,kk)
803 - VAT3( oN, ip1, j,kp1) * VAT3(uPNW, iip1,jj,kk)
804 - VAT3( oE, ip1, j,kp1) * VAT3( uPC, iip1,jj,kk))
805
806 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oN, ip1, j,km1) * VAT3( dPW, iip1,jj,kk)
807 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
808 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNW, iip1,jj,kk)
809 - VAT3( oE, ip1,jp1,km1) * VAT3( dPN, iip1,jj,kk));
810
811 TMP4_XOE =
812 - VAT3(oPNE, ii,jj,kk) * (- VAT3( oN, ip1, j, k) * VAT3( oPW, iip1,jj,kk)
813 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
814 + VAT3( oC, ip1,jp1, k) * VAT3(oPNW, iip1,jj,kk)
815 - VAT3( uC, ip1,jp1, k) * VAT3(uPNW, iip1,jj,kk)
816 - VAT3( oE, ip1,jp1, k) * VAT3( oPN, iip1,jj,kk))
817
818 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oN, ip1, j,kp1) * VAT3( uPW, iip1,jj,kk)
819 - VAT3( uC, ip1,jp1, k) * VAT3(oPNW, iip1,jj,kk)
820 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNW, iip1,jj,kk)
821 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPN, iip1,jj,kk));
822
823 VAT3( XoE, ii,jj,kk) = TMP1_XOE + TMP2_XOE + TMP3_XOE + TMP4_XOE;
824
825 //fprintf(data, "%19.12E\n", VAT3(XoE, ii, jj, kk));
826
828
829 /* *************************************************************
830 * *** > oN;
831 * *************************************************************/
832
833 // VAT3( XoN, ii,jj,kk) =
834 TMP1_XON =
835 VAT3( dPW, ii,jj,kk) * VAT3( oN, im1,j,km1) * VAT3(dPSW, ii,jjp1,kk)
836 + VAT3( oPW, ii,jj,kk) * VAT3( oN, im1,j,k) * VAT3(oPSW, ii,jjp1,kk)
837 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(uPSW, ii,jjp1,kk)
838
839 - VAT3(dPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
840 - VAT3( uC, im1,jp1,km1) * VAT3(oPSW, ii,jjp1,kk)
841 - VAT3( oN, im1,jp1,km1) * VAT3( dPW, ii,jjp1,kk)
842 - VAT3( oE, im1,jp1,km1) * VAT3( dPS, ii,jjp1,kk))
843
844 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
845 + VAT3( oC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
846 - VAT3( uC, im1,jp1,k) * VAT3(uPSW, ii,jjp1,kk)
847 - VAT3( oN, im1,jp1,k) * VAT3( oPW, ii,jjp1,kk)
848 - VAT3( oE, im1,jp1,k) * VAT3( oPS, ii,jjp1,kk));
849
850 TMP2_XON =
851 - VAT3(uPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
852 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
853 - VAT3( oN, im1,jp1,kp1) * VAT3( uPW, ii,jjp1,kk)
854 - VAT3( oE, im1,jp1,kp1) * VAT3( uPS, ii,jjp1,kk))
855
856 + VAT3( dPC, ii,jj,kk) * VAT3( oN, i,j,km1) * VAT3( dPS, ii,jjp1,kk)
857 + VAT3( oPC, ii,jj,kk) * VAT3( oN, i,j,k) * VAT3( oPS, ii,jjp1,kk)
858 + VAT3( uPC, ii,jj,kk) * VAT3( oN, i,j,kp1) * VAT3( uPS, ii,jjp1,kk)
859
860 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
861 + VAT3( oC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
862 - VAT3( uC, i,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
863 - VAT3( oN, i,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
864 - VAT3( oE, i,jp1,km1) * VAT3(dPSE, ii,jjp1,kk));
865
866 TMP3_XON =
867 - VAT3( oPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
868 - VAT3( uC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
869 + VAT3( oC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
870 - VAT3( uC, i,jp1,k) * VAT3( uPS, ii,jjp1,kk)
871 - VAT3( oN, i,jp1,k) * VAT3( oPC, ii,jjp1,kk)
872 - VAT3( oE, i,jp1,k) * VAT3(oPSE, ii,jjp1,kk))
873
874 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
875 - VAT3( uC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
876 + VAT3( oC, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
877 - VAT3( oN, i,jp1,kp1) * VAT3( uPC, ii,jjp1,kk)
878 - VAT3( oE, i,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk))
879
880 + VAT3( dPE, ii,jj,kk) * VAT3( oN, ip1,j,km1) * VAT3(dPSE, ii,jjp1,kk)
881 + VAT3( oPE, ii,jj,kk) * VAT3( oN, ip1,j,k) * VAT3(oPSE, ii,jjp1,kk)
882 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(uPSE, ii,jjp1,kk);
883
884 TMP4_XON =
885 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
886 + VAT3( oC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
887 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSE, ii,jjp1,kk)
888 - VAT3( oN, ip1,jp1,km1) * VAT3( dPE, ii,jjp1,kk))
889
890 - VAT3(oPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
891 - VAT3( uC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
892 + VAT3( oC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
893 - VAT3( uC, ip1,jp1,k) * VAT3(uPSE, ii,jjp1,kk)
894 - VAT3( oN, ip1,jp1,k) * VAT3( oPE, ii,jjp1,kk))
895
896 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
897 - VAT3( uC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
898 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk)
899 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPE, ii,jjp1,kk));
900
901 VAT3( XoN, ii,jj,kk) = TMP1_XON + TMP2_XON + TMP3_XON + TMP4_XON;
902
903 //fprintf(data, "%19.12E\n", VAT3(XoN, ii, jj, kk));
904
905 /* *************************************************************
906 * *** > uC;
907 * *************************************************************/
908
909 // VAT3( XuC, ii,jj,kk) =
910 TMP1_XUC =
911 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPSW, ii,jj,kkp1)
912
913 - VAT3(uPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
914 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSW, ii,jj,kkp1)
915 - VAT3( oN, im1,jm1,kp1) * VAT3( dPW, ii,jj,kkp1)
916 - VAT3( oE, im1,jm1,kp1) * VAT3( dPS, ii,jj,kkp1))
917
918 + VAT3( oPW, ii,jj,kk) * VAT3( uC, im1,j,k) * VAT3( dPW, ii,jj,kkp1)
919
920 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
921 + VAT3( oC, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
922 - VAT3( uC, im1,j,kp1) * VAT3( oPW, ii,jj,kkp1)
923 - VAT3( oN, im1,j,kp1) * VAT3(dPNW, ii,jj,kkp1)
924 - VAT3( oE, im1,j,kp1) * VAT3( dPC, ii,jj,kkp1))
925
926 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPNW, ii,jj,kkp1);
927
928 TMP2_XUC =
929 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oN, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
930 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
931 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNW, ii,jj,kkp1)
932 - VAT3( oE, im1,jp1,kp1) * VAT3( dPN, ii,jj,kkp1))
933
934 + VAT3( oPS, ii,jj,kk) * VAT3( uC, i,jm1,k) * VAT3( dPS, ii,jj,kkp1)
935
936 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
937 + VAT3( oC, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
938 - VAT3( uC, i,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
939 - VAT3( oN, i,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
940 - VAT3( oE, i,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1))
941
942 + VAT3( oPC, ii,jj,kk) * VAT3( uC, i,j,k) * VAT3( dPC, ii,jj,kkp1)
943
944 - VAT3( uPC, ii,jj,kk) * (- VAT3( oE, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
945 - VAT3( oN, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
946 + VAT3( oC, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
947 - VAT3( uC, i,j,kp1) * VAT3( oPC, ii,jj,kkp1)
948 - VAT3( oN, i,j,kp1) * VAT3( dPN, ii,jj,kkp1)
949 - VAT3( oE, i,j,kp1) * VAT3( dPE, ii,jj,kkp1));
950
951 TMP3_XUC =
952 + VAT3( oPN, ii,jj,kk) * VAT3( uC, i,jp1,k) * VAT3( dPN, ii,jj,kkp1)
953
954 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
955 - VAT3( oN, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
956 + VAT3( oC, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
957 - VAT3( uC, i,jp1,kp1) * VAT3( oPN, ii,jj,kkp1)
958 - VAT3( oE, i,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1))
959
960 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPSE, ii,jj,kkp1)
961
962 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
963 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
964 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSE, ii,jj,kkp1)
965 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPE, ii,jj,kkp1));
966
967 TMP4_XUC =
968 + VAT3( oPE, ii,jj,kk) * VAT3( uC, ip1,j,k) * VAT3( dPE, ii,jj,kkp1)
969
970 - VAT3( uPE, ii,jj,kk) * (- VAT3( oE, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
971 - VAT3( oN, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
972 + VAT3( oC, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
973 - VAT3( uC, ip1,j,kp1) * VAT3( oPE, ii,jj,kkp1)
974 - VAT3( oN, ip1,j,kp1) * VAT3(dPNE, ii,jj,kkp1))
975
976 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPNE, ii,jj,kkp1)
977
978 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
979 - VAT3( oN, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
980 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1)
981 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNE, ii,jj,kkp1));
982
983 VAT3( XuC, ii,jj,kk) = TMP1_XUC + TMP2_XUC + TMP3_XUC + TMP4_XUC;
984
985 //fprintf(data, "%19.12E\n", VAT3(XuC, ii, jj, kk));
986
987 /* *************************************************************
988 * *** > oNE;
989 * *************************************************************/
990
991 VAT3(XoNE, ii,jj,kk) =
992 VAT3( dPN, ii,jj,kk) * VAT3( oE, i,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
993 + VAT3( oPN, ii,jj,kk) * VAT3( oE, i,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
994 + VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
995 + VAT3( dPE, ii,jj,kk) * VAT3( oN, ip1,j,km1) * VAT3(dPSW, iip1,jjp1,kk)
996 + VAT3( oPE, ii,jj,kk) * VAT3( oN, ip1,j,k) * VAT3(oPSW, iip1,jjp1,kk)
997 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(uPSW, iip1,jjp1,kk)
998
999 - VAT3(dPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
1000 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSW, iip1,jjp1,kk)
1001 - VAT3( oN, ip1,jp1,km1) * VAT3( dPW, iip1,jjp1,kk)
1002 - VAT3( oE, ip1,jp1,km1) * VAT3( dPS, iip1,jjp1,kk))
1003
1004 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
1005 + VAT3( oC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
1006 - VAT3( uC, ip1,jp1,k) * VAT3(uPSW, iip1,jjp1,kk)
1007 - VAT3( oN, ip1,jp1,k) * VAT3( oPW, iip1,jjp1,kk)
1008 - VAT3( oE, ip1,jp1,k) * VAT3( oPS, iip1,jjp1,kk))
1009
1010 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
1011 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
1012 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPW, iip1,jjp1,kk)
1013 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPS, iip1,jjp1,kk));
1014
1015 //fprintf(data, "%19.12E\n", VAT3(XoNE, ii, jj, kk));
1016
1017 /* *************************************************************
1018 * *** > oNW;
1019 * *************************************************************/
1020
1021 VAT3(XoNW, ii,jj,kk) =
1022 VAT3( dPW, ii,jj,kk) * VAT3( oN, im1,j,km1) * VAT3(dPSE, iim1,jjp1,kk)
1023 + VAT3( oPW, ii,jj,kk) * VAT3( oN, im1,j,k) * VAT3(oPSE, iim1,jjp1,kk)
1024 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(uPSE, iim1,jjp1,kk)
1025
1026 - VAT3(dPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
1027 + VAT3( oC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
1028 - VAT3( uC, im1,jp1,km1) * VAT3(oPSE, iim1,jjp1,kk)
1029 - VAT3( oN, im1,jp1,km1) * VAT3( dPE, iim1,jjp1,kk))
1030
1031 - VAT3(oPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
1032 - VAT3( uC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
1033 + VAT3( oC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
1034 - VAT3( uC, im1,jp1,k) * VAT3(uPSE, iim1,jjp1,kk)
1035 - VAT3( oN, im1,jp1,k) * VAT3( oPE, iim1,jjp1,kk))
1036
1037 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( uPS, iim1,jjp1,kk)
1038 - VAT3( uC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
1039 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk)
1040 - VAT3( oN, im1,jp1,kp1) * VAT3( uPE, iim1,jjp1,kk))
1041
1042 + VAT3( dPN, ii,jj,kk) * VAT3( oE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
1043 + VAT3( oPN, ii,jj,kk) * VAT3( oE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
1044 + VAT3( uPN, ii,jj,kk) * VAT3( oE, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk);
1045
1046 //fprintf(data, "%19.12E\n", VAT3(XoNW, ii, jj, kk));
1047
1048 /* *************************************************************
1049 * *** > uE;
1050 * *************************************************************/
1051
1052 VAT3( XuE, ii,jj,kk) =
1053 VAT3( uPS, ii,jj,kk) * VAT3( oE, i,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
1054 + VAT3( uPC, ii,jj,kk) * VAT3( oE, i,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
1055 + VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
1056 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPSW, iip1,jj,kkp1)
1057
1058 - VAT3(uPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
1059 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSW, iip1,jj,kkp1)
1060 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPW, iip1,jj,kkp1)
1061 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPS, iip1,jj,kkp1))
1062
1063 + VAT3( oPE, ii,jj,kk) * VAT3( uC, ip1,j,k) * VAT3( dPW, iip1,jj,kkp1)
1064
1065 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
1066 + VAT3( oC, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
1067 - VAT3( uC, ip1,j,kp1) * VAT3( oPW, iip1,jj,kkp1)
1068 - VAT3( oN, ip1,j,kp1) * VAT3(dPNW, iip1,jj,kkp1)
1069 - VAT3( oE, ip1,j,kp1) * VAT3( dPC, iip1,jj,kkp1))
1070
1071 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPNW, iip1,jj,kkp1)
1072
1073 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oN, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
1074 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
1075 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNW, iip1,jj,kkp1)
1076 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPN, iip1,jj,kkp1));
1077
1078 //fprintf(data, "%19.12E\n", VAT3(XuE, ii, jj, kk));
1079
1080 /* *************************************************************
1081 * *** > uW;
1082 * *************************************************************/
1083
1084 VAT3( XuW, ii,jj,kk) =
1085 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPSE, iim1,jj,kkp1)
1086
1087 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oE, im2,jm1,kp1) * VAT3( dPS, iim1,jj,kkp1)
1088 + VAT3( oC, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
1089 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSE, iim1,jj,kkp1)
1090 - VAT3( oN, im1,jm1,kp1) * VAT3( dPE, iim1,jj,kkp1))
1091
1092 + VAT3( oPW, ii,jj,kk) * VAT3( uC, im1,j,k) * VAT3( dPE, iim1,jj,kkp1)
1093
1094 - VAT3( uPW, ii,jj,kk) * (- VAT3( oE, im2,j,kp1) * VAT3( dPC, iim1,jj,kkp1)
1095 - VAT3( oN, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
1096 + VAT3( oC, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
1097 - VAT3( uC, im1,j,kp1) * VAT3( oPE, iim1,jj,kkp1)
1098 - VAT3( oN, im1,j,kp1) * VAT3(dPNE, iim1,jj,kkp1))
1099
1100 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPNE, iim1,jj,kkp1)
1101
1102 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( dPN, iim1,jj,kkp1)
1103 - VAT3( oN, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
1104 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1)
1105 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNE, iim1,jj,kkp1))
1106
1107 + VAT3( uPS, ii,jj,kk) * VAT3( oE, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
1108 + VAT3( uPC, ii,jj,kk) * VAT3( oE, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
1109 + VAT3( uPN, ii,jj,kk) * VAT3( oE, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1);
1110
1111 //fprintf(data, "%19.12E\n", VAT3(XuW, ii, jj, kk));
1112
1113 /* *************************************************************
1114 * *** > uN;
1115 * *************************************************************/
1116
1117 VAT3( XuN, ii,jj,kk) =
1118 VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
1119 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPSW, ii,jjp1,kkp1)
1120
1121 - VAT3(uPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
1122 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
1123 - VAT3( oN, im1,jp1,kp1) * VAT3( dPW, ii,jjp1,kkp1)
1124 - VAT3( oE, im1,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1))
1125
1126 + VAT3( uPC, ii,jj,kk) * VAT3( oN, i,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
1127 + VAT3( oPN, ii,jj,kk) * VAT3( uC, i,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
1128
1129 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
1130 + VAT3( oC, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
1131 - VAT3( uC, i,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
1132 - VAT3( oN, i,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
1133 - VAT3( oE, i,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1))
1134
1135 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
1136 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPSE, ii,jjp1,kkp1)
1137
1138 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
1139 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
1140 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSE, ii,jjp1,kkp1)
1141 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPE, ii,jjp1,kkp1));
1142
1143 //fprintf(data, "%19.12E\n", VAT3(XuN, ii, jj, kk));
1144
1145 /* *************************************************************
1146 * *** > uS;
1147 * *************************************************************/
1148
1149 VAT3( XuS, ii,jj,kk) =
1150 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPNW, ii,jjm1,kkp1)
1151
1152 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oN, im1,jm2,kp1) * VAT3( dPW, ii,jjm1,kkp1)
1153 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
1154 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
1155 - VAT3( oE, im1,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1))
1156
1157 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
1158 + VAT3( oPS, ii,jj,kk) * VAT3( uC, i,jm1,k) * VAT3( dPN, ii,jjm1,kkp1)
1159
1160 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
1161 - VAT3( oN, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
1162 + VAT3( oC, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
1163 - VAT3( uC, i,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1)
1164 - VAT3( oE, i,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1))
1165
1166 + VAT3( uPC, ii,jj,kk) * VAT3( oN, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
1167 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPNE, ii,jjm1,kkp1)
1168
1169 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
1170 - VAT3( oN, ip1,jm2,kp1) * VAT3( dPE, ii,jjm1,kkp1)
1171 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
1172 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNE, ii,jjm1,kkp1))
1173
1174 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1);
1175
1176 //fprintf(data, "%19.12E\n", VAT3(XuS, ii, jj, kk));
1177
1178 /* *************************************************************
1179 * *** > uNE;
1180 * *************************************************************/
1181
1182 VAT3(XuNE, ii,jj,kk) =
1183 VAT3( uPN, ii,jj,kk) * VAT3( oE, i,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
1184 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,j,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
1185 + VAT3(oPNE, ii,jj,kk) * VAT3( uC, ip1,jp1,k) * VAT3(dPSW, iip1,jjp1,kkp1)
1186
1187 - VAT3(uPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
1188 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
1189 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPW, iip1,jjp1,kkp1)
1190 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPS, iip1,jjp1,kkp1));
1191
1192 //fprintf(data, "%19.12E\n", VAT3(XuNE, ii, jj, kk));
1193
1194 /* *************************************************************
1195 * *** > uNW;
1196 * *************************************************************/
1197
1198 VAT3(XuNW, ii,jj,kk) =
1199 VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,j,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
1200 + VAT3(oPNW, ii,jj,kk) * VAT3( uC, im1,jp1,k) * VAT3(dPSE, iim1,jjp1,kkp1)
1201
1202 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( dPS, iim1,jjp1,kkp1)
1203 + VAT3( oC, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
1204 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSE, iim1,jjp1,kkp1)
1205 - VAT3( oN, im1,jp1,kp1) * VAT3( dPE, iim1,jjp1,kkp1))
1206
1207 + VAT3( uPN, ii,jj,kk) * VAT3( oE, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1);
1208
1209 //fprintf(data, "%19.12E\n", VAT3(XuNW, ii, jj, kk));
1210
1211 /* *************************************************************
1212 * *** > uSE;
1213 * *************************************************************/
1214
1215 VAT3(XuSE, ii,jj,kk) =
1216 VAT3( uPS, ii,jj,kk) * VAT3( oE, i,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
1217 + VAT3(oPSE, ii,jj,kk) * VAT3( uC, ip1,jm1,k) * VAT3(dPNW, iip1,jjm1,kkp1)
1218
1219 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oN, ip1,jm2,kp1) * VAT3( dPW, iip1,jjm1,kkp1)
1220 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
1221 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNW, iip1,jjm1,kkp1)
1222 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPN, iip1,jjm1,kkp1))
1223
1224 + VAT3( uPE, ii,jj,kk) * VAT3( oN, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1);
1225
1226 //fprintf(data, "%19.12E\n", VAT3(XuSE, ii, jj, kk));
1227
1228 /* *************************************************************
1229 * *** > uSW;
1230 * *************************************************************/
1231
1232 VAT3(XuSW, ii,jj,kk) =
1233 VAT3(oPSW, ii,jj,kk) * VAT3( uC, im1,jm1,k) * VAT3(dPNE, iim1,jjm1,kkp1)
1234
1235 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oE, im2,jm1,kp1) * VAT3( dPN, iim1,jjm1,kkp1)
1236 - VAT3( oN, im1,jm2,kp1) * VAT3( dPE, iim1,jjm1,kkp1)
1237 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
1238 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
1239
1240 + VAT3( uPW, ii,jj,kk) * VAT3( oN, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
1241 + VAT3( uPS, ii,jj,kk) * VAT3( oE, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1);
1242
1243 //fprintf(data, "%19.12E\n", VAT3(XuSW, ii, jj, kk));
1244
1245 }
1246 }
1247 }
1248}
1249
1250
1251
1252VPUBLIC void VbuildG_27(int *nxf, int *nyf, int *nzf,
1253 int *nx, int *ny, int *nz,
1254 double *oPC, double *oPN, double *oPS, double *oPE, double *oPW,
1255 double *oPNE, double *oPNW, double *oPSE, double *oPSW,
1256 double *uPC, double *uPN, double *uPS, double *uPE, double *uPW,
1257 double *uPNE, double *uPNW, double *uPSE, double *uPSW,
1258 double *dPC, double *dPN, double *dPS, double *dPE, double *dPW,
1259 double *dPNE, double *dPNW, double *dPSE, double *dPSW,
1260 double *oC, double *oE, double *oN, double *uC,
1261 double *oNE, double *oNW, double *uE, double *uW, double *uN,
1262 double *uS, double *uNE, double *uNW, double *uSE, double *uSW,
1263 double *XoC, double *XoE, double *XoN,
1264 double *XuC,
1265 double *XoNE, double *XoNW,
1266 double *XuE, double *XuW, double *XuN, double *XuS,
1267 double *XuNE, double *XuNW, double *XuSE, double *XuSW) {
1268
1269 int i, j, k;
1270 int ii, jj, kk;
1271 int im1, ip1, im2, ip2;
1272 int jm1, jp1, jm2, jp2;
1273 int km1, kp1, km2, kp2;
1274 int iim1, iip1;
1275 int jjm1, jjp1;
1276 int kkm1, kkp1;
1277 int nxm1, nym1, nzm1;
1278
1279 double TMP1_XOC, TMP2_XOC, TMP3_XOC, TMP4_XOC;
1280 double TMP5_XOC, TMP6_XOC, TMP7_XOC, TMP8_XOC;
1281 double TMP9_XOC, TMP10_XOC, TMP11_XOC, TMP12_XOC;
1282 double TMP13_XOC, TMP14_XOC, TMP15_XOC, TMP16_XOC;
1283 double TMP17_XOC, TMP18_XOC, TMP19_XOC, TMP20_XOC;
1284 double TMP21_XOC, TMP22_XOC, TMP23_XOC, TMP24_XOC;
1285 double TMP25_XOC, TMP26_XOC, TMP27_XOC;
1286
1287 double TMP1_XOE, TMP2_XOE, TMP3_XOE, TMP4_XOE;
1288 double TMP5_XOE, TMP6_XOE, TMP7_XOE, TMP8_XOE;
1289 double TMP9_XOE, TMP10_XOE, TMP11_XOE, TMP12_XOE;
1290
1291 double TMP1_XON, TMP2_XON, TMP3_XON, TMP4_XON;
1292 double TMP5_XON, TMP6_XON, TMP7_XON, TMP8_XON;
1293 double TMP9_XON, TMP10_XON, TMP11_XON, TMP12_XON;
1294
1295 double TMP1_XUC, TMP2_XUC, TMP3_XUC, TMP4_XUC;
1296 double TMP5_XUC, TMP6_XUC, TMP7_XUC, TMP8_XUC;
1297 double TMP9_XUC, TMP10_XUC, TMP11_XUC, TMP12_XUC;
1298
1299 double TMP1_XONE, TMP2_XONE, TMP3_XONE, TMP4_XONE;
1300 double TMP5_XONE, TMP6_XONE;
1301
1302 double TMP1_XONW, TMP2_XONW, TMP3_XONW, TMP4_XONW;
1303 double TMP5_XONW, TMP6_XONW;
1304
1305 double TMP1_XUE, TMP2_XUE, TMP3_XUE, TMP4_XUE;
1306 double TMP5_XUE, TMP6_XUE;
1307
1308 double TMP1_XUW, TMP2_XUW, TMP3_XUW, TMP4_XUW;
1309 double TMP5_XUW, TMP6_XUW;
1310
1311 double TMP1_XUN, TMP2_XUN, TMP3_XUN, TMP4_XUN;
1312 double TMP5_XUN, TMP6_XUN;
1313
1314 double TMP1_XUS, TMP2_XUS, TMP3_XUS, TMP4_XUS;
1315 double TMP5_XUS, TMP6_XUS;
1316
1317 double TMP1_XUNE, TMP2_XUNE, TMP1_XUNW, TMP2_XUNW;
1318 double TMP1_XUSE, TMP2_XUSE, TMP1_XUSW, TMP2_XUSW;
1319
1320 MAT3( oC, *nxf, *nyf, *nzf);
1321 MAT3( oE, *nxf, *nyf, *nzf);
1322 MAT3( oN, *nxf, *nyf, *nzf);
1323 MAT3( oNE, *nxf, *nyf, *nzf);
1324 MAT3( oNW, *nxf, *nyf, *nzf);
1325
1326 MAT3( uC, *nxf, *nyf, *nzf);
1327 MAT3( uE, *nxf, *nyf, *nzf);
1328 MAT3( uW, *nxf, *nyf, *nzf);
1329 MAT3( uN, *nxf, *nyf, *nzf);
1330 MAT3( uS, *nxf, *nyf, *nzf);
1331 MAT3( uNE, *nxf, *nyf, *nzf);
1332 MAT3( uNW, *nxf, *nyf, *nzf);
1333 MAT3( uSE, *nxf, *nyf, *nzf);
1334 MAT3( uSW, *nxf, *nyf, *nzf);
1335
1336 MAT3( XoC, *nx, *ny, *nz);
1337 MAT3( XoE, *nx, *ny, *nz);
1338 MAT3( XoN, *nx, *ny, *nz);
1339 MAT3(XoNE, *nx, *ny, *nz);
1340 MAT3(XoNW, *nx, *ny, *nz);
1341
1342 MAT3( XuC, *nx, *ny, *nz);
1343 MAT3( XuE, *nx, *ny, *nz);
1344 MAT3( XuW, *nx, *ny, *nz);
1345 MAT3( XuN, *nx, *ny, *nz);
1346 MAT3( XuS, *nx, *ny, *nz);
1347 MAT3(XuNE, *nx, *ny, *nz);
1348 MAT3(XuNW, *nx, *ny, *nz);
1349 MAT3(XuSE, *nx, *ny, *nz);
1350 MAT3(XuSW, *nx, *ny, *nz);
1351
1352 MAT3( oPC, *nx, *ny, *nz);
1353 MAT3( oPN, *nx, *ny, *nz);
1354 MAT3( oPS, *nx, *ny, *nz);
1355 MAT3( oPE, *nx, *ny, *nz);
1356 MAT3( oPW, *nx, *ny, *nz);
1357 MAT3(oPNE, *nx, *ny, *nz);
1358 MAT3(oPNW, *nx, *ny, *nz);
1359 MAT3(oPSE, *nx, *ny, *nz);
1360 MAT3(oPSW, *nx, *ny, *nz);
1361
1362 MAT3( uPC, *nx, *ny, *nz);
1363 MAT3( uPN, *nx, *ny, *nz);
1364 MAT3( uPS, *nx, *ny, *nz);
1365 MAT3( uPE, *nx, *ny, *nz);
1366 MAT3( uPW, *nx, *ny, *nz);
1367 MAT3(uPNE, *nx, *ny, *nz);
1368 MAT3(uPNW, *nx, *ny, *nz);
1369 MAT3(uPSE, *nx, *ny, *nz);
1370 MAT3(uPSW, *nx, *ny, *nz);
1371
1372 MAT3( dPC, *nx, *ny, *nz);
1373 MAT3( dPN, *nx, *ny, *nz);
1374 MAT3( dPS, *nx, *ny, *nz);
1375 MAT3( dPE, *nx, *ny, *nz);
1376 MAT3( dPW, *nx, *ny, *nz);
1377 MAT3(dPNE, *nx, *ny, *nz);
1378 MAT3(dPNW, *nx, *ny, *nz);
1379 MAT3(dPSE, *nx, *ny, *nz);
1380 MAT3(dPSW, *nx, *ny, *nz);
1381
1382 // define n and determine number of mesh points ***
1383 nxm1 = *nx - 1;
1384 nym1 = *ny - 1;
1385 nzm1 = *nz - 1;
1386
1387 //fprintf(data, "%s\n", PRINT_FUNC);
1388
1389 // Build the operator ***
1390 for(kk=2; kk<=*nz-1; kk++) {
1391 k = 2 * kk - 1;
1392
1393 for(jj=2; jj<=*ny-1; jj++) {
1394 j = 2 * jj - 1;
1395
1396 for(ii=2; ii<=*nx-1; ii++) {
1397 i = 2 * ii - 1;
1398
1399 // Index computations
1400 im1 = i - 1;
1401 ip1 = i + 1;
1402 im2 = i - 2;
1403 ip2 = i + 2;
1404 jm1 = j - 1;
1405 jp1 = j + 1;
1406 jm2 = j - 2;
1407 jp2 = j + 2;
1408 km1 = k - 1;
1409 kp1 = k + 1;
1410 km2 = k - 2;
1411 kp2 = k + 2;
1412 iim1 = ii - 1;
1413 iip1 = ii + 1;
1414 jjm1 = jj - 1;
1415 jjp1 = jj + 1;
1416 kkm1 = kk - 1;
1417 kkp1 = kk + 1;
1418
1419
1420
1421 /* *************************************************************
1422 * *** > oC;
1423 * *************************************************************/
1424
1425 // VAT3( XoC, ii,jj,kk) =
1426 TMP1_XOC =
1427 VAT3( oPN, ii,jj,kk) * (- VAT3( uNE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1428 - VAT3( oNE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1429 - VAT3( uSW, i,jp1,k) * VAT3( uPW, ii,jj,kk)
1430 - VAT3( uE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1431 - VAT3( oE, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1432 - VAT3( uW, i,jp1,k) * VAT3(uPNW, ii,jj,kk)
1433 - VAT3( uN, i,j,km1) * VAT3( dPC, ii,jj,kk)
1434 - VAT3( oN, i,j,k) * VAT3( oPC, ii,jj,kk)
1435 - VAT3( uS, i,jp1,k) * VAT3( uPC, ii,jj,kk)
1436 - VAT3( uC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1437 + VAT3( oC, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1438 - VAT3( uC, i,jp1,k) * VAT3( uPN, ii,jj,kk)
1439 - VAT3( uNW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1440 - VAT3( oNW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1441 - VAT3( uSE, i,jp1,k) * VAT3( uPE, ii,jj,kk)
1442 - VAT3( uW, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1443 - VAT3( oE, i,jp1,k) * VAT3(oPNE, ii,jj,kk)
1444 - VAT3( uE, i,jp1,k) * VAT3(uPNE, ii,jj,kk));
1445
1446 TMP2_XOC =
1447 + VAT3( dPN, ii,jj,kk) * (- VAT3( oNE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1448 - VAT3( uSW, i,jp1,km1) * VAT3( oPW, ii,jj,kk)
1449 - VAT3( oE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1450 - VAT3( uW, i,jp1,km1) * VAT3(oPNW, ii,jj,kk)
1451 - VAT3( oN, i,j,km1) * VAT3( dPC, ii,jj,kk)
1452 - VAT3( uS, i,jp1,km1) * VAT3( oPC, ii,jj,kk)
1453 + VAT3( oC, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1454 - VAT3( uC, i,jp1,km1) * VAT3( oPN, ii,jj,kk)
1455 - VAT3( oNW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1456 - VAT3( uSE, i,jp1,km1) * VAT3( oPE, ii,jj,kk)
1457 - VAT3( oE, i,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1458 - VAT3( uE, i,jp1,km1) * VAT3(oPNE, ii,jj,kk));
1459
1460 TMP3_XOC =
1461 + VAT3( dPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1462 - VAT3( uSW, i,j,km1) * VAT3(oPSW, ii,jj,kk)
1463 - VAT3( oE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1464 - VAT3( uW, i,j,km1) * VAT3( oPW, ii,jj,kk)
1465 - VAT3( oNW, i,j,km1) * VAT3(dPNW, ii,jj,kk)
1466 - VAT3( uNW, i,j,km1) * VAT3(oPNW, ii,jj,kk)
1467 - VAT3( oN, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1468 - VAT3( uS, i,j,km1) * VAT3( oPS, ii,jj,kk)
1469 + VAT3( oC, i,j,km1) * VAT3( dPC, ii,jj,kk)
1470 - VAT3( uC, i,j,km1) * VAT3( oPC, ii,jj,kk)
1471 - VAT3( oN, i,j,km1) * VAT3( dPN, ii,jj,kk)
1472 - VAT3( uN, i,j,km1) * VAT3( oPN, ii,jj,kk)
1473 - VAT3( oNW, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1474 - VAT3( uSE, i,j,km1) * VAT3(oPSE, ii,jj,kk)
1475 - VAT3( oE, i,j,km1) * VAT3( dPE, ii,jj,kk)
1476 - VAT3( uE, i,j,km1) * VAT3( oPE, ii,jj,kk)
1477 - VAT3( oNE, i,j,km1) * VAT3(dPNE, ii,jj,kk)
1478 - VAT3( uNE, i,j,km1) * VAT3(oPNE, ii,jj,kk));
1479
1480 TMP4_XOC =
1481 + VAT3( uPC, ii,jj,kk) * (- VAT3( uNE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1482 - VAT3( oNE, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1483 - VAT3( uE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1484 - VAT3( oE, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1485 - VAT3( uSE, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1486 - VAT3( oNW, i,j,kp1) * VAT3(uPNW, ii,jj,kk)
1487 - VAT3( uN, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1488 - VAT3( oN, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1489 - VAT3( uC, i,j,k) * VAT3( oPC, ii,jj,kk)
1490 + VAT3( oC, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1491 - VAT3( uS, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1492 - VAT3( oN, i,j,kp1) * VAT3( uPN, ii,jj,kk)
1493 - VAT3( uNW, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1494 - VAT3( oNW, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1495 - VAT3( uW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1496 - VAT3( oE, i,j,kp1) * VAT3( uPE, ii,jj,kk)
1497 - VAT3( uSW, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1498 - VAT3( oNE, i,j,kp1) * VAT3(uPNE, ii,jj,kk));
1499
1500 TMP5_XOC =
1501 + VAT3( oPC, ii,jj,kk) * (- VAT3( uW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1502 - VAT3( oE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1503 - VAT3( uSE, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1504 - VAT3( uNE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1505 - VAT3( uN, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1506 - VAT3( oNE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1507 - VAT3( uE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1508 - VAT3( oNW, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1509 - VAT3( uC, i,j,km1) * VAT3( dPC, ii,jj,kk)
1510 - VAT3( uNW, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1511 - VAT3( uSW, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1512 - VAT3( uS, i,jp1,km1) * VAT3( dPN, ii,jj,kk) - VAT3( oN, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1513 - VAT3( uNE, i,j,k) * VAT3(uPNE, ii,jj,kk) - VAT3( oNE, i,j,k) * VAT3(oPNE, ii,jj,kk)
1514 - VAT3( uE, i,j,k) * VAT3( uPE, ii,jj,kk) - VAT3( uSE, i,j,k) * VAT3(uPSE, ii,jj,kk)
1515 - VAT3( oN, i,j,k) * VAT3( oPN, ii,jj,kk) - VAT3( oE, i,j,k) * VAT3( oPE, ii,jj,kk)
1516 - VAT3( uS, i,j,k) * VAT3( uPS, ii,jj,kk) + VAT3( oC, i,j,k) * VAT3( oPC, ii,jj,kk)
1517 - VAT3( uSW, i,j,k) * VAT3(uPSW, ii,jj,kk) - VAT3( uN, i,j,k) * VAT3( uPN, ii,jj,kk)
1518 - VAT3( uC, i,j,k) * VAT3( uPC, ii,jj,kk) - VAT3( uW, i,j,k) * VAT3( uPW, ii,jj,kk)
1519 - VAT3( oNW, i,j,k) * VAT3(oPNW, ii,jj,kk) - VAT3( uNW, i,j,k) * VAT3(uPNW, ii,jj,kk));
1520
1521 TMP6_XOC =
1522 + VAT3( uPS, ii,jj,kk) * (- VAT3( uE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1523 - VAT3( oE, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1524 - VAT3( uSE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1525 - VAT3( oNW, i,jm1,kp1) * VAT3( uPW, ii,jj,kk)
1526 - VAT3( uC, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1527 + VAT3( oC, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1528 - VAT3( uS, i,j,k) * VAT3( oPC, ii,jj,kk)
1529 - VAT3( oN, i,jm1,kp1) * VAT3( uPC, ii,jj,kk)
1530 - VAT3( uW, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1531 - VAT3( oE, i,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1532 - VAT3( uSW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1533 - VAT3( oNE, i,jm1,kp1) * VAT3( uPE, ii,jj,kk));
1534
1535 TMP7_XOC =
1536 + VAT3( oPS, ii,jj,kk) * (- VAT3( uE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1537 - VAT3( oE, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1538 - VAT3( uW, i,jm1,k) * VAT3(uPSW, ii,jj,kk)
1539 - VAT3( uSE, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1540 - VAT3( oNW, i,jm1,k) * VAT3( oPW, ii,jj,kk)
1541 - VAT3( uNW, i,jm1,k) * VAT3( uPW, ii,jj,kk)
1542 - VAT3( uC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1543 + VAT3( oC, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1544 - VAT3( uC, i,jm1,k) * VAT3( uPS, ii,jj,kk)
1545 - VAT3( uS, i,j,km1) * VAT3( dPC, ii,jj,kk)
1546 - VAT3( oN, i,jm1,k) * VAT3( oPC, ii,jj,kk)
1547 - VAT3( uN, i,jm1,k) * VAT3( uPC, ii,jj,kk)
1548 - VAT3( uW, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1549 - VAT3( oE, i,jm1,k) * VAT3(oPSE, ii,jj,kk)
1550 - VAT3( uE, i,jm1,k) * VAT3(uPSE, ii,jj,kk)
1551 - VAT3( uSW, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1552 - VAT3( oNE, i,jm1,k) * VAT3( oPE, ii,jj,kk)
1553 - VAT3( uNE, i,jm1,k) * VAT3( uPE, ii,jj,kk));
1554
1555 TMP8_XOC =
1556 + VAT3( dPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1557 - VAT3( uW, i,jm1,km1) * VAT3(oPSW, ii,jj,kk)
1558 - VAT3( oNW, i,jm1,km1) * VAT3( dPW, ii,jj,kk)
1559 - VAT3( uNW, i,jm1,km1) * VAT3( oPW, ii,jj,kk)
1560 + VAT3( oC, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1561 - VAT3( uC, i,jm1,km1) * VAT3( oPS, ii,jj,kk)
1562 - VAT3( oN, i,jm1,km1) * VAT3( dPC, ii,jj,kk)
1563 - VAT3( uN, i,jm1,km1) * VAT3( oPC, ii,jj,kk)
1564 - VAT3( oE, i,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1565 - VAT3( uE, i,jm1,km1) * VAT3(oPSE, ii,jj,kk)
1566 - VAT3( oNE, i,jm1,km1) * VAT3( dPE, ii,jj,kk)
1567 - VAT3( uNE, i,jm1,km1) * VAT3( oPE, ii,jj,kk));
1568
1569 TMP9_XOC =
1570 + VAT3(uPNW, ii,jj,kk) * (- VAT3( uN, im1,j,k) * VAT3( oPW, ii,jj,kk)
1571 - VAT3( oN, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1572 - VAT3( uC, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1573 + VAT3( oC, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
1574 - VAT3( uNW, i,j,k) * VAT3( oPC, ii,jj,kk)
1575 - VAT3( oNW, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1576 - VAT3( uW, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1577 - VAT3( oE, im1,jp1,kp1) * VAT3( uPN, ii,jj,kk));
1578
1579 TMP10_XOC =
1580 + VAT3(oPNW, ii,jj,kk) * (- VAT3( uN, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1581 - VAT3( oN, im1,j,k) * VAT3( oPW, ii,jj,kk)
1582 - VAT3( uS, im1,jp1,k) * VAT3( uPW, ii,jj,kk)
1583 - VAT3( uC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1584 + VAT3( oC, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1585 - VAT3( uC, im1,jp1,k) * VAT3(uPNW, ii,jj,kk)
1586 - VAT3( uNW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1587 - VAT3( oNW, i,j,k) * VAT3( oPC, ii,jj,kk)
1588 - VAT3( uSE, im1,jp1,k) * VAT3( uPC, ii,jj,kk)
1589 - VAT3( uW, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1590 - VAT3( oE, im1,jp1,k) * VAT3( oPN, ii,jj,kk)
1591 - VAT3( uE, im1,jp1,k) * VAT3( uPN, ii,jj,kk));
1592
1593 TMP11_XOC =
1594 + VAT3( uPW, ii,jj,kk) * (- VAT3( uN, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1595 - VAT3( oN, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1596 - VAT3( uC, im1,j,k) * VAT3( oPW, ii,jj,kk)
1597 + VAT3( oC, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1598 - VAT3( uS, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1599 - VAT3( oN, im1,j,kp1) * VAT3(uPNW, ii,jj,kk)
1600 - VAT3( uNW, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1601 - VAT3( oNW, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1602 - VAT3( uW, i,j,k) * VAT3( oPC, ii,jj,kk)
1603 - VAT3( oE, im1,j,kp1) * VAT3( uPC, ii,jj,kk)
1604 - VAT3( uSW, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1605 - VAT3( oNE, im1,j,kp1) * VAT3( uPN, ii,jj,kk));
1606
1607 TMP12_XOC =
1608 + VAT3(dPNW, ii,jj,kk) * (- VAT3( oN, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1609 - VAT3( uS, im1,jp1,km1) * VAT3( oPW, ii,jj,kk)
1610 + VAT3( oC, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1611 - VAT3( uC, im1,jp1,km1) * VAT3(oPNW, ii,jj,kk)
1612 - VAT3( oNW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1613 - VAT3( uSE, im1,jp1,km1) * VAT3( oPC, ii,jj,kk)
1614 - VAT3( oE, im1,jp1,km1) * VAT3( dPN, ii,jj,kk)
1615 - VAT3( uE, im1,jp1,km1) * VAT3( oPN, ii,jj,kk));
1616
1617 TMP13_XOC =
1618 + VAT3( oPW, ii,jj,kk) * (- VAT3( uN, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1619 - VAT3( oN, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1620 - VAT3( uS, im1,j,k) * VAT3(uPSW, ii,jj,kk)
1621 - VAT3( uC, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1622 + VAT3( oC, im1,j,k) * VAT3( oPW, ii,jj,kk)
1623 - VAT3( uC, im1,j,k) * VAT3( uPW, ii,jj,kk)
1624 - VAT3( uS, im1,jp1,km1) * VAT3(dPNW, ii,jj,kk)
1625 - VAT3( oN, im1,j,k) * VAT3(oPNW, ii,jj,kk)
1626 - VAT3( uN, im1,j,k) * VAT3(uPNW, ii,jj,kk)
1627 - VAT3( uNW, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1628 - VAT3( oNW, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1629 - VAT3( uSE, im1,j,k) * VAT3( uPS, ii,jj,kk)
1630 - VAT3( uW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1631 - VAT3( oE, im1,j,k) * VAT3( oPC, ii,jj,kk)
1632 - VAT3( uE, im1,j,k) * VAT3( uPC, ii,jj,kk)
1633 - VAT3( uSW, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1634 - VAT3( oNE, im1,j,k) * VAT3( oPN, ii,jj,kk)
1635 - VAT3( uNE, im1,j,k) * VAT3( uPN, ii,jj,kk));
1636
1637 TMP14_XOC =
1638 + VAT3(uPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1639 + VAT3( oC, im1,jm1,kp1) * VAT3(uPSW, ii,jj,kk)
1640 - VAT3( uS, im1,j,k) * VAT3( oPW, ii,jj,kk)
1641 - VAT3( oN, im1,jm1,kp1) * VAT3( uPW, ii,jj,kk)
1642 - VAT3( uW, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1643 - VAT3( oE, im1,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1644 - VAT3( uSW, i,j,k) * VAT3( oPC, ii,jj,kk)
1645 - VAT3( oNE, im1,jm1,kp1) * VAT3( uPC, ii,jj,kk));
1646
1647 TMP15_XOC =
1648 + VAT3(oPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1649 + VAT3( oC, im1,jm1,k) * VAT3(oPSW, ii,jj,kk)
1650 - VAT3( uC, im1,jm1,k) * VAT3(uPSW, ii,jj,kk)
1651 - VAT3( uS, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1652 - VAT3( oN, im1,jm1,k) * VAT3( oPW, ii,jj,kk)
1653 - VAT3( uN, im1,jm1,k) * VAT3( uPW, ii,jj,kk)
1654 - VAT3( uW, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1655 - VAT3( oE, im1,jm1,k) * VAT3( oPS, ii,jj,kk)
1656 - VAT3( uE, im1,jm1,k) * VAT3( uPS, ii,jj,kk)
1657 - VAT3( uSW, i,j,km1) * VAT3( dPC, ii,jj,kk)
1658 - VAT3( oNE, im1,jm1,k) * VAT3( oPC, ii,jj,kk)
1659 - VAT3( uNE, im1,jm1,k) * VAT3( uPC, ii,jj,kk));
1660
1661 TMP16_XOC =
1662 + VAT3( dPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1663 - VAT3( uS, im1,j,km1) * VAT3(oPSW, ii,jj,kk)
1664 + VAT3( oC, im1,j,km1) * VAT3( dPW, ii,jj,kk)
1665 - VAT3( uC, im1,j,km1) * VAT3( oPW, ii,jj,kk)
1666 - VAT3( oN, im1,j,km1) * VAT3(dPNW, ii,jj,kk)
1667 - VAT3( uN, im1,j,km1) * VAT3(oPNW, ii,jj,kk)
1668 - VAT3( oNW, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1669 - VAT3( uSE, im1,j,km1) * VAT3( oPS, ii,jj,kk)
1670 - VAT3( oE, im1,j,km1) * VAT3( dPC, ii,jj,kk)
1671 - VAT3( uE, im1,j,km1) * VAT3( oPC, ii,jj,kk)
1672 - VAT3( oNE, im1,j,km1) * VAT3( dPN, ii,jj,kk)
1673 - VAT3( uNE, im1,j,km1) * VAT3( oPN, ii,jj,kk));
1674
1675 TMP17_XOC =
1676 + VAT3(uPNE, ii,jj,kk) * (- VAT3( uNE, i,j,k) * VAT3( oPC, ii,jj,kk)
1677 - VAT3( oNE, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1678 - VAT3( uE, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1679 - VAT3( oE, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
1680 - VAT3( uN, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1681 - VAT3( oN, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
1682 - VAT3( uC, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1683 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNE, ii,jj,kk));
1684
1685 TMP18_XOC =
1686 + VAT3( uPE, ii,jj,kk) * (- VAT3( uNE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1687 - VAT3( oNE, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1688 - VAT3( uE, i,j,k) * VAT3( oPC, ii,jj,kk)
1689 - VAT3( oE, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1690 - VAT3( uSE, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1691 - VAT3( oNW, ip1,j,kp1) * VAT3( uPN, ii,jj,kk)
1692 - VAT3( uN, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1693 - VAT3( oN, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1694 - VAT3( uC, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1695 + VAT3( oC, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
1696 - VAT3( uS, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1697 - VAT3( oN, ip1,j,kp1) * VAT3(uPNE, ii,jj,kk));
1698
1699 TMP19_XOC =
1700 + VAT3(dPNE, ii,jj,kk) * (- VAT3( oNE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1701 - VAT3( uSW, ip1,jp1,km1) * VAT3( oPC, ii,jj,kk)
1702 - VAT3( oE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1703 - VAT3( uW, ip1,jp1,km1) * VAT3( oPN, ii,jj,kk)
1704 - VAT3( oN, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1705 - VAT3( uS, ip1,jp1,km1) * VAT3( oPE, ii,jj,kk)
1706 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1707 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNE, ii,jj,kk));
1708
1709 TMP20_XOC =
1710 + VAT3(oPNE, ii,jj,kk) * (- VAT3( uNE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1711 - VAT3( oNE, i,j,k) * VAT3( oPC, ii,jj,kk)
1712 - VAT3( uSW, ip1,jp1,k) * VAT3( uPC, ii,jj,kk)
1713 - VAT3( uE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1714 - VAT3( oE, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1715 - VAT3( uW, ip1,jp1,k) * VAT3( uPN, ii,jj,kk)
1716 - VAT3( uN, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1717 - VAT3( oN, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1718 - VAT3( uS, ip1,jp1,k) * VAT3( uPE, ii,jj,kk)
1719 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1720 + VAT3( oC, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1721 - VAT3( uC, ip1,jp1,k) * VAT3(uPNE, ii,jj,kk));
1722
1723 TMP21_XOC =
1724 + VAT3(oPSE, ii,jj,kk) * (- VAT3( uE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1725 - VAT3( oE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1726 - VAT3( uW, ip1,jm1,k) * VAT3( uPS, ii,jj,kk)
1727 - VAT3( uSE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1728 - VAT3( oNW, ip1,jm1,k) * VAT3( oPC, ii,jj,kk)
1729 - VAT3( uNW, ip1,jm1,k) * VAT3( uPC, ii,jj,kk)
1730 - VAT3( uC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1731 + VAT3( oC, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1732 - VAT3( uC, ip1,jm1,k) * VAT3(uPSE, ii,jj,kk)
1733 - VAT3( uS, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1734 - VAT3( oN, ip1,jm1,k) * VAT3( oPE, ii,jj,kk)
1735 - VAT3( uN, ip1,jm1,k) * VAT3( uPE, ii,jj,kk));
1736
1737 TMP22_XOC =
1738 + VAT3(dPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1739 - VAT3( uW, ip1,jm1,km1) * VAT3( oPS, ii,jj,kk)
1740 - VAT3( oNW, ip1,jm1,km1) * VAT3( dPC, ii,jj,kk)
1741 - VAT3( uNW, ip1,jm1,km1) * VAT3( oPC, ii,jj,kk)
1742 + VAT3( oC, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1743 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSE, ii,jj,kk)
1744 - VAT3( oN, ip1,jm1,km1) * VAT3( dPE, ii,jj,kk)
1745 - VAT3( uN, ip1,jm1,km1) * VAT3( oPE, ii,jj,kk));
1746
1747 TMP23_XOC =
1748 + VAT3(uPSE, ii,jj,kk) * (- VAT3( uE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1749 - VAT3( oE, i,jm1,kp1) * VAT3( uPS, ii,jj,kk)
1750 - VAT3( uSE, i,j,k) * VAT3( oPC, ii,jj,kk)
1751 - VAT3( oNW, ip1,jm1,kp1) * VAT3( uPC, ii,jj,kk)
1752 - VAT3( uC, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1753 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSE, ii,jj,kk)
1754 - VAT3( uS, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1755 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPE, ii,jj,kk));
1756
1757 TMP24_XOC =
1758 + VAT3( oPE, ii,jj,kk) * (- VAT3( uNE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1759 - VAT3( oNE, i,jm1,k) * VAT3( oPS, ii,jj,kk)
1760 - VAT3( uSW, ip1,j,k) * VAT3( uPS, ii,jj,kk)
1761 - VAT3( uE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1762 - VAT3( oE, i,j,k) * VAT3( oPC, ii,jj,kk)
1763 - VAT3( uW, ip1,j,k) * VAT3( uPC, ii,jj,kk)
1764 - VAT3( uSE, i,jp1,km1) * VAT3( dPN, ii,jj,kk)
1765 - VAT3( oNW, ip1,j,k) * VAT3( oPN, ii,jj,kk)
1766 - VAT3( uNW, ip1,j,k) * VAT3( uPN, ii,jj,kk)
1767 - VAT3( uN, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1768 - VAT3( oN, ip1,jm1,k) * VAT3(oPSE, ii,jj,kk)
1769 - VAT3( uS, ip1,j,k) * VAT3(uPSE, ii,jj,kk)
1770 - VAT3( uC, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1771 + VAT3( oC, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1772 - VAT3( uC, ip1,j,k) * VAT3( uPE, ii,jj,kk)
1773 - VAT3( uS, ip1,jp1,km1) * VAT3(dPNE, ii,jj,kk)
1774 - VAT3( oN, ip1,j,k) * VAT3(oPNE, ii,jj,kk)
1775 - VAT3( uN, ip1,j,k) * VAT3(uPNE, ii,jj,kk));
1776
1777 TMP25_XOC =
1778 + VAT3( dPE, ii,jj,kk) * (- VAT3( oNE, i,jm1,km1) * VAT3( dPS, ii,jj,kk)
1779 - VAT3( uSW, ip1,j,km1) * VAT3( oPS, ii,jj,kk)
1780 - VAT3( oE, i,j,km1) * VAT3( dPC, ii,jj,kk)
1781 - VAT3( uW, ip1,j,km1) * VAT3( oPC, ii,jj,kk)
1782 - VAT3( oNW, ip1,j,km1) * VAT3( dPN, ii,jj,kk)
1783 - VAT3( uNW, ip1,j,km1) * VAT3( oPN, ii,jj,kk)
1784 - VAT3( oN, ip1,jm1,km1) * VAT3(dPSE, ii,jj,kk)
1785 - VAT3( uS, ip1,j,km1) * VAT3(oPSE, ii,jj,kk)
1786 + VAT3( oC, ip1,j,km1) * VAT3( dPE, ii,jj,kk)
1787 - VAT3( uC, ip1,j,km1) * VAT3( oPE, ii,jj,kk)
1788 - VAT3( oN, ip1,j,km1) * VAT3(dPNE, ii,jj,kk)
1789 - VAT3( uN, ip1,j,km1) * VAT3(oPNE, ii,jj,kk));
1790
1791 TMP26_XOC =
1792 + VAT3( uPN, ii,jj,kk) * (- VAT3( uNE, im1,j,k) * VAT3( oPW, ii,jj,kk)
1793 - VAT3( oNE, im1,j,kp1) * VAT3( uPW, ii,jj,kk)
1794 - VAT3( uE, im1,jp1,k) * VAT3(oPNW, ii,jj,kk)
1795 - VAT3( oE, im1,jp1,kp1) * VAT3(uPNW, ii,jj,kk)
1796 - VAT3( uN, i,j,k) * VAT3( oPC, ii,jj,kk)
1797 - VAT3( oN, i,j,kp1) * VAT3( uPC, ii,jj,kk)
1798 - VAT3( uC, i,jp1,k) * VAT3( oPN, ii,jj,kk)
1799 + VAT3( oC, i,jp1,kp1) * VAT3( uPN, ii,jj,kk)
1800 - VAT3( uNW, ip1,j,k) * VAT3( oPE, ii,jj,kk)
1801 - VAT3( oNW, ip1,j,kp1) * VAT3( uPE, ii,jj,kk)
1802 - VAT3( uW, ip1,jp1,k) * VAT3(oPNE, ii,jj,kk)
1803 - VAT3( oE, i,jp1,kp1) * VAT3(uPNE, ii,jj,kk));
1804
1805 TMP27_XOC =
1806 + VAT3(dPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,km1) * VAT3(dPSW, ii,jj,kk)
1807 - VAT3( uC, im1,jm1,km1) * VAT3(oPSW, ii,jj,kk)
1808 - VAT3( oN, im1,jm1,km1) * VAT3( dPW, ii,jj,kk)
1809 - VAT3( uN, im1,jm1,km1) * VAT3( oPW, ii,jj,kk)
1810 - VAT3( oE, im1,jm1,km1) * VAT3( dPS, ii,jj,kk)
1811 - VAT3( uE, im1,jm1,km1) * VAT3( oPS, ii,jj,kk)
1812 - VAT3( oNE, im1,jm1,km1) * VAT3( dPC, ii,jj,kk)
1813 - VAT3( uNE, im1,jm1,km1) * VAT3( oPC, ii,jj,kk));
1814
1815 VAT3( XoC, ii,jj,kk) = TMP1_XOC + TMP2_XOC + TMP3_XOC + TMP4_XOC
1816 + TMP5_XOC + TMP6_XOC + TMP7_XOC + TMP8_XOC + TMP9_XOC
1817 + TMP10_XOC + TMP11_XOC + TMP12_XOC + TMP13_XOC + TMP14_XOC
1818 + TMP15_XOC + TMP16_XOC + TMP17_XOC + TMP18_XOC + TMP19_XOC
1819 + TMP20_XOC + TMP21_XOC + TMP22_XOC + TMP23_XOC + TMP24_XOC
1820 + TMP25_XOC + TMP26_XOC + TMP27_XOC;
1821
1822 //fprintf(data, "%19.12E\n", VAT3(XoC, ii, jj, kk));
1823
1824 /* *************************************************************
1825 * *** > OE;
1826 * *************************************************************/
1827
1828 // VAT3( XoE, ii,jj,kk) =
1829 TMP1_XOE =
1830 - VAT3( dPS, ii,jj,kk) * (- VAT3( oE, i,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1831 - VAT3( uE, i,jm1,km1) * VAT3(oPSW, iip1,jj,kk)
1832 - VAT3( oNE, i,jm1,km1) * VAT3( dPW, iip1,jj,kk)
1833 - VAT3( uNE, i,jm1,km1) * VAT3( oPW, iip1,jj,kk))
1834
1835 - VAT3( oPS, ii,jj,kk) * (- VAT3( uW, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1836 - VAT3( oE, i,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1837 - VAT3( uE, i,jm1,k) * VAT3(uPSW, iip1,jj,kk)
1838 - VAT3( uSW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1839 - VAT3( oNE, i,jm1,k) * VAT3( oPW, iip1,jj,kk)
1840 - VAT3( uNE, i,jm1,k) * VAT3( uPW, iip1,jj,kk))
1841
1842 - VAT3( uPS, ii,jj,kk) * (- VAT3( uW, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1843 - VAT3( oE, i,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1844 - VAT3( uSW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1845 - VAT3( oNE, i,jm1,kp1) * VAT3( uPW, iip1,jj,kk));
1846
1847 TMP2_XOE =
1848 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNW, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1849 - VAT3( uSE, i,j,km1) * VAT3(oPSW, iip1,jj,kk)
1850 - VAT3( oE, i,j,km1) * VAT3( dPW, iip1,jj,kk)
1851 - VAT3( uE, i,j,km1) * VAT3( oPW, iip1,jj,kk)
1852 - VAT3( oNE, i,j,km1) * VAT3(dPNW, iip1,jj,kk)
1853 - VAT3( uNE, i,j,km1) * VAT3(oPNW, iip1,jj,kk))
1854
1855 - VAT3( oPC, ii,jj,kk) * (- VAT3( uNW, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1856 - VAT3( oNW, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1857 - VAT3( uSE, i,j,k) * VAT3(uPSW, iip1,jj,kk)
1858 - VAT3( uW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1859 - VAT3( oE, i,j,k) * VAT3( oPW, iip1,jj,kk)
1860 - VAT3( uE, i,j,k) * VAT3( uPW, iip1,jj,kk)
1861 - VAT3( uSW, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1862 - VAT3( oNE, i,j,k) * VAT3(oPNW, iip1,jj,kk)
1863 - VAT3( uNE, i,j,k) * VAT3(uPNW, iip1,jj,kk));
1864
1865 TMP3_XOE =
1866 - VAT3( uPC, ii,jj,kk) * (- VAT3( uNW, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1867 - VAT3( oNW, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1868 - VAT3( uW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1869 - VAT3( oE, i,j,kp1) * VAT3( uPW, iip1,jj,kk)
1870 - VAT3( uSW, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1871 - VAT3( oNE, i,j,kp1) * VAT3(uPNW, iip1,jj,kk))
1872
1873 - VAT3( dPN, ii,jj,kk) * (- VAT3( oNW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1874 - VAT3( uSE, i,jp1,km1) * VAT3( oPW, iip1,jj,kk)
1875 - VAT3( oE, i,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1876 - VAT3( uE, i,jp1,km1) * VAT3(oPNW, iip1,jj,kk))
1877
1878 - VAT3( oPN, ii,jj,kk) * (- VAT3( uNW, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1879 - VAT3( oNW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1880 - VAT3( uSE, i,jp1,k) * VAT3( uPW, iip1,jj,kk)
1881 - VAT3( uW, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1882 - VAT3( oE, i,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1883 - VAT3( uE, i,jp1,k) * VAT3(uPNW, iip1,jj,kk));
1884
1885 TMP4_XOE =
1886 - VAT3( uPN, ii,jj,kk) * (- VAT3( uNW, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1887 - VAT3( oNW, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
1888 - VAT3( uW, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1889 - VAT3( oE, i,jp1,kp1) * VAT3(uPNW, iip1,jj,kk))
1890
1891 - VAT3(dPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1892 - VAT3( uC, ip1,jm1,km1) * VAT3(oPSW, iip1,jj,kk)
1893 - VAT3( oN, ip1,jm1,km1) * VAT3( dPW, iip1,jj,kk)
1894 - VAT3( uN, ip1,jm1,km1) * VAT3( oPW, iip1,jj,kk)
1895 - VAT3( oE, ip1,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1896 - VAT3( uE, ip1,jm1,km1) * VAT3( oPS, iip1,jj,kk)
1897 - VAT3( oNE, ip1,jm1,km1) * VAT3( dPC, iip1,jj,kk)
1898 - VAT3( uNE, ip1,jm1,km1) * VAT3( oPC, iip1,jj,kk));
1899
1900 TMP5_XOE =
1901 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1902 + VAT3( oC, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1903 - VAT3( uC, ip1,jm1,k) * VAT3(uPSW, iip1,jj,kk)
1904 - VAT3( uS, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1905 - VAT3( oN, ip1,jm1,k) * VAT3( oPW, iip1,jj,kk)
1906 - VAT3( uN, ip1,jm1,k) * VAT3( uPW, iip1,jj,kk)
1907 - VAT3( uW, ip2,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1908 - VAT3( oE, ip1,jm1,k) * VAT3( oPS, iip1,jj,kk)
1909 - VAT3( uE, ip1,jm1,k) * VAT3( uPS, iip1,jj,kk)
1910 - VAT3( uSW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1911 - VAT3( oNE, ip1,jm1,k) * VAT3( oPC, iip1,jj,kk)
1912 - VAT3( uNE, ip1,jm1,k) * VAT3( uPC, iip1,jj,kk));
1913
1914 TMP6_XOE =
1915 - VAT3(uPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1916 + VAT3( oC, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1917 - VAT3( uS, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1918 - VAT3( oN, ip1,jm1,kp1) * VAT3( uPW, iip1,jj,kk)
1919 - VAT3( uW, ip2,jm1,k) * VAT3( oPS, iip1,jj,kk)
1920 - VAT3( oE, ip1,jm1,kp1) * VAT3( uPS, iip1,jj,kk)
1921 - VAT3( uSW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
1922 - VAT3( oNE, ip1,jm1,kp1) * VAT3( uPC, iip1,jj,kk));
1923
1924 TMP7_XOE =
1925 - VAT3( dPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1926 - VAT3( uS, ip1,j,km1) * VAT3(oPSW, iip1,jj,kk)
1927 + VAT3( oC, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1928 - VAT3( uC, ip1,j,km1) * VAT3( oPW, iip1,jj,kk)
1929 - VAT3( oN, ip1,j,km1) * VAT3(dPNW, iip1,jj,kk)
1930 - VAT3( uN, ip1,j,km1) * VAT3(oPNW, iip1,jj,kk)
1931 - VAT3( oNW, ip2,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1932 - VAT3( uSE, ip1,j,km1) * VAT3( oPS, iip1,jj,kk)
1933 - VAT3( oE, ip1,j,km1) * VAT3( dPC, iip1,jj,kk)
1934 - VAT3( uE, ip1,j,km1) * VAT3( oPC, iip1,jj,kk)
1935 - VAT3( oNE, ip1,j,km1) * VAT3( dPN, iip1,jj,kk)
1936 - VAT3( uNE, ip1,j,km1) * VAT3( oPN, iip1,jj,kk));
1937
1938 TMP8_XOE =
1939 - VAT3( oPE, ii,jj,kk) * (- VAT3( uN, ip1,jm1,km1) * VAT3(dPSW, iip1,jj,kk)
1940 - VAT3( oN, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1941 - VAT3( uS, ip1,j,k) * VAT3(uPSW, iip1,jj,kk)
1942 - VAT3( uC, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1943 + VAT3( oC, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1944 - VAT3( uC, ip1,j,k) * VAT3( uPW, iip1,jj,kk)
1945 - VAT3( uS, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1946 - VAT3( oN, ip1,j,k) * VAT3(oPNW, iip1,jj,kk)
1947 - VAT3( uN, ip1,j,k) * VAT3(uPNW, iip1,jj,kk)
1948 - VAT3( uNW, ip2,jm1,km1) * VAT3( dPS, iip1,jj,kk)
1949 - VAT3( oNW, ip2,jm1,k) * VAT3( oPS, iip1,jj,kk)
1950 - VAT3( uSE, ip1,j,k) * VAT3( uPS, iip1,jj,kk)
1951 - VAT3( uW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1952 - VAT3( oE, ip1,j,k) * VAT3( oPC, iip1,jj,kk)
1953 - VAT3( uE, ip1,j,k) * VAT3( uPC, iip1,jj,kk)
1954 - VAT3( uSW, ip2,jp1,km1) * VAT3( dPN, iip1,jj,kk)
1955 - VAT3( oNE, ip1,j,k) * VAT3( oPN, iip1,jj,kk)
1956 - VAT3( uNE, ip1,j,k) * VAT3( uPN, iip1,jj,kk));
1957
1958 TMP9_XOE =
1959 - VAT3( uPE, ii,jj,kk) * (- VAT3( uN, ip1,jm1,k) * VAT3(oPSW, iip1,jj,kk)
1960 - VAT3( oN, ip1,jm1,kp1) * VAT3(uPSW, iip1,jj,kk)
1961 - VAT3( uC, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1962 + VAT3( oC, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
1963 - VAT3( uS, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1964 - VAT3( oN, ip1,j,kp1) * VAT3(uPNW, iip1,jj,kk)
1965 - VAT3( uNW, ip2,jm1,k) * VAT3( oPS, iip1,jj,kk)
1966 - VAT3( oNW, ip2,jm1,kp1) * VAT3( uPS, iip1,jj,kk)
1967 - VAT3( uW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
1968 - VAT3( oE, ip1,j,kp1) * VAT3( uPC, iip1,jj,kk)
1969 - VAT3( uSW, ip2,jp1,k) * VAT3( oPN, iip1,jj,kk)
1970 - VAT3( oNE, ip1,j,kp1) * VAT3( uPN, iip1,jj,kk));
1971
1972 TMP10_XOE =
1973 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oN, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1974 - VAT3( uS, ip1,jp1,km1) * VAT3( oPW, iip1,jj,kk)
1975 + VAT3( oC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1976 - VAT3( uC, ip1,jp1,km1) * VAT3(oPNW, iip1,jj,kk)
1977 - VAT3( oNW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1978 - VAT3( uSE, ip1,jp1,km1) * VAT3( oPC, iip1,jj,kk)
1979 - VAT3( oE, ip1,jp1,km1) * VAT3( dPN, iip1,jj,kk)
1980 - VAT3( uE, ip1,jp1,km1) * VAT3( oPN, iip1,jj,kk));
1981
1982 TMP11_XOE =
1983 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uN, ip1,j,km1) * VAT3( dPW, iip1,jj,kk)
1984 - VAT3( oN, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1985 - VAT3( uS, ip1,jp1,k) * VAT3( uPW, iip1,jj,kk)
1986 - VAT3( uC, ip1,jp1,km1) * VAT3(dPNW, iip1,jj,kk)
1987 + VAT3( oC, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
1988 - VAT3( uC, ip1,jp1,k) * VAT3(uPNW, iip1,jj,kk)
1989 - VAT3( uNW, ip2,j,km1) * VAT3( dPC, iip1,jj,kk)
1990 - VAT3( oNW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
1991 - VAT3( uSE, ip1,jp1,k) * VAT3( uPC, iip1,jj,kk)
1992 - VAT3( uW, ip2,jp1,km1) * VAT3( dPN, iip1,jj,kk)
1993 - VAT3( oE, ip1,jp1,k) * VAT3( oPN, iip1,jj,kk)
1994 - VAT3( uE, ip1,jp1,k) * VAT3( uPN, iip1,jj,kk));
1995
1996 TMP12_XOE =
1997 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uN, ip1,j,k) * VAT3( oPW, iip1,jj,kk)
1998 - VAT3( oN, ip1,j,kp1) * VAT3( uPW, iip1,jj,kk)
1999 - VAT3( uC, ip1,jp1,k) * VAT3(oPNW, iip1,jj,kk)
2000 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPNW, iip1,jj,kk)
2001 - VAT3( uNW, ip2,j,k) * VAT3( oPC, iip1,jj,kk)
2002 - VAT3( oNW, ip2,j,kp1) * VAT3( uPC, iip1,jj,kk)
2003 - VAT3( uW, ip2,jp1,k) * VAT3( oPN, iip1,jj,kk)
2004 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPN, iip1,jj,kk));
2005
2006 VAT3( XoE, ii,jj,kk) = TMP1_XOE + TMP2_XOE + TMP3_XOE + TMP4_XOE
2007 + TMP5_XOE + TMP6_XOE + TMP7_XOE + TMP8_XOE + TMP9_XOE
2008 + TMP10_XOE + TMP11_XOE + TMP12_XOE;
2009
2010 //fprintf(data, "%19.12E\n", VAT3(XoE, ii, jj, kk));
2011
2012 /* *************************************************************
2013 * *** > ON;
2014 * *************************************************************/
2015
2016 // VAT3( XoN, ii,jj,kk) =
2017 TMP1_XON =
2018 - VAT3( dPW, ii,jj,kk) * (- VAT3( oN, im1,j,km1) * VAT3(dPSW, ii,jjp1,kk)
2019 - VAT3( uN, im1,j,km1) * VAT3(oPSW, ii,jjp1,kk)
2020 - VAT3( oNE, im1,j,km1) * VAT3( dPS, ii,jjp1,kk)
2021 - VAT3( uNE, im1,j,km1) * VAT3( oPS, ii,jjp1,kk))
2022
2023 - VAT3( oPW, ii,jj,kk) * (- VAT3( uS, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2024 - VAT3( oN, im1,j,k) * VAT3(oPSW, ii,jjp1,kk)
2025 - VAT3( uN, im1,j,k) * VAT3(uPSW, ii,jjp1,kk)
2026 - VAT3( uSW, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2027 - VAT3( oNE, im1,j,k) * VAT3( oPS, ii,jjp1,kk)
2028 - VAT3( uNE, im1,j,k) * VAT3( uPS, ii,jjp1,kk))
2029
2030 - VAT3( uPW, ii,jj,kk) * (- VAT3( uS, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2031 - VAT3( oN, im1,j,kp1) * VAT3(uPSW, ii,jjp1,kk)
2032 - VAT3( uSW, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2033 - VAT3( oNE, im1,j,kp1) * VAT3( uPS, ii,jjp1,kk));
2034
2035 TMP2_XON =
2036 - VAT3(dPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2037 - VAT3( uC, im1,jp1,km1) * VAT3(oPSW, ii,jjp1,kk)
2038 - VAT3( oN, im1,jp1,km1) * VAT3( dPW, ii,jjp1,kk)
2039 - VAT3( uN, im1,jp1,km1) * VAT3( oPW, ii,jjp1,kk)
2040 - VAT3( oE, im1,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2041 - VAT3( uE, im1,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
2042 - VAT3( oNE, im1,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
2043 - VAT3( uNE, im1,jp1,km1) * VAT3( oPC, ii,jjp1,kk));
2044
2045 TMP3_XON =
2046 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2047 + VAT3( oC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2048 - VAT3( uC, im1,jp1,k) * VAT3(uPSW, ii,jjp1,kk)
2049 - VAT3( uS, im1,jp2,km1) * VAT3( dPW, ii,jjp1,kk)
2050 - VAT3( oN, im1,jp1,k) * VAT3( oPW, ii,jjp1,kk)
2051 - VAT3( uN, im1,jp1,k) * VAT3( uPW, ii,jjp1,kk)
2052 - VAT3( uW, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2053 - VAT3( oE, im1,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2054 - VAT3( uE, im1,jp1,k) * VAT3( uPS, ii,jjp1,kk)
2055 - VAT3( uSW, i,jp2,km1) * VAT3( dPC, ii,jjp1,kk)
2056 - VAT3( oNE, im1,jp1,k) * VAT3( oPC, ii,jjp1,kk)
2057 - VAT3( uNE, im1,jp1,k) * VAT3( uPC, ii,jjp1,kk));
2058
2059 TMP4_XON =
2060 - VAT3(uPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2061 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
2062 - VAT3( uS, im1,jp2,k) * VAT3( oPW, ii,jjp1,kk)
2063 - VAT3( oN, im1,jp1,kp1) * VAT3( uPW, ii,jjp1,kk)
2064 - VAT3( uW, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2065 - VAT3( oE, im1,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
2066 - VAT3( uSW, i,jp2,k) * VAT3( oPC, ii,jjp1,kk)
2067 - VAT3( oNE, im1,jp1,kp1) * VAT3( uPC, ii,jjp1,kk))
2068
2069 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNW, i,j,km1) * VAT3(dPSW, ii,jjp1,kk)
2070 - VAT3( uNW, i,j,km1) * VAT3(oPSW, ii,jjp1,kk)
2071 - VAT3( oN, i,j,km1) * VAT3( dPS, ii,jjp1,kk)
2072 - VAT3( uN, i,j,km1) * VAT3( oPS, ii,jjp1,kk)
2073 - VAT3( oNE, i,j,km1) * VAT3(dPSE, ii,jjp1,kk)
2074 - VAT3( uNE, i,j,km1) * VAT3(oPSE, ii,jjp1,kk));
2075
2076 TMP5_XON =
2077 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2078 - VAT3( oNW, i,j,k) * VAT3(oPSW, ii,jjp1,kk)
2079 - VAT3( uNW, i,j,k) * VAT3(uPSW, ii,jjp1,kk)
2080 - VAT3( uS, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2081 - VAT3( oN, i,j,k) * VAT3( oPS, ii,jjp1,kk)
2082 - VAT3( uN, i,j,k) * VAT3( uPS, ii,jjp1,kk)
2083 - VAT3( uSW, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2084 - VAT3( oNE, i,j,k) * VAT3(oPSE, ii,jjp1,kk)
2085 - VAT3( uNE, i,j,k) * VAT3(uPSE, ii,jjp1,kk))
2086
2087 - VAT3( uPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2088 - VAT3( oNW, i,j,kp1) * VAT3(uPSW, ii,jjp1,kk)
2089 - VAT3( uS, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2090 - VAT3( oN, i,j,kp1) * VAT3( uPS, ii,jjp1,kk)
2091 - VAT3( uSW, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2092 - VAT3( oNE, i,j,kp1) * VAT3(uPSE, ii,jjp1,kk));
2093
2094 TMP6_XON =
2095 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2096 - VAT3( uW, i,jp1,km1) * VAT3(oPSW, ii,jjp1,kk)
2097 - VAT3( oNW, i,jp1,km1) * VAT3( dPW, ii,jjp1,kk)
2098 - VAT3( uNW, i,jp1,km1) * VAT3( oPW, ii,jjp1,kk)
2099 + VAT3( oC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2100 - VAT3( uC, i,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
2101 - VAT3( oN, i,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
2102 - VAT3( uN, i,jp1,km1) * VAT3( oPC, ii,jjp1,kk)
2103 - VAT3( oE, i,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2104 - VAT3( uE, i,jp1,km1) * VAT3(oPSE, ii,jjp1,kk)
2105 - VAT3( oNE, i,jp1,km1) * VAT3( dPE, ii,jjp1,kk)
2106 - VAT3( uNE, i,jp1,km1) * VAT3( oPE, ii,jjp1,kk));
2107
2108 TMP7_XON =
2109 - VAT3( oPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,km1) * VAT3(dPSW, ii,jjp1,kk)
2110 - VAT3( oE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2111 - VAT3( uW, i,jp1,k) * VAT3(uPSW, ii,jjp1,kk)
2112 - VAT3( uSE, im1,jp2,km1) * VAT3( dPW, ii,jjp1,kk)
2113 - VAT3( oNW, i,jp1,k) * VAT3( oPW, ii,jjp1,kk)
2114 - VAT3( uNW, i,jp1,k) * VAT3( uPW, ii,jjp1,kk)
2115 - VAT3( uC, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2116 + VAT3( oC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2117 - VAT3( uC, i,jp1,k) * VAT3( uPS, ii,jjp1,kk)
2118 - VAT3( uS, i,jp2,km1) * VAT3( dPC, ii,jjp1,kk)
2119 - VAT3( oN, i,jp1,k) * VAT3( oPC, ii,jjp1,kk)
2120 - VAT3( uN, i,jp1,k) * VAT3( uPC, ii,jjp1,kk)
2121 - VAT3( uW, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2122 - VAT3( oE, i,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2123 - VAT3( uE, i,jp1,k) * VAT3(uPSE, ii,jjp1,kk)
2124 - VAT3( uSW, ip1,jp2,km1) * VAT3( dPE, ii,jjp1,kk)
2125 - VAT3( oNE, i,jp1,k) * VAT3( oPE, ii,jjp1,kk)
2126 - VAT3( uNE, i,jp1,k) * VAT3( uPE, ii,jjp1,kk));
2127
2128 TMP8_XON =
2129 - VAT3( uPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,k) * VAT3(oPSW, ii,jjp1,kk)
2130 - VAT3( oE, im1,jp1,kp1) * VAT3(uPSW, ii,jjp1,kk)
2131 - VAT3( uSE, im1,jp2,k) * VAT3( oPW, ii,jjp1,kk)
2132 - VAT3( oNW, i,jp1,kp1) * VAT3( uPW, ii,jjp1,kk)
2133 - VAT3( uC, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2134 + VAT3( oC, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
2135 - VAT3( uS, i,jp2,k) * VAT3( oPC, ii,jjp1,kk)
2136 - VAT3( oN, i,jp1,kp1) * VAT3( uPC, ii,jjp1,kk)
2137 - VAT3( uW, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2138 - VAT3( oE, i,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk)
2139 - VAT3( uSW, ip1,jp2,k) * VAT3( oPE, ii,jjp1,kk)
2140 - VAT3( oNE, i,jp1,kp1) * VAT3( uPE, ii,jjp1,kk))
2141
2142 - VAT3( dPE, ii,jj,kk) * (- VAT3( oNW, ip1,j,km1) * VAT3( dPS, ii,jjp1,kk)
2143 - VAT3( uNW, ip1,j,km1) * VAT3( oPS, ii,jjp1,kk)
2144 - VAT3( oN, ip1,j,km1) * VAT3(dPSE, ii,jjp1,kk)
2145 - VAT3( uN, ip1,j,km1) * VAT3(oPSE, ii,jjp1,kk));
2146
2147 TMP9_XON =
2148 - VAT3( oPE, ii,jj,kk) * (- VAT3( uSE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2149 - VAT3( oNW, ip1,j,k) * VAT3( oPS, ii,jjp1,kk)
2150 - VAT3( uNW, ip1,j,k) * VAT3( uPS, ii,jjp1,kk)
2151 - VAT3( uS, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2152 - VAT3( oN, ip1,j,k) * VAT3(oPSE, ii,jjp1,kk)
2153 - VAT3( uN, ip1,j,k) * VAT3(uPSE, ii,jjp1,kk))
2154
2155 - VAT3( uPE, ii,jj,kk) * (- VAT3( uSE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2156 - VAT3( oNW, ip1,j,kp1) * VAT3( uPS, ii,jjp1,kk)
2157 - VAT3( uS, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2158 - VAT3( oN, ip1,j,kp1) * VAT3(uPSE, ii,jjp1,kk));
2159
2160 TMP10_XON =
2161 - VAT3(dPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2162 - VAT3( uW, ip1,jp1,km1) * VAT3( oPS, ii,jjp1,kk)
2163 - VAT3( oNW, ip1,jp1,km1) * VAT3( dPC, ii,jjp1,kk)
2164 - VAT3( uNW, ip1,jp1,km1) * VAT3( oPC, ii,jjp1,kk)
2165 + VAT3( oC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2166 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSE, ii,jjp1,kk)
2167 - VAT3( oN, ip1,jp1,km1) * VAT3( dPE, ii,jjp1,kk)
2168 - VAT3( uN, ip1,jp1,km1) * VAT3( oPE, ii,jjp1,kk));
2169
2170 TMP11_XON =
2171 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uE, i,jp1,km1) * VAT3( dPS, ii,jjp1,kk)
2172 - VAT3( oE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2173 - VAT3( uW, ip1,jp1,k) * VAT3( uPS, ii,jjp1,kk)
2174 - VAT3( uSE, i,jp2,km1) * VAT3( dPC, ii,jjp1,kk)
2175 - VAT3( oNW, ip1,jp1,k) * VAT3( oPC, ii,jjp1,kk)
2176 - VAT3( uNW, ip1,jp1,k) * VAT3( uPC, ii,jjp1,kk)
2177 - VAT3( uC, ip1,jp1,km1) * VAT3(dPSE, ii,jjp1,kk)
2178 + VAT3( oC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2179 - VAT3( uC, ip1,jp1,k) * VAT3(uPSE, ii,jjp1,kk)
2180 - VAT3( uS, ip1,jp2,km1) * VAT3( dPE, ii,jjp1,kk)
2181 - VAT3( oN, ip1,jp1,k) * VAT3( oPE, ii,jjp1,kk)
2182 - VAT3( uN, ip1,jp1,k) * VAT3( uPE, ii,jjp1,kk));
2183
2184 TMP12_XON =
2185 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uE, i,jp1,k) * VAT3( oPS, ii,jjp1,kk)
2186 - VAT3( oE, i,jp1,kp1) * VAT3( uPS, ii,jjp1,kk)
2187 - VAT3( uSE, i,jp2,k) * VAT3( oPC, ii,jjp1,kk)
2188 - VAT3( oNW, ip1,jp1,kp1) * VAT3( uPC, ii,jjp1,kk)
2189 - VAT3( uC, ip1,jp1,k) * VAT3(oPSE, ii,jjp1,kk)
2190 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSE, ii,jjp1,kk)
2191 - VAT3( uS, ip1,jp2,k) * VAT3( oPE, ii,jjp1,kk)
2192 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPE, ii,jjp1,kk));
2193
2194 VAT3( XoN, ii,jj,kk) = TMP1_XON + TMP2_XON + TMP3_XON + TMP4_XON
2195 + TMP5_XON + TMP6_XON + TMP7_XON + TMP8_XON + TMP9_XON
2196 + TMP10_XON + TMP11_XON + TMP12_XON;
2197
2198 //fprintf(data, "%19.12E\n", VAT3(XoN, ii, jj, kk));
2199
2200 /* *************************************************************
2201 * *** > UC;
2202 * *************************************************************/
2203
2204
2205 // VAT3( XuC, ii,jj,kk) =
2206 TMP1_XUC =
2207 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uC, im1,jm1,k) * VAT3(dPSW, ii,jj,kkp1)
2208 - VAT3( uN, im1,jm1,k) * VAT3( dPW, ii,jj,kkp1)
2209 - VAT3( uE, im1,jm1,k) * VAT3( dPS, ii,jj,kkp1)
2210 - VAT3( uNE, im1,jm1,k) * VAT3( dPC, ii,jj,kkp1))
2211
2212 - VAT3(uPSW, ii,jj,kk) * ( VAT3( oC, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2213 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSW, ii,jj,kkp1)
2214 - VAT3( oN, im1,jm1,kp1) * VAT3( dPW, ii,jj,kkp1)
2215 - VAT3( uN, im1,jm1,kp1) * VAT3( oPW, ii,jj,kkp1)
2216 - VAT3( oE, im1,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2217 - VAT3( uE, im1,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
2218 - VAT3( oNE, im1,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
2219 - VAT3( uNE, im1,jm1,kp1) * VAT3( oPC, ii,jj,kkp1));
2220
2221 TMP2_XUC =
2222 - VAT3( oPW, ii,jj,kk) * (- VAT3( uS, im1,j,k) * VAT3(dPSW, ii,jj,kkp1)
2223 - VAT3( uC, im1,j,k) * VAT3( dPW, ii,jj,kkp1)
2224 - VAT3( uN, im1,j,k) * VAT3(dPNW, ii,jj,kkp1)
2225 - VAT3( uSE, im1,j,k) * VAT3( dPS, ii,jj,kkp1)
2226 - VAT3( uE, im1,j,k) * VAT3( dPC, ii,jj,kkp1)
2227 - VAT3( uNE, im1,j,k) * VAT3( dPN, ii,jj,kkp1));
2228
2229 TMP3_XUC =
2230 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2231 - VAT3( uS, im1,j,kp1) * VAT3(oPSW, ii,jj,kkp1)
2232 + VAT3( oC, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2233 - VAT3( uC, im1,j,kp1) * VAT3( oPW, ii,jj,kkp1)
2234 - VAT3( oN, im1,j,kp1) * VAT3(dPNW, ii,jj,kkp1)
2235 - VAT3( uN, im1,j,kp1) * VAT3(oPNW, ii,jj,kkp1)
2236 - VAT3( oNW, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2237 - VAT3( uSE, im1,j,kp1) * VAT3( oPS, ii,jj,kkp1)
2238 - VAT3( oE, im1,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2239 - VAT3( uE, im1,j,kp1) * VAT3( oPC, ii,jj,kkp1)
2240 - VAT3( oNE, im1,j,kp1) * VAT3( dPN, ii,jj,kkp1)
2241 - VAT3( uNE, im1,j,kp1) * VAT3( oPN, ii,jj,kkp1))
2242
2243 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uS, im1,jp1,k) * VAT3( dPW, ii,jj,kkp1)
2244 - VAT3( uC, im1,jp1,k) * VAT3(dPNW, ii,jj,kkp1)
2245 - VAT3( uSE, im1,jp1,k) * VAT3( dPC, ii,jj,kkp1)
2246 - VAT3( uE, im1,jp1,k) * VAT3( dPN, ii,jj,kkp1));
2247
2248 TMP4_XUC =
2249 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oN, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2250 - VAT3( uS, im1,jp1,kp1) * VAT3( oPW, ii,jj,kkp1)
2251 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
2252 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNW, ii,jj,kkp1)
2253 - VAT3( oNW, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2254 - VAT3( uSE, im1,jp1,kp1) * VAT3( oPC, ii,jj,kkp1)
2255 - VAT3( oE, im1,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
2256 - VAT3( uE, im1,jp1,kp1) * VAT3( oPN, ii,jj,kkp1))
2257
2258 - VAT3( oPS, ii,jj,kk) * (- VAT3( uW, i,jm1,k) * VAT3(dPSW, ii,jj,kkp1)
2259 - VAT3( uNW, i,jm1,k) * VAT3( dPW, ii,jj,kkp1)
2260 - VAT3( uC, i,jm1,k) * VAT3( dPS, ii,jj,kkp1)
2261 - VAT3( uN, i,jm1,k) * VAT3( dPC, ii,jj,kkp1)
2262 - VAT3( uE, i,jm1,k) * VAT3(dPSE, ii,jj,kkp1)
2263 - VAT3( uNE, i,jm1,k) * VAT3( dPE, ii,jj,kkp1));
2264
2265 TMP5_XUC =
2266 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2267 - VAT3( uW, i,jm1,kp1) * VAT3(oPSW, ii,jj,kkp1)
2268 - VAT3( oNW, i,jm1,kp1) * VAT3( dPW, ii,jj,kkp1)
2269 - VAT3( uNW, i,jm1,kp1) * VAT3( oPW, ii,jj,kkp1)
2270 + VAT3( oC, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2271 - VAT3( uC, i,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
2272 - VAT3( oN, i,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
2273 - VAT3( uN, i,jm1,kp1) * VAT3( oPC, ii,jj,kkp1)
2274 - VAT3( oE, i,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2275 - VAT3( uE, i,jm1,kp1) * VAT3(oPSE, ii,jj,kkp1)
2276 - VAT3( oNE, i,jm1,kp1) * VAT3( dPE, ii,jj,kkp1)
2277 - VAT3( uNE, i,jm1,kp1) * VAT3( oPE, ii,jj,kkp1));
2278
2279 TMP6_XUC =
2280 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, i,j,k) * VAT3(dPSW, ii,jj,kkp1)
2281 - VAT3( uW, i,j,k) * VAT3( dPW, ii,jj,kkp1)
2282 - VAT3( uNW, i,j,k) * VAT3(dPNW, ii,jj,kkp1)
2283 - VAT3( uS, i,j,k) * VAT3( dPS, ii,jj,kkp1)
2284 - VAT3( uC, i,j,k) * VAT3( dPC, ii,jj,kkp1)
2285 - VAT3( uN, i,j,k) * VAT3( dPN, ii,jj,kkp1)
2286 - VAT3( uSE, i,j,k) * VAT3(dPSE, ii,jj,kkp1)
2287 - VAT3( uE, i,j,k) * VAT3( dPE, ii,jj,kkp1)
2288 - VAT3( uNE, i,j,k) * VAT3(dPNE, ii,jj,kkp1));
2289
2290 TMP7_XUC =
2291 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPSW, ii,jj,kkp1)
2292 - VAT3( uSW, i,j,kp1) * VAT3(oPSW, ii,jj,kkp1)
2293 - VAT3( oE, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2294 - VAT3( uW, i,j,kp1) * VAT3( oPW, ii,jj,kkp1)
2295 - VAT3( oNW, i,j,kp1) * VAT3(dPNW, ii,jj,kkp1)
2296 - VAT3( uNW, i,j,kp1) * VAT3(oPNW, ii,jj,kkp1)
2297 - VAT3( oN, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2298 - VAT3( uS, i,j,kp1) * VAT3( oPS, ii,jj,kkp1)
2299 + VAT3( oC, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2300 - VAT3( uC, i,j,kp1) * VAT3( oPC, ii,jj,kkp1)
2301 - VAT3( oN, i,j,kp1) * VAT3( dPN, ii,jj,kkp1)
2302 - VAT3( uN, i,j,kp1) * VAT3( oPN, ii,jj,kkp1)
2303 - VAT3( oNW, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2304 - VAT3( uSE, i,j,kp1) * VAT3(oPSE, ii,jj,kkp1)
2305 - VAT3( oE, i,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2306 - VAT3( uE, i,j,kp1) * VAT3( oPE, ii,jj,kkp1)
2307 - VAT3( oNE, i,j,kp1) * VAT3(dPNE, ii,jj,kkp1)
2308 - VAT3( uNE, i,j,kp1) * VAT3(oPNE, ii,jj,kkp1));
2309
2310 TMP8_XUC =
2311 - VAT3( oPN, ii,jj,kk) * (- VAT3( uSW, i,jp1,k) * VAT3( dPW, ii,jj,kkp1)
2312 - VAT3( uW, i,jp1,k) * VAT3(dPNW, ii,jj,kkp1)
2313 - VAT3( uS, i,jp1,k) * VAT3( dPC, ii,jj,kkp1)
2314 - VAT3( uC, i,jp1,k) * VAT3( dPN, ii,jj,kkp1)
2315 - VAT3( uSE, i,jp1,k) * VAT3( dPE, ii,jj,kkp1)
2316 - VAT3( uE, i,jp1,k) * VAT3(dPNE, ii,jj,kkp1));
2317
2318 TMP9_XUC =
2319 - VAT3( uPN, ii,jj,kk) * (- VAT3( oNE, im1,j,kp1) * VAT3( dPW, ii,jj,kkp1)
2320 - VAT3( uSW, i,jp1,kp1) * VAT3( oPW, ii,jj,kkp1)
2321 - VAT3( oE, im1,jp1,kp1) * VAT3(dPNW, ii,jj,kkp1)
2322 - VAT3( uW, i,jp1,kp1) * VAT3(oPNW, ii,jj,kkp1)
2323 - VAT3( oN, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2324 - VAT3( uS, i,jp1,kp1) * VAT3( oPC, ii,jj,kkp1)
2325 + VAT3( oC, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
2326 - VAT3( uC, i,jp1,kp1) * VAT3( oPN, ii,jj,kkp1)
2327 - VAT3( oNW, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2328 - VAT3( uSE, i,jp1,kp1) * VAT3( oPE, ii,jj,kkp1)
2329 - VAT3( oE, i,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1)
2330 - VAT3( uE, i,jp1,kp1) * VAT3(oPNE, ii,jj,kkp1))
2331
2332 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uW, ip1,jm1,k) * VAT3( dPS, ii,jj,kkp1)
2333 - VAT3( uNW, ip1,jm1,k) * VAT3( dPC, ii,jj,kkp1)
2334 - VAT3( uC, ip1,jm1,k) * VAT3(dPSE, ii,jj,kkp1)
2335 - VAT3( uN, ip1,jm1,k) * VAT3( dPE, ii,jj,kkp1));
2336
2337 TMP10_XUC =
2338 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2339 - VAT3( uW, ip1,jm1,kp1) * VAT3( oPS, ii,jj,kkp1)
2340 - VAT3( oNW, ip1,jm1,kp1) * VAT3( dPC, ii,jj,kkp1)
2341 - VAT3( uNW, ip1,jm1,kp1) * VAT3( oPC, ii,jj,kkp1)
2342 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2343 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSE, ii,jj,kkp1)
2344 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPE, ii,jj,kkp1)
2345 - VAT3( uN, ip1,jm1,kp1) * VAT3( oPE, ii,jj,kkp1))
2346
2347 - VAT3( oPE, ii,jj,kk) * (- VAT3( uSW, ip1,j,k) * VAT3( dPS, ii,jj,kkp1)
2348 - VAT3( uW, ip1,j,k) * VAT3( dPC, ii,jj,kkp1)
2349 - VAT3( uNW, ip1,j,k) * VAT3( dPN, ii,jj,kkp1)
2350 - VAT3( uS, ip1,j,k) * VAT3(dPSE, ii,jj,kkp1)
2351 - VAT3( uC, ip1,j,k) * VAT3( dPE, ii,jj,kkp1)
2352 - VAT3( uN, ip1,j,k) * VAT3(dPNE, ii,jj,kkp1));
2353
2354 TMP11_XUC =
2355 - VAT3( uPE, ii,jj,kk) * (- VAT3( oNE, i,jm1,kp1) * VAT3( dPS, ii,jj,kkp1)
2356 - VAT3( uSW, ip1,j,kp1) * VAT3( oPS, ii,jj,kkp1)
2357 - VAT3( oE, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2358 - VAT3( uW, ip1,j,kp1) * VAT3( oPC, ii,jj,kkp1)
2359 - VAT3( oNW, ip1,j,kp1) * VAT3( dPN, ii,jj,kkp1)
2360 - VAT3( uNW, ip1,j,kp1) * VAT3( oPN, ii,jj,kkp1)
2361 - VAT3( oN, ip1,jm1,kp1) * VAT3(dPSE, ii,jj,kkp1)
2362 - VAT3( uS, ip1,j,kp1) * VAT3(oPSE, ii,jj,kkp1)
2363 + VAT3( oC, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2364 - VAT3( uC, ip1,j,kp1) * VAT3( oPE, ii,jj,kkp1)
2365 - VAT3( oN, ip1,j,kp1) * VAT3(dPNE, ii,jj,kkp1)
2366 - VAT3( uN, ip1,j,kp1) * VAT3(oPNE, ii,jj,kkp1));
2367
2368 TMP12_XUC =
2369 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uSW, ip1,jp1,k) * VAT3( dPC, ii,jj,kkp1)
2370 - VAT3( uW, ip1,jp1,k) * VAT3( dPN, ii,jj,kkp1)
2371 - VAT3( uS, ip1,jp1,k) * VAT3( dPE, ii,jj,kkp1)
2372 - VAT3( uC, ip1,jp1,k) * VAT3(dPNE, ii,jj,kkp1))
2373
2374 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oNE, i,j,kp1) * VAT3( dPC, ii,jj,kkp1)
2375 - VAT3( uSW, ip1,jp1,kp1) * VAT3( oPC, ii,jj,kkp1)
2376 - VAT3( oE, i,jp1,kp1) * VAT3( dPN, ii,jj,kkp1)
2377 - VAT3( uW, ip1,jp1,kp1) * VAT3( oPN, ii,jj,kkp1)
2378 - VAT3( oN, ip1,j,kp1) * VAT3( dPE, ii,jj,kkp1)
2379 - VAT3( uS, ip1,jp1,kp1) * VAT3( oPE, ii,jj,kkp1)
2380 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNE, ii,jj,kkp1)
2381 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNE, ii,jj,kkp1));
2382
2383 VAT3( XuC, ii,jj,kk) = TMP1_XUC + TMP2_XUC + TMP3_XUC + TMP4_XUC
2384 + TMP5_XUC + TMP6_XUC + TMP7_XUC + TMP8_XUC + TMP9_XUC
2385 + TMP10_XUC + TMP11_XUC + TMP12_XUC;
2386
2387 //fprintf(data, "%19.12E\n", VAT3(XuC, ii, jj, kk));
2388
2389 /* *************************************************************
2390 * *** > ONE;
2391 * *************************************************************/
2392
2393 // VAT3(XoNE, ii,jj,kk) =
2394 TMP1_XONE =
2395 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNE, i,j,km1) * VAT3(dPSW, iip1,jjp1,kk)
2396 - VAT3( uNE, i,j,km1) * VAT3(oPSW, iip1,jjp1,kk))
2397
2398 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2399 - VAT3( oNE, i,j,k) * VAT3(oPSW, iip1,jjp1,kk)
2400 - VAT3( uNE, i,j,k) * VAT3(uPSW, iip1,jjp1,kk))
2401
2402 - VAT3( uPC, ii,jj,kk) * (- VAT3( uSW, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2403 - VAT3( oNE, i,j,kp1) * VAT3(uPSW, iip1,jjp1,kk))
2404
2405 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, i,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2406 - VAT3( uE, i,jp1,km1) * VAT3(oPSW, iip1,jjp1,kk)
2407 - VAT3( oNE, i,jp1,km1) * VAT3( dPW, iip1,jjp1,kk)
2408 - VAT3( uNE, i,jp1,km1) * VAT3( oPW, iip1,jjp1,kk));
2409
2410 TMP2_XONE =
2411 - VAT3( oPN, ii,jj,kk) * (- VAT3( uW, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2412 - VAT3( oE, i,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2413 - VAT3( uE, i,jp1,k) * VAT3(uPSW, iip1,jjp1,kk)
2414 - VAT3( uSW, ip1,jp2,km1) * VAT3( dPW, iip1,jjp1,kk)
2415 - VAT3( oNE, i,jp1,k) * VAT3( oPW, iip1,jjp1,kk)
2416 - VAT3( uNE, i,jp1,k) * VAT3( uPW, iip1,jjp1,kk))
2417
2418 - VAT3( uPN, ii,jj,kk) * (- VAT3( uW, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2419 - VAT3( oE, i,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
2420 - VAT3( uSW, ip1,jp2,k) * VAT3( oPW, iip1,jjp1,kk)
2421 - VAT3( oNE, i,jp1,kp1) * VAT3( uPW, iip1,jjp1,kk))
2422
2423 - VAT3( dPE, ii,jj,kk) * (- VAT3( oN, ip1,j,km1) * VAT3(dPSW, iip1,jjp1,kk)
2424 - VAT3( uN, ip1,j,km1) * VAT3(oPSW, iip1,jjp1,kk)
2425 - VAT3( oNE, ip1,j,km1) * VAT3( dPS, iip1,jjp1,kk)
2426 - VAT3( uNE, ip1,j,km1) * VAT3( oPS, iip1,jjp1,kk));
2427
2428 TMP3_XONE =
2429 - VAT3( oPE, ii,jj,kk) * (- VAT3( uS, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2430 - VAT3( oN, ip1,j,k) * VAT3(oPSW, iip1,jjp1,kk)
2431 - VAT3( uN, ip1,j,k) * VAT3(uPSW, iip1,jjp1,kk)
2432 - VAT3( uSW, ip2,jp1,km1) * VAT3( dPS, iip1,jjp1,kk)
2433 - VAT3( oNE, ip1,j,k) * VAT3( oPS, iip1,jjp1,kk)
2434 - VAT3( uNE, ip1,j,k) * VAT3( uPS, iip1,jjp1,kk))
2435
2436 - VAT3( uPE, ii,jj,kk) * (- VAT3( uS, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2437 - VAT3( oN, ip1,j,kp1) * VAT3(uPSW, iip1,jjp1,kk)
2438 - VAT3( uSW, ip2,jp1,k) * VAT3( oPS, iip1,jjp1,kk)
2439 - VAT3( oNE, ip1,j,kp1) * VAT3( uPS, iip1,jjp1,kk));
2440
2441 TMP4_XONE =
2442 - VAT3(dPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2443 - VAT3( uC, ip1,jp1,km1) * VAT3(oPSW, iip1,jjp1,kk)
2444 - VAT3( oN, ip1,jp1,km1) * VAT3( dPW, iip1,jjp1,kk)
2445 - VAT3( uN, ip1,jp1,km1) * VAT3( oPW, iip1,jjp1,kk)
2446 - VAT3( oE, ip1,jp1,km1) * VAT3( dPS, iip1,jjp1,kk)
2447 - VAT3( uE, ip1,jp1,km1) * VAT3( oPS, iip1,jjp1,kk)
2448 - VAT3( oNE, ip1,jp1,km1) * VAT3( dPC, iip1,jjp1,kk)
2449 - VAT3( uNE, ip1,jp1,km1) * VAT3( oPC, iip1,jjp1,kk));
2450
2451 TMP5_XONE =
2452 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,km1) * VAT3(dPSW, iip1,jjp1,kk)
2453 + VAT3( oC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2454 - VAT3( uC, ip1,jp1,k) * VAT3(uPSW, iip1,jjp1,kk)
2455 - VAT3( uS, ip1,jp2,km1) * VAT3( dPW, iip1,jjp1,kk)
2456 - VAT3( oN, ip1,jp1,k) * VAT3( oPW, iip1,jjp1,kk)
2457 - VAT3( uN, ip1,jp1,k) * VAT3( uPW, iip1,jjp1,kk)
2458 - VAT3( uW, ip2,jp1,km1) * VAT3( dPS, iip1,jjp1,kk)
2459 - VAT3( oE, ip1,jp1,k) * VAT3( oPS, iip1,jjp1,kk)
2460 - VAT3( uE, ip1,jp1,k) * VAT3( uPS, iip1,jjp1,kk)
2461 - VAT3( uSW, ip2,jp2,km1) * VAT3( dPC, iip1,jjp1,kk)
2462 - VAT3( oNE, ip1,jp1,k) * VAT3( oPC, iip1,jjp1,kk)
2463 - VAT3( uNE, ip1,jp1,k) * VAT3( uPC, iip1,jjp1,kk));
2464
2465 TMP6_XONE =
2466 - VAT3(uPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,k) * VAT3(oPSW, iip1,jjp1,kk)
2467 + VAT3( oC, ip1,jp1,kp1) * VAT3(uPSW, iip1,jjp1,kk)
2468 - VAT3( uS, ip1,jp2,k) * VAT3( oPW, iip1,jjp1,kk)
2469 - VAT3( oN, ip1,jp1,kp1) * VAT3( uPW, iip1,jjp1,kk)
2470 - VAT3( uW, ip2,jp1,k) * VAT3( oPS, iip1,jjp1,kk)
2471 - VAT3( oE, ip1,jp1,kp1) * VAT3( uPS, iip1,jjp1,kk)
2472 - VAT3( uSW, ip2,jp2,k) * VAT3( oPC, iip1,jjp1,kk)
2473 - VAT3( oNE, ip1,jp1,kp1) * VAT3( uPC, iip1,jjp1,kk));
2474
2475 VAT3(XoNE, ii,jj,kk) = TMP1_XONE + TMP2_XONE + TMP3_XONE + TMP4_XONE
2476 + TMP5_XONE + TMP6_XONE;
2477
2478 //fprintf(data, "%19.12E\n", VAT3(XoNE, ii, jj, kk));
2479
2480 /* *************************************************************
2481 * *** > ONW;
2482 * *************************************************************/
2483
2484 // VAT3(XoNW, ii,jj,kk) =
2485 TMP1_XONW =
2486 - VAT3( dPW, ii,jj,kk) * (- VAT3( oNW, im1,j,km1) * VAT3( dPS, iim1,jjp1,kk)
2487 - VAT3( uNW, im1,j,km1) * VAT3( oPS, iim1,jjp1,kk)
2488 - VAT3( oN, im1,j,km1) * VAT3(dPSE, iim1,jjp1,kk)
2489 - VAT3( uN, im1,j,km1) * VAT3(oPSE, iim1,jjp1,kk))
2490
2491 - VAT3( oPW, ii,jj,kk) * (- VAT3( uSE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
2492 - VAT3( oNW, im1,j,k) * VAT3( oPS, iim1,jjp1,kk)
2493 - VAT3( uNW, im1,j,k) * VAT3( uPS, iim1,jjp1,kk)
2494 - VAT3( uS, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2495 - VAT3( oN, im1,j,k) * VAT3(oPSE, iim1,jjp1,kk)
2496 - VAT3( uN, im1,j,k) * VAT3(uPSE, iim1,jjp1,kk));
2497
2498 TMP2_XONW =
2499 - VAT3( uPW, ii,jj,kk) * (- VAT3( uSE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
2500 - VAT3( oNW, im1,j,kp1) * VAT3( uPS, iim1,jjp1,kk)
2501 - VAT3( uS, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2502 - VAT3( oN, im1,j,kp1) * VAT3(uPSE, iim1,jjp1,kk))
2503
2504 - VAT3(dPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
2505 - VAT3( uW, im1,jp1,km1) * VAT3( oPS, iim1,jjp1,kk)
2506 - VAT3( oNW, im1,jp1,km1) * VAT3( dPC, iim1,jjp1,kk)
2507 - VAT3( uNW, im1,jp1,km1) * VAT3( oPC, iim1,jjp1,kk)
2508 + VAT3( oC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2509 - VAT3( uC, im1,jp1,km1) * VAT3(oPSE, iim1,jjp1,kk)
2510 - VAT3( oN, im1,jp1,km1) * VAT3( dPE, iim1,jjp1,kk)
2511 - VAT3( uN, im1,jp1,km1) * VAT3( oPE, iim1,jjp1,kk));
2512
2513 TMP3_XONW =
2514 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uE, im2,jp1,km1) * VAT3( dPS, iim1,jjp1,kk)
2515 - VAT3( oE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
2516 - VAT3( uW, im1,jp1,k) * VAT3( uPS, iim1,jjp1,kk)
2517 - VAT3( uSE, im2,jp2,km1) * VAT3( dPC, iim1,jjp1,kk)
2518 - VAT3( oNW, im1,jp1,k) * VAT3( oPC, iim1,jjp1,kk)
2519 - VAT3( uNW, im1,jp1,k) * VAT3( uPC, iim1,jjp1,kk)
2520 - VAT3( uC, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2521 + VAT3( oC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2522 - VAT3( uC, im1,jp1,k) * VAT3(uPSE, iim1,jjp1,kk)
2523 - VAT3( uS, im1,jp2,km1) * VAT3( dPE, iim1,jjp1,kk)
2524 - VAT3( oN, im1,jp1,k) * VAT3( oPE, iim1,jjp1,kk)
2525 - VAT3( uN, im1,jp1,k) * VAT3( uPE, iim1,jjp1,kk));
2526
2527 TMP4_XONW =
2528 - VAT3(uPNW, ii,jj,kk) * (- VAT3( uE, im2,jp1,k) * VAT3( oPS, iim1,jjp1,kk)
2529 - VAT3( oE, im2,jp1,kp1) * VAT3( uPS, iim1,jjp1,kk)
2530 - VAT3( uSE, im2,jp2,k) * VAT3( oPC, iim1,jjp1,kk)
2531 - VAT3( oNW, im1,jp1,kp1) * VAT3( uPC, iim1,jjp1,kk)
2532 - VAT3( uC, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2533 + VAT3( oC, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk)
2534 - VAT3( uS, im1,jp2,k) * VAT3( oPE, iim1,jjp1,kk)
2535 - VAT3( oN, im1,jp1,kp1) * VAT3( uPE, iim1,jjp1,kk))
2536
2537 - VAT3( dPC, ii,jj,kk) * (- VAT3( oNW, i,j,km1) * VAT3(dPSE, iim1,jjp1,kk)
2538 - VAT3( uNW, i,j,km1) * VAT3(oPSE, iim1,jjp1,kk));
2539
2540 TMP5_XONW =
2541 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2542 - VAT3( oNW, i,j,k) * VAT3(oPSE, iim1,jjp1,kk)
2543 - VAT3( uNW, i,j,k) * VAT3(uPSE, iim1,jjp1,kk))
2544
2545 - VAT3( uPC, ii,jj,kk) * (- VAT3( uSE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2546 - VAT3( oNW, i,j,kp1) * VAT3(uPSE, iim1,jjp1,kk))
2547
2548 - VAT3( dPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2549 - VAT3( uW, i,jp1,km1) * VAT3(oPSE, iim1,jjp1,kk)
2550 - VAT3( oNW, i,jp1,km1) * VAT3( dPE, iim1,jjp1,kk)
2551 - VAT3( uNW, i,jp1,km1) * VAT3( oPE, iim1,jjp1,kk));
2552
2553 TMP6_XONW =
2554 - VAT3( oPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,km1) * VAT3(dPSE, iim1,jjp1,kk)
2555 - VAT3( oE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2556 - VAT3( uW, i,jp1,k) * VAT3(uPSE, iim1,jjp1,kk)
2557 - VAT3( uSE, im1,jp2,km1) * VAT3( dPE, iim1,jjp1,kk)
2558 - VAT3( oNW, i,jp1,k) * VAT3( oPE, iim1,jjp1,kk)
2559 - VAT3( uNW, i,jp1,k) * VAT3( uPE, iim1,jjp1,kk))
2560
2561 - VAT3( uPN, ii,jj,kk) * (- VAT3( uE, im1,jp1,k) * VAT3(oPSE, iim1,jjp1,kk)
2562 - VAT3( oE, im1,jp1,kp1) * VAT3(uPSE, iim1,jjp1,kk)
2563 - VAT3( uSE, im1,jp2,k) * VAT3( oPE, iim1,jjp1,kk)
2564 - VAT3( oNW, i,jp1,kp1) * VAT3( uPE, iim1,jjp1,kk));
2565
2566 VAT3(XoNW, ii,jj,kk) = TMP1_XONW + TMP2_XONW + TMP3_XONW + TMP4_XONW
2567 + TMP5_XONW + TMP6_XONW;
2568
2569 //fprintf(data, "%19.12E\n", VAT3(XoNW, ii, jj, kk));
2570
2571 /* *************************************************************
2572 * *** > UE;
2573 * *************************************************************/
2574 // VAT3( XuE, ii,jj,kk) =
2575 TMP1_XUE =
2576 - VAT3( oPS, ii,jj,kk) * (- VAT3( uE, i,jm1,k) * VAT3(dPSW, iip1,jj,kkp1)
2577 - VAT3( uNE, i,jm1,k) * VAT3( dPW, iip1,jj,kkp1))
2578
2579 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, i,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2580 - VAT3( uE, i,jm1,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2581 - VAT3( oNE, i,jm1,kp1) * VAT3( dPW, iip1,jj,kkp1)
2582 - VAT3( uNE, i,jm1,kp1) * VAT3( oPW, iip1,jj,kkp1))
2583
2584 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSE, i,j,k) * VAT3(dPSW, iip1,jj,kkp1)
2585 - VAT3( uE, i,j,k) * VAT3( dPW, iip1,jj,kkp1)
2586 - VAT3( uNE, i,j,k) * VAT3(dPNW, iip1,jj,kkp1));
2587
2588 TMP2_XUE =
2589 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2590 - VAT3( uSE, i,j,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2591 - VAT3( oE, i,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2592 - VAT3( uE, i,j,kp1) * VAT3( oPW, iip1,jj,kkp1)
2593 - VAT3( oNE, i,j,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2594 - VAT3( uNE, i,j,kp1) * VAT3(oPNW, iip1,jj,kkp1))
2595
2596 - VAT3( oPN, ii,jj,kk) * (- VAT3( uSE, i,jp1,k) * VAT3( dPW, iip1,jj,kkp1)
2597 - VAT3( uE, i,jp1,k) * VAT3(dPNW, iip1,jj,kkp1))
2598
2599 - VAT3( uPN, ii,jj,kk) * (- VAT3( oNW, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2600 - VAT3( uSE, i,jp1,kp1) * VAT3( oPW, iip1,jj,kkp1)
2601 - VAT3( oE, i,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2602 - VAT3( uE, i,jp1,kp1) * VAT3(oPNW, iip1,jj,kkp1));
2603
2604 TMP3_XUE =
2605 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uC, ip1,jm1,k) * VAT3(dPSW, iip1,jj,kkp1)
2606 - VAT3( uN, ip1,jm1,k) * VAT3( dPW, iip1,jj,kkp1)
2607 - VAT3( uE, ip1,jm1,k) * VAT3( dPS, iip1,jj,kkp1)
2608 - VAT3( uNE, ip1,jm1,k) * VAT3( dPC, iip1,jj,kkp1))
2609
2610 - VAT3(uPSE, ii,jj,kk) * ( VAT3( oC, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2611 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2612 - VAT3( oN, ip1,jm1,kp1) * VAT3( dPW, iip1,jj,kkp1)
2613 - VAT3( uN, ip1,jm1,kp1) * VAT3( oPW, iip1,jj,kkp1)
2614 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPS, iip1,jj,kkp1)
2615 - VAT3( uE, ip1,jm1,kp1) * VAT3( oPS, iip1,jj,kkp1)
2616 - VAT3( oNE, ip1,jm1,kp1) * VAT3( dPC, iip1,jj,kkp1)
2617 - VAT3( uNE, ip1,jm1,kp1) * VAT3( oPC, iip1,jj,kkp1));
2618
2619 TMP4_XUE =
2620 - VAT3( oPE, ii,jj,kk) * (- VAT3( uS, ip1,j,k) * VAT3(dPSW, iip1,jj,kkp1)
2621 - VAT3( uC, ip1,j,k) * VAT3( dPW, iip1,jj,kkp1)
2622 - VAT3( uN, ip1,j,k) * VAT3(dPNW, iip1,jj,kkp1)
2623 - VAT3( uSE, ip1,j,k) * VAT3( dPS, iip1,jj,kkp1)
2624 - VAT3( uE, ip1,j,k) * VAT3( dPC, iip1,jj,kkp1)
2625 - VAT3( uNE, ip1,j,k) * VAT3( dPN, iip1,jj,kkp1));
2626
2627 TMP5_XUE =
2628 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(dPSW, iip1,jj,kkp1)
2629 - VAT3( uS, ip1,j,kp1) * VAT3(oPSW, iip1,jj,kkp1)
2630 + VAT3( oC, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2631 - VAT3( uC, ip1,j,kp1) * VAT3( oPW, iip1,jj,kkp1)
2632 - VAT3( oN, ip1,j,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2633 - VAT3( uN, ip1,j,kp1) * VAT3(oPNW, iip1,jj,kkp1)
2634 - VAT3( oNW, ip2,jm1,kp1) * VAT3( dPS, iip1,jj,kkp1)
2635 - VAT3( uSE, ip1,j,kp1) * VAT3( oPS, iip1,jj,kkp1)
2636 - VAT3( oE, ip1,j,kp1) * VAT3( dPC, iip1,jj,kkp1)
2637 - VAT3( uE, ip1,j,kp1) * VAT3( oPC, iip1,jj,kkp1)
2638 - VAT3( oNE, ip1,j,kp1) * VAT3( dPN, iip1,jj,kkp1)
2639 - VAT3( uNE, ip1,j,kp1) * VAT3( oPN, iip1,jj,kkp1));
2640
2641 TMP6_XUE =
2642 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uS, ip1,jp1,k) * VAT3( dPW, iip1,jj,kkp1)
2643 - VAT3( uC, ip1,jp1,k) * VAT3(dPNW, iip1,jj,kkp1)
2644 - VAT3( uSE, ip1,jp1,k) * VAT3( dPC, iip1,jj,kkp1)
2645 - VAT3( uE, ip1,jp1,k) * VAT3( dPN, iip1,jj,kkp1))
2646
2647 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oN, ip1,j,kp1) * VAT3( dPW, iip1,jj,kkp1)
2648 - VAT3( uS, ip1,jp1,kp1) * VAT3( oPW, iip1,jj,kkp1)
2649 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPNW, iip1,jj,kkp1)
2650 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPNW, iip1,jj,kkp1)
2651 - VAT3( oNW, ip2,j,kp1) * VAT3( dPC, iip1,jj,kkp1)
2652 - VAT3( uSE, ip1,jp1,kp1) * VAT3( oPC, iip1,jj,kkp1)
2653 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPN, iip1,jj,kkp1)
2654 - VAT3( uE, ip1,jp1,kp1) * VAT3( oPN, iip1,jj,kkp1));
2655
2656 VAT3( XuE, ii,jj,kk) = TMP1_XUE + TMP2_XUE + TMP3_XUE + TMP4_XUE
2657 + TMP5_XUE + TMP6_XUE;
2658
2659 //fprintf(data, "%19.12E\n", VAT3(XuE, ii, jj, kk));
2660
2661 /* *************************************************************
2662 * *** > UW;
2663 * *************************************************************/
2664 // VAT3( XuW, ii,jj,kk) =
2665 TMP1_XUW =
2666 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uW, im1,jm1,k) * VAT3( dPS, iim1,jj,kkp1)
2667 - VAT3( uNW, im1,jm1,k) * VAT3( dPC, iim1,jj,kkp1)
2668 - VAT3( uC, im1,jm1,k) * VAT3(dPSE, iim1,jj,kkp1)
2669 - VAT3( uN, im1,jm1,k) * VAT3( dPE, iim1,jj,kkp1))
2670
2671 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oE, im2,jm1,kp1) * VAT3( dPS, iim1,jj,kkp1)
2672 - VAT3( uW, im1,jm1,kp1) * VAT3( oPS, iim1,jj,kkp1)
2673 - VAT3( oNW, im1,jm1,kp1) * VAT3( dPC, iim1,jj,kkp1)
2674 - VAT3( uNW, im1,jm1,kp1) * VAT3( oPC, iim1,jj,kkp1)
2675 + VAT3( oC, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2676 - VAT3( uC, im1,jm1,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2677 - VAT3( oN, im1,jm1,kp1) * VAT3( dPE, iim1,jj,kkp1)
2678 - VAT3( uN, im1,jm1,kp1) * VAT3( oPE, iim1,jj,kkp1));
2679
2680 TMP2_XUW =
2681 - VAT3( oPW, ii,jj,kk) * (- VAT3( uSW, im1,j,k) * VAT3( dPS, iim1,jj,kkp1)
2682 - VAT3( uW, im1,j,k) * VAT3( dPC, iim1,jj,kkp1)
2683 - VAT3( uNW, im1,j,k) * VAT3( dPN, iim1,jj,kkp1)
2684 - VAT3( uS, im1,j,k) * VAT3(dPSE, iim1,jj,kkp1)
2685 - VAT3( uC, im1,j,k) * VAT3( dPE, iim1,jj,kkp1)
2686 - VAT3( uN, im1,j,k) * VAT3(dPNE, iim1,jj,kkp1));
2687
2688 TMP3_XUW =
2689 - VAT3( uPW, ii,jj,kk) * (- VAT3( oNE, im2,jm1,kp1) * VAT3( dPS, iim1,jj,kkp1)
2690 - VAT3( uSW, im1,j,kp1) * VAT3( oPS, iim1,jj,kkp1)
2691 - VAT3( oE, im2,j,kp1) * VAT3( dPC, iim1,jj,kkp1)
2692 - VAT3( uW, im1,j,kp1) * VAT3( oPC, iim1,jj,kkp1)
2693 - VAT3( oNW, im1,j,kp1) * VAT3( dPN, iim1,jj,kkp1)
2694 - VAT3( uNW, im1,j,kp1) * VAT3( oPN, iim1,jj,kkp1)
2695 - VAT3( oN, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2696 - VAT3( uS, im1,j,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2697 + VAT3( oC, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2698 - VAT3( uC, im1,j,kp1) * VAT3( oPE, iim1,jj,kkp1)
2699 - VAT3( oN, im1,j,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2700 - VAT3( uN, im1,j,kp1) * VAT3(oPNE, iim1,jj,kkp1));
2701
2702 TMP4_XUW =
2703 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uSW, im1,jp1,k) * VAT3( dPC, iim1,jj,kkp1)
2704 - VAT3( uW, im1,jp1,k) * VAT3( dPN, iim1,jj,kkp1)
2705 - VAT3( uS, im1,jp1,k) * VAT3( dPE, iim1,jj,kkp1)
2706 - VAT3( uC, im1,jp1,k) * VAT3(dPNE, iim1,jj,kkp1))
2707
2708 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oNE, im2,j,kp1) * VAT3( dPC, iim1,jj,kkp1)
2709 - VAT3( uSW, im1,jp1,kp1) * VAT3( oPC, iim1,jj,kkp1)
2710 - VAT3( oE, im2,jp1,kp1) * VAT3( dPN, iim1,jj,kkp1)
2711 - VAT3( uW, im1,jp1,kp1) * VAT3( oPN, iim1,jj,kkp1)
2712 - VAT3( oN, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2713 - VAT3( uS, im1,jp1,kp1) * VAT3( oPE, iim1,jj,kkp1)
2714 + VAT3( oC, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2715 - VAT3( uC, im1,jp1,kp1) * VAT3(oPNE, iim1,jj,kkp1));
2716
2717 TMP5_XUW =
2718 - VAT3( oPS, ii,jj,kk) * (- VAT3( uW, i,jm1,k) * VAT3(dPSE, iim1,jj,kkp1)
2719 - VAT3( uNW, i,jm1,k) * VAT3( dPE, iim1,jj,kkp1))
2720
2721 - VAT3( uPS, ii,jj,kk) * (- VAT3( oE, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2722 - VAT3( uW, i,jm1,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2723 - VAT3( oNW, i,jm1,kp1) * VAT3( dPE, iim1,jj,kkp1)
2724 - VAT3( uNW, i,jm1,kp1) * VAT3( oPE, iim1,jj,kkp1))
2725
2726 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, i,j,k) * VAT3(dPSE, iim1,jj,kkp1)
2727 - VAT3( uW, i,j,k) * VAT3( dPE, iim1,jj,kkp1)
2728 - VAT3( uNW, i,j,k) * VAT3(dPNE, iim1,jj,kkp1));
2729
2730 TMP6_XUW =
2731 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPSE, iim1,jj,kkp1)
2732 - VAT3( uSW, i,j,kp1) * VAT3(oPSE, iim1,jj,kkp1)
2733 - VAT3( oE, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2734 - VAT3( uW, i,j,kp1) * VAT3( oPE, iim1,jj,kkp1)
2735 - VAT3( oNW, i,j,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2736 - VAT3( uNW, i,j,kp1) * VAT3(oPNE, iim1,jj,kkp1))
2737
2738 - VAT3( oPN, ii,jj,kk) * (- VAT3( uSW, i,jp1,k) * VAT3( dPE, iim1,jj,kkp1)
2739 - VAT3( uW, i,jp1,k) * VAT3(dPNE, iim1,jj,kkp1))
2740
2741 - VAT3( uPN, ii,jj,kk) * (- VAT3( oNE, im1,j,kp1) * VAT3( dPE, iim1,jj,kkp1)
2742 - VAT3( uSW, i,jp1,kp1) * VAT3( oPE, iim1,jj,kkp1)
2743 - VAT3( oE, im1,jp1,kp1) * VAT3(dPNE, iim1,jj,kkp1)
2744 - VAT3( uW, i,jp1,kp1) * VAT3(oPNE, iim1,jj,kkp1));
2745
2746 VAT3( XuW, ii,jj,kk) = TMP1_XUW + TMP2_XUW + TMP3_XUW + TMP4_XUW
2747 + TMP5_XUW + TMP6_XUW;
2748
2749 //fprintf(data, "%19.12E\n", VAT3(XuW, ii, jj, kk));
2750
2751 /* *************************************************************
2752 * *** > UN;
2753 * *************************************************************/
2754 // VAT3( XuN, ii,jj,kk) =
2755 TMP1_XUN =
2756 - VAT3( oPW, ii,jj,kk) * (- VAT3( uN, im1,j,k) * VAT3(dPSW, ii,jjp1,kkp1)
2757 - VAT3( uNE, im1,j,k) * VAT3( dPS, ii,jjp1,kkp1))
2758
2759 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,j,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2760 - VAT3( uN, im1,j,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2761 - VAT3( oNE, im1,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2762 - VAT3( uNE, im1,j,kp1) * VAT3( oPS, ii,jjp1,kkp1))
2763
2764 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uC, im1,jp1,k) * VAT3(dPSW, ii,jjp1,kkp1)
2765 - VAT3( uN, im1,jp1,k) * VAT3( dPW, ii,jjp1,kkp1)
2766 - VAT3( uE, im1,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
2767 - VAT3( uNE, im1,jp1,k) * VAT3( dPC, ii,jjp1,kkp1));
2768
2769 TMP2_XUN =
2770 - VAT3(uPNW, ii,jj,kk) * ( VAT3( oC, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2771 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2772 - VAT3( oN, im1,jp1,kp1) * VAT3( dPW, ii,jjp1,kkp1)
2773 - VAT3( uN, im1,jp1,kp1) * VAT3( oPW, ii,jjp1,kkp1)
2774 - VAT3( oE, im1,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2775 - VAT3( uE, im1,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2776 - VAT3( oNE, im1,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
2777 - VAT3( uNE, im1,jp1,kp1) * VAT3( oPC, ii,jjp1,kkp1))
2778
2779 - VAT3( oPC, ii,jj,kk) * (- VAT3( uNW, i,j,k) * VAT3(dPSW, ii,jjp1,kkp1)
2780 - VAT3( uN, i,j,k) * VAT3( dPS, ii,jjp1,kkp1)
2781 - VAT3( uNE, i,j,k) * VAT3(dPSE, ii,jjp1,kkp1));
2782
2783 TMP3_XUN =
2784 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, i,j,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2785 - VAT3( uNW, i,j,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2786 - VAT3( oN, i,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2787 - VAT3( uN, i,j,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2788 - VAT3( oNE, i,j,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2789 - VAT3( uNE, i,j,kp1) * VAT3(oPSE, ii,jjp1,kkp1))
2790
2791 - VAT3( oPN, ii,jj,kk) * (- VAT3( uW, i,jp1,k) * VAT3(dPSW, ii,jjp1,kkp1)
2792 - VAT3( uNW, i,jp1,k) * VAT3( dPW, ii,jjp1,kkp1)
2793 - VAT3( uC, i,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
2794 - VAT3( uN, i,jp1,k) * VAT3( dPC, ii,jjp1,kkp1)
2795 - VAT3( uE, i,jp1,k) * VAT3(dPSE, ii,jjp1,kkp1)
2796 - VAT3( uNE, i,jp1,k) * VAT3( dPE, ii,jjp1,kkp1));
2797
2798 TMP4_XUN =
2799 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPSW, ii,jjp1,kkp1)
2800 - VAT3( uW, i,jp1,kp1) * VAT3(oPSW, ii,jjp1,kkp1)
2801 - VAT3( oNW, i,jp1,kp1) * VAT3( dPW, ii,jjp1,kkp1)
2802 - VAT3( uNW, i,jp1,kp1) * VAT3( oPW, ii,jjp1,kkp1)
2803 + VAT3( oC, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2804 - VAT3( uC, i,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2805 - VAT3( oN, i,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
2806 - VAT3( uN, i,jp1,kp1) * VAT3( oPC, ii,jjp1,kkp1)
2807 - VAT3( oE, i,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2808 - VAT3( uE, i,jp1,kp1) * VAT3(oPSE, ii,jjp1,kkp1)
2809 - VAT3( oNE, i,jp1,kp1) * VAT3( dPE, ii,jjp1,kkp1)
2810 - VAT3( uNE, i,jp1,kp1) * VAT3( oPE, ii,jjp1,kkp1));
2811
2812 TMP5_XUN =
2813 - VAT3( oPE, ii,jj,kk) * (- VAT3( uNW, ip1,j,k) * VAT3( dPS, ii,jjp1,kkp1)
2814 - VAT3( uN, ip1,j,k) * VAT3(dPSE, ii,jjp1,kkp1))
2815
2816 - VAT3( uPE, ii,jj,kk) * (- VAT3( oNW, ip1,j,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2817 - VAT3( uNW, ip1,j,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2818 - VAT3( oN, ip1,j,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2819 - VAT3( uN, ip1,j,kp1) * VAT3(oPSE, ii,jjp1,kkp1))
2820
2821 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uW, ip1,jp1,k) * VAT3( dPS, ii,jjp1,kkp1)
2822 - VAT3( uNW, ip1,jp1,k) * VAT3( dPC, ii,jjp1,kkp1)
2823 - VAT3( uC, ip1,jp1,k) * VAT3(dPSE, ii,jjp1,kkp1)
2824 - VAT3( uN, ip1,jp1,k) * VAT3( dPE, ii,jjp1,kkp1));
2825
2826 TMP6_XUN =
2827 - VAT3(uPNE, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3( dPS, ii,jjp1,kkp1)
2828 - VAT3( uW, ip1,jp1,kp1) * VAT3( oPS, ii,jjp1,kkp1)
2829 - VAT3( oNW, ip1,jp1,kp1) * VAT3( dPC, ii,jjp1,kkp1)
2830 - VAT3( uNW, ip1,jp1,kp1) * VAT3( oPC, ii,jjp1,kkp1)
2831 + VAT3( oC, ip1,jp1,kp1) * VAT3(dPSE, ii,jjp1,kkp1)
2832 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSE, ii,jjp1,kkp1)
2833 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPE, ii,jjp1,kkp1)
2834 - VAT3( uN, ip1,jp1,kp1) * VAT3( oPE, ii,jjp1,kkp1));
2835
2836 VAT3( XuN, ii,jj,kk) = TMP1_XUN + TMP2_XUN + TMP3_XUN + TMP4_XUN
2837 + TMP5_XUN + TMP6_XUN;
2838
2839 //fprintf(data, "%19.12E\n", VAT3(XuN, ii, jj, kk));
2840
2841 /* *************************************************************
2842 * *** > US;
2843 * *************************************************************/
2844
2845 // VAT3( XuS, ii,jj,kk) =
2846 TMP1_XUS =
2847 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uS, im1,jm1,k) * VAT3( dPW, ii,jjm1,kkp1)
2848 - VAT3( uC, im1,jm1,k) * VAT3(dPNW, ii,jjm1,kkp1)
2849 - VAT3( uSE, im1,jm1,k) * VAT3( dPC, ii,jjm1,kkp1)
2850 - VAT3( uE, im1,jm1,k) * VAT3( dPN, ii,jjm1,kkp1))
2851
2852 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oN, im1,jm2,kp1) * VAT3( dPW, ii,jjm1,kkp1)
2853 - VAT3( uS, im1,jm1,kp1) * VAT3( oPW, ii,jjm1,kkp1)
2854 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2855 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2856 - VAT3( oNW, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
2857 - VAT3( uSE, im1,jm1,kp1) * VAT3( oPC, ii,jjm1,kkp1)
2858 - VAT3( oE, im1,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2859 - VAT3( uE, im1,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1));
2860
2861 TMP2_XUS =
2862 - VAT3( oPW, ii,jj,kk) * (- VAT3( uS, im1,j,k) * VAT3(dPNW, ii,jjm1,kkp1)
2863 - VAT3( uSE, im1,j,k) * VAT3( dPN, ii,jjm1,kkp1))
2864
2865 - VAT3( uPW, ii,jj,kk) * (- VAT3( oN, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2866 - VAT3( uS, im1,j,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2867 - VAT3( oNW, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2868 - VAT3( uSE, im1,j,kp1) * VAT3( oPN, ii,jjm1,kkp1))
2869
2870 - VAT3( oPS, ii,jj,kk) * (- VAT3( uSW, i,jm1,k) * VAT3( dPW, ii,jjm1,kkp1)
2871 - VAT3( uW, i,jm1,k) * VAT3(dPNW, ii,jjm1,kkp1)
2872 - VAT3( uS, i,jm1,k) * VAT3( dPC, ii,jjm1,kkp1)
2873 - VAT3( uC, i,jm1,k) * VAT3( dPN, ii,jjm1,kkp1)
2874 - VAT3( uSE, i,jm1,k) * VAT3( dPE, ii,jjm1,kkp1)
2875 - VAT3( uE, i,jm1,k) * VAT3(dPNE, ii,jjm1,kkp1));
2876
2877 TMP3_XUS =
2878 - VAT3( uPS, ii,jj,kk) * (- VAT3( oNE, im1,jm2,kp1) * VAT3( dPW, ii,jjm1,kkp1)
2879 - VAT3( uSW, i,jm1,kp1) * VAT3( oPW, ii,jjm1,kkp1)
2880 - VAT3( oE, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2881 - VAT3( uW, i,jm1,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2882 - VAT3( oN, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
2883 - VAT3( uS, i,jm1,kp1) * VAT3( oPC, ii,jjm1,kkp1)
2884 + VAT3( oC, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2885 - VAT3( uC, i,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2886 - VAT3( oNW, ip1,jm2,kp1) * VAT3( dPE, ii,jjm1,kkp1)
2887 - VAT3( uSE, i,jm1,kp1) * VAT3( oPE, ii,jjm1,kkp1)
2888 - VAT3( oE, i,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2889 - VAT3( uE, i,jm1,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2890
2891 TMP4_XUS =
2892 - VAT3( oPC, ii,jj,kk) * (- VAT3( uSW, i,j,k) * VAT3(dPNW, ii,jjm1,kkp1)
2893 - VAT3( uS, i,j,k) * VAT3( dPN, ii,jjm1,kkp1)
2894 - VAT3( uSE, i,j,k) * VAT3(dPNE, ii,jjm1,kkp1))
2895
2896 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPNW, ii,jjm1,kkp1)
2897 - VAT3( uSW, i,j,kp1) * VAT3(oPNW, ii,jjm1,kkp1)
2898 - VAT3( oN, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2899 - VAT3( uS, i,j,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2900 - VAT3( oNW, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2901 - VAT3( uSE, i,j,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2902
2903 TMP5_XUS =
2904 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uSW, ip1,jm1,k) * VAT3( dPC, ii,jjm1,kkp1)
2905 - VAT3( uW, ip1,jm1,k) * VAT3( dPN, ii,jjm1,kkp1)
2906 - VAT3( uS, ip1,jm1,k) * VAT3( dPE, ii,jjm1,kkp1)
2907 - VAT3( uC, ip1,jm1,k) * VAT3(dPNE, ii,jjm1,kkp1))
2908
2909 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oNE, i,jm2,kp1) * VAT3( dPC, ii,jjm1,kkp1)
2910 - VAT3( uSW, ip1,jm1,kp1) * VAT3( oPC, ii,jjm1,kkp1)
2911 - VAT3( oE, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2912 - VAT3( uW, ip1,jm1,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2913 - VAT3( oN, ip1,jm2,kp1) * VAT3( dPE, ii,jjm1,kkp1)
2914 - VAT3( uS, ip1,jm1,kp1) * VAT3( oPE, ii,jjm1,kkp1)
2915 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2916 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2917
2918 TMP6_XUS =
2919 - VAT3( oPE, ii,jj,kk) * (- VAT3( uSW, ip1,j,k) * VAT3( dPN, ii,jjm1,kkp1)
2920 - VAT3( uS, ip1,j,k) * VAT3(dPNE, ii,jjm1,kkp1))
2921
2922 - VAT3( uPE, ii,jj,kk) * (- VAT3( oNE, i,jm1,kp1) * VAT3( dPN, ii,jjm1,kkp1)
2923 - VAT3( uSW, ip1,j,kp1) * VAT3( oPN, ii,jjm1,kkp1)
2924 - VAT3( oN, ip1,jm1,kp1) * VAT3(dPNE, ii,jjm1,kkp1)
2925 - VAT3( uS, ip1,j,kp1) * VAT3(oPNE, ii,jjm1,kkp1));
2926
2927 VAT3( XuS, ii,jj,kk) = TMP1_XUS + TMP2_XUS + TMP3_XUS + TMP4_XUS
2928 + TMP5_XUS + TMP6_XUS;
2929
2930 //fprintf(data, "%19.12E\n", VAT3(XuS, ii, jj, kk));
2931
2932 /* *************************************************************
2933 * *** > UNE;
2934 * *************************************************************/
2935 // VAT3(XuNE, ii,jj,kk) =
2936 TMP1_XUNE =
2937 VAT3( oPC, ii,jj,kk) * VAT3( uNE, i,j,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2938
2939 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, i,j,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2940 - VAT3( uNE, i,j,kp1) * VAT3(oPSW, iip1,jjp1,kkp1))
2941
2942 - VAT3( oPN, ii,jj,kk) * (- VAT3( uE, i,jp1,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2943 - VAT3( uNE, i,jp1,k) * VAT3( dPW, iip1,jjp1,kkp1))
2944
2945 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, i,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2946 - VAT3( uE, i,jp1,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
2947 - VAT3( oNE, i,jp1,kp1) * VAT3( dPW, iip1,jjp1,kkp1)
2948 - VAT3( uNE, i,jp1,kp1) * VAT3( oPW, iip1,jjp1,kkp1))
2949
2950 - VAT3( oPE, ii,jj,kk) * (- VAT3( uN, ip1,j,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2951 - VAT3( uNE, ip1,j,k) * VAT3( dPS, iip1,jjp1,kkp1));
2952
2953 TMP2_XUNE =
2954 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,j,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2955 - VAT3( uN, ip1,j,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
2956 - VAT3( oNE, ip1,j,kp1) * VAT3( dPS, iip1,jjp1,kkp1)
2957 - VAT3( uNE, ip1,j,kp1) * VAT3( oPS, iip1,jjp1,kkp1))
2958
2959 - VAT3(oPNE, ii,jj,kk) * (- VAT3( uC, ip1,jp1,k) * VAT3(dPSW, iip1,jjp1,kkp1)
2960 - VAT3( uN, ip1,jp1,k) * VAT3( dPW, iip1,jjp1,kkp1)
2961 - VAT3( uE, ip1,jp1,k) * VAT3( dPS, iip1,jjp1,kkp1)
2962 - VAT3( uNE, ip1,jp1,k) * VAT3( dPC, iip1,jjp1,kkp1))
2963
2964 - VAT3(uPNE, ii,jj,kk) * ( VAT3( oC, ip1,jp1,kp1) * VAT3(dPSW, iip1,jjp1,kkp1)
2965 - VAT3( uC, ip1,jp1,kp1) * VAT3(oPSW, iip1,jjp1,kkp1)
2966 - VAT3( oN, ip1,jp1,kp1) * VAT3( dPW, iip1,jjp1,kkp1)
2967 - VAT3( uN, ip1,jp1,kp1) * VAT3( oPW, iip1,jjp1,kkp1)
2968 - VAT3( oE, ip1,jp1,kp1) * VAT3( dPS, iip1,jjp1,kkp1)
2969 - VAT3( uE, ip1,jp1,kp1) * VAT3( oPS, iip1,jjp1,kkp1)
2970 - VAT3( oNE, ip1,jp1,kp1) * VAT3( dPC, iip1,jjp1,kkp1)
2971 - VAT3( uNE, ip1,jp1,kp1) * VAT3( oPC, iip1,jjp1,kkp1));
2972 VAT3(XuNE, ii,jj,kk) = TMP1_XUNE + TMP2_XUNE;
2973
2974 //fprintf(data, "%19.12E\n", VAT3(XuNE, ii, jj, kk));
2975
2976 /* *************************************************************
2977 * *** > UNW;
2978 * *************************************************************/
2979
2980 // VAT3(XuNW, ii,jj,kk) =
2981 TMP1_XUNW =
2982 - VAT3( oPW, ii,jj,kk) * (- VAT3( uNW, im1,j,k) * VAT3( dPS, iim1,jjp1,kkp1)
2983 - VAT3( uN, im1,j,k) * VAT3(dPSE, iim1,jjp1,kkp1))
2984
2985 - VAT3( uPW, ii,jj,kk) * (- VAT3( oNW, im1,j,kp1) * VAT3( dPS, iim1,jjp1,kkp1)
2986 - VAT3( uNW, im1,j,kp1) * VAT3( oPS, iim1,jjp1,kkp1)
2987 - VAT3( oN, im1,j,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
2988 - VAT3( uN, im1,j,kp1) * VAT3(oPSE, iim1,jjp1,kkp1))
2989
2990 - VAT3(oPNW, ii,jj,kk) * (- VAT3( uW, im1,jp1,k) * VAT3( dPS, iim1,jjp1,kkp1)
2991 - VAT3( uNW, im1,jp1,k) * VAT3( dPC, iim1,jjp1,kkp1)
2992 - VAT3( uC, im1,jp1,k) * VAT3(dPSE, iim1,jjp1,kkp1)
2993 - VAT3( uN, im1,jp1,k) * VAT3( dPE, iim1,jjp1,kkp1));
2994
2995 TMP2_XUNW =
2996 - VAT3(uPNW, ii,jj,kk) * (- VAT3( oE, im2,jp1,kp1) * VAT3( dPS, iim1,jjp1,kkp1)
2997 - VAT3( uW, im1,jp1,kp1) * VAT3( oPS, iim1,jjp1,kkp1)
2998 - VAT3( oNW, im1,jp1,kp1) * VAT3( dPC, iim1,jjp1,kkp1)
2999 - VAT3( uNW, im1,jp1,kp1) * VAT3( oPC, iim1,jjp1,kkp1)
3000 + VAT3( oC, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
3001 - VAT3( uC, im1,jp1,kp1) * VAT3(oPSE, iim1,jjp1,kkp1)
3002 - VAT3( oN, im1,jp1,kp1) * VAT3( dPE, iim1,jjp1,kkp1)
3003 - VAT3( uN, im1,jp1,kp1) * VAT3( oPE, iim1,jjp1,kkp1))
3004
3005 + VAT3( oPC, ii,jj,kk) * VAT3( uNW, i,j,k) * VAT3(dPSE, iim1,jjp1,kkp1)
3006
3007 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, i,j,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
3008 - VAT3( uNW, i,j,kp1) * VAT3(oPSE, iim1,jjp1,kkp1))
3009
3010 - VAT3( oPN, ii,jj,kk) * (- VAT3( uW, i,jp1,k) * VAT3(dPSE, iim1,jjp1,kkp1)
3011 - VAT3( uNW, i,jp1,k) * VAT3( dPE, iim1,jjp1,kkp1))
3012
3013 - VAT3( uPN, ii,jj,kk) * (- VAT3( oE, im1,jp1,kp1) * VAT3(dPSE, iim1,jjp1,kkp1)
3014 - VAT3( uW, i,jp1,kp1) * VAT3(oPSE, iim1,jjp1,kkp1)
3015 - VAT3( oNW, i,jp1,kp1) * VAT3( dPE, iim1,jjp1,kkp1)
3016 - VAT3( uNW, i,jp1,kp1) * VAT3( oPE, iim1,jjp1,kkp1));
3017
3018 VAT3(XuNW, ii,jj,kk) = TMP1_XUNW + TMP2_XUNW;
3019
3020 //fprintf(data, "%19.12E\n", VAT3(XuNW, ii, jj, kk));
3021
3022 /* *************************************************************
3023 * *** > USE;
3024 * *************************************************************/
3025
3026 // VAT3(XuSE, ii,jj,kk) =
3027 TMP1_XUSE =
3028 - VAT3( oPS, ii,jj,kk) * (- VAT3( uSE, i,jm1,k) * VAT3( dPW, iip1,jjm1,kkp1)
3029 - VAT3( uE, i,jm1,k) * VAT3(dPNW, iip1,jjm1,kkp1))
3030
3031 - VAT3( uPS, ii,jj,kk) * (- VAT3( oNW, ip1,jm2,kp1) * VAT3( dPW, iip1,jjm1,kkp1)
3032 - VAT3( uSE, i,jm1,kp1) * VAT3( oPW, iip1,jjm1,kkp1)
3033 - VAT3( oE, i,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3034 - VAT3( uE, i,jm1,kp1) * VAT3(oPNW, iip1,jjm1,kkp1))
3035
3036 + VAT3( oPC, ii,jj,kk) * VAT3( uSE, i,j,k) * VAT3(dPNW, iip1,jjm1,kkp1)
3037
3038 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNW, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3039 - VAT3( uSE, i,j,kp1) * VAT3(oPNW, iip1,jjm1,kkp1));
3040
3041 TMP2_XUSE =
3042 - VAT3(oPSE, ii,jj,kk) * (- VAT3( uS, ip1,jm1,k) * VAT3( dPW, iip1,jjm1,kkp1)
3043 - VAT3( uC, ip1,jm1,k) * VAT3(dPNW, iip1,jjm1,kkp1)
3044 - VAT3( uSE, ip1,jm1,k) * VAT3( dPC, iip1,jjm1,kkp1)
3045 - VAT3( uE, ip1,jm1,k) * VAT3( dPN, iip1,jjm1,kkp1))
3046
3047 - VAT3(uPSE, ii,jj,kk) * (- VAT3( oN, ip1,jm2,kp1) * VAT3( dPW, iip1,jjm1,kkp1)
3048 - VAT3( uS, ip1,jm1,kp1) * VAT3( oPW, iip1,jjm1,kkp1)
3049 + VAT3( oC, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3050 - VAT3( uC, ip1,jm1,kp1) * VAT3(oPNW, iip1,jjm1,kkp1)
3051 - VAT3( oNW, ip2,jm2,kp1) * VAT3( dPC, iip1,jjm1,kkp1)
3052 - VAT3( uSE, ip1,jm1,kp1) * VAT3( oPC, iip1,jjm1,kkp1)
3053 - VAT3( oE, ip1,jm1,kp1) * VAT3( dPN, iip1,jjm1,kkp1)
3054 - VAT3( uE, ip1,jm1,kp1) * VAT3( oPN, iip1,jjm1,kkp1))
3055
3056 - VAT3( oPE, ii,jj,kk) * (- VAT3( uS, ip1,j,k) * VAT3(dPNW, iip1,jjm1,kkp1)
3057 - VAT3( uSE, ip1,j,k) * VAT3( dPN, iip1,jjm1,kkp1))
3058
3059 - VAT3( uPE, ii,jj,kk) * (- VAT3( oN, ip1,jm1,kp1) * VAT3(dPNW, iip1,jjm1,kkp1)
3060 - VAT3( uS, ip1,j,kp1) * VAT3(oPNW, iip1,jjm1,kkp1)
3061 - VAT3( oNW, ip2,jm1,kp1) * VAT3( dPN, iip1,jjm1,kkp1)
3062 - VAT3( uSE, ip1,j,kp1) * VAT3( oPN, iip1,jjm1,kkp1));
3063 VAT3(XuSE, ii,jj,kk) = TMP1_XUSE + TMP2_XUSE;
3064
3065 //fprintf(data, "%19.12E\n", VAT3(XuSE, ii, jj, kk));
3066
3067 /* *************************************************************
3068 * *** > USW;
3069 * *************************************************************/
3070
3071 // VAT3(XuSW, ii,jj,kk) =
3072 TMP1_XUSW =
3073 - VAT3(oPSW, ii,jj,kk) * (- VAT3( uSW, im1,jm1,k) * VAT3( dPC, iim1,jjm1,kkp1)
3074 - VAT3( uW, im1,jm1,k) * VAT3( dPN, iim1,jjm1,kkp1)
3075 - VAT3( uS, im1,jm1,k) * VAT3( dPE, iim1,jjm1,kkp1)
3076 - VAT3( uC, im1,jm1,k) * VAT3(dPNE, iim1,jjm1,kkp1))
3077
3078 - VAT3(uPSW, ii,jj,kk) * (- VAT3( oNE, im2,jm2,kp1) * VAT3( dPC, iim1,jjm1,kkp1)
3079 - VAT3( uSW, im1,jm1,kp1) * VAT3( oPC, iim1,jjm1,kkp1)
3080 - VAT3( oE, im2,jm1,kp1) * VAT3( dPN, iim1,jjm1,kkp1)
3081 - VAT3( uW, im1,jm1,kp1) * VAT3( oPN, iim1,jjm1,kkp1)
3082 - VAT3( oN, im1,jm2,kp1) * VAT3( dPE, iim1,jjm1,kkp1)
3083 - VAT3( uS, im1,jm1,kp1) * VAT3( oPE, iim1,jjm1,kkp1)
3084 + VAT3( oC, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3085 - VAT3( uC, im1,jm1,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
3086
3087 - VAT3( oPW, ii,jj,kk) * (- VAT3( uSW, im1,j,k) * VAT3( dPN, iim1,jjm1,kkp1)
3088 - VAT3( uS, im1,j,k) * VAT3(dPNE, iim1,jjm1,kkp1));
3089
3090 TMP2_XUSW =
3091 - VAT3( uPW, ii,jj,kk) * (- VAT3( oNE, im2,jm1,kp1) * VAT3( dPN, iim1,jjm1,kkp1)
3092 - VAT3( uSW, im1,j,kp1) * VAT3( oPN, iim1,jjm1,kkp1)
3093 - VAT3( oN, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3094 - VAT3( uS, im1,j,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
3095
3096 - VAT3( oPS, ii,jj,kk) * (- VAT3( uSW, i,jm1,k) * VAT3( dPE, iim1,jjm1,kkp1)
3097 - VAT3( uW, i,jm1,k) * VAT3(dPNE, iim1,jjm1,kkp1))
3098
3099 - VAT3( uPS, ii,jj,kk) * (- VAT3( oNE, im1,jm2,kp1) * VAT3( dPE, iim1,jjm1,kkp1)
3100 - VAT3( uSW, i,jm1,kp1) * VAT3( oPE, iim1,jjm1,kkp1)
3101 - VAT3( oE, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3102 - VAT3( uW, i,jm1,kp1) * VAT3(oPNE, iim1,jjm1,kkp1))
3103
3104 + VAT3( oPC, ii,jj,kk) * VAT3( uSW, i,j,k) * VAT3(dPNE, iim1,jjm1,kkp1)
3105
3106 - VAT3( uPC, ii,jj,kk) * (- VAT3( oNE, im1,jm1,kp1) * VAT3(dPNE, iim1,jjm1,kkp1)
3107 - VAT3( uSW, i,j,kp1) * VAT3(oPNE, iim1,jjm1,kkp1));
3108 VAT3(XuSW, ii,jj,kk) = TMP1_XUSW + TMP2_XUSW;
3109
3110 //fprintf(data, "%19.12E\n", VAT3(XuSW, ii, jj, kk));
3111
3112 }
3113 }
3114 }
3115}
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.
Definition buildGd.c:450
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.
Definition buildGd.c:1252
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.
Definition buildGd.c:145
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.
Definition buildGd.c:57