anchor.diff:  
         deallocate(triangle)
         deallocate(vertex_xyz)
  +
  +      call enter_exit(sub_name,2)
   
         end subroutine make_data_grid
   ! 
  
positive.code:      !#This subroutine is called by make_data_grid subroutine.
       !#Application: to grow a grid in 2D surface.
       
  
       use arrays,only: dp,num_elems_2d,elem_nodes_2d
       integer,intent(in) :: surface_elems(:)
       integer,allocatable :: triangle(:,:)
  
negative.code: 
       use arrays,only: dp,data_xyz,data_weight,num_data
       use mesh_utilities,only: volume_internal_to_surface,point_internal_to_surface
       use diagnostics,only: enter_exit
       ! Parameters
       integer,intent(in) :: surface_elems(:)
       real(dp),intent(in) :: spacing
  