diff --git a/src/input_xml.F90 b/src/input_xml.F90
index 1c8bcc22a3..8db0b1ee45 100644
--- a/src/input_xml.F90
+++ b/src/input_xml.F90
@@ -27,7 +27,7 @@ module input_xml
   use plot_header
   use random_lcg,       only: prn
   use surface_header
-  use set_header,       only: SetChar
+  use set_header,       only: SetChar, SetInt
   use settings
   use stl_vector,       only: VectorInt, VectorReal, VectorChar
   use string,           only: to_lower, to_str, str_to_int, str_to_real, &
diff --git a/src/settings.F90 b/src/settings.F90
index a64d60b704..19ec3f7714 100644
--- a/src/settings.F90
+++ b/src/settings.F90
@@ -3,7 +3,6 @@ module settings
   use, intrinsic :: ISO_C_BINDING
 
   use constants
-  use set_header,    only: SetInt
 
   implicit none
 
@@ -67,9 +66,6 @@ module settings
   integer(C_INT32_T), bind(C) :: index_ufs_mesh
 
   ! Write source at end of simulation
-  logical(C_BOOL), bind(C) :: source_separate
-  logical(C_BOOL), bind(C) :: source_write
-  logical(C_BOOL), bind(C) :: source_latest
 
   ! Variance reduction settins
   logical(C_BOOL), bind(C) :: survival_biasing
@@ -100,12 +96,6 @@ module settings
   ! Particle tracks
   logical(C_BOOL), bind(C) :: write_all_tracks
 
-  ! Particle restart run
-  logical(C_BOOL), bind(C) :: particle_restart_run
-
-  ! Write out initial source
-  logical(C_BOOL), bind(C) :: write_initial_source
-
   ! Whether create fission neutrons or not. Only applied for MODE_FIXEDSOURCE
   logical(C_BOOL), bind(C) :: create_fission_neutrons
 
@@ -119,7 +109,6 @@ module settings
 
   ! Various output options
   logical(C_BOOL), bind(C) :: output_summary
-  logical(C_BOOL), bind(C) :: output_tallies
 
   ! Resonance scattering settings
   logical(C_BOOL), bind(C) :: res_scat_on ! is resonance scattering treated?
diff --git a/src/simulation_header.F90 b/src/simulation_header.F90
index e16b934322..1df4b0a4e3 100644
--- a/src/simulation_header.F90
+++ b/src/simulation_header.F90
@@ -23,7 +23,6 @@ module simulation_header
   integer(C_INT), bind(C) :: current_batch     ! current batch
   integer(C_INT), bind(C) :: current_gen       ! current generation within a batch
   integer(C_INT), bind(C) :: total_gen         ! total number of generations simulated
-  logical(C_BOOL), bind(C, name='initialized') :: simulation_initialized
   logical(C_BOOL), bind(C) :: need_depletion_rx ! need to calculate depletion reaction rx?
 
   ! ============================================================================
@@ -32,7 +31,6 @@ module simulation_header
   logical(C_BOOL), bind(C) :: satisfy_triggers  ! whether triggers are satisfied
 
   integer(C_INT64_T), bind(C) :: work         ! number of particles per processor
-  integer(C_INT64_T), bind(C) :: current_work ! index in source bank of current history simulated
 
   ! ============================================================================
   ! K-EIGENVALUE SIMULATION VARIABLES
@@ -59,7 +57,7 @@ module simulation_header
 
   logical(C_BOOL), bind(C) :: trace
 
-!$omp threadprivate(trace, thread_id, current_work)
+!$omp threadprivate(trace, thread_id)
 
   interface
     subroutine entropy_clear() bind(C)