 						flexIdentifierD,
 						flexMapEachMarkerOnAllBeamNodesD,
 						flexSPH_MeshPos_LRF_D,
+						flexSPH_MeshSlope_Initial_D,
 						flexParametricDistD,
 						ANCF_Beam_LengthD,
 						referenceArray,


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

		const thrust::device_vector<int> & flexIdentifierD,
		const thrust::device_vector<int2> & flexMapEachMarkerOnAllBeamNodesD,
		const thrust::device_vector<real3> & flexSPH_MeshPos_LRF_D,
		const thrust::device_vector<real3> & flexSPH_MeshSlope_Initial_D,
		const thrust::device_vector<real_> & flexParametricDistD,
		const thrust::device_vector<real_> & ANCF_Beam_LengthD,
		const thrust::host_vector<int3> & referenceArray,
