anchor.diff:          iteration = iteration + 1
   
           if (mod(iteration, output_freq) == 0 ) then
  -           if (arg == 'load' .and. iteration > stp_value) then
  +           if ((arg == 'load' .or. arg == 'restart') &
  +                .and. iteration-start_iteration > stp_value) then
                 stop
              endif
              call snapshot(s, iteration, t, 13)
  
positive.code: 
          T_converged = .false.
  
          t_steps = 0
          do while (.not. T_converged)
             t_steps = t_steps + 1
             if (t_steps >= 20 .and. mod(t_steps, 20) == 0) then
                print*, 'W: T is taking time to converge, already', t_steps, 'iterations'
             end if
  
             call compute_variables(s)
  
negative.code:      s%S = s%S / state_0%S_0 * x_state%x
       s%Mdot = s%Mdot / state_0%Mdot_0
       close(21)
  
       if (maxval(s%T - T_c) < 0) then
         ! …or right from its critical S…
         unstable = maxval(s%S - S_c) > 0
  