anchor.diff:    end subroutine curve
   
   
  -
  +  
     !-------------------------------------------------------------------------
     ! Subroutine in order to find the first critical point
     !-------------------------------------------------------------------------
  
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:     integer                                                :: index_scp
      real(kind = x_precision)                               :: sigma_c_thin
      real(kind = x_precision)                               :: temp_c_thin
      
      !------------------------------------------------------------------------
  
      call display_parameters()
  