anchor.diff:      !-------------------------------------------------------------------------
       ! Subroutine in order to find the second critical point
       !-------------------------------------------------------------------------
  -    subroutine second_critical_point(sigma_real_1, sigma_real_0, temp_real_1,&
  +    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)
       implicit none
       integer,intent(in)                                   :: nb_it
   
  -    real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_1
  -    real(kind = x_precision),dimension(nb_it),intent(in) :: temp_real_1
  -    real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_0
  -
  -    real(kind = x_precision)                 ,intent(in):: eps2
  +    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
  +
  +    real(kind = x_precision)                 ,intent(in) :: eps2
   
       integer                                  ,intent(in) :: index_fcp
       integer                                  ,intent(out):: index_scp
  
positive.code:       integer::i = 0
  
        do i = 1,index_scp
           sigma_real(i) = sigma_real_1(i)
           temp_real(i) = temp_real_1(i)
        enddo
        do i = index_scp + 1, nb_it
           sigma_real(i) = sigma_real_0(i)
           temp_real(i) = temp_real_1(i)
        enddo
  
      end subroutine build_s_curve
  
negative.code:             ! write(*,*)index_fcp
  
      ! do i = index_fcp, nb_it
      !   write(*,*)sigma_real_thin(i), sigma_real_thick(i)
      !   write(*,*)(dabs(sigma_real_thick(i) - sigma_real_thin(i+1)))
      ! end do
  
      do i = index_fcp, nb_it - 1
  
         if(dabs(sigma_real_thick(i) - sigma_real_thin(i+1)) .lt. eps2)then
            index_scp = i
            ! sigma_c_thin = (sigma_real_thick(i) + sigma_real_thin(i+1))/2._x_precision
            sigma_c_thin =  sigma_real_thin(i+1)
  
            temp_c_thin = temp_real_thick(i)
         end if
      enddo
  
  