anchor.diff:      endif
   
       if(diagnostics_on)then
  -       write(*,'( ''  '',i6,'' seeds for element'',i7,'' and'',i6,'' for element'',i7)') &
  -            dat1,ne+1,dat2,ne+2
  +       write(*,'( ''Parent '',i6,'' with '',i6,'' seeds for element'',i7,'' and'',i6,'' for element'',i7)') &
  +            ne1,dat1,ne+1,dat2,ne+2
       endif
   
       call enter_exit(sub_name,2)
  
positive.code:                 ne = local_parent(noelem)
                  np = elem_nodes(2,ne)
                  dist = distance_between_points(data_xyz(1,nd),node_xyz(1,np))
                  if(DIST.lt.MIN_DIST)then
                     ne_min = ne
                     MIN_DIST=DIST
                  endif
               enddo
               if(min_dist.lt.distance_limit/real(elem_ordrs(no_gen,ne_min),kind=dp))then !keep seed points
                  map_array(nd)=ne_min
                  if(ne_min.eq.2404)then
                     write(*,*) nd,min_dist
                  endif
               else
                  map_array(nd)=0 !too far from branch ends, so discard
               endif
  
negative.code:     deallocate(map_seed_to_space)
      deallocate(num_seeds_from_elem)
  
  !    call smooth_1d_tree(ne_start,length_limit)
      
      call enter_exit(sub_name,2)
      
  