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