anchor.diff:  
              ! Do a snapshot
              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:         iteration = iteration + 1
  
          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
  