SUBROUTINE ORIENT(T,NOEL,NPT,LAYER,KSPT,COORDS,BASIS,ANAME, 1 NNODES,CNODES,JNNUM) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 ANAME DIMENSION T(3,3),COORDS(3),BASIS(3,3),CNODES(3,NNODES) DIMENSION JNNUM(NNODES) C PARAMETER(ZERO=0.D0,ONE=1.D0,P1=0.70711D0) PARAMETER (MSGDEV=7) C IF (JNNUM(2).EQ.8) THEN C C LOCAL DIRECTION 3 OF COUPLING NODE 8 IS C IN THE RADIAL DIRECTION C T(1,1)= P1 T(2,1)= P1 T(3,1)=ZERO T(1,2)=ZERO T(2,2)=ZERO T(3,2)= ONE T(1,3)= P1 T(2,3)= -P1 T(3,3)= ONE C ELSE IF (JNNUM(2).EQ.1) THEN C C LOCAL DIRECTION 2 OF COUPLING NODE 1 IS C IN THE RADIAL DIRECTION C T(1,1)=ZERO T(2,1)=-ONE T(3,1)=ZERO T(1,2)= ONE T(2,2)=ZERO T(3,2)=ZERO T(1,3)=ZERO T(2,3)=ZERO T(3,3)= ONE C ELSE C WRITE (MSGDEV,*) 'ERROR. ORIENTATION REQUESTED FOR ', $ 'INVALID NODE' CALL XIT C END IF RETURN END