 	paramsH_B.bodyForce4 = R4(0);
 	paramsH_B.gravity = R3(0);
 
+	SimParams currentParamsH = paramsH;
+
 	//for (int tStep = 0; tStep < 0; tStep ++) {
 	for (int tStep = 0; tStep < stepEnd + 1; tStep++) {
 		//edit  since yu deleted cyliderRotOmegaJD

//		if (tStep < 1000) delT = 0.25 * delTOrig; else delT = delTOrig;

		if (tStep <= numPause) 	{
			setParameters(&paramsH_B); 														// sets paramsD in SDKCollisionSystem
			cutilSafeCall( cudaMemcpyToSymbolAsync(paramsD, &paramsH_B, sizeof(SimParams))); 	//sets paramsD for this file
		} else {
			setParameters(&paramsH); 														// sets paramsD in SDKCollisionSystem
			cutilSafeCall( cudaMemcpyToSymbolAsync(paramsD, &paramsH, sizeof(SimParams))); 	//sets paramsD for this file
		}

		//computations
				//markers

							ANCF_IsCantilever,
							referenceArray,

							flexParams,
							numObjects,
							float(tStep)/stepEnd,
