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