anchor.diff:        sigma_c_thick = sigma_real_thick(i)
         i             = i + 1
       enddo
  -    temp_c_thick  = temp_real_thick(i)
  +    temp_c_thick  = temperature(i)
   
     end subroutine first_critical_point
   
  
positive.code:     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
  
      ! For variables
      real(kind = x_precision)                    :: f              ! Q+ - Q-
  
negative.code:     ! For the separate curves given the thickness
      real(kind = x_precision), dimension(nb_it)  :: sigma_t_thick
      real(kind = x_precision), dimension(nb_it)  :: sigma_t_thin
      real(kind = x_precision), dimension(nb_it)  :: temp_real      ! Dimensioned temperature for the S curve
      real(kind = x_precision), dimension(nb_it)  :: sigma          ! Surface density for the S curve
      real(kind = x_precision), dimension(nb_it)  :: sigma_real     ! Dimensioned surface density for the S curve
  
      ! For the output file
      character(len = 8)                          :: number_of_cell
  