32
33
34
35#include "implicit_f.inc"
36
37
38
39#include "com01_c.inc"
40#include "com04_c.inc"
41#include "param_c.inc"
42
43
44
45 INTEGER IXC(NIXC,*),IXTG(NIXTG,*),IPARG(NPARG,*),
46 . IXQ(NIXQ,*),EL2FA(*),INUM(*),
47 . NBF,NELCUT,NBPART,IDCMAX
48
49
50
51
53 . off
54 INTEGER IE, NG, , LFT, LLT, , I, J,
55 . NEL, IAD, NFT,
56 . NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9
57
58
59 nn1 = 1
60 nn2 = 1
61 nn3 = 1
62 nn4 = nn3 + numelq
63 nn5 = nn4 + numelc
64 nn6 = nn5 + numeltg
65 nn7 = nn6
66 nn8 = nn7
67 nn9 = nn8
68
69 DO ng=1,ngroup
70 nel =iparg(2,ng)
71 nft =iparg(3,ng)
72 iad =iparg(4,ng)
73 ity =iparg(5,ng)
74 lft=1
75 llt=nel
76
77
78
79 IF(ity==2)THEN
80 DO i=lft,llt
81 n = i + nft
82 inum(el2fa(nn3+n)) = ixq(nixq,n)
83 ENDDO
84
85
86
87 ELSEIF(ity==3)THEN
88 DO i=lft,llt
89 n = i + nft
90 inum(el2fa(nn4+n)) = ixc(nixc,n)
91 ENDDO
92
93
94
95 ELSEIF(ity==7)THEN
96 DO i=lft,llt
97 n = i + nft
98 inum(el2fa(nn5+n)) = ixtg(nixtg,n)
99 ENDDO
100 ELSE
101 ENDIF
102 ENDDO
103
104
105 idcmax=0
106 DO i=1,nbf
107 idcmax=
max(idcmax,inum(i))
108 ENDDO
110 DO i=1,nelcut
112 ENDDO
113
114 RETURN
void write_i_c(int *w, int *len)