2      &                   u3,           dxidx,    rLui,     argument
32      &            dxidx(npro,nsd,nsd), 
44       call e3gijd( dxidx, gijd )
329       subroutine e3uBar (rho,          ui,         dxidx,       argument
335      &           dxidx(npro,nsd,nsd),  rLui(npro,nsd),  local
343       call e3gijd( dxidx, gijd )
396       subroutine e3gijd( dxidx,  gijd )  argument
400       real*8  dxidx(npro,nsd,nsd),  gijd(npro,6),  local
410          gijd(:,1) = dxidx(:,1,1) * dxidx(:,1,1)
411      &             + dxidx(:,2,1) * dxidx(:,2,1)
412      &             + dxidx(:,3,1) * dxidx(:,3,1)
414          gijd(:,4) = dxidx(:,1,1) * dxidx(:,1,2)
415      &             + dxidx(:,2,1) * dxidx(:,2,2)
416      &             + dxidx(:,3,1) * dxidx(:,3,2)
418          gijd(:,2) = dxidx(:,1,2) * dxidx(:,1,2)
419      &             + dxidx(:,2,2) * dxidx(:,2,2)
420      &             + dxidx(:,3,2) * dxidx(:,3,2)
422          gijd(:,5) = dxidx(:,1,2) * dxidx(:,1,3)
423      &             + dxidx(:,2,2) * dxidx(:,2,3)
424      &             + dxidx(:,3,2) * dxidx(:,3,3)
426          gijd(:,6) = dxidx(:,1,1) * dxidx(:,1,3)
427      &             + dxidx(:,2,1) * dxidx(:,2,3)
428      &             + dxidx(:,3,1) * dxidx(:,3,3)
430          gijd(:,3) = dxidx(:,1,3) * dxidx(:,1,3)
431      &             + dxidx(:,2,3) * dxidx(:,2,3)
432      &             + dxidx(:,3,3) * dxidx(:,3,3)
443          tmp1(:) = c1 * dxidx(:,1,1)+c2 *(dxidx(:,2,1)+dxidx(:,3,1))
444          tmp2(:) = c1 * dxidx(:,2,1)+c2 *(dxidx(:,1,1)+dxidx(:,3,1))
445          tmp3(:) = c1 * dxidx(:,3,1)+c2 *(dxidx(:,1,1)+dxidx(:,2,1))
446          gijd(:,1) = dxidx(:,1,1) * tmp1
447      1              + dxidx(:,2,1) * tmp2
448      2              + dxidx(:,3,1) * tmp3
450          tmp1(:) = c1 * dxidx(:,1,2)+c2 *(dxidx(:,2,2)+dxidx(:,3,2))
451          tmp2(:) = c1 * dxidx(:,2,2)+c2 *(dxidx(:,1,2)+dxidx(:,3,2))
452          tmp3(:) = c1 * dxidx(:,3,2)+c2 *(dxidx(:,1,2)+dxidx(:,2,2))
453          gijd(:,2) = dxidx(:,1,2) * tmp1
454      1             + dxidx(:,2,2) * tmp2
455      2             + dxidx(:,3,2) * tmp3
457          gijd(:,4) = dxidx(:,1,1) * tmp1
458      1             + dxidx(:,2,1) * tmp2
459      2             + dxidx(:,3,1) * tmp3
461          tmp1(:) = c1 * dxidx(:,1,3)+c2 *(dxidx(:,2,3)+dxidx(:,3,3))
462          tmp2(:) = c1 * dxidx(:,2,3)+c2 *(dxidx(:,1,3)+dxidx(:,3,3))
463          tmp3(:) = c1 * dxidx(:,3,3)+c2 *(dxidx(:,1,3)+dxidx(:,2,3))
464          gijd(:,3) = dxidx(:,1,3) * tmp1
465      1             + dxidx(:,2,3) * tmp2
466      2             + dxidx(:,3,3) * tmp3
468          gijd(:,5) = dxidx(:,1,2) * tmp1
469      1             + dxidx(:,2,2) * tmp2
470      2             + dxidx(:,3,2) * tmp3
472          gijd(:,6) = dxidx(:,1,1) * tmp1
473      1             + dxidx(:,2,1) * tmp2
474      2             + dxidx(:,3,1) * tmp3
489       subroutine e3StabSclr (uMod,  dxidx,  tauT, diffus, srcR, giju,  argument
496      &            dxidx(npro,nsd,nsd),       diffus(npro),  local
513       call e3gijd( dxidx, gijd )