 
 		//************************************************
 		//edit PrintToFile since yu deleted cyliderRotOmegaJD
-		PrintToFile(posRadD, velMasD, rhoPresMuD, referenceArray, rigidIdentifierD, posRigidD, posRigidCumulativeD, velMassRigidD, qD1, omegaLRF_D, cMax, cMin, paramsH,
+		PrintToFile(posRadD, velMasD, rhoPresMuD, referenceArray, rigidIdentifierD, posRigidD, posRigidCumulativeD, velMassRigidD, qD1, AD1, AD2, AD3, omegaLRF_D, cMax, cMin, paramsH,
 				delT, tStep, channelRadius);
 //		PrintToFileDistribution(distributionD, channelRadius, numberOfSections, tStep);
 		//************

		thrust::device_vector<float3> & posRigidCumulativeD,
		thrust::device_vector<float4> & velMassRigidD,
		thrust::device_vector<float4> & qD1,
		thrust::device_vector<float3> & omegaLRF_D,
		float3 cMax,
		float3 cMin,

				fileRigidParticlesDataForTecplot.open("dataRigidParticlesDataForTecplot.txt");
				fileRigidParticlesDataForTecplot<<"PipeRadius, PipeLength\n";
				fileRigidParticlesDataForTecplot<<	channelRadius <<", "<< cMax.x - cMin.x<<", "<< posRigidH.size()<<endl;
				fileRigidParticlesDataForTecplot<<"variables = \"t(s)\", \"x\", \"y\", \"z\", \"r\", \"CumulX\", \"vX\", \"vY\", \"vZ\", \"velMagnitude\", \"angleAxisXWithPipeAxis\", \"angleAxisYWithPipeAxis\", \"angleAxisZWithPipeAxis\"\n";
			} else {
				fileRigidParticlesDataForTecplot.open("dataRigidParticlesDataForTecplot.txt", ios::app);
			}
