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 IPARG(NPARG,*),
46 . IXT(NIXT,*),IXP(NIXP,*),
47 . IXR(NIXR,*),EL2FA(*),INUM(*),
48
49
50 INTEGER II(4),IE,NG, ITY, LFT, LLT, N, I, J,
51 . IPID, NEL, IAD, , IMID,IALEL,MLW,
52 . NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,
53
54
55
56 nn1 = 1
57 nn2 = 1
58 nn3 = 1
59 nn4 = nn3
60 nn5 = nn4
61 nn6 = nn5
62 nn7 = nn6 + numelt
63 nn8 = nn7 + numelp
64 nn9 = nn8 + numelr
65 nn10= nn9
66
67 DO 490 ng=1,ngroup
68 mlw =iparg(1,ng)
69 nel =iparg(2,ng)
70 nft =iparg(3,ng)
71 iad =iparg(4,ng)
72 ity =iparg(5,ng)
73 lft=1
74 llt=nel
75
76
77
78 IF(ity==4)THEN
79 DO i=lft,llt
80 n = i + nft
81 inum(el2fa(nn6+n)) = ixt(nixt,n)
82 ENDDO
83
84
85
86 ELSEIF(ity==5)THEN
87 DO i=lft,llt
88 n = i + nft
89 inum(el2fa(nn7+n)) = ixp(nixp,n)
90 ENDDO
91
92
93
94 ELSEIF(ity==6)THEN
95 DO i=lft,llt
96 n = i + nft
97 inum(el2fa(nn8+n)) = ixr(nixr,n)
98 IF(mlw==3)inum(el2fa(nn8+n)+1) = ixr(nixr,n)
99 ENDDO
100
101 ELSE
102 ENDIF
103 490 CONTINUE
104
107
108
109 RETURN
void write_i_c(int *w, int *len)