1 subroutine evtricub(xget,yget,zget,x,nx,y,ny,z,nz,
3 > f,inf2,inf3,ict,fval,ier)
25 real f(0:7,inf2,inf3,nz)
109 integer i(1),j(1),k(1)
116 real xparam(1),yparam(1),zparam(1)
121 real hx(1),hy(1),hz(1)
122 real hxi(1),hyi(1),hzi(1)
134 call herm3xyz(xget,yget,zget,x,nx,y,ny,z,nz,ilinx,iliny,ilinz,
135 > i(1),j(1),k(1),xparam(1),yparam(1),zparam(1),
136 > hx(1),hxi(1),hy(1),hyi(1),hz(1),hzi(1),ier)
139 call fvtricub(ict,1,1,
140 > fval,i,j,k,xparam,yparam,zparam,
141 > hx,hxi,hy,hyi,hz,hzi,
150 subroutine fvtricub(ict,ivec,ivecd,
151 > fval,ii,jj,kk,xparam,yparam,zparam,
152 > hx,hxi,hy,hyi,hz,hzi,
161 integer ii(ivec),jj(ivec),kk(ivec)
162 real xparam(ivec),yparam(ivec),zparam(ivec)
165 real hx(ivec),hy(ivec),hz(ivec)
166 real hxi(ivec),hyi(ivec),hzi(ivec)
169 real fin(0:7,inf2,inf3,nz)
189 real,
parameter :: sixth = 0.166666666666666667
194 z216th=sixth*sixth*sixth
197 if(abs(ict(1)).le.2)
then
248 > xpi*(ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k))+
249 > xp*(ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
251 > xpi*(ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1))+
252 > xp*(ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
256 > cxi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
257 > cx*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
259 > cxi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
260 > cx*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
264 > xpi*(cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k))+
265 > xp*(cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
267 > xpi*(cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1))+
268 > xp*(cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
272 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
273 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
275 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
276 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
278 sum=sum+z36th*hx2*hy2*(
280 > cxi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
281 > cx*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
283 > cxi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
284 > cx*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
286 sum=sum+z36th*hx2*hz2*(
288 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
289 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
291 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
292 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
294 sum=sum+z36th*hy2*hz2*(
296 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
297 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
299 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
300 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
302 sum=sum+z216th*hx2*hy2*hz2*(
304 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
305 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
307 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
308 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
358 > -(ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k))
359 > +(ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
361 > -(ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1))
362 > +(ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
364 sum=sum+sixth*hx(v)*(
366 > cxdi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
367 > cxd*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
369 > cxdi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
370 > cxd*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
372 sum=sum+sixth*hxi(v)*hy2*(
374 > -(cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k))
375 > +(cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
377 > -(cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1))
378 > +(cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
380 sum=sum+sixth*hxi(v)*hz2*(
382 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
383 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
385 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
386 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
388 sum=sum+z36th*hx(v)*hy2*(
390 > cxdi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
391 > cxd*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
393 > cxdi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
394 > cxd*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
396 sum=sum+z36th*hx(v)*hz2*(
398 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
399 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
401 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
402 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
404 sum=sum+z36th*hxi(v)*hy2*hz2*(
406 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
407 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
409 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
410 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
412 sum=sum+z216th*hx(v)*hy2*hz2*(
414 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
415 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
417 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
418 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
468 > xpi*(-fin(0,i,j,k) +fin(0,i,j+1,k))+
469 > xp*(-fin(0,i+1,j,k)+fin(0,i+1,j+1,k)))
471 > xpi*(-fin(0,i,j,k+1) +fin(0,i,j+1,k+1))+
472 > xp*(-fin(0,i+1,j,k+1)+fin(0,i+1,j+1,k+1))))
474 sum=sum+sixth*hyi(v)*hx2*(
476 > cxi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
477 > cx*(-fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
479 > cxi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
480 > cx*(-fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
482 sum=sum+sixth*hy(v)*(
484 > xpi*(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))+
485 > xp*(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
487 > xpi*(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))+
488 > xp*(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
490 sum=sum+sixth*hyi(v)*hz2*(
492 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
493 > xp*(-fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
495 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
496 > xp*(-fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
498 sum=sum+z36th*hx2*hy(v)*(
500 > cxi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
501 > cx*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
503 > cxi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
504 > cx*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
506 sum=sum+z36th*hyi(v)*hx2*hz2*(
508 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
509 > cx*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
511 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
512 > cx*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
514 sum=sum+z36th*hy(v)*hz2*(
516 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
517 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
519 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
520 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
522 sum=sum+z216th*hx2*hy(v)*hz2*(
524 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
525 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
527 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
528 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
578 > xpi*(ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k))+
579 > xp*(ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
581 > xpi*(ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1))+
582 > xp*(ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
584 sum=sum+sixth*hx2*hzi(v)*(
586 > cxi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
587 > cx*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
589 > cxi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
590 > cx*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
592 sum=sum+sixth*hy2*hzi(v)*(
594 > xpi*(cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k))+
595 > xp*(cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
597 > xpi*(cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1))+
598 > xp*(cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
600 sum=sum+sixth*hz(v)*(
602 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
603 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
605 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
606 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
608 sum=sum+z36th*hx2*hy2*hzi(v)*(
610 > cxi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
611 > cx*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
613 > cxi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
614 > cx*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
616 sum=sum+z36th*hx2*hz(v)*(
618 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
619 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
621 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
622 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
624 sum=sum+z36th*hy2*hz(v)*(
626 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
627 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
629 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
630 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
632 sum=sum+z216th*hx2*hy2*hz(v)*(
634 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
635 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
637 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
638 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
683 > xpi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
684 > xp*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
686 > xpi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
687 > xp*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
691 > xpi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
692 > xp*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
694 > xpi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
695 > xp*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
699 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
700 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
702 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
703 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
705 sum=sum+z36th*hy2*hz2*(
707 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
708 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
710 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
711 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
756 > xpi*(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))+
757 > xp*(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
759 > xpi*(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))+
760 > xp*(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
764 > cxi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
765 > cx*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
767 > cxi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
768 > cx*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
772 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
773 > xp*(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
775 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
776 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
778 sum=sum+z36th*hx2*hz2*(
780 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
781 > cx*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
783 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
784 > cx*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
829 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
830 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
832 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
833 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
837 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
838 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
840 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
841 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
845 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
846 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
848 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
849 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
851 sum=sum+z36th*hx2*hy2*(
853 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
854 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
856 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
857 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
906 > (fin(0,i,j,k) -fin(0,i,j+1,k))-
907 > (fin(0,i+1,j,k)-fin(0,i+1,j+1,k)))
909 > (fin(0,i,j,k+1) -fin(0,i,j+1,k+1))-
910 > (fin(0,i+1,j,k+1)-fin(0,i+1,j+1,k+1))))
912 sum=sum+sixth*hyi(v)*hx(v)*(
914 > cxdi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
915 > cxd*(-fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
917 > cxdi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
918 > cxd*(-fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
920 sum=sum+sixth*hxi(v)*hy(v)*(
922 > -(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))
923 > +(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
925 > -(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))
926 > +(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
928 sum=sum+sixth*hxi(v)*hyi(v)*hz2*(
930 > (fin(3,i,j,k) -fin(3,i,j+1,k))-
931 > (fin(3,i+1,j,k)-fin(3,i+1,j+1,k)))
933 > (fin(3,i,j,k+1) -fin(3,i,j+1,k+1))-
934 > (fin(3,i+1,j,k+1)-fin(3,i+1,j+1,k+1))))
936 sum=sum+z36th*hx(v)*hy(v)*(
938 > cxdi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
939 > cxd*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
941 > cxdi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
942 > cxd*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
944 sum=sum+z36th*hyi(v)*hx(v)*hz2*(
946 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
947 > cxd*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
949 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
950 > cxd*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
952 sum=sum+z36th*hxi(v)*hy(v)*hz2*(
954 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
955 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
957 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
958 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
960 sum=sum+z216th*hx(v)*hy(v)*hz2*(
962 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
963 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
965 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
966 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
1015 > (ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k)) -
1016 > (ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
1018 > (ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1)) -
1019 > (ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
1021 sum=sum+sixth*hx(v)*hzi(v)*(
1023 > cxdi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
1024 > cxd*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
1026 > cxdi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
1027 > cxd*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
1029 sum=sum+sixth*hxi(v)*hy2*hzi(v)*(
1031 > (cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k)) -
1032 > (cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
1034 > (cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1)) -
1035 > (cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
1037 sum=sum+sixth*hxi(v)*hz(v)*(
1039 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
1040 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
1042 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
1043 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
1045 sum=sum+z36th*hx(v)*hy2*hzi(v)*(
1047 > cxdi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
1048 > cxd*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
1050 > cxdi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
1051 > cxd*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
1053 sum=sum+z36th*hx(v)*hz(v)*(
1055 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
1056 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
1058 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
1059 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
1061 sum=sum+z36th*hxi(v)*hy2*hz(v)*(
1063 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
1064 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
1066 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
1067 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
1069 sum=sum+z216th*hx(v)*hy2*hz(v)*(
1071 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
1072 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
1074 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
1075 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
1081 if(ict(10).eq.1)
then
1124 > xpi*(fin(0,i,j,k) -fin(0,i,j+1,k))+
1125 > xp*(fin(0,i+1,j,k)-fin(0,i+1,j+1,k)))
1127 > xpi*(fin(0,i,j,k+1) -fin(0,i,j+1,k+1))+
1128 > xp*(fin(0,i+1,j,k+1)-fin(0,i+1,j+1,k+1))))
1130 sum=sum+sixth*hyi(v)*hx2*hzi(v)*(
1132 > cxi*(fin(1,i,j,k) -fin(1,i,j+1,k))+
1133 > cx*(fin(1,i+1,j,k)-fin(1,i+1,j+1,k)))
1135 > cxi*(fin(1,i,j,k+1) -fin(1,i,j+1,k+1))+
1136 > cx*(fin(1,i+1,j,k+1)-fin(1,i+1,j+1,k+1))))
1138 sum=sum+sixth*hy(v)*hzi(v)*(
1140 > xpi*(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))+
1141 > xp*(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
1143 > xpi*(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))+
1144 > xp*(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
1146 sum=sum+sixth*hyi(v)*hz(v)*(
1148 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
1149 > xp*(-fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
1151 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
1152 > xp*(-fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
1154 sum=sum+z36th*hx2*hy(v)*hzi(v)*(
1156 > cxi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
1157 > cx*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
1159 > cxi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
1160 > cx*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
1162 sum=sum+z36th*hyi(v)*hx2*hz(v)*(
1164 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
1165 > cx*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
1167 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
1168 > cx*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
1170 sum=sum+z36th*hy(v)*hz(v)*(
1172 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
1173 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
1175 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
1176 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
1178 sum=sum+z216th*hx2*hy(v)*hz(v)*(
1180 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
1181 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
1183 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
1184 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
1193 else if(ict(1).eq.3)
then
1194 if(ict(2).eq.1)
then
1230 > xpi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
1231 > xp*( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
1233 > xpi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
1234 > xp*( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
1236 sum=sum+sixth*hy(v)*(
1238 > xpi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
1239 > xp*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
1241 > xpi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
1242 > xp*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
1244 sum=sum+sixth*hz2*hyi(v)*(
1246 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
1247 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
1249 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
1250 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
1252 sum=sum+z36th*hy(v)*hz2*(
1254 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
1255 > xp*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
1257 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
1258 > xp*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
1264 if(ict(3).eq.1)
then
1300 > xpi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
1301 > xp*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
1303 > xpi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
1304 > xp*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
1306 sum=sum+sixth*hy2*hzi(v)*(
1308 > xpi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
1309 > xp*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
1311 > xpi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
1312 > xp*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
1314 sum=sum+sixth*hz(v)*(
1316 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
1317 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
1319 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
1320 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
1322 sum=sum+z36th*hy2*hz(v)*(
1324 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
1325 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
1327 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
1328 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
1334 if(ict(4).eq.1)
then
1370 > -(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))
1371 > +(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
1373 > -(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))
1374 > +(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
1376 sum=sum+sixth*hx(v)*(
1378 > cxdi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
1379 > cxd*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
1381 > cxdi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
1382 > cxd*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
1384 sum=sum+sixth*hz2*hxi(v)*(
1386 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
1387 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
1389 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
1390 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
1392 sum=sum+z36th*hx(v)*hz2*(
1394 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
1395 > cxd*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
1397 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
1398 > cxd*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
1404 if(ict(5).eq.1)
then
1442 sum=hxi(v)*hyi(v)*hzi(v)*(
1444 > (fin(0,i,j,k) -fin(0,i,j+1,k))-
1445 > (fin(0,i+1,j,k)-fin(0,i+1,j+1,k)))
1447 > (fin(0,i,j,k+1) -fin(0,i,j+1,k+1))-
1448 > (fin(0,i+1,j,k+1)-fin(0,i+1,j+1,k+1))))
1450 sum=sum+sixth*hyi(v)*hx(v)*hzi(v)*(
1452 > cxdi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
1453 > cxd*(-fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
1455 > cxdi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
1456 > cxd*(-fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
1458 sum=sum+sixth*hxi(v)*hy(v)*hzi(v)*(
1460 > -(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))
1461 > +(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
1463 > -(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))
1464 > +(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
1466 sum=sum+sixth*hxi(v)*hyi(v)*hz(v)*(
1468 > (fin(3,i,j,k) -fin(3,i,j+1,k))-
1469 > (fin(3,i+1,j,k)-fin(3,i+1,j+1,k)))
1471 > (fin(3,i,j,k+1) -fin(3,i,j+1,k+1))-
1472 > (fin(3,i+1,j,k+1)-fin(3,i+1,j+1,k+1))))
1474 sum=sum+z36th*hx(v)*hy(v)*hzi(v)*(
1476 > cxdi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
1477 > cxd*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
1479 > cxdi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
1480 > cxd*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
1482 sum=sum+z36th*hyi(v)*hx(v)*hz(v)*(
1484 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
1485 > cxd*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
1487 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
1488 > cxd*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
1490 sum=sum+z36th*hxi(v)*hy(v)*hz(v)*(
1492 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
1493 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
1495 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
1496 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
1498 sum=sum+z216th*hx(v)*hy(v)*hz(v)*(
1500 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
1501 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
1503 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
1504 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
1510 if(ict(6).eq.1)
then
1547 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
1548 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
1550 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
1551 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
1553 sum=sum+sixth*hx(v)*(
1555 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
1556 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
1558 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
1559 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
1561 sum=sum+sixth*hy2*hxi(v)*(
1563 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
1564 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
1566 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
1567 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
1569 sum=sum+z36th*hx(v)*hy2*(
1571 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
1572 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
1574 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
1575 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
1581 if(ict(7).eq.1)
then
1617 > xpi*(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))+
1618 > xp*(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
1620 > xpi*(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))+
1621 > xp*(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
1623 sum=sum+sixth*hx2*hzi(v)*(
1625 > cxi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
1626 > cx*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
1628 > cxi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
1629 > cx*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
1631 sum=sum+sixth*hz(v)*(
1633 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
1634 > xp*(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
1636 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
1637 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
1639 sum=sum+z36th*hx2*hz(v)*(
1641 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
1642 > cx*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
1644 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
1645 > cx*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
1651 if(ict(8).eq.1)
then
1687 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
1688 > xp*( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
1690 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
1691 > xp*( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
1693 sum=sum+sixth*hx2*hyi(v)*(
1695 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
1696 > cx*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
1698 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
1699 > cx*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
1701 sum=sum+sixth*hy(v)*(
1703 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
1704 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
1706 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
1707 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
1709 sum=sum+z36th*hx2*hy(v)*(
1711 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
1712 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
1714 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
1715 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
1723 else if(ict(1).eq.-3)
then
1724 if(ict(2).eq.1)
then
1756 > -(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))
1757 > +(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
1759 > -(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))
1760 > +(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
1762 sum=sum+sixth*hy2*hxi(v)*(
1764 > -(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))
1765 > +(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
1767 > -(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))
1768 > +(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
1770 sum=sum+sixth*hz2*hxi(v)*(
1772 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
1773 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
1775 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
1776 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
1778 sum=sum+z36th*hy2*hz2*hxi(v)*(
1780 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
1781 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
1783 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
1784 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
1790 if(ict(3).eq.1)
then
1822 > xpi*(-fin(2,i,j,k) +fin(2,i,j+1,k))+
1823 > xp*( -fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
1825 > xpi*(-fin(2,i,j,k+1) +fin(2,i,j+1,k+1))+
1826 > xp*( -fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
1828 sum=sum+sixth*hx2*hyi(v)*(
1830 > cxi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
1831 > cx*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
1833 > cxi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
1834 > cx*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
1836 sum=sum+sixth*hz2*hyi(v)*(
1838 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
1839 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
1841 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
1842 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
1844 sum=sum+z36th*hx2*hz2*hyi(v)*(
1846 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
1847 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
1849 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
1850 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
1856 if(ict(4).eq.1)
then
1888 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
1889 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
1891 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
1892 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
1894 sum=sum+sixth*hx2*hzi(v)*(
1896 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
1897 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
1899 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
1900 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
1902 sum=sum+sixth*hy2*hzi(v)*(
1904 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
1905 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
1907 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
1908 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
1910 sum=sum+z36th*hx2*hy2*hzi(v)*(
1912 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
1913 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
1915 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
1916 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
1925 else if(ict(1).eq.4)
then
1926 if(ict(2).eq.1)
then
1957 > xpi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
1958 > xp*( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
1960 > xpi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
1961 > xp*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
1965 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
1966 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
1968 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
1969 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
1975 if(ict(3).eq.1)
then
2010 > xpi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
2011 > xp*( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
2013 > xpi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
2014 > xp*( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
2016 sum=sum+sixth*hy(v)*hzi(v)*(
2018 > xpi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
2019 > xp*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
2021 > xpi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
2022 > xp*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
2024 sum=sum+sixth*hz(v)*hyi(v)*(
2026 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
2027 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
2029 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
2030 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
2032 sum=sum+z36th*hy(v)*hz(v)*(
2034 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
2035 > xp*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
2037 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
2038 > xp*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
2044 if(ict(4).eq.1)
then
2075 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
2076 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
2078 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
2079 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
2083 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
2084 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
2086 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
2087 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
2093 if(ict(5).eq.1)
then
2128 > -(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))
2129 > +(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
2131 > -(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))
2132 > +(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
2134 sum=sum+sixth*hx(v)*hzi(v)*(
2136 > cxdi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
2137 > cxd*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
2139 > cxdi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
2140 > cxd*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
2142 sum=sum+sixth*hz(v)*hxi(v)*(
2144 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
2145 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
2147 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
2148 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
2150 sum=sum+z36th*hx(v)*hz(v)*(
2152 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
2153 > cxd*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
2155 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
2156 > cxd*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
2162 if(ict(6).eq.1)
then
2197 > ( +fin(3,i,j,k) -fin(3,i,j+1,k))
2198 > +( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
2200 > ( +fin(3,i,j,k+1) -fin(3,i,j+1,k+1))
2201 > +( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
2203 sum=sum+sixth*hx(v)*hyi(v)*(
2205 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
2206 > cxd*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
2208 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
2209 > cxd*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
2211 sum=sum+sixth*hy(v)*hxi(v)*(
2213 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
2214 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
2216 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
2217 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
2219 sum=sum+z36th*hx(v)*hy(v)*(
2221 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
2222 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
2224 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
2225 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
2231 if(ict(7).eq.1)
then
2262 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
2263 > xp*( ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
2265 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
2266 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
2270 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
2271 > cx*( ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
2273 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
2274 > cx*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
2283 else if(ict(1).eq.-4)
then
2284 if(ict(2).eq.1)
then
2315 > ( fin(1,i,j,k) -fin(1,i,j+1,k))+
2316 > ( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
2318 > ( fin(1,i,j,k+1) -fin(1,i,j+1,k+1))+
2319 > ( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
2321 sum=sum+sixth*hy(v)*hxi(v)*(
2323 > -(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
2324 > (cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
2326 > -(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
2327 > (cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
2329 sum=sum+sixth*hz2*hyi(v)*hxi(v)*(
2331 > ( fin(5,i,j,k) -fin(5,i,j+1,k))+
2332 > ( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
2334 > ( fin(5,i,j,k+1) -fin(5,i,j+1,k+1))+
2335 > ( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
2337 sum=sum+z36th*hy(v)*hz2*hxi(v)*(
2339 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
2340 > (cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
2342 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
2343 > (cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
2349 if(ict(3).eq.1)
then
2380 > +(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))
2381 > -(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
2383 > -(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))
2384 > +(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
2386 sum=sum+sixth*hy2*hzi(v)*hxi(v)*(
2388 > +(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))
2389 > -(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
2391 > -(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))
2392 > +(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
2394 sum=sum+sixth*hz(v)*hxi(v)*(
2396 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
2397 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
2399 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
2400 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
2402 sum=sum+z36th*hy2*hz(v)*hxi(v)*(
2404 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
2405 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
2407 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
2408 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
2414 if(ict(4).eq.1)
then
2445 > ( fin(2,i,j,k) -fin(2,i,j+1,k))
2446 > +(-fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
2448 > ( fin(2,i,j,k+1) -fin(2,i,j+1,k+1))
2449 > +(-fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
2451 sum=sum+sixth*hx(v)*hyi(v)*(
2453 > cxdi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
2454 > cxd*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
2456 > cxdi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
2457 > cxd*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
2459 sum=sum+sixth*hz2*hxi(v)*hyi(v)*(
2461 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
2462 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
2464 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
2465 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
2467 sum=sum+z36th*hx(v)*hz2*hyi(v)*(
2469 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
2470 > cxd*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
2472 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
2473 > cxd*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
2479 if(ict(5).eq.1)
then
2510 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
2511 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
2513 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
2514 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
2516 sum=sum+sixth*hx(v)*hzi(v)*(
2518 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
2519 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
2521 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
2522 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
2524 sum=sum+sixth*hy2*hxi(v)*hzi(v)*(
2526 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
2527 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
2529 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
2530 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
2532 sum=sum+z36th*hx(v)*hy2*hzi(v)*(
2534 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
2535 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
2537 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
2538 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
2544 if(ict(6).eq.1)
then
2575 > xpi*(-fin(2,i,j,k) +fin(2,i,j+1,k))+
2576 > xp*( -fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
2578 > xpi*(-fin(2,i,j,k+1) +fin(2,i,j+1,k+1))+
2579 > xp*( -fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
2581 sum=sum+sixth*hx2*hyi(v)*hzi(v)*(
2583 > cxi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
2584 > cx*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
2586 > cxi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
2587 > cx*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
2589 sum=sum+sixth*hz(v)*hyi(v)*(
2591 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
2592 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
2594 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
2595 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
2597 sum=sum+z36th*hx2*hz(v)*hyi(v)*(
2599 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
2600 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
2602 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
2603 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
2609 if(ict(7).eq.1)
then
2640 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
2641 > xp*( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
2643 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
2644 > xp*( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
2646 sum=sum+sixth*hx2*hyi(v)*hzi(v)*(
2648 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
2649 > cx*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
2651 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
2652 > cx*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
2654 sum=sum+sixth*hy(v)*hzi(v)*(
2656 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
2657 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
2659 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
2660 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
2662 sum=sum+z36th*hx2*hy(v)*hzi(v)*(
2664 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
2665 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
2667 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
2668 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
2677 else if(ict(1).eq.5)
then
2678 if(ict(2).eq.1)
then
2708 > xpi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
2709 > xp*( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
2711 > xpi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
2712 > xp*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
2714 sum=sum+sixth*hz(v)*(
2716 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
2717 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
2719 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
2720 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
2726 if(ict(3).eq.1)
then
2756 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
2757 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
2759 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
2760 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
2762 sum=sum+sixth*hy(v)*(
2764 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
2765 > xp*( cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
2767 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
2768 > xp*( cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
2774 if(ict(4).eq.1)
then
2804 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
2805 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
2807 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
2808 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
2810 sum=sum+sixth*hx(v)*(
2812 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
2813 > cxd*(ypi*fin(7,i+1,j,k) +yp*fin(7,i+1,j+1,k)))
2815 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
2816 > cxd*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
2825 else if(ict(1).eq.-5)
then
2826 if(ict(2).eq.1)
then
2852 > -(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))
2853 > +( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
2855 > -(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))
2856 > +(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
2858 sum=sum+sixth*hz2*hxi(v)*(
2860 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
2861 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
2863 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
2864 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
2870 if(ict(3).eq.1)
then
2898 sum=hyi(v)*hzi(v)*hxi(v)*(
2900 > -(-fin(1,i,j,k) +fin(1,i,j+1,k))
2901 > +( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
2903 > -(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))
2904 > +( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
2906 sum=sum+sixth*hy(v)*hzi(v)*hxi(v)*(
2908 > -(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))
2909 > +(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
2911 > -(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))
2912 > +(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
2914 sum=sum+sixth*hz(v)*hyi(v)*hxi(v)*(
2916 > -(-fin(5,i,j,k) +fin(5,i,j+1,k))
2917 > +( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
2919 > -(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))
2920 > +( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
2922 sum=sum+z36th*hy(v)*hz(v)*hxi(v)*(
2924 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))
2925 > +(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
2927 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))
2928 > +(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
2934 if(ict(4).eq.1)
then
2960 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
2961 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
2963 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
2964 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
2966 sum=sum+sixth*hy2*hxi(v)*(
2968 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
2969 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
2971 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
2972 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
2978 if(ict(5).eq.1)
then
3005 > xpi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
3006 > xp*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
3008 > xpi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
3009 > xp*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
3011 sum=sum+sixth*hz2*hyi(v)*(
3013 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3014 > xp*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3016 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3017 > xp*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3023 if(ict(6).eq.1)
then
3049 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
3050 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
3052 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
3053 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
3055 sum=sum+sixth*hy2*hzi(v)*(
3057 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
3058 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
3060 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
3061 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
3067 if(ict(7).eq.1)
then
3095 sum=hxi(v)*hzi(v)*hyi(v)*(
3097 > ( fin(2,i,j,k) -fin(2,i,j+1,k))
3098 > +(-fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
3100 > ( fin(2,i,j,k+1) -fin(2,i,j+1,k+1))
3101 > +(-fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
3103 sum=sum+sixth*hx(v)*hzi(v)*hyi(v)*(
3105 > cxdi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
3106 > cxd*(-fin(4,i+1,j,k) +fin(4,i+1,j+1,k)))
3108 > cxdi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
3109 > cxd*(-fin(4,i+1,j,k+1) +fin(4,i+1,j+1,k+1))))
3111 sum=sum+sixth*hz(v)*hxi(v)*hyi(v)*(
3113 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
3114 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
3116 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
3117 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
3119 sum=sum+z36th*hx(v)*hz(v)*hyi(v)*(
3121 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3122 > cxd*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
3124 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3125 > cxd*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
3131 if(ict(8).eq.1)
then
3159 sum=hyi(v)*hxi(v)*hzi(v)*(
3161 > ( +fin(3,i,j,k) -fin(3,i,j+1,k))
3162 > +( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
3164 > ( +fin(3,i,j,k+1) -fin(3,i,j+1,k+1))
3165 > +( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
3167 sum=sum+sixth*hx(v)*hyi(v)*hzi(v)*(
3169 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
3170 > cxd*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
3172 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
3173 > cxd*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
3175 sum=sum+sixth*hy(v)*hxi(v)*hzi(v)*(
3177 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
3178 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
3180 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
3181 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
3183 sum=sum+z36th*hx(v)*hy(v)*hzi(v)*(
3185 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
3186 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
3188 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
3189 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
3195 if(ict(9).eq.1)
then
3221 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
3222 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
3224 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
3225 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
3227 sum=sum+sixth*hx2*hyi(v)*(
3229 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3230 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3232 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3233 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3239 if(ict(10).eq.1)
then
3265 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
3266 > xp*( ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
3268 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
3269 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
3271 sum=sum+sixth*hx2*hzi(v)*(
3273 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
3274 > cx*( ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
3276 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
3277 > cx*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
3286 else if(ict(1).eq.6)
then
3311 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
3312 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
3314 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
3315 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
3324 if(ict(1).eq.-6)
then
3325 if(ict(2).eq.1)
then
3346 > ( fin(4,i,j,k) -fin(4,i,j+1,k))
3347 > +(-fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
3349 > ( fin(4,i,j,k+1) -fin(4,i,j+1,k+1))
3350 > +(-fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
3352 sum=sum+sixth*hz2*hyi(v)*hxi(v)*(
3354 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
3355 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3357 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
3358 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3364 if(ict(3).eq.1)
then
3389 > -(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))
3390 > +( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
3392 > -(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))
3393 > +(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
3395 sum=sum+sixth*hz(v)*hxi(v)*(
3397 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
3398 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
3400 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
3401 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
3407 if(ict(4).eq.1)
then
3432 > ( fin(5,i,j,k) -fin(5,i,j+1,k))
3433 > +(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
3435 > ( fin(5,i,j,k+1) -fin(5,i,j+1,k+1))
3436 > +(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
3438 sum=sum+sixth*hy(v)*hxi(v)*(
3440 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))
3441 > +(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
3443 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))
3444 > +(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
3450 if(ict(5).eq.1)
then
3471 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
3472 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
3474 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
3475 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
3477 sum=sum+sixth*hy2*hxi(v)*hzi(v)*(
3479 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
3480 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
3482 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
3483 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
3489 if(ict(6).eq.1)
then
3514 > xpi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
3515 > xp*(-fin(4,i+1,j,k) +fin(4,i+1,j+1,k)))
3517 > xpi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
3518 > xp*(-fin(4,i+1,j,k+1) +fin(4,i+1,j+1,k+1))))
3520 sum=sum+sixth*hz(v)*hyi(v)*(
3522 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3523 > xp*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
3525 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3526 > xp*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
3532 if(ict(7).eq.1)
then
3557 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
3558 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
3560 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
3561 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
3563 sum=sum+sixth*hy(v)*hzi(v)*(
3565 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
3566 > xp*( cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
3568 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
3569 > xp*( cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
3575 if(ict(8).eq.1)
then
3600 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
3601 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
3603 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
3604 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
3606 sum=sum+sixth*hx(v)*hyi(v)*(
3608 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3609 > cxd*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
3611 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3612 > cxd*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
3618 if(ict(9).eq.1)
then
3643 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
3644 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
3646 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
3647 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
3649 sum=sum+sixth*hx(v)*hzi(v)*(
3651 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
3652 > cxd*(ypi*fin(7,i+1,j,k) +yp*fin(7,i+1,j+1,k)))
3654 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
3655 > cxd*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
3661 if(ict(10).eq.1)
then
3682 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
3683 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
3685 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
3686 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
3688 sum=sum+sixth*hx2*hyi(v)*hzi(v)*(
3690 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3691 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3693 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3694 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3703 else if(abs(ict(1)).eq.7)
then
3704 if(ict(2).eq.1)
then
3722 sum=hyi(v)*hxi(v)*hzi(v)*(
3724 > ( fin(4,i,j,k) -fin(4,i,j+1,k))
3725 > +(-fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
3727 > ( fin(4,i,j,k+1) -fin(4,i,j+1,k+1))
3728 > +(-fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
3730 sum=sum+sixth*hz(v)*hyi(v)*hxi(v)*(
3732 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
3733 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3735 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
3736 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3742 if(ict(3).eq.1)
then
3762 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
3763 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
3765 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
3766 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
3772 if(ict(4).eq.1)
then
3790 sum=hxi(v)*hyi(v)*hzi(v)*(
3792 > ( fin(5,i,j,k) -fin(5,i,j+1,k))
3793 > +(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
3795 > ( fin(5,i,j,k+1) -fin(5,i,j+1,k+1))
3796 > +(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
3798 sum=sum+sixth*hy(v)*hxi(v)*hzi(v)*(
3800 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))
3801 > +(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
3803 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))
3804 > +(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
3810 if(ict(5).eq.1)
then
3830 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3831 > xp*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3833 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3834 > xp*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3840 if(ict(6).eq.1)
then
3860 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
3861 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
3863 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
3864 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
3870 if(ict(7).eq.1)
then
3888 sum=hxi(v)*hyi(v)*hzi(v)*(
3890 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
3891 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
3893 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
3894 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
3896 sum=sum+sixth*hx(v)*hyi(v)*hzi(v)*(
3898 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3899 > cxd*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
3901 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3902 > cxd*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
3911 else if(abs(ict(1)).eq.8)
then
3912 if(ict(2).eq.1)
then
3927 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
3928 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3930 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
3931 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3937 if(ict(3).eq.1)
then
3952 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
3953 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
3955 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
3956 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
3962 if(ict(4).eq.1)
then
3977 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
3978 > xp*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
3980 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
3981 > xp*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
3990 else if(abs(ict(1)).eq.9)
then
3998 sum=hyi(v)*hxi(v)*hzi(v)*(
4000 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
4001 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
4003 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
4004 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
4017 subroutine fvtricubx(ict,ivec,ivecd,
4018 > fval,ii,jj,kk,xparam,yparam,zparam,
4019 > hx,hxi,hy,hyi,hz,hzi,
4028 integer ii(ivec),jj,kk
4029 real xparam(ivec),yparam,zparam
4033 real hxi(ivec),hyi,hzi
4036 real fin(0:7,inf2,inf3,nz)
4056 real,
parameter :: sixth = 0.166666666666666667
4061 z216th=sixth*sixth*sixth
4064 if(abs(ict(1)).le.2)
then
4070 if(ict(1).eq.1)
then
4116 > xpi*(ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k))+
4117 > xp*(ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
4119 > xpi*(ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1))+
4120 > xp*(ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
4124 > cxi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
4125 > cx*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
4127 > cxi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
4128 > cx*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
4132 > xpi*(cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k))+
4133 > xp*(cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
4135 > xpi*(cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1))+
4136 > xp*(cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
4140 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
4141 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
4143 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
4144 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
4146 sum=sum+z36th*hx2*hy2*(
4148 > cxi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
4149 > cx*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
4151 > cxi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
4152 > cx*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
4154 sum=sum+z36th*hx2*hz2*(
4156 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
4157 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
4159 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
4160 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
4162 sum=sum+z36th*hy2*hz2*(
4164 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
4165 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
4167 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
4168 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
4170 sum=sum+z216th*hx2*hy2*hz2*(
4172 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
4173 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
4175 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
4176 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
4182 if(ict(2).eq.1)
then
4227 > -(ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k))
4228 > +(ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
4230 > -(ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1))
4231 > +(ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
4233 sum=sum+sixth*hx(v)*(
4235 > cxdi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
4236 > cxd*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
4238 > cxdi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
4239 > cxd*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
4241 sum=sum+sixth*hxi(v)*hy2*(
4243 > -(cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k))
4244 > +(cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
4246 > -(cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1))
4247 > +(cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
4249 sum=sum+sixth*hxi(v)*hz2*(
4251 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
4252 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
4254 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
4255 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
4257 sum=sum+z36th*hx(v)*hy2*(
4259 > cxdi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
4260 > cxd*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
4262 > cxdi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
4263 > cxd*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
4265 sum=sum+z36th*hx(v)*hz2*(
4267 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
4268 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
4270 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
4271 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
4273 sum=sum+z36th*hxi(v)*hy2*hz2*(
4275 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
4276 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
4278 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
4279 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
4281 sum=sum+z216th*hx(v)*hy2*hz2*(
4283 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
4284 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
4286 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
4287 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
4293 if(ict(3).eq.1)
then
4338 > xpi*(-fin(0,i,j,k) +fin(0,i,j+1,k))+
4339 > xp*(-fin(0,i+1,j,k)+fin(0,i+1,j+1,k)))
4341 > xpi*(-fin(0,i,j,k+1) +fin(0,i,j+1,k+1))+
4342 > xp*(-fin(0,i+1,j,k+1)+fin(0,i+1,j+1,k+1))))
4344 sum=sum+sixth*hyi*hx2*(
4346 > cxi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
4347 > cx*(-fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
4349 > cxi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
4350 > cx*(-fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
4354 > xpi*(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))+
4355 > xp*(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
4357 > xpi*(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))+
4358 > xp*(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
4360 sum=sum+sixth*hyi*hz2*(
4362 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
4363 > xp*(-fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
4365 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
4366 > xp*(-fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
4368 sum=sum+z36th*hx2*hy*(
4370 > cxi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
4371 > cx*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
4373 > cxi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
4374 > cx*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
4376 sum=sum+z36th*hyi*hx2*hz2*(
4378 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
4379 > cx*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
4381 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
4382 > cx*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
4384 sum=sum+z36th*hy*hz2*(
4386 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
4387 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
4389 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
4390 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
4392 sum=sum+z216th*hx2*hy*hz2*(
4394 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
4395 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
4397 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
4398 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
4404 if(ict(4).eq.1)
then
4449 > xpi*(ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k))+
4450 > xp*(ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
4452 > xpi*(ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1))+
4453 > xp*(ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
4455 sum=sum+sixth*hx2*hzi*(
4457 > cxi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
4458 > cx*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
4460 > cxi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
4461 > cx*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
4463 sum=sum+sixth*hy2*hzi*(
4465 > xpi*(cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k))+
4466 > xp*(cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
4468 > xpi*(cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1))+
4469 > xp*(cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
4473 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
4474 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
4476 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
4477 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
4479 sum=sum+z36th*hx2*hy2*hzi*(
4481 > cxi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
4482 > cx*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
4484 > cxi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
4485 > cx*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
4487 sum=sum+z36th*hx2*hz*(
4489 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
4490 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
4492 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
4493 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
4495 sum=sum+z36th*hy2*hz*(
4497 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
4498 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
4500 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
4501 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
4503 sum=sum+z216th*hx2*hy2*hz*(
4505 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
4506 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
4508 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
4509 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
4515 if(ict(5).eq.1)
then
4555 > xpi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
4556 > xp*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
4558 > xpi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
4559 > xp*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
4563 > xpi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
4564 > xp*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
4566 > xpi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
4567 > xp*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
4571 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
4572 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
4574 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
4575 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
4577 sum=sum+z36th*hy2*hz2*(
4579 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
4580 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
4582 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
4583 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
4589 if(ict(6).eq.1)
then
4629 > xpi*(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))+
4630 > xp*(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
4632 > xpi*(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))+
4633 > xp*(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
4637 > cxi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
4638 > cx*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
4640 > cxi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
4641 > cx*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
4645 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
4646 > xp*(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
4648 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
4649 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
4651 sum=sum+z36th*hx2*hz2*(
4653 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
4654 > cx*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
4656 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
4657 > cx*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
4663 if(ict(7).eq.1)
then
4703 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
4704 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
4706 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
4707 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
4711 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
4712 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
4714 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
4715 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
4719 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
4720 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
4722 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
4723 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
4725 sum=sum+z36th*hx2*hy2*(
4727 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
4728 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
4730 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
4731 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
4737 if(ict(8).eq.1)
then
4781 > (fin(0,i,j,k) -fin(0,i,j+1,k))-
4782 > (fin(0,i+1,j,k)-fin(0,i+1,j+1,k)))
4784 > (fin(0,i,j,k+1) -fin(0,i,j+1,k+1))-
4785 > (fin(0,i+1,j,k+1)-fin(0,i+1,j+1,k+1))))
4787 sum=sum+sixth*hyi*hx(v)*(
4789 > cxdi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
4790 > cxd*(-fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
4792 > cxdi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
4793 > cxd*(-fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
4795 sum=sum+sixth*hxi(v)*hy*(
4797 > -(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))
4798 > +(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
4800 > -(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))
4801 > +(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
4803 sum=sum+sixth*hxi(v)*hyi*hz2*(
4805 > (fin(3,i,j,k) -fin(3,i,j+1,k))-
4806 > (fin(3,i+1,j,k)-fin(3,i+1,j+1,k)))
4808 > (fin(3,i,j,k+1) -fin(3,i,j+1,k+1))-
4809 > (fin(3,i+1,j,k+1)-fin(3,i+1,j+1,k+1))))
4811 sum=sum+z36th*hx(v)*hy*(
4813 > cxdi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
4814 > cxd*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
4816 > cxdi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
4817 > cxd*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
4819 sum=sum+z36th*hyi*hx(v)*hz2*(
4821 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
4822 > cxd*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
4824 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
4825 > cxd*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
4827 sum=sum+z36th*hxi(v)*hy*hz2*(
4829 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
4830 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
4832 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
4833 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
4835 sum=sum+z216th*hx(v)*hy*hz2*(
4837 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
4838 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
4840 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
4841 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
4847 if(ict(9).eq.1)
then
4891 > (ypi*fin(0,i,j,k) +yp*fin(0,i,j+1,k)) -
4892 > (ypi*fin(0,i+1,j,k)+yp*fin(0,i+1,j+1,k)))
4894 > (ypi*fin(0,i,j,k+1) +yp*fin(0,i,j+1,k+1)) -
4895 > (ypi*fin(0,i+1,j,k+1)+yp*fin(0,i+1,j+1,k+1))))
4897 sum=sum+sixth*hx(v)*hzi*(
4899 > cxdi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
4900 > cxd*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
4902 > cxdi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
4903 > cxd*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
4905 sum=sum+sixth*hxi(v)*hy2*hzi*(
4907 > (cyi*fin(2,i,j,k) +cy*fin(2,i,j+1,k)) -
4908 > (cyi*fin(2,i+1,j,k)+cy*fin(2,i+1,j+1,k)))
4910 > (cyi*fin(2,i,j,k+1) +cy*fin(2,i,j+1,k+1)) -
4911 > (cyi*fin(2,i+1,j,k+1)+cy*fin(2,i+1,j+1,k+1))))
4913 sum=sum+sixth*hxi(v)*hz*(
4915 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
4916 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
4918 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
4919 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
4921 sum=sum+z36th*hx(v)*hy2*hzi*(
4923 > cxdi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
4924 > cxd*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
4926 > cxdi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
4927 > cxd*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
4929 sum=sum+z36th*hx(v)*hz*(
4931 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
4932 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
4934 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
4935 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
4937 sum=sum+z36th*hxi(v)*hy2*hz*(
4939 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
4940 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
4942 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
4943 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
4945 sum=sum+z216th*hx(v)*hy2*hz*(
4947 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
4948 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
4950 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
4951 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
4957 if(ict(10).eq.1)
then
5001 > xpi*(fin(0,i,j,k) -fin(0,i,j+1,k))+
5002 > xp*(fin(0,i+1,j,k)-fin(0,i+1,j+1,k)))
5004 > xpi*(fin(0,i,j,k+1) -fin(0,i,j+1,k+1))+
5005 > xp*(fin(0,i+1,j,k+1)-fin(0,i+1,j+1,k+1))))
5007 sum=sum+sixth*hyi*hx2*hzi*(
5009 > cxi*(fin(1,i,j,k) -fin(1,i,j+1,k))+
5010 > cx*(fin(1,i+1,j,k)-fin(1,i+1,j+1,k)))
5012 > cxi*(fin(1,i,j,k+1) -fin(1,i,j+1,k+1))+
5013 > cx*(fin(1,i+1,j,k+1)-fin(1,i+1,j+1,k+1))))
5015 sum=sum+sixth*hy*hzi*(
5017 > xpi*(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))+
5018 > xp*(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
5020 > xpi*(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))+
5021 > xp*(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
5023 sum=sum+sixth*hyi*hz*(
5025 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
5026 > xp*(-fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
5028 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
5029 > xp*(-fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
5031 sum=sum+z36th*hx2*hy*hzi*(
5033 > cxi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
5034 > cx*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
5036 > cxi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
5037 > cx*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
5039 sum=sum+z36th*hyi*hx2*hz*(
5041 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
5042 > cx*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
5044 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
5045 > cx*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
5047 sum=sum+z36th*hy*hz*(
5049 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
5050 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
5052 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
5053 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
5055 sum=sum+z216th*hx2*hy*hz*(
5057 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
5058 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
5060 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
5061 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
5070 else if(ict(1).eq.3)
then
5071 if(ict(2).eq.1)
then
5108 > xpi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
5109 > xp*( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
5111 > xpi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
5112 > xp*( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
5116 > xpi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
5117 > xp*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
5119 > xpi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
5120 > xp*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
5122 sum=sum+sixth*hz2*hyi*(
5124 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
5125 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
5127 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
5128 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
5130 sum=sum+z36th*hy*hz2*(
5132 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
5133 > xp*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
5135 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
5136 > xp*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
5142 if(ict(3).eq.1)
then
5179 > xpi*(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))+
5180 > xp*(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
5182 > xpi*(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))+
5183 > xp*(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
5185 sum=sum+sixth*hy2*hzi*(
5187 > xpi*(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))+
5188 > xp*(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
5190 > xpi*(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))+
5191 > xp*(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
5195 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
5196 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
5198 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
5199 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
5201 sum=sum+z36th*hy2*hz*(
5203 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
5204 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
5206 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
5207 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
5213 if(ict(4).eq.1)
then
5250 > -(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))
5251 > +(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
5253 > -(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))
5254 > +(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
5256 sum=sum+sixth*hx(v)*(
5258 > cxdi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
5259 > cxd*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
5261 > cxdi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
5262 > cxd*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
5264 sum=sum+sixth*hz2*hxi(v)*(
5266 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
5267 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
5269 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
5270 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
5272 sum=sum+z36th*hx(v)*hz2*(
5274 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
5275 > cxd*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
5277 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
5278 > cxd*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
5284 if(ict(5).eq.1)
then
5323 sum=hxi(v)*hyi*hzi*(
5325 > (fin(0,i,j,k) -fin(0,i,j+1,k))-
5326 > (fin(0,i+1,j,k)-fin(0,i+1,j+1,k)))
5328 > (fin(0,i,j,k+1) -fin(0,i,j+1,k+1))-
5329 > (fin(0,i+1,j,k+1)-fin(0,i+1,j+1,k+1))))
5331 sum=sum+sixth*hyi*hx(v)*hzi*(
5333 > cxdi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
5334 > cxd*(-fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
5336 > cxdi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
5337 > cxd*(-fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
5339 sum=sum+sixth*hxi(v)*hy*hzi*(
5341 > -(cydi*fin(2,i,j,k) +cyd*fin(2,i,j+1,k))
5342 > +(cydi*fin(2,i+1,j,k)+cyd*fin(2,i+1,j+1,k)))
5344 > -(cydi*fin(2,i,j,k+1) +cyd*fin(2,i,j+1,k+1))
5345 > +(cydi*fin(2,i+1,j,k+1)+cyd*fin(2,i+1,j+1,k+1))))
5347 sum=sum+sixth*hxi(v)*hyi*hz*(
5349 > (fin(3,i,j,k) -fin(3,i,j+1,k))-
5350 > (fin(3,i+1,j,k)-fin(3,i+1,j+1,k)))
5352 > (fin(3,i,j,k+1) -fin(3,i,j+1,k+1))-
5353 > (fin(3,i+1,j,k+1)-fin(3,i+1,j+1,k+1))))
5355 sum=sum+z36th*hx(v)*hy*hzi*(
5357 > cxdi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
5358 > cxd*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
5360 > cxdi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
5361 > cxd*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
5363 sum=sum+z36th*hyi*hx(v)*hz*(
5365 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
5366 > cxd*(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
5368 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
5369 > cxd*(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
5371 sum=sum+z36th*hxi(v)*hy*hz*(
5373 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
5374 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
5376 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
5377 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
5379 sum=sum+z216th*hx(v)*hy*hz*(
5381 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
5382 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
5384 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
5385 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
5391 if(ict(6).eq.1)
then
5429 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
5430 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
5432 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
5433 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
5435 sum=sum+sixth*hx(v)*(
5437 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
5438 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
5440 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
5441 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
5443 sum=sum+sixth*hy2*hxi(v)*(
5445 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
5446 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
5448 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
5449 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
5451 sum=sum+z36th*hx(v)*hy2*(
5453 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
5454 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
5456 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
5457 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
5463 if(ict(7).eq.1)
then
5501 > xpi*(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))+
5502 > xp*(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
5504 > xpi*(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))+
5505 > xp*(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
5507 sum=sum+sixth*hx2*hzi*(
5509 > cxi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
5510 > cx*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
5512 > cxi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
5513 > cx*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
5517 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
5518 > xp*(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
5520 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
5521 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
5523 sum=sum+z36th*hx2*hz*(
5525 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
5526 > cx*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
5528 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
5529 > cx*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
5535 if(ict(8).eq.1)
then
5572 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
5573 > xp*( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
5575 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
5576 > xp*( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
5578 sum=sum+sixth*hx2*hyi*(
5580 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
5581 > cx*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
5583 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
5584 > cx*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
5588 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
5589 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
5591 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
5592 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
5594 sum=sum+z36th*hx2*hy*(
5596 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
5597 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
5599 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
5600 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
5608 else if(ict(1).eq.-3)
then
5609 if(ict(2).eq.1)
then
5641 > -(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))
5642 > +(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
5644 > -(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))
5645 > +(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
5647 sum=sum+sixth*hy2*hxi(v)*(
5649 > -(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))
5650 > +(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
5652 > -(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))
5653 > +(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
5655 sum=sum+sixth*hz2*hxi(v)*(
5657 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
5658 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
5660 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
5661 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
5663 sum=sum+z36th*hy2*hz2*hxi(v)*(
5665 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
5666 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
5668 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
5669 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
5675 if(ict(3).eq.1)
then
5708 > xpi*(-fin(2,i,j,k) +fin(2,i,j+1,k))+
5709 > xp*( -fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
5711 > xpi*(-fin(2,i,j,k+1) +fin(2,i,j+1,k+1))+
5712 > xp*( -fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
5714 sum=sum+sixth*hx2*hyi*(
5716 > cxi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
5717 > cx*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
5719 > cxi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
5720 > cx*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
5722 sum=sum+sixth*hz2*hyi*(
5724 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
5725 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
5727 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
5728 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
5730 sum=sum+z36th*hx2*hz2*hyi*(
5732 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
5733 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
5735 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
5736 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
5742 if(ict(4).eq.1)
then
5775 > xpi*(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))+
5776 > xp*(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
5778 > xpi*(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))+
5779 > xp*(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
5781 sum=sum+sixth*hx2*hzi*(
5783 > cxi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
5784 > cx*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
5786 > cxi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
5787 > cx*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
5789 sum=sum+sixth*hy2*hzi*(
5791 > xpi*(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))+
5792 > xp*(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
5794 > xpi*(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))+
5795 > xp*(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
5797 sum=sum+z36th*hx2*hy2*hzi*(
5799 > cxi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
5800 > cx*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
5802 > cxi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
5803 > cx*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
5812 else if(ict(1).eq.4)
then
5813 if(ict(2).eq.1)
then
5845 > xpi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
5846 > xp*( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
5848 > xpi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
5849 > xp*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
5853 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
5854 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
5856 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
5857 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
5863 if(ict(3).eq.1)
then
5899 > xpi*(-fin(1,i,j,k) +fin(1,i,j+1,k))+
5900 > xp*( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
5902 > xpi*(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))+
5903 > xp*( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
5905 sum=sum+sixth*hy*hzi*(
5907 > xpi*(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
5908 > xp*(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
5910 > xpi*(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
5911 > xp*(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
5913 sum=sum+sixth*hz*hyi*(
5915 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
5916 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
5918 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
5919 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
5921 sum=sum+z36th*hy*hz*(
5923 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
5924 > xp*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
5926 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
5927 > xp*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
5933 if(ict(4).eq.1)
then
5965 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
5966 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
5968 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
5969 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
5973 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
5974 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
5976 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
5977 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
5983 if(ict(5).eq.1)
then
6019 > -(ypi*fin(2,i,j,k) +yp*fin(2,i,j+1,k))
6020 > +(ypi*fin(2,i+1,j,k)+yp*fin(2,i+1,j+1,k)))
6022 > -(ypi*fin(2,i,j,k+1) +yp*fin(2,i,j+1,k+1))
6023 > +(ypi*fin(2,i+1,j,k+1)+yp*fin(2,i+1,j+1,k+1))))
6025 sum=sum+sixth*hx(v)*hzi*(
6027 > cxdi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
6028 > cxd*(ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
6030 > cxdi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
6031 > cxd*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
6033 sum=sum+sixth*hz*hxi(v)*(
6035 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
6036 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
6038 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
6039 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
6041 sum=sum+z36th*hx(v)*hz*(
6043 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
6044 > cxd*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
6046 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
6047 > cxd*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
6053 if(ict(6).eq.1)
then
6089 > ( +fin(3,i,j,k) -fin(3,i,j+1,k))
6090 > +( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
6092 > ( +fin(3,i,j,k+1) -fin(3,i,j+1,k+1))
6093 > +( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
6095 sum=sum+sixth*hx(v)*hyi*(
6097 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
6098 > cxd*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
6100 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
6101 > cxd*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
6103 sum=sum+sixth*hy*hxi(v)*(
6105 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
6106 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
6108 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
6109 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
6111 sum=sum+z36th*hx(v)*hy*(
6113 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
6114 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
6116 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
6117 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
6123 if(ict(7).eq.1)
then
6155 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
6156 > xp*( ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
6158 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
6159 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
6163 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
6164 > cx*( ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
6166 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
6167 > cx*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
6176 else if(ict(1).eq.-4)
then
6177 if(ict(2).eq.1)
then
6209 > ( fin(1,i,j,k) -fin(1,i,j+1,k))+
6210 > ( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
6212 > ( fin(1,i,j,k+1) -fin(1,i,j+1,k+1))+
6213 > ( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
6215 sum=sum+sixth*hy*hxi(v)*(
6217 > -(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))+
6218 > (cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
6220 > -(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))+
6221 > (cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
6223 sum=sum+sixth*hz2*hyi*hxi(v)*(
6225 > ( fin(5,i,j,k) -fin(5,i,j+1,k))+
6226 > ( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
6228 > ( fin(5,i,j,k+1) -fin(5,i,j+1,k+1))+
6229 > ( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
6231 sum=sum+z36th*hy*hz2*hxi(v)*(
6233 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
6234 > (cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
6236 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
6237 > (cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
6243 if(ict(3).eq.1)
then
6275 > +(ypi*fin(1,i,j,k) +yp*fin(1,i,j+1,k))
6276 > -(ypi*fin(1,i+1,j,k)+yp*fin(1,i+1,j+1,k)))
6278 > -(ypi*fin(1,i,j,k+1) +yp*fin(1,i,j+1,k+1))
6279 > +(ypi*fin(1,i+1,j,k+1)+yp*fin(1,i+1,j+1,k+1))))
6281 sum=sum+sixth*hy2*hzi*hxi(v)*(
6283 > +(cyi*fin(4,i,j,k) +cy*fin(4,i,j+1,k))
6284 > -(cyi*fin(4,i+1,j,k)+cy*fin(4,i+1,j+1,k)))
6286 > -(cyi*fin(4,i,j,k+1) +cy*fin(4,i,j+1,k+1))
6287 > +(cyi*fin(4,i+1,j,k+1)+cy*fin(4,i+1,j+1,k+1))))
6289 sum=sum+sixth*hz*hxi(v)*(
6291 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
6292 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
6294 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
6295 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
6297 sum=sum+z36th*hy2*hz*hxi(v)*(
6299 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
6300 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
6302 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
6303 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
6309 if(ict(4).eq.1)
then
6341 > ( fin(2,i,j,k) -fin(2,i,j+1,k))
6342 > +(-fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
6344 > ( fin(2,i,j,k+1) -fin(2,i,j+1,k+1))
6345 > +(-fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
6347 sum=sum+sixth*hx(v)*hyi*(
6349 > cxdi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
6350 > cxd*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
6352 > cxdi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
6353 > cxd*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
6355 sum=sum+sixth*hz2*hxi(v)*hyi*(
6357 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
6358 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
6360 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
6361 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
6363 sum=sum+z36th*hx(v)*hz2*hyi*(
6365 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
6366 > cxd*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
6368 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
6369 > cxd*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
6375 if(ict(5).eq.1)
then
6407 > -(ypi*fin(3,i,j,k) +yp*fin(3,i,j+1,k))
6408 > +(ypi*fin(3,i+1,j,k)+yp*fin(3,i+1,j+1,k)))
6410 > -(ypi*fin(3,i,j,k+1) +yp*fin(3,i,j+1,k+1))
6411 > +(ypi*fin(3,i+1,j,k+1)+yp*fin(3,i+1,j+1,k+1))))
6413 sum=sum+sixth*hx(v)*hzi*(
6415 > cxdi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
6416 > cxd*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
6418 > cxdi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
6419 > cxd*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
6421 sum=sum+sixth*hy2*hxi(v)*hzi*(
6423 > -(cyi*fin(6,i,j,k) +cy*fin(6,i,j+1,k))
6424 > +(cyi*fin(6,i+1,j,k)+cy*fin(6,i+1,j+1,k)))
6426 > -(cyi*fin(6,i,j,k+1) +cy*fin(6,i,j+1,k+1))
6427 > +(cyi*fin(6,i+1,j,k+1)+cy*fin(6,i+1,j+1,k+1))))
6429 sum=sum+z36th*hx(v)*hy2*hzi*(
6431 > cxdi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
6432 > cxd*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
6434 > cxdi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
6435 > cxd*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
6441 if(ict(6).eq.1)
then
6473 > xpi*(-fin(2,i,j,k) +fin(2,i,j+1,k))+
6474 > xp*( -fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
6476 > xpi*(-fin(2,i,j,k+1) +fin(2,i,j+1,k+1))+
6477 > xp*( -fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
6479 sum=sum+sixth*hx2*hyi*hzi*(
6481 > cxi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
6482 > cx*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
6484 > cxi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
6485 > cx*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
6487 sum=sum+sixth*hz*hyi*(
6489 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
6490 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
6492 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
6493 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
6495 sum=sum+z36th*hx2*hz*hyi*(
6497 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
6498 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
6500 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
6501 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
6507 if(ict(7).eq.1)
then
6539 > xpi*(-fin(3,i,j,k) +fin(3,i,j+1,k))+
6540 > xp*( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
6542 > xpi*(-fin(3,i,j,k+1) +fin(3,i,j+1,k+1))+
6543 > xp*( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
6545 sum=sum+sixth*hx2*hyi*hzi*(
6547 > cxi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
6548 > cx*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
6550 > cxi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
6551 > cx*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
6553 sum=sum+sixth*hy*hzi*(
6555 > xpi*(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))+
6556 > xp*(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
6558 > xpi*(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))+
6559 > xp*(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
6561 sum=sum+z36th*hx2*hy*hzi*(
6563 > cxi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
6564 > cx*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
6566 > cxi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
6567 > cx*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
6576 else if(ict(1).eq.5)
then
6577 if(ict(2).eq.1)
then
6608 > xpi*(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))+
6609 > xp*( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
6611 > xpi*(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))+
6612 > xp*(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
6616 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
6617 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
6619 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
6620 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
6626 if(ict(3).eq.1)
then
6657 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
6658 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
6660 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
6661 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
6665 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
6666 > xp*( cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
6668 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
6669 > xp*( cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
6675 if(ict(4).eq.1)
then
6706 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
6707 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
6709 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
6710 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
6712 sum=sum+sixth*hx(v)*(
6714 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
6715 > cxd*(ypi*fin(7,i+1,j,k) +yp*fin(7,i+1,j+1,k)))
6717 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
6718 > cxd*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
6727 else if(ict(1).eq.-5)
then
6728 if(ict(2).eq.1)
then
6755 > -(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))
6756 > +( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
6758 > -(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))
6759 > +(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
6761 sum=sum+sixth*hz2*hxi(v)*(
6763 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
6764 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
6766 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
6767 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
6773 if(ict(3).eq.1)
then
6802 sum=hyi*hzi*hxi(v)*(
6804 > -(-fin(1,i,j,k) +fin(1,i,j+1,k))
6805 > +( -fin(1,i+1,j,k)+fin(1,i+1,j+1,k)))
6807 > -(-fin(1,i,j,k+1) +fin(1,i,j+1,k+1))
6808 > +( -fin(1,i+1,j,k+1)+fin(1,i+1,j+1,k+1))))
6810 sum=sum+sixth*hy*hzi*hxi(v)*(
6812 > -(cydi*fin(4,i,j,k) +cyd*fin(4,i,j+1,k))
6813 > +(cydi*fin(4,i+1,j,k)+cyd*fin(4,i+1,j+1,k)))
6815 > -(cydi*fin(4,i,j,k+1) +cyd*fin(4,i,j+1,k+1))
6816 > +(cydi*fin(4,i+1,j,k+1)+cyd*fin(4,i+1,j+1,k+1))))
6818 sum=sum+sixth*hz*hyi*hxi(v)*(
6820 > -(-fin(5,i,j,k) +fin(5,i,j+1,k))
6821 > +( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
6823 > -(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))
6824 > +( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
6826 sum=sum+z36th*hy*hz*hxi(v)*(
6828 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))
6829 > +(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
6831 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))
6832 > +(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
6838 if(ict(4).eq.1)
then
6865 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
6866 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
6868 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
6869 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
6871 sum=sum+sixth*hy2*hxi(v)*(
6873 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
6874 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
6876 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
6877 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
6883 if(ict(5).eq.1)
then
6910 > xpi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
6911 > xp*( -fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
6913 > xpi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
6914 > xp*( -fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
6916 sum=sum+sixth*hz2*hyi*(
6918 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
6919 > xp*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
6921 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
6922 > xp*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
6928 if(ict(6).eq.1)
then
6955 > xpi*(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))+
6956 > xp*(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
6958 > xpi*(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))+
6959 > xp*(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
6961 sum=sum+sixth*hy2*hzi*(
6963 > xpi*(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))+
6964 > xp*(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
6966 > xpi*(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))+
6967 > xp*(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
6973 if(ict(7).eq.1)
then
7002 sum=hxi(v)*hzi*hyi*(
7004 > ( fin(2,i,j,k) -fin(2,i,j+1,k))
7005 > +(-fin(2,i+1,j,k)+fin(2,i+1,j+1,k)))
7007 > ( fin(2,i,j,k+1) -fin(2,i,j+1,k+1))
7008 > +(-fin(2,i+1,j,k+1)+fin(2,i+1,j+1,k+1))))
7010 sum=sum+sixth*hx(v)*hzi*hyi*(
7012 > cxdi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
7013 > cxd*(-fin(4,i+1,j,k) +fin(4,i+1,j+1,k)))
7015 > cxdi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
7016 > cxd*(-fin(4,i+1,j,k+1) +fin(4,i+1,j+1,k+1))))
7018 sum=sum+sixth*hz*hxi(v)*hyi*(
7020 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
7021 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
7023 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
7024 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
7026 sum=sum+z36th*hx(v)*hz*hyi*(
7028 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7029 > cxd*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
7031 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7032 > cxd*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
7038 if(ict(8).eq.1)
then
7067 sum=hyi*hxi(v)*hzi*(
7069 > ( +fin(3,i,j,k) -fin(3,i,j+1,k))
7070 > +( -fin(3,i+1,j,k)+fin(3,i+1,j+1,k)))
7072 > ( +fin(3,i,j,k+1) -fin(3,i,j+1,k+1))
7073 > +( -fin(3,i+1,j,k+1)+fin(3,i+1,j+1,k+1))))
7075 sum=sum+sixth*hx(v)*hyi*hzi*(
7077 > cxdi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
7078 > cxd*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
7080 > cxdi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
7081 > cxd*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
7083 sum=sum+sixth*hy*hxi(v)*hzi*(
7085 > -(cydi*fin(6,i,j,k) +cyd*fin(6,i,j+1,k))
7086 > +(cydi*fin(6,i+1,j,k)+cyd*fin(6,i+1,j+1,k)))
7088 > -(cydi*fin(6,i,j,k+1) +cyd*fin(6,i,j+1,k+1))
7089 > +(cydi*fin(6,i+1,j,k+1)+cyd*fin(6,i+1,j+1,k+1))))
7091 sum=sum+z36th*hx(v)*hy*hzi*(
7093 > cxdi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
7094 > cxd*(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
7096 > cxdi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
7097 > cxd*(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
7103 if(ict(9).eq.1)
then
7130 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
7131 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
7133 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
7134 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
7136 sum=sum+sixth*hx2*hyi*(
7138 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7139 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7141 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7142 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
7148 if(ict(10).eq.1)
then
7175 > xpi*(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))+
7176 > xp*( ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
7178 > xpi*(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))+
7179 > xp*(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
7181 sum=sum+sixth*hx2*hzi*(
7183 > cxi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
7184 > cx*( ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
7186 > cxi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
7187 > cx*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
7196 else if(ict(1).eq.6)
then
7222 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
7223 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
7225 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
7226 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
7235 if(ict(1).eq.-6)
then
7236 if(ict(2).eq.1)
then
7258 > ( fin(4,i,j,k) -fin(4,i,j+1,k))
7259 > +(-fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
7261 > ( fin(4,i,j,k+1) -fin(4,i,j+1,k+1))
7262 > +(-fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
7264 sum=sum+sixth*hz2*hyi*hxi(v)*(
7266 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
7267 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7269 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
7270 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
7276 if(ict(3).eq.1)
then
7302 > -(ypi*fin(4,i,j,k) +yp*fin(4,i,j+1,k))
7303 > +( ypi*fin(4,i+1,j,k)+yp*fin(4,i+1,j+1,k)))
7305 > -(ypi*fin(4,i,j,k+1) +yp*fin(4,i,j+1,k+1))
7306 > +(ypi*fin(4,i+1,j,k+1)+yp*fin(4,i+1,j+1,k+1))))
7308 sum=sum+sixth*hz*hxi(v)*(
7310 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
7311 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
7313 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
7314 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
7320 if(ict(4).eq.1)
then
7346 > ( fin(5,i,j,k) -fin(5,i,j+1,k))
7347 > +(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
7349 > ( fin(5,i,j,k+1) -fin(5,i,j+1,k+1))
7350 > +(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
7352 sum=sum+sixth*hy*hxi(v)*(
7354 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))
7355 > +(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
7357 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))
7358 > +(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
7364 if(ict(5).eq.1)
then
7386 > -(ypi*fin(5,i,j,k) +yp*fin(5,i,j+1,k))
7387 > +(ypi*fin(5,i+1,j,k)+yp*fin(5,i+1,j+1,k)))
7389 > -(ypi*fin(5,i,j,k+1) +yp*fin(5,i,j+1,k+1))
7390 > +(ypi*fin(5,i+1,j,k+1)+yp*fin(5,i+1,j+1,k+1))))
7392 sum=sum+sixth*hy2*hxi(v)*hzi*(
7394 > -(cyi*fin(7,i,j,k) +cy*fin(7,i,j+1,k))
7395 > +(cyi*fin(7,i+1,j,k)+cy*fin(7,i+1,j+1,k)))
7397 > -(cyi*fin(7,i,j,k+1) +cy*fin(7,i,j+1,k+1))
7398 > +(cyi*fin(7,i+1,j,k+1)+cy*fin(7,i+1,j+1,k+1))))
7404 if(ict(6).eq.1)
then
7430 > xpi*(-fin(4,i,j,k) +fin(4,i,j+1,k))+
7431 > xp*(-fin(4,i+1,j,k) +fin(4,i+1,j+1,k)))
7433 > xpi*(-fin(4,i,j,k+1) +fin(4,i,j+1,k+1))+
7434 > xp*(-fin(4,i+1,j,k+1) +fin(4,i+1,j+1,k+1))))
7436 sum=sum+sixth*hz*hyi*(
7438 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7439 > xp*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
7441 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7442 > xp*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
7448 if(ict(7).eq.1)
then
7474 > xpi*(-fin(5,i,j,k) +fin(5,i,j+1,k))+
7475 > xp*( -fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
7477 > xpi*(-fin(5,i,j,k+1) +fin(5,i,j+1,k+1))+
7478 > xp*( -fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
7480 sum=sum+sixth*hy*hzi*(
7482 > xpi*(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))+
7483 > xp*( cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
7485 > xpi*(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))+
7486 > xp*( cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
7492 if(ict(8).eq.1)
then
7518 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
7519 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
7521 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
7522 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
7524 sum=sum+sixth*hx(v)*hyi*(
7526 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7527 > cxd*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
7529 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7530 > cxd*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
7536 if(ict(9).eq.1)
then
7562 > -(ypi*fin(6,i,j,k) +yp*fin(6,i,j+1,k))
7563 > +(ypi*fin(6,i+1,j,k)+yp*fin(6,i+1,j+1,k)))
7565 > -(ypi*fin(6,i,j,k+1) +yp*fin(6,i,j+1,k+1))
7566 > +(ypi*fin(6,i+1,j,k+1)+yp*fin(6,i+1,j+1,k+1))))
7568 sum=sum+sixth*hx(v)*hzi*(
7570 > cxdi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
7571 > cxd*(ypi*fin(7,i+1,j,k) +yp*fin(7,i+1,j+1,k)))
7573 > cxdi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
7574 > cxd*( ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
7580 if(ict(10).eq.1)
then
7602 > xpi*(-fin(6,i,j,k) +fin(6,i,j+1,k))+
7603 > xp*( -fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
7605 > xpi*(-fin(6,i,j,k+1) +fin(6,i,j+1,k+1))+
7606 > xp*( -fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
7608 sum=sum+sixth*hx2*hyi*hzi*(
7610 > cxi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7611 > cx*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7613 > cxi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7614 > cx*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
7623 else if(abs(ict(1)).eq.7)
then
7624 if(ict(2).eq.1)
then
7643 sum=hyi*hxi(v)*hzi*(
7645 > ( fin(4,i,j,k) -fin(4,i,j+1,k))
7646 > +(-fin(4,i+1,j,k)+fin(4,i+1,j+1,k)))
7648 > ( fin(4,i,j,k+1) -fin(4,i,j+1,k+1))
7649 > +(-fin(4,i+1,j,k+1)+fin(4,i+1,j+1,k+1))))
7651 sum=sum+sixth*hz*hyi*hxi(v)*(
7653 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
7654 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7656 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
7657 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
7663 if(ict(3).eq.1)
then
7684 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
7685 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
7687 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
7688 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
7694 if(ict(4).eq.1)
then
7713 sum=hxi(v)*hyi*hzi*(
7715 > ( fin(5,i,j,k) -fin(5,i,j+1,k))
7716 > +(-fin(5,i+1,j,k)+fin(5,i+1,j+1,k)))
7718 > ( fin(5,i,j,k+1) -fin(5,i,j+1,k+1))
7719 > +(-fin(5,i+1,j,k+1)+fin(5,i+1,j+1,k+1))))
7721 sum=sum+sixth*hy*hxi(v)*hzi*(
7723 > -(cydi*fin(7,i,j,k) +cyd*fin(7,i,j+1,k))
7724 > +(cydi*fin(7,i+1,j,k)+cyd*fin(7,i+1,j+1,k)))
7726 > -(cydi*fin(7,i,j,k+1) +cyd*fin(7,i,j+1,k+1))
7727 > +(cydi*fin(7,i+1,j,k+1)+cyd*fin(7,i+1,j+1,k+1))))
7733 if(ict(5).eq.1)
then
7754 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7755 > xp*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7757 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7758 > xp*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
7764 if(ict(6).eq.1)
then
7785 > xpi*(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))+
7786 > xp*(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
7788 > xpi*(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))+
7789 > xp*(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
7795 if(ict(7).eq.1)
then
7814 sum=hxi(v)*hyi*hzi*(
7816 > ( fin(6,i,j,k) -fin(6,i,j+1,k))
7817 > +(-fin(6,i+1,j,k)+fin(6,i+1,j+1,k)))
7819 > ( fin(6,i,j,k+1) -fin(6,i,j+1,k+1))
7820 > +(-fin(6,i+1,j,k+1)+fin(6,i+1,j+1,k+1))))
7822 sum=sum+sixth*hx(v)*hyi*hzi*(
7824 > cxdi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7825 > cxd*(-fin(7,i+1,j,k) +fin(7,i+1,j+1,k)))
7827 > cxdi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7828 > cxd*(-fin(7,i+1,j,k+1) +fin(7,i+1,j+1,k+1))))
7837 else if(abs(ict(1)).eq.8)
then
7838 if(ict(2).eq.1)
then
7854 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
7855 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7857 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
7858 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
7864 if(ict(3).eq.1)
then
7880 > -(ypi*fin(7,i,j,k) +yp*fin(7,i,j+1,k))
7881 > +(ypi*fin(7,i+1,j,k)+yp*fin(7,i+1,j+1,k)))
7883 > -(ypi*fin(7,i,j,k+1) +yp*fin(7,i,j+1,k+1))
7884 > +(ypi*fin(7,i+1,j,k+1)+yp*fin(7,i+1,j+1,k+1))))
7890 if(ict(4).eq.1)
then
7906 > xpi*(-fin(7,i,j,k) +fin(7,i,j+1,k))+
7907 > xp*( -fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7909 > xpi*(-fin(7,i,j,k+1) +fin(7,i,j+1,k+1))+
7910 > xp*( -fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))
7919 else if(abs(ict(1)).eq.9)
then
7928 sum=hyi*hxi(v)*hzi*(
7930 > ( fin(7,i,j,k) -fin(7,i,j+1,k))
7931 > +(-fin(7,i+1,j,k)+fin(7,i+1,j+1,k)))
7933 > ( fin(7,i,j,k+1) -fin(7,i,j+1,k+1))
7934 > +(-fin(7,i+1,j,k+1)+fin(7,i+1,j+1,k+1))))