anchor.diff:          iteration = iteration + 1
   
           if (mod(iteration, output_freq) == 0 ) then
  +           if (arg == 'start' .and. iteration > stp_value) then
  +              stop
  +           endif
              call snapshot(s, iteration, t, 13)
              print*,'snapshot', iteration, t, 'S', dt_nu
           end if
  
positive.code: 
    real(x_precision), dimension(n_cell) :: dist_crit
  
    !----------------------------------------------
    ! Convergence criteria for S and T
    !----------------------------------------------
  
negative.code: 
             ! Do a snapshot
             if (mod(iteration, output_freq) == 0) then
                if (arg == 'start' .and. iteration > stp_value) then
                   stop
                endif
                call snapshot(s, iteration, t, 13)
                print*,'snapshot', iteration, t, 'T', dt_T
             end if
  