anchor.diff:    if(main_process()) then
        write(*,*) repeat('-',60)
        write(*,'(" Total CPU time elapsed: ",F16.2)') time
  +     call hdf5_write_keyword(handle_out, '/', 'cpu_time', time)
        datetime = now()
        write(*,'(" Ended on ",A)') trim(datetime)
  +     call hdf5_write_keyword(handle_out, '/', 'date_ended', trim(datetime))
        write(*,*) repeat('-',60)
     end if
   
  
positive.code:   ! Prepare output directory
    if(main_process()) then
       handle_out = hdf5_open_new(output_file)
    end if
  
    ! Wait for all threads
  
negative.code:   ! Start up multi-processing if needed
    call mp_initialize()
  
    ! SETUP
  
    call check_file_exists(input_file)
  