anchor.diff:             ! Recompute variables
              call compute_variables(s)
   
  -           ! FIXME : increment time
  -           print *, "Before integration of T", j, maxval(abs(prev_T - s%T)/s%T)
              ! Iterate while T hasn't converged
              j = 0
              do while (maxval(abs(prev_T - s%T)/s%T) > delta_T_max)
  
positive.code:   ! we add initially 2*delta to prevent the code from thinking it converged
    prev_S = s%S + 2*delta_S_max
    prev_T = s%T + 2*delta_T_max
  
    call snapshot(s, iteration, t, 13)
    
  
negative.code:   delta_T_max = 1e-2
  
    ! Initial time = 0
    t = 0._x_precision
  
    ! Read the parameters 
    call get_parameters()
  