32
33
34
36
37
38
39#include "implicit_f.inc"
40
41
42
43#include "param_c.inc"
44#include "tablen_c.inc"
45
46
47
48 INTEGER OFF,NUMELP,
49 . NUMMAT,NUMGEO, IXP(NIXP,*), IGEO(NPROPGI,*),
50 . IPM(NPROPMI,*)
51
53 . pm(npropm,*), geo(npropg,*)
54 REAL (*)
55
56 INTEGER NPN, MID, PID, JHBE, IGT, MLN,
57 . ISTRAIN, ITHK, IHBE, IPLA, ISSN, MTN, I, J, K,L,
58 . NFUNC,MPT,NPTS,NPTT,NPTR,NPTOT,IFLAG,JSROT,
59 . ,I_PID,I_MID_OLD,I_PID_OLD,PUID,MUID,
60 . ELM_TYP,ELM_TYP_OLD,ILAW,ILAW_OLD,TEST_MAT,
61 . I_PRO,ISOL2,MID_OLD,PID_OLD,MUID_OLD,PUID_OLD,
62 . TEST,NFUNC1,NFUNC2
63
64 real
65 . wtype(9),fwihbe,fac8,
66 . tabmat(3),tabx(3),timmat,npt,telt,poids,w,
67 . batozmult,tmat
68
70 . cc, invtref,a,b,a1,a2
71 DATA wtype /1.6 ,1. ,1. ,.9 ,1.1 ,1.4 ,0.65 ,.9 ,2.0/
72
73 DO i = 1, numelp
74 mid= ixp(1,i)
75 mln= pm(19,mid)
76 pid= ixp(5,i)
77 issn = 0
78 IF(geo(5,pid)/=zero) issn=1
79
80 IF (wd(i+off)==0.)
81 . wd(i+off) = wtype(5)
82 ENDDO
83 RETURN