anchor.diff:    call mp_join()
   
     if(main_process()) then
  +     write(*,*) repeat('-',60)
        write(*,'(" Total CPU time elapsed: ",F16.2)') time
  +     write(*,'(" Ended on ",A)') trim(now())
  +     write(*,*) repeat('-',60)
     end if
   
     ! Stop multi-processing
  
positive.code:   if(trim(input_file)=="".or.trim(output_file)=="") then
       stop "Usage: bin/rt input_file output_file"
    end if
  
    ! Start up multi-processing if needed
    call mp_initialize()
  
    ! SETUP
  
    if(main_process()) call message_section('Initial setup')
  
    handle_in = hdf5_open_read(input_file)
  
  
negative.code:   character(len=1000) :: input_file, output_file
    integer(hid_t) :: handle_in, handle_out, g_peeled, g_binned
    real(dp) :: time1, time2, time
    
    call hdf5_set_compression(.true.)
  
    ! Retrieve command-line arguments
  