anchor.diff:                prev_T = s%T
   
                 ! Integrate T
  -              call do_timestep_T(s)
  +              call do_timestep_T(s, dt_T)
  +              ! Increment time
  +              t = t + dt_V
                 ! Recompute variables
                 call compute_variables(s)
  -              print *, "After integration of T", j, maxval(abs(prev_T - s%T)/s%T)
                 call snapshot(s, iteration, t, 13)
                 j = j+1
   
                 ! FIXME : increment time
              end do
              ! Output things here
  -           write(13, *) "YOLO DUDE"
  -           call snapshot(s, iteration, t, 13)
           end if
        end do
        
  
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, prev_T, S_crit
  
    ! FIXME
    S_crit = 1.e99_x_precision
    ! FIXME
    delta_S_max = 1e-2
    ! FIXME I love it
    delta_T_max = 1e-2
    
  
    ! Initial time = 0
    t = 0
  
negative.code:            prev_S = s%S
  
             ! Integrate S
             call do_timestep_S(s, dt_V)
             ! Increment time
             t = t + dt_V
             
             ! Recompute variables
             call compute_variables(s)
  