i12rot3.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i12rot3 (xs, xm, x, skew, iskwn, ilev, iref, icen, nsn, nmn, tets, tetm)

Function/Subroutine Documentation

i12rot3()

subroutine i12rot3	(		xs,
			xm,
			x,
			skew,
		integer, dimension(liskn,*)	iskwn,
		integer	ilev,
		integer	iref,
		integer	icen,
		integer	nsn,
		integer	nmn,
			tets,
			tetm )

Definition at line 28 of file i12rot3.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 ILEV,IREF,ICEN,ISKWN(LISKN,*),NSN,NMN
      my_real
     .   xs(3,*),xm(3,*),x(3,*),skew(lskew,*),tets(*),tetm(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,ORIG
      my_real
     . xx(3),xo(3),xr(3),rr,tet,p(9)
C
      IF(ilev==1)THEN
C        ORIG=ISKWN(1,IREF+1)
        IF(icen>0)THEN
          xo(1)=x(1,icen)
          xo(2)=x(2,icen)
          xo(3)=x(3,icen)
        ELSE
          xo(1)=zero
          xo(2)=zero
          xo(3)=zero
        ENDIF
        DO i=1,9
          p(i)=skew(i,iref+1)
        ENDDO
        DO i=1,nsn
          xx(1)=xs(1,i)-xo(1)
          xx(2)=xs(2,i)-xo(2)
          xx(3)=xs(3,i)-xo(3)
          xr(1)=p(1)*xx(1)+p(2)*xx(2)+p(3)*xx(3)
          xr(2)=p(4)*xx(1)+p(5)*xx(2)+p(6)*xx(3)
          xr(3)=p(7)*xx(1)+p(8)*xx(2)+p(9)*xx(3)
          rr=sqrt(xr(2)**2+xr(3)**2)
          IF(xr(3)/=zero.AND.xr(2)/=zero) THEN
            tet=atan2(xr(3),xr(2))
          ELSE
            tet=zero
          ENDIF
          xs(1,i)=xr(1)
          xs(2,i)=rr
          xs(3,i)=tet
          tets(i)=tet
        ENDDO
        DO i=1,nmn
          xx(1)=xm(1,i)-xo(1)
          xx(2)=xm(2,i)-xo(2)
          xx(3)=xm(3,i)-xo(3)
          xr(1)=p(1)*xx(1)+p(2)*xx(2)+p(3)*xx(3)
          xr(2)=p(4)*xx(1)+p(5)*xx(2)+p(6)*xx(3)
          xr(3)=p(7)*xx(1)+p(8)*xx(2)+p(9)*xx(3)
          rr=sqrt(xr(2)**2+xr(3)**2)
          IF(xr(3)/=zero.AND.xr(2)/=zero) THEN
            tet=atan2(xr(3),xr(2))
          ELSE
            tet=zero
          ENDIF
          xm(1,i)=xr(1)
          xm(2,i)=rr
          xm(3,i)=tet
          tetm(i)=tet
        ENDDO
      ENDIF
 
      RETURN