32
33
34
35 use element_mod , only : nixr
36
37
38
39#include "implicit_f.inc"
40
41
42
43#include "com04_c.inc"
44#include "spmd_c.inc"
45#include "com01_c.inc"
46#include "param_c.inc"
47#include "sphcom.inc"
48
49
50
51 INTEGER IGEO(NPROPGI,*),IXR(NIXR,*),
52 . NSTRF(*)
53
54
55
56 INTEGER NLOCAL
58
59
60
61 INTEGER I,P,NNOD,NSEGC,NSEGTG,NSEGS,K0
62
63
64 numnodg = numnod
65
66
67 numelqg = numelq
68 numelcg = numelc
69 numeltgg = numeltg
70 numelsg = numels
71 numelpg = numelp
72 numelrg = numelr
73 numeltrg = numelt
74 numels10g = numels10
75 numels16g = numels16
76 numels20g = numels20
77 numsphg = numsph
78 numelxg = numelx
79
80 numnodm = 0
81
82 nrbe2g = nrbe2
83
84 nrbe3g = nrbe3
85
86 DO p = 1,nspmd
87 nnod=0
88 DO i = 1, numnod
90 nnod = nnod + 1
91 ENDIF
92 ENDDO
93 numnodm =
max(numnodm,nnod)
94 ENDDO
95
96
97
98 nb1dg = numelt + numelr + numelp
99 nskewag = numelp + numelt + numskw
100
101 DO i=1,numelr
102 IF(igeo(11,ixr(1,i))==4 .OR.
103 . igeo(11,ixr(1,i))==26.OR.
104 . igeo(11,ixr(1,i))==13.OR.
105 . igeo(11,ixr(1,i))==25.OR.
106 . igeo(11,ixr(1,i))==45.OR.
107 . igeo(11,ixr(1,i))>=29.AND.
108 . igeo(11,ixr(1,i))<=33 .OR.
109 . igeo(11,ixr(1,i))==23)THEN
110
111 nskewag=nskewag+1
112
113 ELSEIF(igeo(11,ixr(1,i))==12)THEN
114
115 nskewag=nskewag+2
116 nb1dg =nb1dg+1
117 ENDIF
118
119 ENDDO
120
121
122
123 sectiom = 0
124 IF (nsect>0) THEN
125 k0 = nstrf(25)
126
127 DO i=1,nsect
128
129
130 nsegc = nstrf(k0+9)
131 nsegtg = nstrf(k0+13)
132 nsegs=nstrf(k0+7)
133 sectiom =
max(sectiom,nsegc,nsegtg,nsegs)
134 k0 = nstrf(k0+24)
135 ENDDO
136 ENDIF
137
138 RETURN