diff --git a/compute.py b/compute.py
index 48bf7c4..8ae3cfb 100644
--- a/compute.py
+++ b/compute.py
@@ -1445,6 +1445,7 @@ def run_simulation(save_trajectory=0):
         header_fields={"DT": str(DT), "NSTEP": str(NSTEP), "method": str(METHOD),
                        "warmup_steps": str(warmup_steps or 0),
                        "dynamic_steps": str(record_steps),
+                       "T_total": str(NT * DT),
                        "X_MAX": str(X_MAX), "X_MIN": str(X_MIN),
                        "Y_MAX": str(Y_MAX), "Y_MIN": str(Y_MIN),
                        "Z_MAX": str(Z_MAX), "Z_MIN": str(Z_MIN),
diff --git a/dynamics.py b/dynamics.py
index a0724b7..c5ee746 100644
--- a/dynamics.py
+++ b/dynamics.py
@@ -228,6 +228,7 @@ def run_case(config_path, runtime_base, input_dir="input", output_dir="output",
             hf = {"DT": str(DT), "NSTEP": str(NSTEP), "method": str(data.get("method", "")),
                   "warmup_steps": str(data.get("warmup_steps", 0)),
                   "dynamic_steps": str(record_steps),
+                  "T_total": str(NT * DT),
                   "X_MAX": str(data.get("X_MAX", 10)), "X_MIN": str(data.get("X_MIN", -10)),
                   "Y_MAX": str(data.get("Y_MAX", 10)), "Y_MIN": str(data.get("Y_MIN", -10)),
                   "Z_MAX": str(data.get("Z_MAX", 10)), "Z_MIN": str(data.get("Z_MIN", -10)),