anchor.diff:      end if
       if(toLower(dir)/='t')then
         allocate(coords(n))
  -      !/!\ need a correction
         do i = 1,n
  -        coords(i) = minCoord + real(i,cp)*delta
  +        coords(i) = minCoord + real(i-1,cp)*delta
         end do
       end if
       select case(toLower(dir))
  
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 select
      select case (toLower(domainType))
        case('d')!discrete information
          domainSize = (nPoint-1)*delta_space
          positive_size = (count(delta_space>epsilon(1.0))==nDimension)
        case('c')!continuous information
          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
  