anchor.diff:          iteration = iteration + 1
   
           if (mod(iteration, output_freq) == 0 ) then
  -           if (arg == 'load' .and. iteration > stp_value) then
  +           if (arg == 'load' .and. iteration > n_iterations) then
                 stop
              endif
              call snapshot(s, iteration, t, 13)
  
positive.code: 
             ! Do a snapshot
             if (mod(iteration, output_freq) == 0) then
                if (arg == 'load' .and. iteration > stp_value) then
                   stop
                endif
                call snapshot(s, iteration, t, 13)
  
negative.code:   implicit none
  
    character(len=10)                      :: arg
    !  don't kill the job without stp_value
    real(x_precision), dimension(n_cell)   :: jnk1, jnk3
    character(len=100)                     :: line
  