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