subroutine orient(mu,nu,kcu,dorien)
c
c CONTACT: ap@marc.de
c      
c* * * * * *
c
c     user subroutine for input of orientation for anisotropic option.
c
c     mu           element number
c     nu           integration point number
c     kc           layer number
c     dorien       orientation matrix
c
c* * * * * *
c
c  preferred system is cylindrical for 2D-planar problems
c  dorien(1,ijk) = 1st preferred direction is radial
c  dorien(2,ijk) = 2nd preferred direction is circumferential
c  dorien(3,ijk) = 3rd preferred direction is normal to plane
c
c* * * * * *
c
      implicit real*8 (a-h,o-z)                                               dp
      dimension dorien(3,3),m(2),xip(12)
c
      include '../common/space'
      include '../common/heat'
      include '../common/lass'
      include '../common/dimen'
      include '../common/arrays'
      include '../common/array4'
c
c use el. integr. point coordinates to define the transformations
      la1 = icrxpt + (n-1)*nelstr + (nn-1)*ncrd
      do ijk=1,ncrd
       xip(ijk) = vars(la1+ijk-1)
      enddo
c
      anorm = 0d0
      do ijk=1,ncrd
       anorm =anorm + xip(ijk)**2
      enddo
      anorm = dsqrt(anorm)
c
      do ijk=1,ncrd
       dorien(1,ijk) = xip(ijk)/anorm
      enddo
c
      dorien(2,1) = -dorien(1,2)
      dorien(2,2) =  dorien(1,1)
c
      dorien(1,3) = 0d0
      dorien(3,1) = 0d0
      dorien(2,3) = 0d0
      dorien(3,2) = 0d0
      dorien(3,3) = 1d0
c
      return
      end