1 | include(dom.inc)
2 |
3 | SUBROUTINE GATHER_FACES(vect, vect_face, ndata)
4 |
5 | USE mod_pmm
6 | USE mod_slave
7 |
8 | IMPLICIT NONE
9 |
10 | ! IN
11 | DOM_INT :: ndim, ndata
12 | DOM_REAL, DIMENSION(ndata,is_nnodes) :: vect
13 |
14 | ! LOCAL
15 | DOM_INT :: ibnd, k, iface, inode
16 |
17 | ! OUT
18 | DOM_REAL,DIMENSION(ndata,is_nbfaces) :: vect_face
19 |
20 | vect_face = 0.
21 |
22 | DO ibnd = 1, is_nbfaces
23 | iface = is_bface(ibnd)
24 |
25 | DO k=1, is_facenodes(1,iface)
26 | inode = is_facenodes(k+1,iface)
27 | vect_face(:,ibnd) = vect_face(:,ibnd) + vect(:,inode)
28 | ENDDO
29 |
30 | vect_face(:,ibnd) = vect_face(:,ibnd) / &
31 | & real(is_facenodes(1,iface))
32 |
33 | ENDDO
34 |
35 | END SUBROUTINE GATHER_FACES
gather_faces.F could be called by: