diff --git a/src/output.F90 b/src/output.F90
index f188874d26..883165556c 100644
--- a/src/output.F90
+++ b/src/output.F90
@@ -778,21 +778,26 @@ contains
       endif
 
       ! Write derivative information.
-      !if (allocated(t % deriv)) then
-      !  select case (t % deriv % variable)
-      !  case (DIFF_DENSITY)
-      !    write(unit=unit_tally, fmt="(' Density derivative  Material ',A)") &
-      !         to_str(t % deriv % diff_material)
-      !  case (DIFF_NUCLIDE_DENSITY)
-      !    i_listing = nuclides(t % deriv % diff_nuclide) % listing
-      !    write(unit=unit_tally, fmt="(' Nuclide density derivative  Material '&
-      !         &,A,'  Nuclide ',A)") trim(to_str(t % deriv % diff_material)), &
-      !         trim(xs_listings(i_listing) % alias)
-      !  case default
-      !    call fatal_error("Differential tally dependent variable for tally " &
-      !         // trim(to_str(t % id)) // " not defined in output.F90.")
-      !  end select
-      !end if
+      if (t % deriv /= NONE) then
+        associate(deriv => tally_derivs(t % deriv))
+          select case (deriv % variable)
+          case (DIFF_DENSITY)
+            write(unit=unit_tally, fmt="(' Density derivative  Material ',A)") &
+                 to_str(deriv % diff_material)
+          case (DIFF_NUCLIDE_DENSITY)
+            write(unit=unit_tally, fmt="(' Nuclide density derivative  &
+                 &Material ',A,'  Nuclide ',A)") &
+                 trim(to_str(deriv % diff_material)), &
+                 trim(nuclides(deriv % diff_nuclide) % name)
+          case (DIFF_TEMPERATURE)
+            write(unit=unit_tally, fmt="(' Temperature derivative  Material ',&
+                 &A)") to_str(deriv % diff_material)
+          case default
+            call fatal_error("Differential tally dependent variable for tally "&
+                 // trim(to_str(t % id)) // " not defined in output.F90.")
+          end select
+        end associate
+      end if
 
       ! Handle surface current tallies separately
       if (t % type == TALLY_SURFACE_CURRENT) then