anchor.diff:  
     implicit none
   
  +  character(len=10)   :: arg
  +  integer            :: stp_value = 100000 ! select the stop iteration 
  +  !  don't kill the job without stp_value
  +  real(x_precision), dimension(n_cell)   :: jnk1, jnk3
     !--------------------------- Parameters for s_curve-----------------------
     real(x_precision), parameter         :: eps_in = 1.e-6_x_precision ! Precision required for the dichotomy
     real(x_precision), parameter         :: Tmin   = 10._x_precision**5.3_x_precision
  
positive.code: 
             ! Do a snapshot
             if (mod(iteration, output_freq) == 0) then
                call snapshot(s, iteration, t, 13)
                print*,'snapshot', iteration, t, 'T', dt_T
             end if
  
negative.code: 
    real(x_precision), dimension(n_cell) :: dist_crit
  
    call getarg(1,arg)
  
    !----------------------------------------------
    ! Convergence criteria for S and T
    !----------------------------------------------
  