anchor.diff:    end if
   
     ! Display message
  -  if(main_process()) call message(1,' [main] exiting final iteration')
  +  if(main_process()) write(*,'(" [main] exiting final iteration")')
   
     if(use_raytracing) then
   
        ! Display message
  -     if(main_process()) call message(1,' [main] starting raytracing iteration')
  +     if(main_process()) write(*,'(" [main] starting raytracing iteration")')
   
        ! Do the raytracing
        call do_raytracing(n_raytracing_photons_star,n_raytracing_photons_dust, n_stats)
   
        ! Display message
  -     if(main_process()) call message(1,' [main] exiting raytracing iteration')
  +     if(main_process()) write(*,'(" [main] exiting raytracing iteration")')
   
     end if
   
  
positive.code:   do iter=1,n_lucy_iter
  
       ! Display message
       if(main_process()) call message_number(1,' [main] starting Lucy iteration ',iter,'(I0)','')
  
       ! Do the RT
       call do_lucy(n_lucy_photons, n_stats)
  
negative.code:      call mp_join()
  
       ! Display message
       if(main_process()) write(*,'(" [main] exiting Lucy iteration")')
  
       ! Output files
       if(main_process()) call output_grid(handle_out, iter, n_lucy_iter)
  