anchor.diff:      itplBlock2NodeScal(:)=vrst(1,:)
     end function
     
  +  !> interpolate vector v within grid from block to facet
  +  function itplBlock2FacetVect(v,grid)
  +    use moduleGrid
  +    use moduleBasicDataStruct
  +    double precision,intent(in)::v(:,:) !< block data to be interpolated
  +    type(typeGrid),intent(inout)::grid !< grid on which v is defined
  +    double precision,allocatable::itplBlock2FacetVect(:,:) !< interpolated data on Facet
  +    
  +    call grid%updateFacetNeib()
  +    call reallocArr(itplBlock2FacetVect,size(v,1),grid%nFacet)
  +    forall(i=1:grid%nFacet)
  +      itplBlock2FacetVect(:,i)=v(:,maxval(grid%FacetNeibBlock(:,i)))
  +    end forall
  +  end function
  +  
  +  !> interpolate scalar v within grid from block to facet
  +  function itplBlock2FacetScal(v,grid)
  +    use moduleGrid
  +    use moduleBasicDataStruct
  +    double precision,intent(in)::v(:) !< block data to be interpolated
  +    type(typeGrid),intent(inout)::grid !< grid on which v is defined
  +    double precision,allocatable::itplBlock2FacetScal(:) !< interpolated data on facet
  +    double precision vv(1,size(v))
  +    double precision,allocatable::vrst(:,:)
  +    
  +    vv(1,:)=v(:)
  +    vrst=itplBlock2FacetVect(vv,grid)
  +    call reallocArr(itplBlock2FacetScal,size(vrst,2))
  +    itplBlock2FacetScal(:)=vrst(1,:)
  +  end function
  +  
     !> interpolate matrix v within grid from block to interface using block center direction scheme
     function itplBCDMat(v,grid)
       use moduleGrid
  
positive.code:     module procedure::itplNode2BlockScal
    end interface
    public itplNode2Block
    
    !> interpolate from node to edge
    interface itplNode2Edge
  
negative.code:     module procedure::itplBlock2NodeScal
    end interface
    public itplBlock2Node
    
    !> interpolate from block to facet
    interface itplBlock2Facet
      module procedure::itplBlock2FacetVect
      module procedure::itplBlock2FacetScal
    end interface
    public itplBlock2Facet
    
    !> interpolate from block to interface
    interface itplBlock2Intf
  