anchor.diff:  
       select case(rank)
       case(rank_main)
  +
  +       ! Find the number of photons per chunk
  +       if(n_photons_chunk == 0) then
  +          n_photons_chunk = nint(real(n_photons_tot, dp) / real(nproc, dp) / real(n_steps, dp))
  +       end if
   
          time_curr = time_curr + time2-time1
   
  
positive.code:              if(n_photons_curr > n_photons_tot) stop "n_photons_curr > n_photons_tot"
  
               ! Find how many photons to request and increment counter
               if(n_photons_tot - n_photons_curr <= n_photons_chunk*2_idp) then
                  n_photons_send(ir) = n_photons_chunk / 10_idp
               else
                  n_photons_send(ir) = n_photons_chunk
               end if
               n_photons_curr = n_photons_curr + n_photons_send(ir)
  
               ! Send number of photons and initialize receive for status check
  
negative.code:   public :: mp_set_random_seed
  
    real(dp) :: time_curr
    integer(idp) :: n_completed, n_photons_chunk
    integer(idp) :: n_steps = 10
    integer(idp) :: n_stats_last
  
  contains
  
  