q4coor2.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	q4coor2 (x, ncp, y1, y2, y3, y4, z1, z2, z3, z4, nc1, nc2, nc3, nc4, ngl, mat, ngeo, vd2, vis, v, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, yavg, ay, exx, nel, jhbe)

Function/Subroutine Documentation

q4coor2()

subroutine q4coor2	(		x,
		integer, dimension(7,*)	ncp,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
		integer, dimension(*)	ngl,
		integer, dimension(*)	mat,
		integer, dimension(*)	ngeo,
			vd2,
			vis,
			v,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
			yavg,
			ay,
			exx,
		integer, intent(in)	nel,
		integer, intent(in)	jhbe )

Definition at line 28 of file q4coor2.F.

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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: JHBE
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*), MAT(*), NGL(*), NGEO(*),
     .        NCP(7,*)
C     REAL
      my_real
     .   x(3,*),y1(*),y2(*),y3(*),y4(*),z1(*),z2(*),z3(*),z4(*),vd2(*),
     .    vy1(*),vy2(*),vy3(*),vy4(*),
     .    vz1(*),vz2(*),vz3(*),vz4(*),
     .    v(3,*),yavg(*),exx(*),ay(*),vis(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,II,IV
      my_real
     .   vyg(mvsiz)
C-----------------------------------------------
      DO i=1,nel
      ngeo(i)=ncp(6,i)
      ngl(i)=ncp(7,i)
      mat(i)=ncp(1,i)
      nc1(i)=ncp(2,i)
      nc2(i)=ncp(3,i)
      nc3(i)=ncp(4,i)
      nc4(i)=ncp(5,i)
      END DO
C
C----------------------------
C     COORDONNEES NODALES
C----------------------------
      DO i=1,nel
      y1(i)=x(2,nc1(i))
      z1(i)=x(3,nc1(i))
      y2(i)=x(2,nc2(i))
      z2(i)=x(3,nc2(i))
      y3(i)=x(2,nc3(i))
      z3(i)=x(3,nc3(i))
      y4(i)=x(2,nc4(i))
      z4(i)=x(3,nc4(i))
      END DO
C
      DO i=1,nel
        vy1(i)=v(2,nc1(i))
        vz1(i)=v(3,nc1(i))
        vy2(i)=v(2,nc2(i))
        vz2(i)=v(3,nc2(i))
        vy3(i)=v(2,nc3(i))
        vz3(i)=v(3,nc3(i))
        vy4(i)=v(2,nc4(i))
        vz4(i)=v(3,nc4(i))
      END DO
C-----Exx ia calculated anyway in global system
      IF(n2d==1) THEN
        DO i=1,nel
          yavg(i) = y1(i)+y2(i)+y3(i)+y4(i)
          vyg(i) = vy1(i)+vy2(i)+vy3(i)+vy4(i)
        ENDDO
C-----JHBE=17 -> petrov G, constant Exx, JHBE=22-> Bubnov
         DO i=1,nel
          ay(i) =one/max(em20,yavg(i))
          exx(i) = vyg(i)*ay(i)
         ENDDO
      ENDIF
      DO i=1,nel
        vd2(i)=zero
        vis(i)=zero
      ENDDO
C
      RETURN