APBS  3.0.0
buildGd.c
1 
55 #include "buildGd.h"
56 
57 VPUBLIC 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 
69  VbuildG_1(
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 
92  VbuildG_7(
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 
116  VbuildG_27(
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 
145 VPUBLIC 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 
450 VPUBLIC 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 
1252 VPUBLIC 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