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:     case(1)
  
         Fz = 2._x_precision * c**2 * T**4 /(27._x_precision * sqrt(3._x_precision) &
              * (K_ff + K_e) * Sigma * Sigma_0 * Omega**(1._x_precision / 3._x_precision))
      case (0)
  
         Fz = 4._x_precision * rs * E_ff * H / (Omega_0 * Sigma_0)
  
negative.code:          sigma_real(i) = sigma_real_0(i)
           temp_real(i) = temp_real_1(i)
        enddo
        
      end subroutine build_s_curve
  
  
  