anchor.diff:            prev_S = s%S
             
             ! Integrate S
  -          call do_timestep_S(s, dt_nu)
  +          
  +          call do_timestep_S(s, minval(dt_nu))
             
             ! Increment time, number of iterations
  -          t = t + dt_nu
  +          t = t + minval(dt_nu)
             iteration = iteration + 1
             ! Output things here
             if (mod(iteration, output_freq) == 0) then
  
positive.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)                    :: dt_nu, dt_T
    logical                                     :: T_converged
  
    ! FIXME
    S_crit = 1.e99_x_precision
  
negative.code:   ! we multiply delta to prevent the code from thinking it converged
    prev_S = 1.2*s%S
  
   ! t_T  = params%t_T  !* state_0%temps_0
   ! t_nu = params%t_nu !*  state_0%temps_0
  
   ! dt_nu = t_nu / cst_dt
   ! dt_T  = t_T / cst_dt
  
    write(*,*) 'dt_T, dt_nu:', dt_T, dt_nu
  
  