diff --git a/src/output.F90 b/src/output.F90
index fb2342efb3..5c758b09a2 100644
--- a/src/output.F90
+++ b/src/output.F90
@@ -1318,20 +1318,28 @@ contains
     do i = 1, n_plots
       pl => plots(i)
 
-      ! Write plot id
+      ! Plot id
       write(ou,100) "Plot ID:", trim(to_str(pl % id))
+      
+      ! Plot type
+      if (pl % type == PLOT_TYPE_SLICE) then
+        write(ou,100) "Plot Type:", "Slice"
+      else if (pl % type == PLOT_TYPE_3DRASTER) then
+        write(ou,100) "Plot Type:", "3D Raster"
+      end if
 
-      ! Write plotting origin
+      ! Plot parameters
       write(ou,100) "Origin:", trim(to_str(pl % origin(1))) // &
            " " // trim(to_str(pl % origin(2))) // " " // &
            trim(to_str(pl % origin(3)))
-
-      ! Write plotting width
+      write(ou,100) "Width:", trim(to_str(pl % width(1))) // &
+           " " // trim(to_str(pl % width(2)))
+      if (pl % color_by == PLOT_COLOR_CELLS) then
+        write(ou,100) "Coloring:", "Cells"
+      else if (pl % color_by == PLOT_COLOR_MATS) then
+        write(ou,100) "Coloring:", "Materials"
+      end if
       if (pl % type == PLOT_TYPE_SLICE) then
-
-        write(ou,100) "Width:", trim(to_str(pl % width(1))) // &
-             " " // trim(to_str(pl % width(2)))
-        write(ou,100) "Coloring:", trim(to_str(pl % color_by))
         select case (pl % basis)
         case (PLOT_BASIS_XY)
           write(ou,100) "Basis:", "xy"
@@ -1342,6 +1350,9 @@ contains
         end select
         write(ou,100) "Pixels:", trim(to_str(pl % pixels(1))) // " " // &
              trim(to_str(pl % pixels(2)))
+      else if (pl % type == PLOT_TYPE_3DRASTER) then
+        write(ou,100) "Voxels:", trim(to_str(pl % pixels(1))) // " " // &
+             trim(to_str(pl % pixels(2))) // " " // trim(to_str(pl % pixels(3))) 
       end if
 
       write(ou,*)