33 use element_mod , only : nixt,nixr,nixp
34
35
36
37#include "implicit_f.inc"
38
39
40
41#include "com01_c.inc"
42#include "com04_c.inc"
43#include "param_c.inc"
44
45
46
47 INTEGER IPARG(NPARG,*),
48 . IXT(NIXT,*),IXP(NIXP,*),
49 . IXR(NIXR,*),EL2FA(*),INUM(*),
50 . NBF,INUMX1(*),BUF
51
52 INTEGER II(4),IE,NG, ITY, LFT, LLT, N, I, J,
53 . IPID, NEL, IAD, NFT, IMID,IALEL,MLW,
54 . NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,
55 . NBF2
56
57
58 nn1 = 1
59 nn2 = 1
60 nn3 = 1
61 nn4 = nn3
62 nn5 = nn4
63 nn6 = nn5
64 nn7 = nn6 + numelt
65 nn8 = nn7 + numelp
66 nn9 = nn8 + numelr
67 nn10= nn9
68
69 DO 490 ng=1,ngroup
70 mlw =iparg(1,ng)
71 nel =iparg(2,ng)
72 nft =iparg(3,ng)
73 iad =iparg(4,ng)
74 ity =iparg(5,ng)
75 lft=1
76 llt=nel
77
78
79
80 IF(ity==4)THEN
81 DO i=lft,llt
82 n = i + nft
83 inum(el2fa(nn6+n)) = ixt(nixt,n)
84 ENDDO
85
86
87
88 ELSEIF(ity==5)THEN
89 DO i=lft,llt
90 n = i + nft
91 inum(el2fa(nn7+n)) = ixp(nixp,n)
92 ENDDO
93
94
95
96 ELSEIF(ity==6)THEN
97 DO i=lft,llt
98 n = i + nft
99 inum(el2fa(nn8+n)) = ixr(nixr,n)
100 IF(mlw==3)inum(el2fa(nn8+n)+1) = ixr(nixr,n)
101 ENDDO
102
103 ELSE
104 ENDIF
105 490 CONTINUE
106
109
110
111 RETURN
void write_i_c(int *w, int *len)