diff --git a/src/BOSampleStepper.C b/src/BOSampleStepper.C
index c70a943..a60ecd5 100644
--- a/src/BOSampleStepper.C
+++ b/src/BOSampleStepper.C
@@ -3,7 +3,7 @@
 // BOSampleStepper.C
 //
 ////////////////////////////////////////////////////////////////////////////////
-// $Id: BOSampleStepper.C,v 1.40 2008-04-05 19:11:25 fgygi Exp $
+// $Id: BOSampleStepper.C,v 1.41 2008-04-15 01:36:15 fgygi Exp $
 
 #include "BOSampleStepper.h"
 #include "EnergyFunctional.h"
@@ -217,7 +217,7 @@ void BOSampleStepper::step(int niter)
     }
     // at this point, positions r0, velocities v0 and forces fion are
     // consistent
-    double ekin_ion, temp_ion;
+    double ekin_ion = 0.0, temp_ion = 0.0;
     if ( ionic_stepper )
     {
       ekin_ion = ionic_stepper->ekin();
diff --git a/src/CPSampleStepper.C b/src/CPSampleStepper.C
index e01a711..bfdc2bc 100644
--- a/src/CPSampleStepper.C
+++ b/src/CPSampleStepper.C
@@ -3,7 +3,7 @@
 // CPSampleStepper.C
 //
 ////////////////////////////////////////////////////////////////////////////////
-// $Id: CPSampleStepper.C,v 1.18 2008-04-05 19:13:03 fgygi Exp $
+// $Id: CPSampleStepper.C,v 1.19 2008-04-15 01:36:15 fgygi Exp $
 
 #include "CPSampleStepper.h"
 #include "SlaterDet.h"
@@ -87,7 +87,7 @@ void CPSampleStepper::step(int niter)
   Wavefunction& wf = s_.wf;
 
   const double dt = s_.ctrl.dt;
-  double ekin_ion=0.0,ekin_e, temp_ion, eta;
+  double ekin_ion=0.0,ekin_e, temp_ion=0.0, eta;
 
   const string wf_dyn = s_.ctrl.wf_dyn;
   assert(wf_dyn=="MD");