initwg_ressort.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_ressort (wd, pm, geo, ixr, igeo, numelr, ipm, nummat, numgeo, off)

Function/Subroutine Documentation

initwg_ressort()

subroutine initwg_ressort	(	real, dimension(*)	wd,
			pm,
			geo,
		integer, dimension(nixr,*)	ixr,
		integer, dimension(npropgi,*)	igeo,
		integer	numelr,
		integer, dimension(npropmi,*)	ipm,
		integer	nummat,
		integer	numgeo,
		integer	off )

Definition at line 30 of file initwg_ressort.F.

C-----------------------------------------------
C            M o d u l e s
C-----------------------------------------------
      USE ddweights_mod
      use element_mod , only : nixr
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  NUMELR,OFF,
     .  NUMMAT,NUMGEO,  IXR(NIXR,*), 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,ITY
     
      real
     .   wtype(3),fwihbe,fac8,
     .   tabmat(3),tabx(3),timmat,npt,telt,poids,w,
     .   batozmult,tmat   
      my_real
     .        cc, invtref,a,b,a1,a2
       DATA wtype /0.5 ,1.4 , 1.4/
!      different spring kind : wtype(1) = type 4    /SPRING
!                              wtype(2) = type 13   /SPR_BEAM
!                              wtype(3) = default   
C-----------------------------------------------
      DO i = 1, numelr
        pid= ixr(1,i)
        mtn= nint(geo(8,pid))
        issn = 0
        ity = igeo(11,pid)
         poids = wtype(3)
        IF(geo(5,pid)/=zero) issn=1
C   test elem delete
        IF (wd(i+off)==0.)
     .    wd(i+off) = poids
      ENDDO
      RETURN