 	int indexOfClosestNode = int(s / l) * nNodes;
 	if (indexOfClosestNode == nNodes) indexOfClosestNode--;
 
-	real3_ beamPointPos = ANCF_Point_Pos(ANCF_Nodes, ANCF_Slopes, indexOfClosestNode, s, l); //interpolation using ANCF beam, cubic hermit equation
-
-	real3_ dist3 = posRadD[absMarkerIndex] - beamPointPos;
+	real3 beamPointPos = ANCF_Point_Pos(ANCF_Nodes, ANCF_Slopes, indexOfClosestNode, s, l); //interpolation using ANCF beam, cubic hermit equation
+
+	real3 dist3 = posRadD[absMarkerIndex] - beamPointPos;
 	flexSPH_MeshPos_LRF_D[index] = dist3;
 }
 

		thrust::device_vector<int> & ANCF_NumNodesMultMarkers_Per_Beam_Cumul,

		const thrust::device_vector<int> & flexIdentifierD,
		const thrust::device_vector<real_> & parametricDist,
		const thrust::device_vector<real_> & ANCF_Beam_Length,
		const thrust::host_vector<int2> & referenceArray,

	int indexOfClosestNode = int(s / l) * nNodes;
	if (indexOfClosestNode == nNodes) indexOfClosestNode--;

	real3 beamPointPos = ANCF_Point_Pos(ANCF_Nodes, ANCF_Slopes, indexOfClosestNode, s, l); //interpolation using ANCF beam, cubic hermit equation
	real3 beamPointSlope = ANCF_Point_Slope(ANCF_Nodes, ANCF_Slopes, indexOfClosestNode, s, l); //interpolation using ANCF beam, cubic hermit equation
	real3 beamPointOmega;

	real3 beamPointVel = ANCF_Point_Vel(ANCF_Nodes, ANCF_Slopes, ANCF_VelNodes, ANCF_VelSlopes, indexOfClosestNode, s, l); //interpolation using ANCF beam, cubic hermit equation

	real3 dist3 = flexSPH_MeshPos_LRF_D[index];
}
//--------------------------------------------------------------------------------------------------------------------------------
void MapSPH_ToGrid(
