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:     first = .true.
      time_curr = 0._dp
      n_completed = 0
    end subroutine mp_reset_first
  
    subroutine mp_n_photons(n_photons_tot, n_photons_curr, n_photons_chunk, n_photons)
  
      implicit none
  
      ! Total number of photons, and size of a chunk
      integer(idp),intent(in) :: n_photons_tot, n_photons_chunk
  
      ! Number of photons requested so far
      integer(idp),intent(inout) :: n_photons_curr
  
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
  
  