diff --git a/src/thermo.cpp b/src/thermo.cpp
index 1e83376864..8de74a47fd 100644
--- a/src/thermo.cpp
+++ b/src/thermo.cpp
@@ -41,7 +41,7 @@ using namespace LAMMPS_NS;
 
 // customize a new keyword by adding to this list:
 
-// step, elapsed, elaplong, dt, cpu
+// step, elapsed, elaplong, dt, cpu, tpcpu, spcpu
 // atoms, temp, press, pe, ke, etotal, enthalpy
 // evdwl, ecoul, epair, ebond, eangle, edihed, eimp, emol, elong, etail
 // vol, lx, ly, lz, xlo, xhi, ylo, yhi, zlo, zhi, xy, xz, yz, xlat, ylat, zlat
@@ -592,6 +592,10 @@ void Thermo::parse_fields(char *str)
       addfield("Dt",&Thermo::compute_dt,FLOAT);
     } else if (strcmp(word,"cpu") == 0) {
       addfield("CPU",&Thermo::compute_cpu,FLOAT);
+    } else if (strcmp(word,"tpcpu") == 0) {
+      addfield("T/CPU",&Thermo::compute_tpcpu,FLOAT);
+    } else if (strcmp(word,"spcpu") == 0) {
+      addfield("S/CPU",&Thermo::compute_spcpu,FLOAT);
 
     } else if (strcmp(word,"atoms") == 0) {
       addfield("Atoms",&Thermo::compute_atoms,INT);
@@ -905,9 +909,20 @@ int Thermo::evaluate_keyword(char *word, double *answer)
     compute_dt();
 
   } else if (strcmp(word,"cpu") == 0) {
-    if (update->whichflag == 0) firststep = 0;
+    if (update->whichflag == 0) 
+      error->all("This variable thermo keyword cannot be used between runs");
     compute_cpu();
 
+  } else if (strcmp(word,"tpcpu") == 0) {
+    if (update->whichflag == 0) 
+      error->all("This variable thermo keyword cannot be used between runs");
+    compute_tpcpu();
+
+  } else if (strcmp(word,"spcpu") == 0) {
+    if (update->whichflag == 0) 
+      error->all("This variable thermo keyword cannot be used between runs");
+    compute_spcpu();
+
   } else if (strcmp(word,"atoms") == 0) {
     compute_atoms();
     dvalue = ivalue;
@@ -1328,6 +1343,50 @@ void Thermo::compute_cpu()
 
 /* ---------------------------------------------------------------------- */
 
+void Thermo::compute_tpcpu()
+{
+  double new_cpu;
+  double new_time = update->ntimestep * update->dt;
+
+  if (firststep == 0) {
+    new_cpu = 0.0;
+    dvalue = 0.0;
+  } else {
+    new_cpu = timer->elapsed(TIME_LOOP);
+    double cpu_diff = new_cpu - last_tpcpu;
+    double time_diff = new_time - last_time;
+    if (time_diff > 0.0 && cpu_diff > 0.0) dvalue = time_diff/cpu_diff;
+    else dvalue = 0.0;
+  }
+
+  last_time = new_time;
+  last_tpcpu = new_cpu;
+}
+
+/* ---------------------------------------------------------------------- */
+
+void Thermo::compute_spcpu()
+{
+  double new_cpu;
+  int new_step = update->ntimestep;
+
+  if (firststep == 0) {
+    new_cpu = 0.0;
+    dvalue = 0.0;
+  } else {
+    new_cpu = timer->elapsed(TIME_LOOP);
+    double cpu_diff = new_cpu - last_spcpu;
+    int step_diff = new_step - last_step;
+    if (cpu_diff > 0.0) dvalue = step_diff/cpu_diff;
+    else dvalue = 0.0;
+  }
+
+  last_step = new_step;
+  last_spcpu = new_cpu;
+}
+
+/* ---------------------------------------------------------------------- */
+
 void Thermo::compute_atoms()
 {
   ivalue = static_cast<int> (natoms);
diff --git a/src/thermo.h b/src/thermo.h
index 0f9d06732e..f261d3c13c 100644
--- a/src/thermo.h
+++ b/src/thermo.h
@@ -57,6 +57,10 @@ class Thermo : protected Pointers {
   int lostflag,lostbefore;
   int flushflag,lineflag;
 
+  double last_tpcpu,last_spcpu;
+  double last_time;
+  int last_step;
+
                          // data used by routines that compute single values
   int ivalue;            // integer value to print
   double dvalue,natoms;  // dvalue = double value to print
@@ -111,6 +115,8 @@ class Thermo : protected Pointers {
   void compute_elapsed_long();
   void compute_dt();
   void compute_cpu();
+  void compute_tpcpu();
+  void compute_spcpu();
 
   void compute_atoms();
   void compute_temp();