 __constant__ int2 portionD;
 __constant__ int flagD;
 __constant__ int numRigidBodiesD;
+__constant__ int numFlexBodiesD;
 __constant__ int startRigidMarkersD;
 __constant__ int startFlexMarkersD;
 __constant__ int numRigid_SphMarkersD;

								float(tStep)/stepEnd,
								0.5 * delT);

		ApplyBoundary(posRadD2, rhoPresMuD2, numAllMarkers, posRigidD2, velMassRigidD2, numRigidBodies);
		//*****
		ForceSPH(posRadD2, velMasD2, vel_XSPH_D, rhoPresMuD2, bodyIndexD, derivVelRhoD, referenceArray, numAllMarkers, paramsH, delT);
		UpdateFluid(posRadD, velMasD, vel_XSPH_D, rhoPresMuD, derivVelRhoD, referenceArray, delT);

//		///ff1
//	 	 real_ hh = paramsD.HSML;
//	 	 printf("dist %f \n", length(dist3)/hh);
	real3 beamPointPos = Calc_ANCF_Point_Pos(ANCF_NodesD, ANCF_SlopesD, indexOfClosestNode, sE, lE); //interpolation using ANCF beam, cubic hermit equation
	posRadD[absMarkerIndex] = beamPointPos + R3(dot(A1, dist3), dot(A2, dist3), dot(A3, dist3));

//	//ask Radu
