anchor.diff:            end do
             ! Recompute variables when the system is stable
             call compute_variables(s)
  +          call timestep (s,dt_T, dt_nu)
          end if
       else
          iteration = n_iterations
  
positive.code:   !use mod_s_curve
    use mod_integrator
    use mod_output
  
    implicit none
  
  
negative.code:   type(state)                                 :: s
    real(kind = x_precision)                    :: delta_S_max, delta_T_max, t
    real(kind = x_precision), dimension(n_cell) :: prev_S, S_crit
   ! real(kind = x_precision)                    :: t_nu, t_T
    real(kind = x_precision), dimension(n_cell) :: dt_nu, dt_T
    logical                                     :: T_converged
    
  
    ! FIXME
    S_crit = 1.e99_x_precision
  