       //cudaThreadSynchronize();
       //potent2Timer.stop();
 
+#ifdef DEBUG_TRACE
+      std::cout << "Field" << std::endl;
+#endif
       //fieldTimer.start();
       // calculate E field at Grid points
       field(dev_ex,dev_ey,dev_phi);


   if(options.getRestartPoint() > 0)
   {
      loadPrevSimState(options.getRestartPoint(), options.getRestartDir(),
         d_eleHotLoc, d_eleHotVel, d_eleColdLoc, d_eleColdVel,
         d_ionHotLoc, d_ionHotVel, d_ionColdLoc, d_ionColdVel,

      // logger.flush();
      // END DEBUG

#ifdef DEBUG_TRACE
      std::cout << "Move" << std::endl;
#endif
      //movepTimer.start();
      // move ions
      cudaStream_t movepStreams[4];
