Get Surface points on boundary elements. For boundary state variable which are evaluated linkwise, extract surface points for each link and for nonlink based variables project barycenter on the boundary surface. Return real coordinates on boundary surface and offset bit which encodes direction.
Type  Intent  Optional  Attributes  Name  

type(boundary_type),  intent(in)  ::  bc  Field boundary type 

type(glob_boundary_type),  intent(in)  ::  globBC  for number of elements in boundary and position in buffer 

type(treelmesh_type),  intent(in)  ::  tree  global treelm mesh 

type(tem_stencilHeader_type),  intent(in)  ::  stencil  for directions 

integer(kind=long_k),  intent(in)  ::  total(:)  global pdf type 

integer,  intent(in)  ::  iLevel  Current level 

integer,  intent(out)  ::  nPoints  Number of points 

real(kind=rk),  intent(out),  allocatable  ::  points(:,:)  3d real coordinates on which boundary variables are evaluated 

character,  intent(out),  allocatable  ::  offset_bit(:)  Offset bit encodes direction of boundary. used by apesmate to translate space coordinate in the offset direction to determine the treeID in remote domain 
Type  Visibility  Attributes  Name  Initial  

real(kind=rk),  private  ::  dx  
real(kind=rk),  private  ::  bary(3)  
real(kind=rk),  private,  allocatable  ::  qVal(:)  
integer,  private  ::  iElem  
integer,  private  ::  iDir  
integer,  private  ::  iPnt  
integer,  private  ::  invDir  
logical,  private  ::  qValIsZero 