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:     end if
      if(toLower(dir)/='t')then
        allocate(coords(n))
        !/!\ need a correction
        do i = 1,n
          coords(i) = minCoord + real(i,cp)*delta
        end do
      end if
      select case(toLower(dir))
  
negative.code: 
  !> @brief Discretization module, define data structures and subroutines for finite difference-type discretization.
  !! @todo
  !! Check the behavior of modules that uses the discretization module. A correction has been made on May 23rd 2014 by Innocent.
  !! The modification should affect function that computes something as a function of the physical coordinates given by the discretization parameter
  !! In the normal case, it should not affects computations that does not used physical coordinates directly.
  !! Now, there is a discretization point at each end of the domain in each direction. The number of discretization points if 1 unit greater than the number of cells.
  !! Previously, there was no discretization point at the min coordinate, The number of discretization points corresponded to the number of cells.
  !!
  !<
  MODULE discretization
  