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#include "scr17_c.inc"
47
48
49
50 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
51 INTEGER IPARG(NPARG,*), IXT(NIXT,*), IXP(NIXP,*), IXR(NIXR,*),
52 . ID_ELEM(*), ITY_ELEM(*),IPART(LIPART1,*), IPARTT(*) ,IPARTP(*) ,IPARTR(*)
54 . oned_scalar(*)
55
56
57
58
59 INTEGER OFFSET,ITY,MLW,NEL,NFT,NG,I
60 TYPE(G_BUFEL_) ,POINTER :: GBUF
61
62
63 DO ng=1,ngroup
64 mlw =iparg(1,ng)
65 nel =iparg(2,ng)
66 ity =iparg(5,ng)
67 nft =iparg(3,ng)
68
69 IF (ity == 4) offset = 0
70 IF (ity == 5) offset = numelt
71 IF (ity == 6) offset = numelt+numelp
72
73 IF(ity == 4 .OR. ity == 5 .OR. ity == 6)THEN
74 DO i=1,nel
75 IF (ity == 4) THEN
76 id_elem(offset+nft+i) = ixt(nixt,nft+i)
77 ity_elem(offset+nft+i) = 4
78 ELSEIF (ity == 5) THEN
79 id_elem(offset+nft+i) = ixp(nixp,nft+i)
80 ity_elem(offset+nft+i) = 5
81 ELSEIF (ity == 6) THEN
82 id_elem(offset+nft+i) = ixr(nixr,nft+i)
83 ity_elem(offset+nft+i) = 6
84 ENDIF
85 ENDDO
86 ENDIF
87
88 IF(ity == 4 )THEN
89 IF(mlw == 0 .OR. mlw == 13)THEN
90 DO i=1,nel
91 oned_scalar(offset+nft+i) = one
92 ENDDO
93 ELSE
94 gbuf => elbuf_tab(ng)%GBUF
95 DO i=1,nel
96 oned_scalar(offset+nft+i) =
min(abs(gbuf%OFF(i)),one)
97 ENDDO
98 ENDIF
99 ELSEIF(ity == 5)THEN
100 IF(mlw == 0 .OR. mlw == 13)THEN
101 DO i=1,nel
102 oned_scalar(offset+nft+i) = one
103 ENDDO
104 ELSE
105 gbuf => elbuf_tab(ng)%GBUF
106 DO i=1,nel
107 oned_scalar(offset+nft+i) =
min(abs(gbuf%OFF(i)),one)
108 ENDDO
109 ENDIF
110 ELSEIF(ity == 6)THEN
111 IF(mlw == 0 .OR. mlw == 13)THEN
112 DO i=1,nel
113 oned_scalar(offset+nft+i) = one
114 ENDDO
115 ELSE
116 gbuf => elbuf_tab(ng)%GBUF
117 DO i=1,nel
118 oned_scalar(offset+nft+i) =
min(abs(gbuf%OFF(i)),one)
119 ENDDO
120 ENDIF
121 ENDIF
122 ENDDO
123 RETURN