 	thrust::device_vector<real3> flex_NodesForcesD1 = flex_FSI_NodesForcesD1;
 	thrust::device_vector<real3> flex_NodesForcesD2 = flex_FSI_NodesForcesD2;
 
-	thrust::host_vector<real3> ANCF_NodesH2 = ANCF_NodesD2;
-	thrust::host_vector<real3> ANCF_SlopesH2 = ANCF_SlopesD2;
-	thrust::host_vector<real3> ANCF_NodesVelH2 = ANCF_NodesVelD2;
-	thrust::host_vector<real3> ANCF_SlopesVelH2 = ANCF_SlopesVelD2;
-
-	thrust::host_vector<real3> ANCF_NodesH = ANCF_NodesD;
-	thrust::host_vector<real3> ANCF_SlopesH = ANCF_SlopesD;
-	thrust::host_vector<real3> ANCF_NodesVelH = ANCF_NodesVelD;
-	thrust::host_vector<real3> ANCF_SlopesVelH = ANCF_SlopesVelD;
-
-	thrust::host_vector<real3> flex_NodesForcesH1 = flex_FSI_NodesForcesD1;
-	thrust::host_vector<real3> flex_NodesForcesH2 = flex_FSI_NodesForcesD2;
-	thrust::host_vector<int2>  ANCF_ReferenceArrayNodesOnBeamsH = ANCF_ReferenceArrayNodesOnBeamsD;
-	thrust::host_vector<real_> ANCF_Beam_LengthH = ANCF_Beam_LengthD;
-	thrust::host_vector<bool> ANCF_IsCantileverH = ANCF_IsCantileverD;
 
 	//---------------------------------------------
 	int totalNumberOfFlexNodes = flex_NodesForcesD1.size();

	ANCF_Beam_LengthH.clear();
	ANCF_IsCantileverH.clear();

	ANCF_NodesAccD.clear();
	ANCF_SlopesAccD.clear();
	ANCF_NodesAccH.clear();
	ANCF_SlopesAccH.clear();

}

	thrust::copy(ANCF_NodesVelH2.begin(), ANCF_NodesVelH2.end(), ANCF_NodesVelD2.begin());
	thrust::copy(ANCF_SlopesVelH2.begin(), ANCF_SlopesVelH2.end(), ANCF_SlopesVelD2.begin());
	//---------------------------------------------
	//------------------ CPU code should work on CPU data
	ANCF_NodesH2.clear();
	ANCF_SlopesH2.clear();
	ANCF_NodesVelH2.clear();
