diff -ruN xf_thermspinxferevolve/xf_thermspinxferevolve.cc xf_thermspinxferevolve-patched/xf_thermspinxferevolve.cc
--- xf_thermspinxferevolve/xf_thermspinxferevolve.cc	2014-03-25 18:12:43.000000000 -0400
+++ xf_thermspinxferevolve-patched/xf_thermspinxferevolve.cc	2014-04-09 18:06:11.000000000 -0400
@@ -76,8 +76,8 @@
     max_step_decrease(0.03125), max_step_increase_limit(4.0),
     max_step_increase_adj_ratio(1.9),
     reject_goal(0.05), reject_ratio(0.05),
-    KBoltzmann(1.38062e-23),tmpStepsize(1e-17),
-    has_tempscript(0),temperature(0.),
+    KBoltzmann(1.38062e-23),temperature(0.),
+    tmpStepsize(1e-17),has_tempscript(0),
     energy_state_id(0),next_timestep(1e-27),
     rkstep_ptr(NULL)
 {
@@ -1573,7 +1573,7 @@
   GetEnergyDensity(next_state_key.GetReadReference(),temp_energy,
                    &vtmpA,NULL,pE_pt);
   UpdateHfluct(cstate->mesh);// Update the random field before calculating slopes
-  Oxs_SimState& workstate = next_state_key.GetWriteReference();
+  // Oxs_SimState& workstate = next_state_key.GetWriteReference();
   tmpStepsize = stepsize;
   Calculate_dm_dt(next_state_key.GetReadReference(),vtmpA,pE_pt,
                   vtmpA,max_dm_dt,dE_dt,timestep_lower_bound);
