anchor.diff:      ! 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, nb_it, eps2)
  +         index_fcp, index_scp, sigma_c_thin, temp_c_thin, nb_it)
       implicit none
       integer,intent(in)                                   :: nb_it
   
  
positive.code:        call first_critical_point(sigma_t_thick, temp_t_thick, index_fcp,sigma_c_thick, temp_c_thick, nb_it)
  
         call second_critical_point(sigma_t_thick, sigma_t_thin, temp_t_thick,&
           index_fcp, index_scp, sigma_c_thin, temp_c_thin, nb_it, eps2)
  
         call build_s_curve(sigma_t_thick, sigma_t_thin, temp_t_thick, nb_it, index_scp, sigma_real, temp_real)
  
  
negative.code:   !-------------------------------------------------------------------------
    ! Subroutine in order to find the first critical point
    !-------------------------------------------------------------------------
    subroutine first_critical_point(sigma_real_thick, temp_real_thick, index_fcp, sigma_c_thick, temp_c_thick, nb_it)
      implicit none
      integer,intent(in)                                   :: nb_it
  
      real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_thick
      real(kind = x_precision),dimension(nb_it),intent(in) :: temp_real_thick
  
      integer                                  ,intent(out):: index_fcp
      real(kind = x_precision)                 ,intent(out):: sigma_c_thick
      real(kind = x_precision)                 ,intent(out):: temp_c_thick
      integer                                              :: i
      i = 1
      do while (sigma_real_thick(i) < sigma_real_thick(i+1) .and. i < nb_it - 1)
         index_fcp     = i
         sigma_c_thick = sigma_real_thick(i)
         i             = i + 1
      enddo
      temp_c_thick  = temp_real_thick(i)
  
      endsubroutine first_critical_point
  
  