anchor.diff:      
       if(allocated(this%NodePos)) deallocate(this%NodePos)
       if(allocated(this%CellPos)) deallocate(this%CellPos)
  +    if(allocated(this%CellWidth)) deallocate(this%CellWidth)
     end subroutine
     
     !> destructor of typeGrid1D
  
positive.code:     h=(boundR-boundL)/dble(nCell)
      allocate(this%NodePos(this%nNode))
      allocate(this%CellPos(this%nCell))
      this%NodePos(1)=boundL
      do i=1,this%nCell
        this%NodePos(i+1)=this%NodePos(i)+h
        this%CellPos(i)=this%NodePos(i)+h/2d0
      end do
    end subroutine
    
    !> cell left node
  
negative.code:     double precision,allocatable::NodePos(:) !< node position
      integer::nCell !< number of cells
      double precision,allocatable::CellPos(:) !< cell position
      double precision,allocatable::CellWidth(:) !< cell width
    contains
      procedure,public::init=>initGrid1D
      procedure,public::clear=>clearGrid1D
  