32
33
34
35 USE elbufdef_mod
36
37
38
39#include "implicit_f.inc"
40
41
42
43#include "com01_c.inc"
44#include "com04_c.inc"
45#include "param_c.inc"
46
47
48
49 INTEGER IPARG(NPARG,*),EL2FA(*),NBF,IOFF(*),
50 . NBPART,IOFFX1(*)
51 TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
52
53
54
55 INTEGER I, NG, NEL, NFT, ITY, LFT, NPT,
56 . N, J, LLT, N1,N2,N3,N4,N_FA,MLW,
57 . NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NBF2,RBUF
58 TYPE(G_BUFEL_) ,POINTER :: GBUF
59
60 nn1 = 1
61 nn2 = 1
62 nn3 = 1
63 nn4 = nn3
64 nn5 = nn4
65 nn6 = nn5
66 nn7 = nn6 + numelt
67 nn8 = nn7 + numelp
68
69 DO ng=1,ngroup
70 mlw =iparg(1,ng)
71 nel =iparg(2,ng)
72 ity =iparg(5,ng)
73 nft =iparg(3,ng)
74 lft=1
75 llt=nel
76
77 gbuf => elbuf_tab(ng)%GBUF
78
79
80
81 IF (ity ==4) THEN
82 n_fa = nn6 + nft
83
84
85
86 ELSEIF (ity == 5) THEN
87 n_fa = nn7 + nft
88
89
90
91 ELSEIF (ity == 6) THEN
92 n_fa = nn8 + nft
93 ELSE
94 ity=0
95 ENDIF
96
97 IF (ity /= 0) THEN
98 IF (mlw == 0) THEN
99
100
101
102 DO i=lft,llt
103 ioff(el2fa(n_fa+i)) = 0
104 ENDDO
105 ELSEIF (ity == 6 .AND. mlw == 3) THEN
106
107
108
109 DO i=lft,llt
110 ioff(el2fa(n_fa+i)) = nint(
min(gbuf%OFF(i),one))
111 ioff(el2fa(n_fa+i)+1) = nint(
min(gbuf%OFF(i),one))
112 ENDDO
113
114 ELSEIF (ity == 1 .or. ity == 2 .or. ity == 51 .or.
115 . ity == 3 .or. ity == 7 .or. ity == 4 .or.
116 . ity == 5 .or. ity == 6) THEN
117 DO i=lft,llt
118 ioff(el2fa(n_fa+i)) = nint(
min(gbuf%OFF(i),one))
119 ENDDO
120
121 ENDIF
122 ENDIF
123
124 ENDDO
125
127
129
130 RETURN
void write_c_c(int *w, int *len)