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:   !-------------------------------------------------------------------------
    ! Subroutine in order to find the second critical point
    !-------------------------------------------------------------------------
    subroutine second_critical_point(sigma_real_thick, sigma_real_thin, temp_real_thick,&
           index_fcp, index_scp, sigma_c_thin, temp_c_thin)
      implicit none
  
      real(kind = x_precision), dimension(nb_it), intent(in)  :: sigma_real_thick
      real(kind = x_precision), dimension(nb_it), intent(in)  :: temp_real_thick
      real(kind = x_precision), dimension(nb_it), intent(in)  :: sigma_real_thin
      integer,                                    intent(in)  :: index_fcp
  
      integer,                                    intent(out) :: index_scp
  
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
  