anchor.diff:    use mod_integrator
     use mod_output
     use mod_timestep
  +  use mod_s_curve
  +  use mod_distance
   
     implicit none
  +
  +  !--------------------------- Parameters for s_curve-----------------------
  +  real(kind = x_precision), parameter         :: eps_in = 1.e-8_x_precision  ! Precision required for the dichotomy
  +  real(kind = x_precision), parameter         :: Tmin   = 2.5e-2_x_precision
  +  real(kind = x_precision), parameter         :: Tmax   = 4.49e0_x_precision
  +  real(kind = x_precision), parameter         :: Smin   = 2.36e1_x_precision
  +  real(kind = x_precision), parameter         :: Smax   = 2.36e3_x_precision
  +
  +  real(kind = x_precision), dimension(n_cell) :: temperature_c
  +  real(kind = x_precision), dimension(n_cell) :: sigma_c
  +  !-------------------------------------------------------------------------
   
     integer                                     :: iteration, ios, i
     type(state)                                 :: s
  
positive.code:           ! Recompute variables when the system is stable
            call compute_variables(s)
            call timestep (s,dt_T, dt_nu)
         end if
      else
         iteration = n_iterations
  
negative.code: 
    ! Read the parameters, generate state_0 and create adim state
    call get_parameters()
    call make_temperature(Tmin, Tmax)
    call set_conditions(eps_in, Smin, Smax)
    call curve(temperature_c, sigma_c)
    ! Initiate the S_curve
    ! call s_curve(foo, bar)
  
  