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 INTEGER IPARG(NPARG,*),
46 . IXS(NIXS,*),EL2FA(*),INUM(*),
47 . NBF,KXSP(NISP,*)
48 INTEGER ISPH3D ,NG, ITY, LFT, LLT, N, I, J, NEL, IAD, NFT,
49 . NN1,NN2,NN3,ISOLNOD
50
51 nn1 = 1
52 nn2 = 1
53 nn3 = nn2 + numels
54
55 DO 490 ng=1,ngroup
56 nel =iparg(2,ng)
57 nft =iparg(3,ng)
58 iad =iparg(4,ng)
59 ity =iparg(5,ng)
60 isolnod = iparg(28,ng)
61 lft=1
62 llt=nel
63
64
65
66 IF(ity==1.AND.isolnod==16)THEN
67 DO i=lft,llt
68 n = i + nft
69 IF(el2fa(nn2+n)/=0)THEN
70 inum(el2fa(nn2+n)) = ixs(nixs,n)
71 inum(el2fa(nn2+n)+1) = ixs(nixs,n)
72 inum(el2fa(nn2+n)+2) = ixs(nixs,n)
73 inum(el2fa(nn2+n)+3) = ixs(nixs,n)
74 ENDIF
75 ENDDO
76
77
78
79 ELSEIF(ity==1)THEN
80 DO i=lft,llt
81 n = i + nft
82 IF(el2fa(nn2+n)/=0)THEN
83 inum(el2fa(nn2+n)) = ixs(nixs,n)
84 ENDIF
85 ENDDO
86 ELSEIF(isph3d==1.AND.ity==51)THEN
87
88
89
90 DO i=lft,llt
91 n = i + nft
92 IF(el2fa(nn3+n)/=0)THEN
93 inum(el2fa(nn3+n)) = kxsp(nisp,n)
94 ENDIF
95 ENDDO
96 ELSE
97 ENDIF
98 490 CONTINUE
99
101
102 RETURN
void write_i_c(int *w, int *len)