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