 //////-+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 	ofstream flexTipPos;
 	int numFlexBodiesInOnePeriod = int(ANCF_ReferenceArrayNodesOnBeams.size() / real_(paramsH.nPeriod) + .5);
-	int tFlexTipPos = 10;
+	int tFlexTipPos = 1000;
 	if (tStep % tFlexTipPos == 0) {
 		if (tStep / tFlexTipPos == 0) {
 			flexTipPos.open("dataFlexTip.txt");

						length(dist2) / channelRadius << ", " << length(dist2) / channelRadius << ", " <<
						p_rigidCumul.x << ", " << p_rigid.y << ", " << p_rigid.z << ", " <<
						v_rigid.x << ", " << v_rigid.y << ", " << v_rigid.z << ", " <<
						omega_rigid.x << ", " << omega_rigid.y << ", " << omega_rigid.z << ", " << aveVel.x << ", " << aveVel.y << ", " << aveVel.z << ", " << length(aveVel) << ", " << channelRadius << ", ";
			}
			ssParticleCenterVsTime<<endl;

using namespace std;


real_ AngleF3F3(real3 a3, real3 b3) {
	real_ cosTheta = dot(a3, b3);
	if(cosTheta>1){
		cosTheta=1;
	} else if(cosTheta<-1){
		cosTheta=-1;
	}
	return acos(cosTheta);
}
//*******************************************************************************************************************************

