setup_owned_faces Subroutine

   subroutine setup_owned_faces(mesh, flow)
      type(mesh_t), intent(in) :: mesh
      type(flow_mpi_t), intent(inout) :: flow
      integer :: f, n

      n = 0
      do f = 1, mesh%nfaces
         if (flow%owned(mesh%faces(f)%owner)) n = n + 1
      end do

      allocate(flow%owned_faces(n))
      n = 0
      do f = 1, mesh%nfaces
         if (.not. flow%owned(mesh%faces(f)%owner)) cycle
         n = n + 1
         flow%owned_faces(n) = f
      end do
   end subroutine setup_owned_faces