anchor.diff:      if(allocated(this%lDmn)) deallocate(this%lDmn)
       if(allocated(this%lPrt)) deallocate(this%lPrt)
       if(allocated(this%NodeNeibBlock)) deallocate(this%NodeNeibBlock)
  +    if(allocated(this%Intf)) deallocate(this%Intf)
  +    if(allocated(this%IntfNeibBlock)) deallocate(this%IntfNeibBlock)
       if(allocated(this%PointPos)) deallocate(this%PointPos)
       if(allocated(this%LinePos)) deallocate(this%LinePos)
       if(allocated(this%FacetPos)) deallocate(this%FacetPos)
  
positive.code:         end do
        end do
        this%isUpNodeNeibBlock=.true.
      end if
    end subroutine
    
  
negative.code:   
    integer,parameter,public::TET_SURF_NUM=4 !< number of surfaces per tet
    integer,parameter,public::HEX_SURF_NUM=6 !< number of surfaces per hex
    
    integer,parameter,public::INTF_NEIB_BLOCK_NUM=2 !< number of neighbour blocks per interface
    
    integer,public::TET_SURF_TAB(TRI_NODE_NUM,TET_SURF_NUM) !< table of surface nodes for tet
    parameter(TET_SURF_TAB=reshape([1,3,2,1,2,4,1,4,3,2,3,4],[TRI_NODE_NUM,TET_SURF_NUM]))
  