 
 
 //*************** parameters of the benchmark
-const unsigned int g_N = 64000;
-const float g_Lx = 55.1294f;
-const float g_Ly = 55.1294f;
-const float g_Lz = 55.1294f;
-const float g_rcut = 3.8f;
+unsigned int g_N;
+float g_Lx;
+float g_Ly;
+float g_Lz;
+float g_rcut;
 const unsigned int g_Nmax = 128;	// Maximum number of particles each cell can hold
 const float tweak_dist = 0.1f;
 

		cudaSetDevice(atoi(getenv("CAC_GPU_ID")));
	#endif
	
	// setup
	allocate_data();
	initialize_data();

	}


int main(int argc, char **argv)
	{
	#ifdef ENABLE_CAC_GPU_ID
	if (!getenv("CAC_GPU_ID"))
