anchor.diff:      real(kind = x_precision), dimension(n_cell), intent(out) :: temperature_c  ! Temperature for the critical point
       real(kind = x_precision), dimension(n_cell), intent(out) :: sigma_c        ! Surface density for the critical point
       !------------------------------------------------------------------------
  -    integer                                     :: i,j,k          ! Iteration counters
  +    integer                                     :: i,k            ! Iteration counters
   
       ! For variables
       real(kind = x_precision)                    :: f              ! Q+ - Q-
  
positive.code:       temperature(i) = dt * (i-1) + Tmin
      enddo
  
    end subroutine make_temperature
  
  
  
negative.code:   implicit none
  
    real(kind = x_precision), dimension(nb_it) :: temperature ! Temperatures for which we need to solve the dichotomy
    real(kind = x_precision), dimension(nb_it) :: temp_real   ! Dimensioned version
    real(kind = x_precision)                   :: eps         ! Precision required for the dichotomy
    real(kind = x_precision)                   :: S_min       ! Minimum surface density limit
    real(kind = x_precision)                   :: S_max       ! Maximum surface density limit
  