diff --git a/src/input_xml.F90 b/src/input_xml.F90
index 61468c898..37d9d7ca5 100644
--- a/src/input_xml.F90
+++ b/src/input_xml.F90
@@ -101,7 +101,7 @@ module input_xml
       import C_PTR
       type(C_PTR) :: node_ptr
     end subroutine read_plots
-    
+
     subroutine set_particle_energy_bounds(particle, E_min, E_max) bind(C)
       import C_INT, C_DOUBLE
       integer(C_INT), value :: particle
@@ -2056,7 +2056,7 @@ contains
     type(XMLNode), allocatable :: node_mask_list(:)
     type(XMLNode), allocatable :: node_meshline_list(:)
 
-    
+
     ! Check if plots.xml exists
     filename = trim(path_input) // "plots.xml"
     inquire(FILE=filename, EXIST=file_exists)
@@ -2074,7 +2074,7 @@ contains
 
     call read_plots(root % ptr)
     return
-    
+
     ! Get list pointer to XML <plot>
     call get_node_list(root, "plot", node_plot_list)
 
diff --git a/src/plot.F90 b/src/plot.F90
index 98b16b45e..faccb317c 100644
--- a/src/plot.F90
+++ b/src/plot.F90
@@ -17,7 +17,7 @@ module plot
       import C_INT
       integer(C_INT) :: err
     end function openmc_plot_geometry_c
-     
+
     subroutine create_ppm(pl) bind(C)
       import ObjectPlot
       type(ObjectPlot), intent(in) :: pl
@@ -40,24 +40,6 @@ contains
 
     err = openmc_plot_geometry_c()
     return
-
-    do i = 1, n_plots
-      associate (pl => plots(i))
-        ! Display output message
-        call write_message("Processing plot " // trim(to_str(pl % id)) &
-             // ": " // trim(pl % path_plot) // " ...", 5)
-
-        if (pl % type == PLOT_TYPE_SLICE) then
-          ! create 2d image
-          call create_ppm(pl)
-        else if (pl % type == PLOT_TYPE_VOXEL) then
-          ! create dump for 3D silomesh utility script
-          call create_voxel(pl)
-        end if
-      end associate
-    end do
-
-    err = 0
   end function openmc_plot_geometry
 
 end module plot