       //logger.flush();
       // END DEBUG
 
+#ifdef DEBUG_TRACE
+      std::cout << "Dens" << std::endl;
+#endif
       //densTimer.start();
       // determine the charge density at the grid points
       dens(dev_rho, dev_rhoe,dev_rhoi, 

   // 6 rands for cold ions
   const int neededRands = neededParticles * 4 * 6;

   // CUDA Variables
   int sharedMemoryBytes;
   // Device Memory Pointers

      //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);
