initwg_truss.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "tablen_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	initwg_truss (wd, pm, geo, ixt, igeo, numelt, ipm, nummat, numgeo, off)

Function/Subroutine Documentation

initwg_truss()

subroutine initwg_truss	(	real, dimension(*)	wd,
			pm,
			geo,
		integer, dimension(nixt,*)	ixt,
		integer, dimension(npropgi,*)	igeo,
		integer	numelt,
		integer, dimension(npropmi,*)	ipm,
		integer	nummat,
		integer	numgeo,
		integer	off )

Definition at line 30 of file initwg_truss.F.

C-----------------------------------------------
C            M o d u l e s
C-----------------------------------------------
      USE ddweights_mod
      use element_mod , only : nixt
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      "param_c.inc"
#include      "tablen_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER  OFF,NUMELT,
     .  NUMMAT,NUMGEO,
     .  IXT(NIXT,*),IGEO(NPROPGI,*),IPM(NPROPMI,*)
 
C     REAL OR REAL*8
      my_real
     .    pm(npropm,*), geo(npropg,*)
      REAL WD(*)
C-----------------------------------------------
      INTEGER  NPN, MID, PID, JHBE, IGT, MLN,
     .    ISTRAIN, ITHK, IHBE, IPLA, ISSN, MTN, I, J, K,L,
     .    NFUNC,MPT,NPTS,NPTT,NPTR,NPTOT,IFLAG,JSROT,
     .    I_MID,I_PID,I_MID_OLD,I_PID_OLD,PUID,MUID,
     .    ELM_TYP,ELM_TYP_OLD,ILAW,ILAW_OLD,TEST_MAT,
     .    I_PRO,ISOL2,MID_OLD,PID_OLD,MUID_OLD,PUID_OLD,
     .    TEST,NFUNC1,NFUNC2
 
      real
     .   wtype(9),fwihbe,fac8,
     .   tabmat(3),tabx(3),timmat,npt,telt,poids,w,
     .   batozmult,tmat
      
      my_real
     .        cc, invtref,a,b,a1,a2
       DATA wtype /1.6 ,1. ,1. ,.9 ,1.1 ,1.4 ,0.65 ,.9 ,2.0/     
C-----------------------------------------------  
      DO i = 1, numelt
        mid= ixt(1,i)
        mln= pm(19,mid)
        pid= ixt(4,i)
        issn = 0
        IF(geo(5,pid)/=zero) issn=1
C   test elem delete
        IF (wd(i+off)==0.)
     .    wd(i+off) = wtype(4)
      ENDDO
      RETURN