anchor.diff:    !  t_V = 1.3e4_x_precision / state_0%temps_0 
     !  t_T = 0.72_x_precision / state_0%temps_0
   
  -  dt_V = t_V / 10._x_precision
  -  dt_T = t_T / 10._x_precision
  +  dt_V = t_V / 2._x_precision
  +  dt_T = t_T / 2._x_precision
   
     write(*,*) 'dt_T, dt_V:', dt_T, dt_V
  -  
  -  ! call snapshot(s, iteration, t, 13)
   
  +  iteration = 0
     ! Start iterating
  -  do iteration = 1, n_iterations
  +  do while iteration < n_iterations
        ! Check that S is at a fixed point
        if (maxval(abs((prev_S - s%S)/s%S)) > delta_S_max) then
           ! Check here that S < S_crit
  
positive.code:            ! Recompute variables
             call compute_variables(s)
  
             ! Iterate while T hasn't converged
             j = 0
             T_converged = .false.
             do while (.not. T_converged)
                ! Integrate T
  
negative.code:   ! FIXME
    S_crit = 1.e99_x_precision
    ! FIXME
    delta_S_max = 1e-5
    ! FIXME I love it
    delta_T_max = 1e-5
  
    ! Initial time = 0
    t = 0._x_precision
  