r2buf3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	r2buf3 (off, geo, x, x0, y0, z0, ix, skew, iposx, iposy, iposz, iposxx, iposyy, iposzz, igeo, skew_id)

Function/Subroutine Documentation

r2buf3()

subroutine r2buf3	(		off,
			geo,
			x,
			x0,
			y0,
			z0,
		integer, dimension(nixr,*)	ix,
			skew,
			iposx,
			iposy,
			iposz,
			iposxx,
			iposyy,
			iposzz,
		integer, dimension(npropgi,*)	igeo,
		integer, dimension(*)	skew_id )

Definition at line 29 of file r2buf3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IX(NIXR,*),IGEO(NPROPGI,*),SKEW_ID(*)
C     REAL
      my_real
     .   off(*), geo(npropg,*), x(3,*), x0(*), y0(*), z0(*),
     .   skew(lskew,*),iposx(5,*) ,iposy(5,*),
     .   iposz(5,*),iposxx(5,*),iposyy(5,*),iposzz(5,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, NG, I1, I2, ISK, K, USENS
C     REAL
      my_real
     .   x1, y1, z1, ex, ey, ez
C-----------------------------------------------
      DO i=lft,llt
        j=i+nft
        usens=igeo(3,ix(1,j))
        IF(usens<=0)THEN
C         no sensor or Isflag=1
          off(i)=one
        ELSE
          off(i)=-ten
        END IF
      ENDDO
C
      DO j=1,5
        DO i=lft,llt
          iposx(j,i)=zero
          iposy(j,i)=zero
          iposz(j,i)=zero
          iposxx(j,i)=zero
          iposyy(j,i)=zero
          iposzz(j,i)=zero
        ENDDO
      ENDDO
C
      DO i=lft,llt
        j=i+nft
        ng=ix(1,j)
        isk=skew_id(i)
        i1=ix(2,j)
        i2=ix(3,j)
        x1=x(1,i2)-x(1,i1)
        y1=x(2,i2)-x(2,i1)
        z1=x(3,i2)-x(3,i1)
        ex=skew(1,isk)
        ey=skew(2,isk)
        ez=skew(3,isk)
        x0(i)=ex*x1+ey*y1+ez*z1
        ex=skew(4,isk)
        ey=skew(5,isk)
        ez=skew(6,isk)
        y0(i)=ex*x1+ey*y1+ez*z1
        ex=skew(7,isk)
        ey=skew(8,isk)
        ez=skew(9,isk)
        z0(i)=ex*x1+ey*y1+ez*z1
      ENDDO
C-----------------------------------------------
      RETURN