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:            add_nodes(cRANK,i) == add_nodes(cRANK,i-1-ndup)) then
          ndup = ndup + 1
        else if (ndup > 0) then
          add_nodes(1:3,i-ndup) = add_nodes(1:3,i)
        endif
      enddo
    end subroutine uniq_add_nodes
  
negative.code:       n_send = n_send + 1
        tag = 3002
        call HECMW_ISEND_INT(exp_cols_item(1,exp_cols_index(idom-1)+1), &
             cNCOL_ITEM * BT_exp(idom)%nnz, irank, tag, hecCOMM%HECMW_COMM, requests(n_send))
        n_send = n_send + 1
        tag = 3003
        call HECMW_ISEND_R(BT_exp(idom)%A, ndof2 * BT_exp(idom)%nnz, &
  