- ofstream datafile(fname.c_str());
-
- if (datafile) {
- datafile << "U: " << FDMExec->GetTranslation()->GetU() << endl;
- datafile << "V: " << FDMExec->GetTranslation()->GetV() << endl;
- datafile << "W: " << FDMExec->GetTranslation()->GetW() << endl;
- datafile << "P: " << FDMExec->GetRotation()->GetP() << endl;
- datafile << "Q: " << FDMExec->GetRotation()->GetQ() << endl;
- datafile << "R: " << FDMExec->GetRotation()->GetR() << endl;
- datafile << "L: " << FDMExec->GetAircraft()->GetL() << endl;
- datafile << "M: " << FDMExec->GetAircraft()->GetM() << endl;
- datafile << "N: " << FDMExec->GetAircraft()->GetN() << endl;
- datafile << "latitude: " << latitude << endl;
- datafile << "longitude: " << longitude << endl;
- datafile << "alpha: " << FDMExec->GetTranslation()->Getalpha() << endl;
- datafile << "beta: " << FDMExec->GetTranslation()->Getbeta() << endl;
- datafile << "gamma: " << FDMExec->GetTranslation()->Getgamma() << endl;
- datafile << "phi: " << FDMExec->GetRotation()->Getphi() << endl;
- datafile << "tht: " << FDMExec->GetRotation()->Gettht() << endl;
- datafile << "psi: " << FDMExec->GetRotation()->Getpsi() << endl;
- datafile << "Pdot: " << FDMExec->GetRotation()->GetPdot() << endl;
- datafile << "Qdot: " << FDMExec->GetRotation()->GetQdot() << endl;
- datafile << "Rdot: " << FDMExec->GetRotation()->GetRdot() << endl;
- datafile << "h: " << h << endl;
- datafile << "a: " << a << endl;
- datafile << "rho: " << FDMExec->GetAtmosphere()->Getrho() << endl;
- datafile << "qbar: " << qbar << endl;
- datafile << "sim_time: " << sim_time << endl;
- datafile << "dt: " << dt << endl;
- datafile << "m: " << FDMExec->GetAircraft()->GetMass() << endl;
- datafile.close();
- return true;
- } else {
- return false;
- }
+ vQdot(1) = -0.5*(vQtrn(2)*vPQR(eP) + vQtrn(3)*vPQR(eQ) + vQtrn(4)*vPQR(eR));
+ vQdot(2) = 0.5*(vQtrn(1)*vPQR(eP) + vQtrn(3)*vPQR(eR) - vQtrn(4)*vPQR(eQ));
+ vQdot(3) = 0.5*(vQtrn(1)*vPQR(eQ) + vQtrn(4)*vPQR(eP) - vQtrn(2)*vPQR(eR));
+ vQdot(4) = 0.5*(vQtrn(1)*vPQR(eR) + vQtrn(2)*vPQR(eQ) - vQtrn(3)*vPQR(eP));
+
+ vQtrn += Integrate(TRAPZ, dt*rate, vQdot, vQdot_prev);
+
+ vQtrn.Normalize();
+}
+
+//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+FGColumnVector3& FGState::CalcEuler(void)
+{
+ if (mTl2b(3,3) == 0.0) mTl2b(3,3) = 0.0000001;
+ if (mTl2b(1,1) == 0.0) mTl2b(1,1) = 0.0000001;
+
+ vEuler(ePhi) = atan2(mTl2b(2,3), mTl2b(3,3));
+ vEuler(eTht) = asin(-mTl2b(1,3));
+ vEuler(ePsi) = atan2(mTl2b(1,2), mTl2b(1,1));
+
+ if (vEuler(ePsi) < 0.0) vEuler(ePsi) += 2*M_PI;
+
+ return vEuler;
+}
+
+//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+FGMatrix33& FGState::GetTs2b(void)
+{
+ double ca, cb, sa, sb;
+
+ double alpha = Translation->Getalpha();
+ double beta = Translation->Getbeta();
+
+ ca = cos(alpha);
+ sa = sin(alpha);
+ cb = cos(beta);
+ sb = sin(beta);
+
+ mTs2b(1,1) = ca*cb;
+ mTs2b(1,2) = -ca*sb;
+ mTs2b(1,3) = -sa;
+ mTs2b(2,1) = sb;
+ mTs2b(2,2) = cb;
+ mTs2b(2,3) = 0.0;
+ mTs2b(3,1) = sa*cb;
+ mTs2b(3,2) = -sa*sb;
+ mTs2b(3,3) = ca;
+
+ return mTs2b;
+}
+
+//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+FGMatrix33& FGState::GetTb2s(void)
+{
+ float alpha,beta;
+ float ca, cb, sa, sb;
+
+ alpha = Translation->Getalpha();
+ beta = Translation->Getbeta();
+
+ ca = cos(alpha);
+ sa = sin(alpha);
+ cb = cos(beta);
+ sb = sin(beta);
+
+ mTb2s(1,1) = ca*cb;
+ mTb2s(1,2) = sb;
+ mTb2s(1,3) = sa*cb;
+ mTb2s(2,1) = -ca*sb;
+ mTb2s(2,2) = cb;
+ mTb2s(2,3) = -sa*sb;
+ mTb2s(3,1) = -sa;
+ mTb2s(3,2) = 0.0;
+ mTb2s(3,3) = ca;
+
+ return mTb2s;
+}
+
+//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+void FGState::ReportState(void)
+{
+#if !defined(__BORLANDCPP__)
+ char out[80], flap[10], gear[12];
+
+ cout << endl << " JSBSim State" << endl;
+ snprintf(out,80," Weight: %7.0f lbs. CG: %5.1f, %5.1f, %5.1f inches\n",
+ FDMExec->GetMassBalance()->GetWeight(),
+ FDMExec->GetMassBalance()->GetXYZcg(1),
+ FDMExec->GetMassBalance()->GetXYZcg(2),
+ FDMExec->GetMassBalance()->GetXYZcg(3));
+ cout << out;
+ if ( FCS->GetDfPos() <= 0.01)
+ snprintf(flap,10,"Up");
+ else
+ snprintf(flap,10,"%2.0f",FCS->GetDfPos());
+
+ if (FCS->GetGearPos() < 0.01)
+ snprintf(gear,12,"Up");
+ else if (FCS->GetGearPos() > 0.99)
+ snprintf(gear,12,"Down");
+ else
+ snprintf(gear,12,"In Transit");
+
+ snprintf(out,80, " Flaps: %3s Gear: %12s\n",flap,gear);
+ cout << out;
+ snprintf(out,80, " Speed: %4.0f KCAS Mach: %5.2f\n",
+ FDMExec->GetAuxiliary()->GetVcalibratedKTS(),
+ Translation->GetMach() );
+ cout << out;
+ snprintf(out,80, " Altitude: %7.0f ft. AGL Altitude: %7.0f ft.\n",
+ Position->Geth(),
+ Position->GetDistanceAGL() );
+ cout << out;
+ snprintf(out,80, " Angle of Attack: %6.2f deg Pitch Angle: %6.2f deg\n",
+ Translation->Getalpha()*radtodeg,
+ Rotation->Gettht()*radtodeg );
+ cout << out;
+ snprintf(out,80, " Flight Path Angle: %6.2f deg Climb Rate: %5.0f ft/min\n",
+ Position->GetGamma()*radtodeg,
+ Position->Gethdot()*60 );
+ cout << out;
+ snprintf(out,80, " Normal Load Factor: %4.2f g's Pitch Rate: %5.2f deg/s\n",
+ Aircraft->GetNlf(),
+ Rotation->GetPQR(2)*radtodeg );
+ cout << out;
+ snprintf(out,80, " Heading: %3.0f deg true Sideslip: %5.2f deg Yaw Rate: %5.2f deg/s\n",
+ Rotation->Getpsi()*radtodeg,
+ Translation->Getbeta()*radtodeg,
+ Rotation->GetPQR(3)*radtodeg );
+ cout << out;
+ snprintf(out,80, " Bank Angle: %5.2f deg Roll Rate: %5.2f deg/s\n",
+ Rotation->Getphi()*radtodeg,
+ Rotation->GetPQR(1)*radtodeg );
+ cout << out;
+ snprintf(out,80, " Elevator: %5.2f deg Left Aileron: %5.2f deg Rudder: %5.2f deg\n",
+ FCS->GetDePos(ofRad)*radtodeg,
+ FCS->GetDaLPos(ofRad)*radtodeg,
+ FCS->GetDrPos(ofRad)*radtodeg );
+ cout << out;
+ snprintf(out,80, " Throttle: %5.2f%c\n",
+ FCS->GetThrottlePos(0)*100,'%' );
+ cout << out;
+
+ snprintf(out,80, " Wind Components: %5.2f kts head wind, %5.2f kts cross wind\n",
+ FDMExec->GetAuxiliary()->GetHeadWind()*fpstokts,
+ FDMExec->GetAuxiliary()->GetCrossWind()*fpstokts );
+ cout << out;
+
+ snprintf(out,80, " Ground Speed: %4.0f knots , Ground Track: %3.0f deg true\n",
+ Position->GetVground()*fpstokts,
+ Position->GetGroundTrack()*radtodeg );
+ cout << out;
+#endif
+}
+
+//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+void FGState::bind(void)
+{
+ PropertyManager->Tie("sim-time-sec",this,
+ &FGState::Getsim_time);