anchor.diff:            enddo !nd
             
          endif !num_seeds_from_elem
  +!       write(*,'('' '',i6,'','')', advance = "no") num_seeds_from_elem(ne_min)
       enddo !N
       
       do N=1,num_next_parents
  
positive.code:     endif
  
      if(diagnostics_on)then
         write(*,'( ''  '',i6,'' seeds for element'',i7,'' and'',i6,'' for element'',i7)') &
              dat1,ne+1,dat2,ne+2
      endif
  
      call enter_exit(sub_name,2)
  
negative.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+zero_tol))then
  !                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)+zero_tol)then !keep seed points
  !             if(min_dist.lt.distance_limit/real(elem_ordrs(no_gen,ne_min),kind=dp))then !keep seed points
                  map_array(nd)=ne_min
               else
                  map_array(nd)=0 !too far from branch ends, so discard
               endif
  