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