h3d_shell_vector.F File Reference

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

Go to the source code of this file.

Functions/Subroutines

subroutine

h3d_shell_vector (elbuf_tab, shell_vector, ifunc, iparg, geo, ixq, ixc, ixtg, pm, el2fa, nbf, iadp, nbf_l, ehour, anim, nbpart, iadg, ipm, igeo, thke, err_thk_sh4, err_thk_sh3, invert, x, v, w, nv46, nercvois, nesdvois, lercvois, lesdvois, stack, id_elem, ity_elem, info1, info2, is_written_shell, ipartc, iparttg, layer_input, ipt_input, ply_input, gauss_input, iuvar_input, h3d_part, keyword, d, multi_fvm)

Function/Subroutine Documentation

h3d_shell_vector()

subroutine h3d_shell_vector	(	type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
			shell_vector,
		integer	ifunc,
		integer, dimension(nparg,ngroup)	iparg,
			geo,
		integer, dimension(nixq,numelq)	ixq,
		integer, dimension(nixc,numelc)	ixc,
		integer, dimension(nixtg,numeltg)	ixtg,
			pm,
		integer, dimension(*)	el2fa,
		integer	nbf,
		integer, dimension(*)	iadp,
		integer	nbf_l,
			ehour,
			anim,
		integer	nbpart,
		integer, dimension(nspmd,*)	iadg,
		integer, dimension(npropmi,nummat)	ipm,
		integer, dimension(npropgi,numgeo)	igeo,
			thke,
			err_thk_sh4,
			err_thk_sh3,
		integer, dimension(*)	invert,
			x,
			v,
			w,
		integer	nv46,
		integer, dimension(*)	nercvois,
		integer, dimension(*)	nesdvois,
		integer, dimension(*)	lercvois,
		integer, dimension(*)	lesdvois,
		type (stack_ply)	stack,
		integer, dimension(*)	id_elem,
		integer, dimension(*)	ity_elem,
		integer	info1,
		integer	info2,
		integer, dimension(*)	is_written_shell,
		integer, dimension(*)	ipartc,
		integer, dimension(*)	iparttg,
		integer	layer_input,
		integer	ipt_input,
		integer	ply_input,
		integer	gauss_input,
		integer	iuvar_input,
		integer, dimension(*)	h3d_part,
		character(len=ncharline100)	keyword,
			d,
		type (multi_fvm_struct), intent(in)	multi_fvm )

Definition at line 38 of file h3d_shell_vector.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE initbuf_mod
      USE elbufdef_mod
      USE schlieren_mod
      USE stack_mod
      USE multi_fvm_mod
      USE names_and_titles_mod, ONLY: ncharline100
      use element_mod , only : nixq,nixc,nixtg
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      "mvsiz_p.inc"
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .   shell_vector(3,*),x(3,numnod),v(3,numnod),w(3,numnod),d(3,numnod),thke(*),ehour(*),geo(npropg,numgeo),
     .   anim(*),pm(npropm,nummat),err_thk_sh4(*), err_thk_sh3(*)
      INTEGER IPARG(NPARG,NGROUP),IXC(NIXC,NUMELC),IXTG(NIXTG,NUMELTG),EL2FA(*),
     .   IXQ(NIXQ,NUMELQ),IFUNC,NBF,
     .   IADP(*),NBF_L, NBPART,IADG(NSPMD,*),IPM(NPROPMI,NUMMAT),
     .   IGEO(NPROPGI,NUMGEO),INVERT(*), NV46,ID_ELEM(*),ITY_ELEM(*),
     .   INFO1,INFO2,IS_WRITTEN_SHELL(*),IPARTC(*),IPARTTG(*),H3D_PART(*),
     .   LAYER_INPUT ,IPT_INPUT,GAUSS_INPUT,PLY_INPUT,IUVAR_INPUT
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
      TYPE (STACK_PLY) :: STACK
      CHARACTER(LEN=NCHARLINE100)::KEYWORD
      TYPE (MULTI_FVM_STRUCT), INTENT(IN) :: MULTI_FVM
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   evar(mvsiz),dam1(mvsiz),dam2(mvsiz),
     .   wpla(mvsiz),dmax(mvsiz),wpmax(mvsiz),fail(mvsiz),
     .   epst1(mvsiz),epst2(mvsiz),epsf1(mvsiz),epsf2(mvsiz),
     .   value(mvsiz),vg(5),vly(5),ve(5),mass(mvsiz)
      my_real
     .   off, p,vonm2,s1,s2,s12,s3,dmgmx,fac,
     .   dir1_1,dir1_2,dir2_1,dir2_2,aa,bb,v1,v2,v3,x21,x32,x34,
     .   x41,y21,y32,y34,y41,z21,z32,z34,z41,suma,vr,vs,x31,y31,
     .   z31,e11,e12,e13,e21,e22,e23,sum,area,x2l,var,
     .   e1x,e1y,e1z,e2x,e2y,e2z,e3x,e3y,e3z,rx,ry,rz,sx,sy,sz,
     .   rho0,thk0,xx1,xx2,xx3,yy1,yy2,yy3,zz1,zz2,zz3,a0
      INTEGER I, NG, 
     .        IR,IS,IT,IL,MLW, NUVAR,IUS,LENF,PTF,PTM,PTS,NFAIL,
     .        N,NN,K,K1,K2,JTURB,MT,IMID,IALEL,IPID,ISH3N,NNI,
     .        NN1,NN2,NN3,NN4,NN5,NN6,NN9,NF,BUF,NVARF,
     .        OFFSET,IHBE,NPTM,NPG, MPT,IPT,IADD,IADR,IPMAT,IFAILT,
     .        IIGEO,IADI,ISUBSTACK,ITHK,NERCVOIS(*),NESDVOIS(*),
     .        LERCVOIS(*),LESDVOIS(*),ID_PLY,NB_PLYOFF,IOK,N1,N2,N3,N4
      INTEGER 
     .        PTE(4),PTP(4),PTMAT(4),PTVAR(4),NPT_ALL,IPLY,
     .        ID_ELEM_TMP(MVSIZ),NIX,IOK_PART(MVSIZ),JJ(5),NPGT,IUVAR,
     .        IS_WRITTEN_VALUE(MVSIZ)
 
C-----------------------------------------------
c      print *,'LAYER_INPUT',LAYER_INPUT
c      print *,'IPT_INPUT',IPT_INPUT
c      print *,'GAUSS_INPUT',GAUSS_INPUT
c      print *,'PLY_INPUT',PLY_INPUT
c      print *,'IUVAR_INPUT',IUVAR_INPUT
c      print *,'KEYWORD',KEYWORD,'KEYWORD'
C
      DO i=1,numelc+numeltg
         is_written_shell(i) = 0
      ENDDO
c a corriger
      nn3 = 0
C
      DO ng=1,ngroup
       CALL h3d_shell_vector_1(
     .                  elbuf_tab   ,shell_vector,ifunc   ,iparg       ,geo        ,
     .                  ixq         ,ixc       ,ixtg      ,pm         ,
     .                  el2fa       ,nbf       ,iadp        ,
     .                  nbf_l       ,ehour     ,anim      ,nbpart      ,iadg       ,
     .                  ipm         ,igeo      ,thke      ,err_thk_sh4 ,err_thk_sh3,
     .                  invert      ,x         ,v         ,w           ,
     .                  nv46        ,nercvois  ,nesdvois  ,lercvois    ,lesdvois,
     .                  stack       ,id_elem   ,ity_elem  ,info1       ,info2      , 
     .                  is_written_shell,ipartc,iparttg   ,layer_input ,ipt_input  ,
     .                  ply_input   ,gauss_input,iuvar_input,h3d_part  ,keyword    ,
     .                  d           ,ng         ,multi_fvm)
      ENDDO
C-----------------------------------------------
      RETURN