anchor.diff:      end select
       select case (toLower(domainType))
         case('d')!discrete information
  -        domainSize = nPoint*delta_space
  -        positive_size = (count(delta_space>0.0)==nDimension)
  +        domainSize = (nPoint-1)*delta_space
  +        positive_size = (count(delta_space>epsilon(1.0))==nDimension)
         case('c')!continuous information
  -        positive_size = (count(domainSize>0.0)==nDimension)
  -        delta_space = domainSize/nPoint
  +        positive_size = (count(domainSize>epsilon(1.0))==nDimension)
  +        delta_space = domainSize/(nPoint-1)
         case default
           call stop_program('Bad type of domain type, should be "c" for continuous or "d" for discrete, check the namelist')
       end select
  
positive.code: 
  !> @brief Discretization module, define data structures and subroutines for finite difference-type discretization.
  !! @todo
  !! Update modules using the discretization module. A correction has been made on May 23rd 2014 by Innocent.
  !! Now, there is a discretization point at each end of the domain in each direction.
  !! Previously, there was no discretization point at the min coordinate, left boundary of the domain.
  !! The change affects coordinates only if the number of discretization points is set in the input file (namelist)
  !!
  !!
  !<
  MODULE discretization
  
negative.code:     end if
      if(toLower(dir)/='t')then
        allocate(coords(n))
        do i = 1,n
          coords(i) = minCoord + real(i-1,cp)*delta
        end do
      end if
      select case(toLower(dir))
  