roto_tens2d_aniso.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	roto_tens2d_aniso (nel, tens, dir_a, dir_b)

Function/Subroutine Documentation

roto_tens2d_aniso()

subroutine roto_tens2d_aniso	(	integer, intent(in)	nel,
		intent(inout)	tens,
		intent(in)	dir_a,
		intent(in)	dir_b )

Definition at line 28 of file roto_tens2d_aniso.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ,INTENT(IN) :: NEL
      my_real ,DIMENSION(NEL,2) ,INTENT(IN)    :: dir_a,dir_b
      my_real ,DIMENSION(NEL,3) ,INTENT(INOUT) :: tens
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I
      my_real :: r1,r2,s1,s2,t1,t2,t3
c=======================================================================
      DO i = 1,nel
        t1 = tens(i,1)
        t2 = tens(i,2)
        t3 = tens(i,3)
        r1 = dir_a(i,1)
        s1 = dir_a(i,2)
        r2 = dir_b(i,1)
        s2 = dir_b(i,2)
        tens(i,1) = r1*r1*t1 + s1*s1*t2 + two*r1*s1*t3        ! eps_x dir1
        tens(i,2) = r2*r2*t1 + s2*s2*t2 + two*r2*s2*t3        ! eps_y dir2           
        tens(i,3) = atan((r1*r2 + s1*s2) / (r1*s2 - r2*s1))   ! angle xy
      ENDDO
c-----------
      RETURN