1*59599516SKenneth E. Jansen 2*59599516SKenneth E. Jansen 3*59599516SKenneth E. Jansen subroutine get_a_not_hex(xc,anot) 4*59599516SKenneth E. Jansen 5*59599516SKenneth E. Jansen include "common.h" 6*59599516SKenneth E. Jansen 7*59599516SKenneth E. Jansen dimension xc(npro,nenl,nsd), anot(npro,nenl,nsd) 8*59599516SKenneth E. Jansen 9*59599516SKenneth E. Jansen 10*59599516SKenneth E. Jansen do i = 1, nsd 11*59599516SKenneth E. Jansen 12*59599516SKenneth E. Jansen anot(:,1,i) = pt125*(xc(:,1,i)+xc(:,2,i)+xc(:,3,i)+xc(:,4,i) 13*59599516SKenneth E. Jansen & +xc(:,5,i)+xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 14*59599516SKenneth E. Jansen 15*59599516SKenneth E. Jansen anot(:,2,i) = pt125*(-xc(:,1,i)+xc(:,2,i)+xc(:,3,i)-xc(:,4,i) 16*59599516SKenneth E. Jansen & -xc(:,5,i)+xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 17*59599516SKenneth E. Jansen 18*59599516SKenneth E. Jansen anot(:,3,i) = pt125*(-xc(:,1,i)-xc(:,2,i)+xc(:,3,i)+xc(:,4,i) 19*59599516SKenneth E. Jansen & -xc(:,5,i)-xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 20*59599516SKenneth E. Jansen 21*59599516SKenneth E. Jansen anot(:,4,i) = pt125*(-xc(:,1,i)-xc(:,2,i)-xc(:,3,i)-xc(:,4,i) 22*59599516SKenneth E. Jansen & +xc(:,5,i)+xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 23*59599516SKenneth E. Jansen 24*59599516SKenneth E. Jansen anot(:,5,i) = pt125*(xc(:,1,i)-xc(:,2,i)+xc(:,3,i)-xc(:,4,i) 25*59599516SKenneth E. Jansen & +xc(:,5,i)-xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 26*59599516SKenneth E. Jansen 27*59599516SKenneth E. Jansen anot(:,6,i) = pt125*(xc(:,1,i)+xc(:,2,i)-xc(:,3,i)-xc(:,4,i) 28*59599516SKenneth E. Jansen & -xc(:,5,i)-xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 29*59599516SKenneth E. Jansen 30*59599516SKenneth E. Jansen anot(:,7,i) = pt125*(xc(:,1,i)-xc(:,2,i)-xc(:,3,i)+xc(:,4,i) 31*59599516SKenneth E. Jansen & -xc(:,5,i)+xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 32*59599516SKenneth E. Jansen 33*59599516SKenneth E. Jansen anot(:,8,i) = pt125*(-xc(:,1,i)+xc(:,2,i)-xc(:,3,i)+xc(:,4,i) 34*59599516SKenneth E. Jansen & +xc(:,5,i)-xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 35*59599516SKenneth E. Jansen 36*59599516SKenneth E. Jansen enddo 37*59599516SKenneth E. Jansen 38*59599516SKenneth E. Jansen return 39*59599516SKenneth E. Jansen end 40*59599516SKenneth E. Jansen 41*59599516SKenneth E. Jansen 42*59599516SKenneth E. Jansen subroutine get_a_not_tet(xc,anot) 43*59599516SKenneth E. Jansen 44*59599516SKenneth E. Jansen include "common.h" 45*59599516SKenneth E. Jansen 46*59599516SKenneth E. Jansen dimension xc(npro,nenl,nsd), anot(npro,nenl,nsd) 47*59599516SKenneth E. Jansen 48*59599516SKenneth E. Jansen 49*59599516SKenneth E. Jansen do i = 1, nsd 50*59599516SKenneth E. Jansen 51*59599516SKenneth E. Jansen anot(:,1,i) = xc(:,4,i) 52*59599516SKenneth E. Jansen anot(:,2,i) = xc(:,1,i)-xc(:,4,i) 53*59599516SKenneth E. Jansen anot(:,3,i) = xc(:,2,i)-xc(:,4,i) 54*59599516SKenneth E. Jansen anot(:,4,i) = xc(:,3,i)-xc(:,4,i) 55*59599516SKenneth E. Jansen 56*59599516SKenneth E. Jansen enddo 57*59599516SKenneth E. Jansen 58*59599516SKenneth E. Jansen return 59*59599516SKenneth E. Jansen end 60*59599516SKenneth E. Jansen 61