anchor.diff:      integer                                                :: index_scp
       real(kind = x_precision)                               :: sigma_c_thin
       real(kind = x_precision)                               :: temp_c_thin
  -
  +    
       !------------------------------------------------------------------------
   
       call display_parameters()
  
positive.code:          sigma_real(i) = sigma_real_0(i)
           temp_real(i) = temp_real_1(i)
        enddo
  
      end subroutine build_s_curve
  
  
  
negative.code:     endsubroutine second_critical_point
  
  
      
      subroutine build_s_curve(sigma_real_1, sigma_real_0, temp_real_1,nb_it, index_scp, sigma_real, temp_real)
        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) :: sigma_real_0
        real(kind = x_precision),dimension(nb_it),intent(in) :: temp_real_1
  