anchor.diff:        endif
       enddo
     end subroutine map_present_nodes
  -
  -  function find_external_gid(hecMESH, gid)
  -    implicit none
  -    integer(kind=kint) :: find_external_gid
  -    type (hecmwST_local_mesh), intent(in) :: hecMESH
  -    integer(kind=kint), intent(in) :: gid
  -    integer(kind=kint) :: i
  -    find_external_gid = -1
  -    do i = hecMESH%nn_internal+1, hecMESH%n_node
  -      if (hecMESH%global_node_ID(i) == gid) then
  -        find_external_gid = i
  -        return
  -      endif
  -    enddo
  -  end function find_external_gid
   
     function find_external_node_ID(hecMESH, lid, rank)
       implicit none
  
positive.code:   subroutine sort_and_uniq_add_nodes(n_add_node, add_nodes)
      implicit none
      integer(kind=kint), intent(inout) :: n_add_node
      integer(kind=kint), intent(inout) :: add_nodes(3,n_add_node)
      integer(kind=kint) :: ndup
      call sort_add_nodes(add_nodes, 1, n_add_node)
      call uniq_add_nodes(add_nodes, n_add_node, ndup)
  
negative.code:       nnz = idx_n - idx_0
        if (nnz == 0) cycle
        tag = 3002
        call HECMW_RECV_INT(imp_cols_item(1, idx_0 + 1), cNCOL_ITEM * nnz, &
             irank, tag, hecCOMM%HECMW_COMM, statuses(:,1))
        tag = 3003
        call HECMW_RECV_R(imp_vals_item(ndof2*idx_0 + 1), ndof2 * nnz, &
  