subroutine orient(m,n,kc,dorien) c* * * * * * c c user subroutine for input of orientation for anisotropic option. c c m element number c n integration point number c kc layer number c dorien orientation matrix c c* * * * * * implicit real*8 (a-h,o-z) dp C include '../common/lass' include '../common/dimen' include '../common/space' include '../common/arrays' include '../common/array4' include '../common/heat' dimension dorien(3,3),m(2),ccint(12) LA1 = (M(1)-1)*NELSTR + ICRXPT + (N-1)*NCRD DO 1 II=1,NCRD CCINT(II)=VARS(LA1) write(72,*) ii,nelstr,icrxpt,n,m,ccint(ii) 1 LA1 = LA1 + 1 do 10 i=1,3 do 20 j=1,2 dorien(i,j)=0. 20 continue 10 continue dorien(1,1)=cosd(25.) dorien(1,2)=-sind(25.) dorien(2,1)=sind(25.) dorien(2,2)=cosd(25.) dorien(3,3)=1. return end