bulkigeo3.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	bulkige3 (volu, nctrl, mat, pm, vnige, bnige, px, py, pz, nctrlmax)

Function/Subroutine Documentation

bulkige3()

subroutine bulkige3	(		volu,
		integer	nctrl,
		integer, dimension(*)	mat,
			pm,
			vnige,
			bnige,
		integer	px,
		integer	py,
		integer	pz,
		integer	nctrlmax )

Definition at line 28 of file bulkigeo3.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 NCTRL, NCTRLMAX, MAT(*), PX,PY,PZ
C     REAL
      my_real
     .   vnige(nctrlmax,*),bnige(nctrlmax,*),volu(*),
     .   pm(npropm,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, N
C     REAL
      my_real
     .   av1,cv1
C-----------------------------------------------------------------------
C
      IF(nctrl==0) RETURN
C 
C     isogeometric elements :
      DO i=lft,llt
        av1 = volu(i)/nctrl
        cv1 = pm(32,mat(i))*av1
        DO n=1,nctrl
          vnige(n,i)=av1
          bnige(n,i)=cv1
        ENDDO
      ENDDO
 
Cnote: the volume can be distributed differently depending on whether the points are at the corners, in the middle
cor on an edge, with PX, PY, PZ you can find the position of each one
CThis is the elements 20 neoud, there are 3/192 * flight to corners and 14/192 * Vol at the backgrounds
C
      RETURN