anchor.diff:       call warn("main", "no dust present, so skipping temperature iterations")
        n_lucy_iter=0
     end if
  -
  +  
     ! Loop over Lucy iterations
     do iter=1,n_lucy_iter
   
  
positive.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
  
negative.code:   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
  