diff --git a/CMakeLists.txt b/CMakeLists.txt
index d956b8d7e1..5a204fb66a 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -323,7 +323,6 @@ add_library(libopenmc SHARED
   src/message_passing.F90
   src/mgxs_interface.F90
   src/nuclide_header.F90
-  src/output.F90
   src/particle_header.F90
   src/particle_restart.F90
   src/photon_header.F90
diff --git a/src/output.F90 b/src/output.F90
deleted file mode 100644
index a9430deb34..0000000000
--- a/src/output.F90
+++ /dev/null
@@ -1,79 +0,0 @@
-module output
-
-  use, intrinsic :: ISO_C_BINDING
-  use, intrinsic :: ISO_FORTRAN_ENV
-
-  use constants
-  use endf,            only: reaction_name
-  use error,           only: fatal_error, warning
-  use geometry_header
-  use math,            only: t_percentile
-  use message_passing, only: master, n_procs
-  use mgxs_interface
-  use nuclide_header
-  use particle_header, only: Particle
-  use settings
-  use simulation_header
-  use surface_header,  only: surfaces
-  use string,          only: to_upper, to_str
-  use tally_header
-  use tally_derivative_header
-  use tally_filter
-
-  implicit none
-
-  ! Short names for output and error units
-  integer :: ou = OUTPUT_UNIT
-  integer :: eu = ERROR_UNIT
-
-  interface
-    subroutine print_particle(p) bind(C)
-      import Particle
-      type(Particle), intent(in) :: p
-    end subroutine
-
-    subroutine write_tallies() bind(C)
-    end subroutine
-  end interface
-
-contains
-
-!===============================================================================
-! HEADER displays a header block according to a specified level. If no level is
-! specified, it is assumed to be a minor header block.
-!===============================================================================
-
-  subroutine header(msg, level, unit)
-    character(*), intent(in)      :: msg   ! header message
-    integer, intent(in)           :: level
-    integer, intent(in), optional :: unit  ! unit to write to
-
-    integer :: n            ! number of = signs on left
-    integer :: m            ! number of = signs on right
-    integer :: unit_        ! unit to write to
-    character(MAX_LINE_LEN) :: line
-
-    ! set default unit
-    if (present(unit)) then
-      unit_ = unit
-    else
-      unit_ = OUTPUT_UNIT
-    end if
-
-    ! determine how many times to repeat '=' character
-    n = (63 - len_trim(msg))/2
-    m = n
-    if (mod(len_trim(msg),2) == 0) m = m + 1
-
-    ! convert line to upper case
-    line = to_upper(msg)
-
-    ! print header based on level
-    if (verbosity >= level) then
-      write(UNIT=unit_, FMT='(/1X,A/)') repeat('=', n) // '>     ' // &
-           trim(line) // '     <' // repeat('=', m)
-    end if
-
-  end subroutine header
-
-end module output
diff --git a/src/particle_restart.F90 b/src/particle_restart.F90
index 9474f58896..dd1c89669d 100644
--- a/src/particle_restart.F90
+++ b/src/particle_restart.F90
@@ -8,7 +8,6 @@ module particle_restart
   use hdf5_interface,   only: file_open, file_close, read_dataset, HID_T
   use mgxs_interface,   only: energy_bin_avg
   use nuclide_header,   only: micro_xs, n_nuclides
-  use output,           only: print_particle
   use particle_header
   use photon_header,    only: micro_photon_xs, n_elements
   use random_lcg,       only: set_particle_seed
@@ -37,6 +36,11 @@ contains
     interface
       subroutine set_micro_xs() bind(C)
       end subroutine
+
+      subroutine print_particle(p) bind(C)
+        import Particle
+        type(Particle), intent(in) :: p
+      end subroutine
     end interface
 
     err = 0