diff --git a/src/solver/impls/pvode/pvode.cxx b/src/solver/impls/pvode/pvode.cxx
index 045d83ed49..74b505e641 100644
--- a/src/solver/impls/pvode/pvode.cxx
+++ b/src/solver/impls/pvode/pvode.cxx
@@ -395,6 +395,7 @@ BoutReal PvodeSolver::run(BoutReal tout) {
 
       for (auto& f : f3d) {
         saveParallel(debug, f.name, *f.var);
+        debug[fmt::format("ddt_{}", f.name)] = *f.F_var;
       }
       bout::OptionsIO::write("BOUT.debug", debug, mesh);
       MPI_Barrier(BoutComm::get());