initwg_quad.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_quad (wd, pm, geo, ixq, igeo, numelq, ipm, off)

Function/Subroutine Documentation

initwg_quad()

subroutine initwg_quad	(	real, dimension(*)	wd,
			pm,
			geo,
		integer, dimension(nixq,*)	ixq,
		integer, dimension(npropgi,*)	igeo,
		integer	numelq,
		integer, dimension(npropmi,*)	ipm,
		integer	off )

Definition at line 30 of file initwg_quad.F.

C-----------------------------------------------
C            M o d u l e s
C-----------------------------------------------
      USE ddweights_mod
      use element_mod , only : nixq
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  
     .  NUMELQ,
     .  NUMMAT,NUMGEO,
     .  IXQ(NIXQ,*),IGEO(NPROPGI,*),
     .  IPM(NPROPMI,*),OFF
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 invtelt_pro,telt_pro
      
      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, numelq
        npn=1
        mid= ixq(1,i)
        pid= ixq(6,i)
        mln = nint(pm(19,abs(mid)))
c        IF(MID<0)THEN
c          IF(MLN==6.AND.JPOR/=2)MLN=17
c          IF(MLN==46)MLN=47
c        ENDIF
C   test elem delete
        IF (wd(i+off)==0.)
     .    wd(i+off) = wtype(2)
      ENDDO
    
      RETURN