anchor.diff:    !
     !-------------------------------------------------------------------------
   
  -  subroutine curve( temperature, s)
  +  subroutine curve(temperature, s)
       implicit none
   
       real(kind = x_precision),dimension(n_cell),intent(out)  :: temperature
  
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:     real(kind = x_precision), parameter                    :: dt    = (t_max-t_min) / (nb_it-1)
  
      real(kind = x_precision)                               :: eps   = 1.0d-4
  
      real(kind = x_precision)                               :: Sigma = 0.0d0
      real(kind = x_precision)                               :: Smin  = 0.0d0
  