 		return;
 	}
 	uint absMarkerIndex = index + startFlexMarkersD;
-	real_ s = flexParametricDist[index];
+	real_ s = flexParametricDistD[index];
 
 	real3 derivVel = F3( derivVelRhoD[absMarkerIndex] );
 	real3 markerForce = markerMass * derivVel;
 
 	int flexBodyIndex = flexIdentifierD[index];
-	int numFlexMarkersPreviousBeamsTotal = ANCF_NumMarkers_Per_Beam_cumul[flexBodyIndex];
-	int numSavedForcesSoFar = ANCF_NumNodesMultMarkers_Per_Beam_Cumul[flexBodyIndex];
+	int numFlexMarkersPreviousBeamsTotal = ANCF_NumMarkers_Per_Beam_cumulD[flexBodyIndex];
+	int numSavedForcesSoFar = ANCF_NumNodesMultMarkers_Per_Beam_CumulD[flexBodyIndex];
 		int markerIndexOnThisBeam = index - numFlexMarkersPreviousBeamsTotal
 
-		int numMarkersOnThisBeam = ANCF_NumMarkers_Per_Beam[flexBodyIndex];
+		int numMarkersOnThisBeam = ANCF_NumMarkers_Per_BeamD[flexBodyIndex];
 
 	//TODO: Map Marker Force to ANCF Nodes, gives you as many forces as the number of nodes per beam
 //	F0, F1, ..., F(m-1) : Forces on nodes 0, 1, 2, ..., m-1

			R3CAST(flexNodesForcesAllMarkers1),
			R3CAST(flexNodesForcesAllMarkers2),
			I1CAST(flexIdentifierD),
			I1CAST(ANCF_NumMarkers_Per_Beam),
			I1CAST(ANCF_NumMarkers_Per_Beam_cumul),
			I1CAST(ANCF_NumNodesMultMarkers_Per_Beam_Cumul),
			R1CAST(flexParametricDist),
			R4CAST(derivVelRhoD),
			markerMass);


	velMasD.clear();
	rhoPresMuD.clear();
	posRigidD.clear();

	ANCF_NodesD.clear();
	ANCF_SlopesD.clear();
	ANCF_VelNodesD.clear();
	ANCF_VelSlopesD.clear();
	ANCF_Beam_LengthD.clear();

	posRigidCumulativeD.clear();
	velMassRigidD.clear();
	omegaLRF_D.clear();
