qcoor2.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "vect01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	qcoor2 (x, ixq, ngl, mxt, pid, ix1, ix2, ix3, ix4, y1, y2, y3, y4, z1, z2, z3, z4, sy, sz, ty, tz)

Function/Subroutine Documentation

qcoor2()

subroutine qcoor2	(		x,
		integer, dimension(nixq,*)	ixq,
		integer, dimension(*)	ngl,
		integer, dimension(*)	mxt,
		integer, dimension(*)	pid,
		integer, dimension(*)	ix1,
		integer, dimension(*)	ix2,
		integer, dimension(*)	ix3,
		integer, dimension(*)	ix4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			sy,
			sz,
			ty,
			tz )

Definition at line 33 of file qcoor2.F.

      use element_mod , only : nixq
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXQ(NIXQ,*), NGL(*),MXT(*)
      INTEGER PID(*), IX1(*), IX2(*), IX3(*), IX4(*)
      my_real
     .   x(3,*), y1(*), y2(*), y3(*), y4(*), z1(*), z2(*), z3(*), z4(*)
      my_real
     .   sy(mvsiz),sz(mvsiz),
     .   ty(mvsiz),tz(mvsiz)
      my_real
     .   suma
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C--------------------------------------------------
C     connectivities and material
C--------------------------------------------------
      DO i=lft,llt
      mxt(i)=ixq(1,i)
      ix1(i)=ixq(2,i)
      ix2(i)=ixq(3,i)
      ix3(i)=ixq(4,i)
      ix4(i)=ixq(5,i)
      pid(i)=ixq(6,i)
        ngl(i)=ixq(nixq,i)
      END DO
 
C----------------------------
C     NODE COORDINATES      |
C----------------------------
      DO i=lft,llt
      y1(i)=x(2,ix1(i))
      z1(i)=x(3,ix1(i))
      y2(i)=x(2,ix2(i))
      z2(i)=x(3,ix2(i))
      y3(i)=x(2,ix3(i))
      z3(i)=x(3,ix3(i))
      y4(i)=x(2,ix4(i))
      z4(i)=x(3,ix4(i))
      END DO
C---------------------------------------
C     LOCAL REFERENCE FRAME (ISOPARAMETRIC)
C---------------------------------------
      DO i=lft,llt
        sy(i)=half*(y2(i)+y3(i)-y1(i)-y4(i))
        sz(i)=half*(z2(i)+z3(i)-z1(i)-z4(i))
        ty(i)=half*(y3(i)+y4(i)-y1(i)-y2(i))
        tz(i)=half*(z3(i)+z4(i)-z1(i)-z2(i))
      END DO
c
      RETURN