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