anchor.diff:    !Interfaces
     private
     public grow_tree,smooth_1d_tree
  +  public grow_tree_wrap
   
   contains
  +
  +!!!#############################################################################
  +  
  +  subroutine grow_tree_wrap(surface_elems,parent_ne,angle_max,angle_min,&
  +       branch_fraction,length_limit,shortest_length,rotation_limit)
  +    !temporary interface to the grow_tree subroutine until bindings sorted out 
  +    !DEC$ ATTRIBUTES DLLEXPORT,ALIAS:"SO_GROW_TREE_WRAP" :: GROW_TREE_WRAP
  +
  +!    use growtree,only: grow_tree
  +    
  +    integer,intent(in)  :: surface_elems(:)         ! list of surface elements defining the host region
  +    integer,intent(in)  :: parent_ne                ! stem branch that supplies 'parents' to grow from
  +    real(dp),intent(in) :: angle_max                ! maximum branch angle with parent; in degrees
  +    real(dp),intent(in) :: angle_min                ! minimum branch angle with parent; in degrees
  +    real(dp),intent(in) :: branch_fraction          ! fraction of distance (to COFM) to branch
  +    real(dp),intent(in) :: length_limit             ! minimum length of a generated branch (shorter == terminal)
  +    real(dp),intent(in) :: shortest_length          ! length that short branches are reset to (shortest in model)
  +    real(dp),intent(in) :: rotation_limit           ! maximum angle of rotation of branching plane
  +
  +    call grow_tree(surface_elems,parent_ne,angle_max,angle_min, &
  +         branch_fraction,length_limit,shortest_length,rotation_limit)
  +    
  +  end subroutine grow_tree_wrap
   
     !###############################################################
     !
  
positive.code:                          map_seed_to_elem(nd_min) = 0 ! remove seed point from list
                           map_seed_to_space(nd_min) = ne ! recording element number
  
                           if(to_export) then
                             write(40,*) nd_min,ne
                           endif
                        endif
                     endif ! .not.internal
  
  
negative.code:     real(dp),intent(in) :: length_limit             ! minimum length of a generated branch (shorter == terminal)
      real(dp),intent(in) :: shortest_length          ! length that short branches are reset to (shortest in model)
      real(dp),intent(in) :: rotation_limit           ! maximum angle of rotation of branching plane
  !    logical,intent(in) :: to_export                 ! option to export terminal element mapping to datapoints
  !    character(len=*),intent(in) :: filename
  
      !Local variables
      integer,allocatable :: local_parent(:)          ! stores current generation of local parent elements
  