anchor.diff:  
     ! 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)
   
  
positive.code:   use mod_integrator
    use mod_output
    use mod_timestep
  
    implicit none
  
    integer                                     :: iteration, ios, i
    type(state)                                 :: s
  
negative.code:  ! real(kind = x_precision)                    :: t_nu, t_T
    real(kind = x_precision), dimension(n_cell) :: dt_nu, dt_T
    logical                                     :: T_converged
    real (kind = x_precision), dimension(n_cell):: dist
    
  
    ! FIXME
  