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: 
             ! Do a snapshot
             if (mod(iteration, output_freq) == 0) then
                if ((arg == 'load' .or. arg == 'restart') &
                     .and. iteration-start_iteration > stp_value) then
                   stop
                endif
                call snapshot(s, iteration, t, 13)
  