anchor.diff:  
   
   
  - !-------------------------------------------------------------------------
  -  ! Subroutine in order to find the second critical point
  -  !-------------------------------------------------------------------------
  +    !-------------------------------------------------------------------------
  +    ! Subroutine in order to find the second critical point
  +    !-------------------------------------------------------------------------
       subroutine second_critical_point(sigma_real_1, sigma_real_0, temp_real_1,&
            index_fcp, index_scp, sigma_c_thin, temp_c_thin, nb_it, eps2)
       implicit none
  
positive.code: 
         call build_s_curve(sigma_real_1, sigma_real_0, temp_real_1,nb_it, index_scp, sigma_real, temp_real)
  
  
            do l          = 1, nb_it
               write(fid_3,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)') sigma_real(l), temp_real(l)
            enddo
         
         
  
negative.code:   subroutine curve()
      implicit none
      integer                                                :: i     = 0
      !integer                                                :: j     = 0
      integer                                                :: k     = 0
      integer                                                :: l     = 0
  
      real(kind = x_precision)                               :: temp  = 0.0d0
      real(kind = x_precision), parameter                    :: t_min = 5.0d-1
      real(kind = x_precision), parameter                    :: t_max = 4.49d0
  
      integer,                  parameter                    :: nb_it = 100
      real(kind = x_precision)                               :: eps   = 1.0d-4
  