anchor.diff:      if(.not.cross_any)then
          point_internal_to_surface = .true.
       else
  -       point_internal_to_surface = .false.
  +       if(ncrossed.eq.2)then
  +          point_internal_to_surface = .true.
  +       else
  +          point_internal_to_surface = .false.
  +       endif
       endif
   
     end function point_internal_to_surface
  
positive.code:     real(dp) :: POINT1(3),POINT2(3),POINT3(3)
      !     Local variables
      real(dp) :: ERR1(3),ERR2(3),LU,LV,U(3),V(3)
      logical :: check_colinear_points
      
      
  
negative.code: !!! Any function that is used by more than one module should appear in here. 
  !!! ALL subroutines and functions in this module are public.
  
    use other_consts
    use arrays, only: dp,zero_tol
    implicit none
  
    private
  
    public  area_between_three_points,area_between_two_vectors,calc_branch_direction,&
         angle_btwn_points,angle_btwn_vectors,calc_scale_factors_2d,check_colinear_points,cross_product,&
         distance_between_points,make_plane_from_3points,mesh_a_x_eq_b,ph3,pl1,&
         point_internal_to_surface,scalar_product_3,scalar_triple_product,scale_mesh,&
         unit_norm_to_plane_two_vectors,unit_norm_to_three_points,unit_vector,&
  