anchor.diff:      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
  
positive.code:     integer, parameter                                     :: nb_it = 100
  
      real(kind = x_precision)                               :: eps   = 1.0d-4
      real(kind = x_precision)                               :: eps2  = 5.0d-1
  
  
      real(kind = x_precision)                               :: sigma = 0.0d0
  
negative.code:           temp_t_0(i)   =  temp 
            sigma_t_0(i)  =  sigma   
  
  
           ! do l          = 1, nb_it
           !    write(fid_1,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)')sigma_real_0(l),temp_real_0(l)
           ! enddo
            
         enddo
  
  
  