s10volnodt3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com04_c.inc"
#include "scr05_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s10volnodt3 (elbuf_tab, volnod6, x, nc, offg, sav, xdp, nel, npt)

Function/Subroutine Documentation

s10volnodt3()

subroutine s10volnodt3	(	type (elbuf_struct_), target	elbuf_tab,
		double precision, dimension(6,2*numnod)	volnod6,
			x,
		integer, dimension(mvsiz,10), intent(in)	nc,
			offg,
		double precision, dimension(nel,30)	sav,
		double precision, dimension(3,sxdp/3), intent(in)	xdp,
		integer, intent(in)	nel,
		integer, intent(in)	npt )

Definition at line 35 of file s10volnodt3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod            
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"
#include      "param_c.inc"
#include      "com04_c.inc"
#include      "scr05_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NPT
      INTEGER, INTENT(IN) ::NEL
      INTEGER, DIMENSION(MVSIZ,10),INTENT(IN) :: NC
      double precision
     .   volnod6(6,2*numnod),sav(nel,30)
      DOUBLE PRECISION , DIMENSION(3,SXDP/3), INTENT(IN) :: XDP     
C     REAL
      my_real
     .   x(3,numnod),offg(nel)
      TYPE (ELBUF_STRUCT_), TARGET  :: ELBUF_TAB
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, K,N1,N2,NN,N
      INTEGER IP,IBID,NPE
      INTEGER IPERM1(10),IPERM2(10)
      DATA iperm1/0,0,0,0,1,2,3,1,2,3/
      DATA iperm2/0,0,0,0,2,3,1,4,4,4/
 
      my_real
     .  mfxx(mvsiz),mfxy(mvsiz),mfxz(mvsiz),
     .  mfyx(mvsiz),mfyy(mvsiz),mfyz(mvsiz),
     .  mfzx(mvsiz),mfzy(mvsiz),mfzz(mvsiz),
     .  px(mvsiz,10),py(mvsiz,10),pz(mvsiz,10),jac(mvsiz),
     .  vx0(mvsiz,10),vy0(mvsiz,10),vz0(mvsiz,10),volp(mvsiz)
      DOUBLE PRECISION 
     .   XX(MVSIZ,10), YY(MVSIZ,10), ZZ(MVSIZ,10),
     .   X0(MVSIZ,10), Y0(MVSIZ,10), Z0(MVSIZ,10),
     .   DET6(6,MVSIZ)
C-----
      TYPE(L_BUFEL_) ,POINTER :: LBUF     
C----------------------------
C     NODAL COORDINATES     
C----------------------------
C----------------------------
      IF(iresp == 1) THEN
        DO n=1,10
          DO i=1,nel               
            nn = max(1,nc(i,n))    
            xx(i,n)=xdp(1,nn)          
            yy(i,n)=xdp(2,nn)          
            zz(i,n)=xdp(3,nn)          
          ENDDO                      
        ENDDO
      ELSE
        DO n=1,10
          DO i=1,nel               
            nn = max(1,nc(i,n))    
            xx(i,n)=x(1,nn)          
            yy(i,n)=x(2,nn)          
            zz(i,n)=x(3,nn)          
          ENDDO                      
        ENDDO
      END IF !(IRESP == 1) THEN
C-------exact for Itetra4=1  
      DO n=5,10
        n1=iperm1(n)
        n2=iperm2(n)
        DO i=1,nel
          IF(nc(i,n)==0)THEN
            xx(i,n) = half*(xx(i,n1)+xx(i,n2))
            yy(i,n) = half*(yy(i,n1)+yy(i,n2))
            zz(i,n) = half*(zz(i,n1)+zz(i,n2))
          ENDIF
        ENDDO
      ENDDO
C----------------------------
C     NODAL COORDINATES INITIALES    |
C----------------------------
      npe = 10
      DO n=1,npe
       DO i=1,nel
        x0(i,n)=sav(i,n)
        y0(i,n)=sav(i,n+npe)
        z0(i,n)=sav(i,n+2*npe)
       ENDDO
      ENDDO 
C----------------------------
C     DISPLACEMENT     |
C----------------------------
        n = npe
        DO i=1,nel
        vx0(i,n)=xx(i,n)-x0(i,n)
        vy0(i,n)=yy(i,n)-y0(i,n)
        vz0(i,n)=zz(i,n)-z0(i,n)
        ENDDO
       DO n=1,npe
        DO i=1,nel
        vx0(i,n)=xx(i,n)-x0(i,n)-vx0(i,npe)
        vy0(i,n)=yy(i,n)-y0(i,n)-vy0(i,npe)
        vz0(i,n)=zz(i,n)-z0(i,n)-vz0(i,npe)
        ENDDO
       ENDDO 
      ibid = 1      
      DO ip=1,npt
         lbuf => elbuf_tab%BUFLY(ibid)%LBUF(ip,ibid,ibid)
         CALL s10pijto3(px,py,pz,lbuf%PIJ,nel)
         CALL s10defot3(
     1   px,      py,      pz,      vx0,
     2   vy0,     vz0,     mfxx,    mfxy,
     3   mfxz,    mfyx,    mfyy,    mfyz,
     4   mfzx,    mfzy,    mfzz,    nel)
         CALL jacob_j33(
     1   jac,     mfxx,    mfxy,    mfxz,
     2   mfyx,    mfyy,    mfyz,    mfzx,
     3   mfzy,    mfzz,    nel)
        DO i=1,nel
         IF (offg(i) == zero .OR. abs(offg(i))>one) THEN
           volp(i)= zero
         ELSE
           volp(i) = jac(i)*lbuf%VOL(i)
         ENDIF
        ENDDO
        
        !Parith-On treatment
        CALL foat_to_6_float(1  ,nel  ,volp ,det6 )
        
        DO i=1,nel
          volnod6(1:6,nc(i,ip)) = volnod6(1:6,nc(i,ip)) + det6(1:6,i)
        ENDDO
C      
      END DO !IP=1,NPT
C
      RETURN