36
37
38
39 USE elbufdef_mod
40 USE my_alloc_mod
41 use element_mod , only : nixq
42
43
44
45#include "implicit_f.inc"
46
47
48
49#include "com01_c.inc"
50#include "com04_c.inc"
51#include "param_c.inc"
52#include "scr17_c.inc"
53#include "spmd_c.inc"
54
55
56
57 INTEGER ITAB(*), IPART(LIPART1,*),
58 . IGEO(NPROPGI,*), IXQ(NIXQ,*),
59 . IPARTQ(*), IPART_STATE(*),
60 . NODTAG(*), STAT_INDXQ(*),
61 . IPARG(NPARG,*),LENGQ
62 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
63
64
65
66 INTEGER I, N, JJ, IPRT
67 INTEGER NG, NEL, NFT, LFT, LLT, ITY, ISOLNOD
68 INTEGER,DIMENSION(:),ALLOCATABLE::IADD
69 INTEGER,DIMENSION(:,:),ALLOCATABLE::IADG
70 INTEGER,DIMENSION(:),ALLOCATABLE:
71INTEGER,DIMENSION(:),ALLOCATABLE::NPGLOB
72 INTEGER,DIMENSION(:,:),ALLOCATABLE:: CLEF
73 TYPE(G_BUFEL_) ,POINTER :: GBUF
74
75
76
77 CALL my_alloc(iadg,nspmd,npart)
78 CALL my_alloc(iadd,npart+1)
79 CALL my_alloc(np,24*numelqg)
80 CALL my_alloc(npglob,24*numelqg)
81 CALL my_alloc(clef,2,numelqg)
82
83 jj = 0
84 DO ng=1,ngroup
85 ity =iparg(5,ng)
86 isolnod = iparg(28,ng)
87 nel =iparg(2,ng)
88 nft =iparg(3,ng)
89 gbuf => elbuf_tab(ng)%GBUF
90 lft=1
91 llt=nel
92 IF(ity == 2) THEN
93 DO i=lft,llt
94 n = i + nft
95
96 iprt=ipartq(n)
97 IF(ipart_state(iprt)==0)cycle
98
99
100 nodtag(ixq(2,n))=1
101 nodtag(ixq(3,n))=1
102 nodtag(ixq(4,n))=1
103 nodtag(ixq(5,n))=1
104
105
106 END DO
107 END IF
108 END DO
109
110
111 DEALLOCATE(iadg)
112 DEALLOCATE(iadd)
113 DEALLOCATE(np)
114 DEALLOCATE(npglob)
115 DEALLOCATE(clef)
116
117 RETURN