35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
54
55
56
57
58
59
60#include "implicit_f.inc"
61
62
63
64#include "com04_c.inc"
65
66
67
68
69 TYPE(GROUP_),INTENT(IN)::IGRNOD(NGRNOD)
70 TYPE(SUBMODEL_DATA),INTENT(IN)::LSUBMODEL(*)
71
72 INTEGER,INTENT(OUT)::NUM
73
74
75
76 INTEGER I,NN,IGU(NUMELX),IGS,STAT
77 INTEGER, DIMENSION(:), ALLOCATABLE :: SUB_XELEM
78
79
80
81 INTEGER NGR2USRN
82
83
84
85
86 ALLOCATE (sub_xelem(numelx),stat=stat)
87 IF (stat /= 0)
CALL ancmsg(msgid=268,anmode=aninfo,
88 . msgtype=msgerror,
89 . c1='SUB_XELEM')
90 sub_xelem(1:numelx) = 0
91
92 num = 0
93 CALL cpp_xelem_preread(igu,sub_xelem)
94
95 DO i=1,numelx
96 igs =
ngr2usrn(igu(i),igrnod,ngrnod,nn)
97 num = num + nn
98 ENDDO
99
100 RETURN
integer function ngr2usrn(iu, igrnod, ngrnod, num)
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)