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,NUMELT,
49 . NUMMAT,NUMGEO,
50 . IXT(NIXT,*),IGEO(NPROPGI,*),IPM(NPROPMI,*)
51
52
54 . pm(npropm,*), geo(npropg,*)
55 REAL WD(*)
56
57 INTEGER NPN, MID, PID, JHBE, IGT, MLN,
58 . ISTRAIN, ITHK, IHBE, IPLA, ISSN, MTN, I, J, K,L,
59 . NFUNC,MPT,NPTS,NPTT,NPTR,NPTOT,IFLAG,JSROT,
60 . I_MID,I_PID,I_MID_OLD,I_PID_OLD,PUID,MUID,
61 . ELM_TYP,ELM_TYP_OLD,ILAW,ILAW_OLD,TEST_MAT,
62 . I_PRO,ISOL2,MID_OLD,PID_OLD,MUID_OLD,PUID_OLD,
63 . TEST,NFUNC1,NFUNC2
64
65 real
66 . wtype(9),fwihbe,fac8,
67 . tabmat(3),tabx(3),timmat,npt,telt,poids,w,
68 . batozmult,tmat
69
71 . cc, invtref,a,b,a1,a2
72 DATA wtype /1.6 ,1. ,1. ,.9 ,1.1 ,1.4 ,0.65 ,.9 ,2.0/
73
74 DO i = 1, numelt
75 mid= ixt(1,i)
76 mln= pm(19,mid)
77 pid= ixt(4,i)
78 issn = 0
79 IF(geo(5,pid)/=zero) issn=1
80
81 IF (wd(i+off)==0.)
82 . wd(i+off) = wtype(4)
83 ENDDO
84 RETURN