anchor.diff:    call set_conditions(eps_in, Smin, Smax)
     call curve(temperature_c, sigma_c)
   
  -  ! Conversion into T*_crit and S*_crit
  -  temperature_c  = temperature_c / state_0%T_0
  -  sigma_c = sigma_c / state_0%S_0 * x_state%x
  +  ! Conversion into S*_crit
  +  S_c = sigma_c * x_state%x
   
     ! Initiate the S_curve
     ! call s_curve(foo, bar)
  
positive.code: 
    real(kind = x_precision), dimension(n_cell) :: temperature_c
    real(kind = x_precision), dimension(n_cell) :: sigma_c
    !-------------------------------------------------------------------------
  
    integer                                     :: iteration, ios, i
  
negative.code:   logical                                     :: T_converged
    real (kind = x_precision), dimension(n_cell):: dist
  
    real (kind = x_precision), dimension(n_cell):: dist_crit
    real (kind = x_precision)                   :: dt_pre_factor
  
    ! FIXME
    S_crit = 1.e99_x_precision
  