anchor.diff:      itplCell2NodeVect(:,grid%nNode)=v(:,NlC(grid%nNode))
     end function
     
  -  !> interpolate scaler v within grid from cell to node (1-D)
  +  !> interpolate scalar v within grid from cell to node (1-D)
     function itplCell2NodeScal(v,grid)
       use moduleGrid1D
       use moduleBasicDataStruct
  
positive.code:     end forall
    end function
    
    !> interpolate scaler v within grid from node to cell (1-D)
    function itplNode2CellScal(v,grid)
      use moduleGrid1D
      use moduleBasicDataStruct
  
negative.code:     end forall
    end function
    
    !> interpolate scalar v within grid from node to edge
    function itplNode2EdgeScal(v,grid)
      use moduleGrid
      use moduleBasicDataStruct
  