31
32
33
34#include "implicit_f.inc"
35
36
37
38#include "com01_c.inc"
39#include "com04_c.inc"
40#include "sphcom.inc"
41#include "param_c.inc"
42
43
44
45
46 INTEGER IADD(*),IPARG(NPARG,*),
47 . MATER(*),EL2FA(*),
48 . KXSP(NISP,*),IPARTSP(*)
49
50
52 . off
53 INTEGER II, IE, NG, ITY, LFT, LLT, KPT, N, ,
54 . NEL, IAD, NPAR, NFT, IPRT,
55 . NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,
56 . JJ, INOD,
57 . NP(NUMSPH),BUF,BUFP
58
59
60 nn1 = 1
61 nn2 = nn1 + numsph+maxpjet
62 ie = 0
63
64
65
66 npar = 0
67 jj = 0
68
69 DO 100 iprt = 1 , npart
70
71 IF(mater(iprt)/=4) GOTO 100
72 npar = npar + 1
73 DO 90 ng=1,ngroup
74 nel =iparg(2,ng)
75 nft =iparg(3,ng)
76 iad =iparg(4,ng)
77 ity =iparg(5,ng)
78 lft=1
79 llt=nel
80 IF(ity==51)THEN
81
82
83
84 DO 20 i=lft,llt
85 n = i + nft
86 IF(ipartsp(n)/=iprt) GOTO 20
87 inod=kxsp(3,n)
88 ii = inod-1
90 jj = jj+1
91 ie = ie + 1
92 el2fa(nn1+n) = ie
93 20 CONTINUE
94 ELSE
95 ENDIF
96 90 CONTINUE
97
98
99
100 iadd(npar) = ie
101 100 CONTINUE
102
103
104 RETURN
void write_i_c(int *w, int *len)