X-Git-Url: https://git.mxchange.org/?a=blobdiff_plain;f=src%2FFDM%2FLaRCsim.cxx;h=47c6bfa93079a73b4cb050749a872e9635185469;hb=7a235abecdb3b17580ed81bd69d20a55b1c8b4e6;hp=6daf0aa6d32bb6d6c40e4b2fb1fbb7ce02ab4383;hpb=c90db01dc8d5462a3da22771ffa7c96f5ea31217;p=flightgear.git diff --git a/src/FDM/LaRCsim.cxx b/src/FDM/LaRCsim.cxx index 6daf0aa6d..47c6bfa93 100644 --- a/src/FDM/LaRCsim.cxx +++ b/src/FDM/LaRCsim.cxx @@ -21,46 +21,173 @@ // $Id$ -#include "LaRCsim.hxx" +#include +#include + +#include #include #include -#include #include #include #include #include +#include + +#include "IO360.hxx" +#include "LaRCsim.hxx" +FGLaRCsim::FGLaRCsim(void) { + ls_toplevel_init( 0.0, + (char *)globals->get_options()->get_aircraft().c_str() ); + lsic=new LaRCsimIC; //this needs to be brought up after LaRCsim is + copy_to_LaRCsim(); // initialize all of LaRCsim's vars + //this should go away someday -- formerly done in fg_init.cxx + Mass = 8.547270E+01; + I_xx = 1.048000E+03; + I_yy = 3.000000E+03; + I_zz = 3.530000E+03; + I_xz = 0.000000E+00; +} + +FGLaRCsim::~FGLaRCsim(void) { + if(lsic != NULL) { + delete lsic; + } +} // Initialize the LaRCsim flight model, dt is the time increment for // each subsequent iteration through the EOM -int fgLaRCsimInit(double dt) { - ls_toplevel_init(dt); +bool FGLaRCsim::init( double dt ) { + + ls_set_model_dt(dt); + // Initialize our little engine that hopefully might + eng.init(dt); + // dcl - in passing dt to init rather than update I am assuming + // that the LaRCsim dt is fixed at one value (yes it is 120hz CLO) + + // update the engines interface + FGEngInterface e; + add_engine( e ); + + + // FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::init()" ); + + double save_alt = 0.0; + + if ( get_Altitude() < -9000.0 ) { + save_alt = get_Altitude(); + set_Altitude( 0.0 ); + } - return(1); + // translate FG to LaRCsim structure + copy_to_LaRCsim(); + + // actual LaRCsim top level init + // ls_toplevel_init( dt, (char *)globals->get_options()->get_aircraft().c_str() ); + + FG_LOG( FG_FLIGHT, FG_INFO, "FG pos = " << + get_Latitude() ); + + // translate LaRCsim back to FG structure + copy_from_LaRCsim(); + + // but lets restore our original bogus altitude when we are done + if ( save_alt < -9000.0 ) { + set_Altitude( save_alt ); + } + + // set valid time for this record + stamp_time(); + + return true; } // Run an iteration of the EOM (equations of motion) -int fgLaRCsimUpdate(FGInterface& f, int multiloop) { +bool FGLaRCsim::update( int multiloop ) { + + if ( globals->get_options()->get_aircraft() == "c172" ) { + // set control inputs + eng.set_IAS( V_calibrated_kts ); + eng.set_Throttle_Lever_Pos( controls.get_throttle( 0 ) * 100.0 ); + eng.set_Propeller_Lever_Pos( 100 ); + if ( controls.get_mixture( 0 ) > 0.60 ) { + eng.set_Mixture_Lever_Pos( controls.get_mixture( 0 ) * 100.0 ); + } else { + eng.set_Mixture_Lever_Pos( 60.0 ); + } + + // update engine model + eng.update(); + + // copy engine state values onto "bus" + FGEngInterface *e = get_engine( 0 ); + e->set_Throttle( controls.get_throttle( 0 ) * 100.0 ); + e->set_Mixture( 80 ); + e->set_Prop_Advance( 100 ); + e->set_RPM( eng.get_RPM() ); + e->set_Manifold_Pressure( eng.get_Manifold_Pressure() ); + e->set_MaxHP( eng.get_MaxHP() ); + e->set_Percentage_Power( eng.get_Percentage_Power() ); + e->set_EGT( eng.get_EGT() ); + e->set_CHT( eng.get_CHT() ); + e->set_prop_thrust( eng.get_prop_thrust_SI() ); + +#if 0 + FG_LOG( FG_FLIGHT, FG_INFO, "Throttle = " + << controls.get_throttle( 0 ) * 100.0); + FG_LOG( FG_FLIGHT, FG_INFO, " Mixture = " << 80); + FG_LOG( FG_FLIGHT, FG_INFO, " RPM = " << eng.get_RPM()); + FG_LOG( FG_FLIGHT, FG_INFO, " MP = " << eng.get_Manifold_Pressure()); + FG_LOG( FG_FLIGHT, FG_INFO, " HP = " + << ( eng.get_MaxHP() * eng.get_Percentage_Power()/ 100.0) ); + FG_LOG( FG_FLIGHT, FG_INFO, " EGT = " << eng.get_EGT()); + FG_LOG( FG_FLIGHT, FG_INFO, " Thrust (N) " + << eng.get_prop_thrust_SI()); // Thrust in Newtons + FG_LOG( FG_FLIGHT, FG_INFO, '\n'); +#endif + // Hmm .. Newtons to lbs is 0.2248 ... + F_X_engine = eng.get_prop_thrust_SI() * 0.07; + } + double save_alt = 0.0; // lets try to avoid really screwing up the LaRCsim model - if ( f.get_Altitude() < -9000.0 ) { - save_alt = f.get_Altitude(); - f.set_Altitude( 0.0 ); + if ( get_Altitude() < -9000.0 ) { + save_alt = get_Altitude(); + set_Altitude( 0.0 ); } // copy control positions into the LaRCsim structure - Lat_control = controls.get_aileron(); + Lat_control = controls.get_aileron() / + globals->get_options()->get_speed_up(); Long_control = controls.get_elevator(); Long_trim = controls.get_elevator_trim(); - Rudder_pedal = controls.get_rudder(); - Throttle_pct = controls.get_throttle( 0 ); - Brake_pct = controls.get_brake( 0 ); + Rudder_pedal = controls.get_rudder() / + globals->get_options()->get_speed_up(); + Flap_handle = 30.0 * controls.get_flaps(); + + if ( globals->get_options()->get_aircraft() == "c172" ) { + Use_External_Engine = 1; + } else { + Use_External_Engine = 0; + } + + Throttle_pct = controls.get_throttle( 0 ) * 1.0; + + Brake_pct[0] = controls.get_brake( 1 ); + Brake_pct[1] = controls.get_brake( 0 ); // Inform LaRCsim of the local terrain altitude - Runway_altitude = f.get_Runway_altitude(); + // Runway_altitude = get_Runway_altitude(); + Runway_altitude = scenery.cur_elev * METER_TO_FEET; + + // Weather + /* V_north_airmass = get_V_north_airmass(); + V_east_airmass = get_V_east_airmass(); + V_down_airmass = get_V_down_airmass(); */ + // old -- FGInterface_2_LaRCsim() not needed except for Init() // translate FG to LaRCsim structure @@ -73,333 +200,532 @@ int fgLaRCsimUpdate(FGInterface& f, int multiloop) { // printf("%d FG_Altitude = %.2f\n", i, FG_Altitude * 0.3048); // printf("%d Altitude = %.2f\n", i, Altitude * 0.3048); - + // translate LaRCsim back to FG structure so that the // autopilot (and the rest of the sim can use the updated // values - fgLaRCsim_2_FGInterface(f); + copy_from_LaRCsim(); // but lets restore our original bogus altitude when we are done if ( save_alt < -9000.0 ) { - f.set_Altitude( save_alt ); + set_Altitude( save_alt ); } - return 1; + return true; } // Convert from the FGInterface struct to the LaRCsim generic_ struct -int FGInterface_2_LaRCsim (FGInterface& f) { - - Mass = f.get_Mass(); - I_xx = f.get_I_xx(); - I_yy = f.get_I_yy(); - I_zz = f.get_I_zz(); - I_xz = f.get_I_xz(); - // Dx_pilot = f.get_Dx_pilot(); - // Dy_pilot = f.get_Dy_pilot(); - // Dz_pilot = f.get_Dz_pilot(); - Dx_cg = f.get_Dx_cg(); - Dy_cg = f.get_Dy_cg(); - Dz_cg = f.get_Dz_cg(); - // F_X = f.get_F_X(); - // F_Y = f.get_F_Y(); - // F_Z = f.get_F_Z(); - // F_north = f.get_F_north(); - // F_east = f.get_F_east(); - // F_down = f.get_F_down(); - // F_X_aero = f.get_F_X_aero(); - // F_Y_aero = f.get_F_Y_aero(); - // F_Z_aero = f.get_F_Z_aero(); - // F_X_engine = f.get_F_X_engine(); - // F_Y_engine = f.get_F_Y_engine(); - // F_Z_engine = f.get_F_Z_engine(); - // F_X_gear = f.get_F_X_gear(); - // F_Y_gear = f.get_F_Y_gear(); - // F_Z_gear = f.get_F_Z_gear(); - // M_l_rp = f.get_M_l_rp(); - // M_m_rp = f.get_M_m_rp(); - // M_n_rp = f.get_M_n_rp(); - // M_l_cg = f.get_M_l_cg(); - // M_m_cg = f.get_M_m_cg(); - // M_n_cg = f.get_M_n_cg(); - // M_l_aero = f.get_M_l_aero(); - // M_m_aero = f.get_M_m_aero(); - // M_n_aero = f.get_M_n_aero(); - // M_l_engine = f.get_M_l_engine(); - // M_m_engine = f.get_M_m_engine(); - // M_n_engine = f.get_M_n_engine(); - // M_l_gear = f.get_M_l_gear(); - // M_m_gear = f.get_M_m_gear(); - // M_n_gear = f.get_M_n_gear(); - // V_dot_north = f.get_V_dot_north(); - // V_dot_east = f.get_V_dot_east(); - // V_dot_down = f.get_V_dot_down(); - // U_dot_body = f.get_U_dot_body(); - // V_dot_body = f.get_V_dot_body(); - // W_dot_body = f.get_W_dot_body(); - // A_X_cg = f.get_A_X_cg(); - // A_Y_cg = f.get_A_Y_cg(); - // A_Z_cg = f.get_A_Z_cg(); - // A_X_pilot = f.get_A_X_pilot(); - // A_Y_pilot = f.get_A_Y_pilot(); - // A_Z_pilot = f.get_A_Z_pilot(); - // N_X_cg = f.get_N_X_cg(); - // N_Y_cg = f.get_N_Y_cg(); - // N_Z_cg = f.get_N_Z_cg(); - // N_X_pilot = f.get_N_X_pilot(); - // N_Y_pilot = f.get_N_Y_pilot(); - // N_Z_pilot = f.get_N_Z_pilot(); - // P_dot_body = f.get_P_dot_body(); - // Q_dot_body = f.get_Q_dot_body(); - // R_dot_body = f.get_R_dot_body(); - V_north = f.get_V_north(); - V_east = f.get_V_east(); - V_down = f.get_V_down(); - // V_north_rel_ground = f.get_V_north_rel_ground(); - // V_east_rel_ground = f.get_V_east_rel_ground(); - // V_down_rel_ground = f.get_V_down_rel_ground(); - // V_north_airmass = f.get_V_north_airmass(); - // V_east_airmass = f.get_V_east_airmass(); - // V_down_airmass = f.get_V_down_airmass(); - // V_north_rel_airmass = f.get_V_north_rel_airmass(); - // V_east_rel_airmass = f.get_V_east_rel_airmass(); - // V_down_rel_airmass = f.get_V_down_rel_airmass(); - // U_gust = f.get_U_gust(); - // V_gust = f.get_V_gust(); - // W_gust = f.get_W_gust(); - // U_body = f.get_U_body(); - // V_body = f.get_V_body(); - // W_body = f.get_W_body(); - // V_rel_wind = f.get_V_rel_wind(); - // V_true_kts = f.get_V_true_kts(); - // V_rel_ground = f.get_V_rel_ground(); - // V_inertial = f.get_V_inertial(); - // V_ground_speed = f.get_V_ground_speed(); - // V_equiv = f.get_V_equiv(); - // V_equiv_kts = f.get_V_equiv_kts(); - // V_calibrated = f.get_V_calibrated(); - // V_calibrated_kts = f.get_V_calibrated_kts(); - P_body = f.get_P_body(); - Q_body = f.get_Q_body(); - R_body = f.get_R_body(); - // P_local = f.get_P_local(); - // Q_local = f.get_Q_local(); - // R_local = f.get_R_local(); - // P_total = f.get_P_total(); - // Q_total = f.get_Q_total(); - // R_total = f.get_R_total(); - // Phi_dot = f.get_Phi_dot(); - // Theta_dot = f.get_Theta_dot(); - // Psi_dot = f.get_Psi_dot(); - // Latitude_dot = f.get_Latitude_dot(); - // Longitude_dot = f.get_Longitude_dot(); - // Radius_dot = f.get_Radius_dot(); - Lat_geocentric = f.get_Lat_geocentric(); - Lon_geocentric = f.get_Lon_geocentric(); - Radius_to_vehicle = f.get_Radius_to_vehicle(); - Latitude = f.get_Latitude(); - Longitude = f.get_Longitude(); - Altitude = f.get_Altitude(); - Phi = f.get_Phi(); - Theta = f.get_Theta(); - Psi = f.get_Psi(); - // T_local_to_body_11 = f.get_T_local_to_body_11(); - // T_local_to_body_12 = f.get_T_local_to_body_12(); - // T_local_to_body_13 = f.get_T_local_to_body_13(); - // T_local_to_body_21 = f.get_T_local_to_body_21(); - // T_local_to_body_22 = f.get_T_local_to_body_22(); - // T_local_to_body_23 = f.get_T_local_to_body_23(); - // T_local_to_body_31 = f.get_T_local_to_body_31(); - // T_local_to_body_32 = f.get_T_local_to_body_32(); - // T_local_to_body_33 = f.get_T_local_to_body_33(); - // Gravity = f.get_Gravity(); - // Centrifugal_relief = f.get_Centrifugal_relief(); - // Alpha = f.get_Alpha(); - // Beta = f.get_Beta(); - // Alpha_dot = f.get_Alpha_dot(); - // Beta_dot = f.get_Beta_dot(); - // Cos_alpha = f.get_Cos_alpha(); - // Sin_alpha = f.get_Sin_alpha(); - // Cos_beta = f.get_Cos_beta(); - // Sin_beta = f.get_Sin_beta(); - // Cos_phi = f.get_Cos_phi(); - // Sin_phi = f.get_Sin_phi(); - // Cos_theta = f.get_Cos_theta(); - // Sin_theta = f.get_Sin_theta(); - // Cos_psi = f.get_Cos_psi(); - // Sin_psi = f.get_Sin_psi(); - // Gamma_vert_rad = f.get_Gamma_vert_rad(); - // Gamma_horiz_rad = f.get_Gamma_horiz_rad(); - // Sigma = f.get_Sigma(); - // Density = f.get_Density(); - // V_sound = f.get_V_sound(); - // Mach_number = f.get_Mach_number(); - // Static_pressure = f.get_Static_pressure(); - // Total_pressure = f.get_Total_pressure(); - // Impact_pressure = f.get_Impact_pressure(); - // Dynamic_pressure = f.get_Dynamic_pressure(); - // Static_temperature = f.get_Static_temperature(); - // Total_temperature = f.get_Total_temperature(); - Sea_level_radius = f.get_Sea_level_radius(); - Earth_position_angle = f.get_Earth_position_angle(); - Runway_altitude = f.get_Runway_altitude(); - // Runway_latitude = f.get_Runway_latitude(); - // Runway_longitude = f.get_Runway_longitude(); - // Runway_heading = f.get_Runway_heading(); - // Radius_to_rwy = f.get_Radius_to_rwy(); - // D_cg_north_of_rwy = f.get_D_cg_north_of_rwy(); - // D_cg_east_of_rwy = f.get_D_cg_east_of_rwy(); - // D_cg_above_rwy = f.get_D_cg_above_rwy(); - // X_cg_rwy = f.get_X_cg_rwy(); - // Y_cg_rwy = f.get_Y_cg_rwy(); - // H_cg_rwy = f.get_H_cg_rwy(); - // D_pilot_north_of_rwy = f.get_D_pilot_north_of_rwy(); - // D_pilot_east_of_rwy = f.get_D_pilot_east_of_rwy(); - // D_pilot_above_rwy = f.get_D_pilot_above_rwy(); - // X_pilot_rwy = f.get_X_pilot_rwy(); - // Y_pilot_rwy = f.get_Y_pilot_rwy(); - // H_pilot_rwy = f.get_H_pilot_rwy(); - - return( 0 ); +bool FGLaRCsim::copy_to_LaRCsim () { + Mass = get_Mass(); + I_xx = get_I_xx(); + I_yy = get_I_yy(); + I_zz = get_I_zz(); + I_xz = get_I_xz(); + // Dx_pilot = get_Dx_pilot(); + // Dy_pilot = get_Dy_pilot(); + // Dz_pilot = get_Dz_pilot(); + Dx_cg = get_Dx_cg(); + Dy_cg = get_Dy_cg(); + Dz_cg = get_Dz_cg(); + // F_X = get_F_X(); + // F_Y = get_F_Y(); + // F_Z = get_F_Z(); + // F_north = get_F_north(); + // F_east = get_F_east(); + // F_down = get_F_down(); + // F_X_aero = get_F_X_aero(); + // F_Y_aero = get_F_Y_aero(); + // F_Z_aero = get_F_Z_aero(); + // F_X_engine = get_F_X_engine(); + // F_Y_engine = get_F_Y_engine(); + // F_Z_engine = get_F_Z_engine(); + // F_X_gear = get_F_X_gear(); + // F_Y_gear = get_F_Y_gear(); + // F_Z_gear = get_F_Z_gear(); + // M_l_rp = get_M_l_rp(); + // M_m_rp = get_M_m_rp(); + // M_n_rp = get_M_n_rp(); + // M_l_cg = get_M_l_cg(); + // M_m_cg = get_M_m_cg(); + // M_n_cg = get_M_n_cg(); + // M_l_aero = get_M_l_aero(); + // M_m_aero = get_M_m_aero(); + // M_n_aero = get_M_n_aero(); + // M_l_engine = get_M_l_engine(); + // M_m_engine = get_M_m_engine(); + // M_n_engine = get_M_n_engine(); + // M_l_gear = get_M_l_gear(); + // M_m_gear = get_M_m_gear(); + // M_n_gear = get_M_n_gear(); + // V_dot_north = get_V_dot_north(); + // V_dot_east = get_V_dot_east(); + // V_dot_down = get_V_dot_down(); + // U_dot_body = get_U_dot_body(); + // V_dot_body = get_V_dot_body(); + // W_dot_body = get_W_dot_body(); + // A_X_cg = get_A_X_cg(); + // A_Y_cg = get_A_Y_cg(); + // A_Z_cg = get_A_Z_cg(); + // A_X_pilot = get_A_X_pilot(); + // A_Y_pilot = get_A_Y_pilot(); + // A_Z_pilot = get_A_Z_pilot(); + // N_X_cg = get_N_X_cg(); + // N_Y_cg = get_N_Y_cg(); + // N_Z_cg = get_N_Z_cg(); + // N_X_pilot = get_N_X_pilot(); + // N_Y_pilot = get_N_Y_pilot(); + // N_Z_pilot = get_N_Z_pilot(); + // P_dot_body = get_P_dot_body(); + // Q_dot_body = get_Q_dot_body(); + // R_dot_body = get_R_dot_body(); + // V_north = get_V_north(); + // V_east = get_V_east(); + // V_down = get_V_down(); + // V_north_rel_ground = get_V_north_rel_ground(); + // V_east_rel_ground = get_V_east_rel_ground(); + // V_down_rel_ground = get_V_down_rel_ground(); + // V_north_airmass = get_V_north_airmass(); + // V_east_airmass = get_V_east_airmass(); + // V_down_airmass = get_V_down_airmass(); + // V_north_rel_airmass = get_V_north_rel_airmass(); + // V_east_rel_airmass = get_V_east_rel_airmass(); + // V_down_rel_airmass = get_V_down_rel_airmass(); + // U_gust = get_U_gust(); + // V_gust = get_V_gust(); + // W_gust = get_W_gust(); + // U_body = get_U_body(); + // V_body = get_V_body(); + // W_body = get_W_body(); + // V_rel_wind = get_V_rel_wind(); + // V_true_kts = get_V_true_kts(); + // V_rel_ground = get_V_rel_ground(); + // V_inertial = get_V_inertial(); + // V_ground_speed = get_V_ground_speed(); + // V_equiv = get_V_equiv(); + // V_equiv_kts = get_V_equiv_kts(); + // V_calibrated = get_V_calibrated(); + // V_calibrated_kts = get_V_calibrated_kts(); + // P_body = get_P_body(); + // Q_body = get_Q_body(); + // R_body = get_R_body(); + // P_local = get_P_local(); + // Q_local = get_Q_local(); + // R_local = get_R_local(); + // P_total = get_P_total(); + // Q_total = get_Q_total(); + // R_total = get_R_total(); + // Phi_dot = get_Phi_dot(); + // Theta_dot = get_Theta_dot(); + // Psi_dot = get_Psi_dot(); + // Latitude_dot = get_Latitude_dot(); + // Longitude_dot = get_Longitude_dot(); + // Radius_dot = get_Radius_dot(); + // Lat_geocentric = get_Lat_geocentric(); + // Lon_geocentric = get_Lon_geocentric(); + // Radius_to_vehicle = get_Radius_to_vehicle(); + // Latitude = get_Latitude(); + // Longitude = get_Longitude(); + // Altitude = get_Altitude(); + // Phi = get_Phi(); + // Theta = get_Theta(); + // Psi = get_Psi(); + // T_local_to_body_11 = get_T_local_to_body_11(); + // T_local_to_body_12 = get_T_local_to_body_12(); + // T_local_to_body_13 = get_T_local_to_body_13(); + // T_local_to_body_21 = get_T_local_to_body_21(); + // T_local_to_body_22 = get_T_local_to_body_22(); + // T_local_to_body_23 = get_T_local_to_body_23(); + // T_local_to_body_31 = get_T_local_to_body_31(); + // T_local_to_body_32 = get_T_local_to_body_32(); + // T_local_to_body_33 = get_T_local_to_body_33(); + // Gravity = get_Gravity(); + // Centrifugal_relief = get_Centrifugal_relief(); + // Alpha = get_Alpha(); + // Beta = get_Beta(); + // Alpha_dot = get_Alpha_dot(); + // Beta_dot = get_Beta_dot(); + // Cos_alpha = get_Cos_alpha(); + // Sin_alpha = get_Sin_alpha(); + // Cos_beta = get_Cos_beta(); + // Sin_beta = get_Sin_beta(); + // Cos_phi = get_Cos_phi(); + // Sin_phi = get_Sin_phi(); + // Cos_theta = get_Cos_theta(); + // Sin_theta = get_Sin_theta(); + // Cos_psi = get_Cos_psi(); + // Sin_psi = get_Sin_psi(); + // Gamma_vert_rad = get_Gamma_vert_rad(); + // Gamma_horiz_rad = get_Gamma_horiz_rad(); + // Sigma = get_Sigma(); + // Density = get_Density(); + // V_sound = get_V_sound(); + // Mach_number = get_Mach_number(); + // Static_pressure = get_Static_pressure(); + // Total_pressure = get_Total_pressure(); + // Impact_pressure = get_Impact_pressure(); + // Dynamic_pressure = get_Dynamic_pressure(); + // Static_temperature = get_Static_temperature(); + // Total_temperature = get_Total_temperature(); + // Sea_level_radius = get_Sea_level_radius(); + // Earth_position_angle = get_Earth_position_angle(); + // Runway_altitude = get_Runway_altitude(); + // Runway_latitude = get_Runway_latitude(); + // Runway_longitude = get_Runway_longitude(); + // Runway_heading = get_Runway_heading(); + // Radius_to_rwy = get_Radius_to_rwy(); + // D_cg_north_of_rwy = get_D_cg_north_of_rwy(); + // D_cg_east_of_rwy = get_D_cg_east_of_rwy(); + // D_cg_above_rwy = get_D_cg_above_rwy(); + // X_cg_rwy = get_X_cg_rwy(); + // Y_cg_rwy = get_Y_cg_rwy(); + // H_cg_rwy = get_H_cg_rwy(); + // D_pilot_north_of_rwy = get_D_pilot_north_of_rwy(); + // D_pilot_east_of_rwy = get_D_pilot_east_of_rwy(); + // D_pilot_above_rwy = get_D_pilot_above_rwy(); + // X_pilot_rwy = get_X_pilot_rwy(); + // Y_pilot_rwy = get_Y_pilot_rwy(); + // H_pilot_rwy = get_H_pilot_rwy(); + + return true; } // Convert from the LaRCsim generic_ struct to the FGInterface struct -int fgLaRCsim_2_FGInterface (FGInterface& f) { +bool FGLaRCsim::copy_from_LaRCsim() { // Mass properties and geometry values - f.set_Inertias( Mass, I_xx, I_yy, I_zz, I_xz ); - // f.set_Pilot_Location( Dx_pilot, Dy_pilot, Dz_pilot ); - f.set_CG_Position( Dx_cg, Dy_cg, Dz_cg ); + _set_Inertias( Mass, I_xx, I_yy, I_zz, I_xz ); + // set_Pilot_Location( Dx_pilot, Dy_pilot, Dz_pilot ); + _set_CG_Position( Dx_cg, Dy_cg, Dz_cg ); // Forces - // f.set_Forces_Body_Total( F_X, F_Y, F_Z ); - // f.set_Forces_Local_Total( F_north, F_east, F_down ); - // f.set_Forces_Aero( F_X_aero, F_Y_aero, F_Z_aero ); - // f.set_Forces_Engine( F_X_engine, F_Y_engine, F_Z_engine ); - // f.set_Forces_Gear( F_X_gear, F_Y_gear, F_Z_gear ); + // set_Forces_Body_Total( F_X, F_Y, F_Z ); + // set_Forces_Local_Total( F_north, F_east, F_down ); + // set_Forces_Aero( F_X_aero, F_Y_aero, F_Z_aero ); + // set_Forces_Engine( F_X_engine, F_Y_engine, F_Z_engine ); + // set_Forces_Gear( F_X_gear, F_Y_gear, F_Z_gear ); // Moments - // f.set_Moments_Total_RP( M_l_rp, M_m_rp, M_n_rp ); - // f.set_Moments_Total_CG( M_l_cg, M_m_cg, M_n_cg ); - // f.set_Moments_Aero( M_l_aero, M_m_aero, M_n_aero ); - // f.set_Moments_Engine( M_l_engine, M_m_engine, M_n_engine ); - // f.set_Moments_Gear( M_l_gear, M_m_gear, M_n_gear ); + // set_Moments_Total_RP( M_l_rp, M_m_rp, M_n_rp ); + // set_Moments_Total_CG( M_l_cg, M_m_cg, M_n_cg ); + // set_Moments_Aero( M_l_aero, M_m_aero, M_n_aero ); + // set_Moments_Engine( M_l_engine, M_m_engine, M_n_engine ); + // set_Moments_Gear( M_l_gear, M_m_gear, M_n_gear ); // Accelerations - // f.set_Accels_Local( V_dot_north, V_dot_east, V_dot_down ); - // f.set_Accels_Body( U_dot_body, V_dot_body, W_dot_body ); - // f.set_Accels_CG_Body( A_X_cg, A_Y_cg, A_Z_cg ); - // f.set_Accels_Pilot_Body( A_X_pilot, A_Y_pilot, A_Z_pilot ); - // f.set_Accels_CG_Body_N( N_X_cg, N_Y_cg, N_Z_cg ); - // f.set_Accels_Pilot_Body_N( N_X_pilot, N_Y_pilot, N_Z_pilot ); - // f.set_Accels_Omega( P_dot_body, Q_dot_body, R_dot_body ); + _set_Accels_Local( V_dot_north, V_dot_east, V_dot_down ); + _set_Accels_Body( U_dot_body, V_dot_body, W_dot_body ); + _set_Accels_CG_Body( A_X_cg, A_Y_cg, A_Z_cg ); + _set_Accels_Pilot_Body( A_X_pilot, A_Y_pilot, A_Z_pilot ); + // set_Accels_CG_Body_N( N_X_cg, N_Y_cg, N_Z_cg ); + // set_Accels_Pilot_Body_N( N_X_pilot, N_Y_pilot, N_Z_pilot ); + // set_Accels_Omega( P_dot_body, Q_dot_body, R_dot_body ); // Velocities - f.set_Velocities_Local( V_north, V_east, V_down ); - // f.set_Velocities_Ground( V_north_rel_ground, V_east_rel_ground, + _set_Velocities_Local( V_north, V_east, V_down ); + // set_Velocities_Ground( V_north_rel_ground, V_east_rel_ground, // V_down_rel_ground ); - // f.set_Velocities_Local_Airmass( V_north_airmass, V_east_airmass, - // V_down_airmass ); - // f.set_Velocities_Local_Rel_Airmass( V_north_rel_airmass, + _set_Velocities_Local_Airmass( V_north_airmass, V_east_airmass, + V_down_airmass ); + // set_Velocities_Local_Rel_Airmass( V_north_rel_airmass, // V_east_rel_airmass, V_down_rel_airmass ); - // f.set_Velocities_Gust( U_gust, V_gust, W_gust ); - // f.set_Velocities_Wind_Body( U_body, V_body, W_body ); - - // f.set_V_rel_wind( V_rel_wind ); - // f.set_V_true_kts( V_true_kts ); - // f.set_V_rel_ground( V_rel_ground ); - // f.set_V_inertial( V_inertial ); - // f.set_V_ground_speed( V_ground_speed ); - // f.set_V_equiv( V_equiv ); - f.set_V_equiv_kts( V_equiv_kts ); - // f.set_V_calibrated( V_calibrated ); - // f.set_V_calibrated_kts( V_calibrated_kts ); - - f.set_Omega_Body( P_body, Q_body, R_body ); - // f.set_Omega_Local( P_local, Q_local, R_local ); - // f.set_Omega_Total( P_total, Q_total, R_total ); - - // f.set_Euler_Rates( Phi_dot, Theta_dot, Psi_dot ); - f.set_Geocentric_Rates( Latitude_dot, Longitude_dot, Radius_dot ); + // set_Velocities_Gust( U_gust, V_gust, W_gust ); + _set_Velocities_Wind_Body( U_body, V_body, W_body ); + + _set_V_rel_wind( V_rel_wind ); + // set_V_true_kts( V_true_kts ); + // set_V_rel_ground( V_rel_ground ); + // set_V_inertial( V_inertial ); + _set_V_ground_speed( V_ground_speed ); + // set_V_equiv( V_equiv ); + _set_V_equiv_kts( V_equiv_kts ); + // set_V_calibrated( V_calibrated ); + _set_V_calibrated_kts( V_calibrated_kts ); + + _set_Omega_Body( P_body, Q_body, R_body ); + // set_Omega_Local( P_local, Q_local, R_local ); + // set_Omega_Total( P_total, Q_total, R_total ); + + _set_Euler_Rates( Phi_dot, Theta_dot, Psi_dot ); + _set_Geocentric_Rates( Latitude_dot, Longitude_dot, Radius_dot ); + + _set_Mach_number( Mach_number ); FG_LOG( FG_FLIGHT, FG_DEBUG, "lon = " << Longitude << " lat_geoc = " << Lat_geocentric << " lat_geod = " << Latitude << " alt = " << Altitude << " sl_radius = " << Sea_level_radius << " radius_to_vehicle = " << Radius_to_vehicle ); - + + double tmp_lon_geoc = Lon_geocentric; + while ( tmp_lon_geoc < -FG_PI ) { tmp_lon_geoc += FG_2PI; } + while ( tmp_lon_geoc > FG_PI ) { tmp_lon_geoc -= FG_2PI; } + + double tmp_lon = Longitude; + while ( tmp_lon < -FG_PI ) { tmp_lon += FG_2PI; } + while ( tmp_lon > FG_PI ) { tmp_lon -= FG_2PI; } + // Positions - f.set_Geocentric_Position( Lat_geocentric, Lon_geocentric, - Radius_to_vehicle ); - f.set_Geodetic_Position( Latitude, Longitude, Altitude ); - f.set_Euler_Angles( Phi, Theta, Psi ); + _set_Geocentric_Position( Lat_geocentric, tmp_lon_geoc, + Radius_to_vehicle ); + _set_Geodetic_Position( Latitude, tmp_lon, Altitude ); + _set_Euler_Angles( Phi, Theta, Psi ); + + _set_Altitude_AGL( Altitude-Runway_altitude ); // Miscellaneous quantities - f.set_T_Local_to_Body(T_local_to_body_m); - // f.set_Gravity( Gravity ); - // f.set_Centrifugal_relief( Centrifugal_relief ); - - f.set_Alpha( Alpha ); - f.set_Beta( Beta ); - // f.set_Alpha_dot( Alpha_dot ); - // f.set_Beta_dot( Beta_dot ); - - // f.set_Cos_alpha( Cos_alpha ); - // f.set_Sin_alpha( Sin_alpha ); - // f.set_Cos_beta( Cos_beta ); - // f.set_Sin_beta( Sin_beta ); - - // f.set_Cos_phi( Cos_phi ); - // f.set_Sin_phi( Sin_phi ); - // f.set_Cos_theta( Cos_theta ); - // f.set_Sin_theta( Sin_theta ); - // f.set_Cos_psi( Cos_psi ); - // f.set_Sin_psi( Sin_psi ); - - f.set_Gamma_vert_rad( Gamma_vert_rad ); - // f.set_Gamma_horiz_rad( Gamma_horiz_rad ); - - // f.set_Sigma( Sigma ); - // f.set_Density( Density ); - // f.set_V_sound( V_sound ); - // f.set_Mach_number( Mach_number ); - - // f.set_Static_pressure( Static_pressure ); - // f.set_Total_pressure( Total_pressure ); - // f.set_Impact_pressure( Impact_pressure ); - // f.set_Dynamic_pressure( Dynamic_pressure ); - - // f.set_Static_temperature( Static_temperature ); - // f.set_Total_temperature( Total_temperature ); - - f.set_Sea_level_radius( Sea_level_radius ); - f.set_Earth_position_angle( Earth_position_angle ); - - f.set_Runway_altitude( Runway_altitude ); - // f.set_Runway_latitude( Runway_latitude ); - // f.set_Runway_longitude( Runway_longitude ); - // f.set_Runway_heading( Runway_heading ); - // f.set_Radius_to_rwy( Radius_to_rwy ); - - // f.set_CG_Rwy_Local( D_cg_north_of_rwy, D_cg_east_of_rwy, D_cg_above_rwy); - // f.set_CG_Rwy_Rwy( X_cg_rwy, Y_cg_rwy, H_cg_rwy ); - // f.set_Pilot_Rwy_Local( D_pilot_north_of_rwy, D_pilot_east_of_rwy, + _set_T_Local_to_Body(T_local_to_body_m); + // set_Gravity( Gravity ); + // set_Centrifugal_relief( Centrifugal_relief ); + + _set_Alpha( Alpha ); + _set_Beta( Beta ); + // set_Alpha_dot( Alpha_dot ); + // set_Beta_dot( Beta_dot ); + + // set_Cos_alpha( Cos_alpha ); + // set_Sin_alpha( Sin_alpha ); + // set_Cos_beta( Cos_beta ); + // set_Sin_beta( Sin_beta ); + + _set_Cos_phi( Cos_phi ); + // set_Sin_phi( Sin_phi ); + _set_Cos_theta( Cos_theta ); + // set_Sin_theta( Sin_theta ); + // set_Cos_psi( Cos_psi ); + // set_Sin_psi( Sin_psi ); + + _set_Gamma_vert_rad( Gamma_vert_rad ); + // set_Gamma_horiz_rad( Gamma_horiz_rad ); + + // set_Sigma( Sigma ); + _set_Density( Density ); + // set_V_sound( V_sound ); + // set_Mach_number( Mach_number ); + + _set_Static_pressure( Static_pressure ); + // set_Total_pressure( Total_pressure ); + // set_Impact_pressure( Impact_pressure ); + // set_Dynamic_pressure( Dynamic_pressure ); + + _set_Static_temperature( Static_temperature ); + // set_Total_temperature( Total_temperature ); + + _set_Sea_level_radius( Sea_level_radius ); + _set_Earth_position_angle( Earth_position_angle ); + + _set_Runway_altitude( Runway_altitude ); + // set_Runway_latitude( Runway_latitude ); + // set_Runway_longitude( Runway_longitude ); + // set_Runway_heading( Runway_heading ); + // set_Radius_to_rwy( Radius_to_rwy ); + + // set_CG_Rwy_Local( D_cg_north_of_rwy, D_cg_east_of_rwy, D_cg_above_rwy); + // set_CG_Rwy_Rwy( X_cg_rwy, Y_cg_rwy, H_cg_rwy ); + // set_Pilot_Rwy_Local( D_pilot_north_of_rwy, D_pilot_east_of_rwy, // D_pilot_above_rwy ); - // f.set_Pilot_Rwy_Rwy( X_pilot_rwy, Y_pilot_rwy, H_pilot_rwy ); + // set_Pilot_Rwy_Rwy( X_pilot_rwy, Y_pilot_rwy, H_pilot_rwy ); - f.set_sin_lat_geocentric(Lat_geocentric); - f.set_cos_lat_geocentric(Lat_geocentric); - f.set_sin_cos_longitude(Longitude); - f.set_sin_cos_latitude(Latitude); + _set_sin_lat_geocentric(Lat_geocentric); + _set_cos_lat_geocentric(Lat_geocentric); + _set_sin_cos_longitude(Longitude); + _set_sin_cos_latitude(Latitude); // printf("sin_lat_geo %f cos_lat_geo %f\n", sin_Lat_geoc, cos_Lat_geoc); // printf("sin_lat %f cos_lat %f\n", - // f.get_sin_latitude(), f.get_cos_latitude()); + // get_sin_latitude(), get_cos_latitude()); // printf("sin_lon %f cos_lon %f\n", - // f.get_sin_longitude(), f.get_cos_longitude()); + // get_sin_longitude(), get_cos_longitude()); - return 0; + _set_Climb_Rate( -1 * V_down ); + + return true; +} + + +void FGLaRCsim::set_ls(void) { + Phi=lsic->GetRollAngleRadIC(); + Theta=lsic->GetPitchAngleRadIC(); + Psi=lsic->GetHeadingRadIC(); + V_north=lsic->GetVnorthFpsIC(); + V_east=lsic->GetVeastFpsIC(); + V_down=lsic->GetVdownFpsIC(); + Altitude=lsic->GetAltitudeFtIC(); + Latitude=lsic->GetLatitudeGDRadIC(); + Longitude=lsic->GetLongitudeRadIC(); + Runway_altitude=lsic->GetRunwayAltitudeFtIC(); + V_north_airmass = lsic->GetVnorthAirmassFpsIC(); + V_east_airmass = lsic->GetVeastAirmassFpsIC(); + V_down_airmass = lsic->GetVdownAirmassFpsIC(); + ls_loop(0.0,-1); + copy_from_LaRCsim(); + FG_LOG( FG_FLIGHT, FG_INFO, " FGLaRCsim::set_ls(): " ); + FG_LOG( FG_FLIGHT, FG_INFO, " Phi: " << Phi ); + FG_LOG( FG_FLIGHT, FG_INFO, " Theta: " << Theta ); + FG_LOG( FG_FLIGHT, FG_INFO, " Psi: " << Psi ); + FG_LOG( FG_FLIGHT, FG_INFO, " V_north: " << V_north ); + FG_LOG( FG_FLIGHT, FG_INFO, " V_east: " << V_east ); + FG_LOG( FG_FLIGHT, FG_INFO, " V_down: " << V_down ); + FG_LOG( FG_FLIGHT, FG_INFO, " Altitude: " << Altitude ); + FG_LOG( FG_FLIGHT, FG_INFO, " Latitude: " << Latitude ); + FG_LOG( FG_FLIGHT, FG_INFO, " Longitude: " << Longitude ); + FG_LOG( FG_FLIGHT, FG_INFO, " Runway_altitude: " << Runway_altitude ); + FG_LOG( FG_FLIGHT, FG_INFO, " V_north_airmass: " << V_north_airmass ); + FG_LOG( FG_FLIGHT, FG_INFO, " V_east_airmass: " << V_east_airmass ); + FG_LOG( FG_FLIGHT, FG_INFO, " V_down_airmass: " << V_down_airmass ); +} + +void FGLaRCsim::snap_shot(void) { + lsic->SetLatitudeGDRadIC( get_Latitude() ); + lsic->SetLongitudeRadIC( get_Longitude() ); + lsic->SetAltitudeFtIC( get_Altitude() ); + lsic->SetRunwayAltitudeFtIC( get_Runway_altitude() ); + lsic->SetVtrueFpsIC( get_V_rel_wind() ); + lsic->SetPitchAngleRadIC( get_Theta() ); + lsic->SetRollAngleRadIC( get_Phi() ); + lsic->SetHeadingRadIC( get_Psi() ); + lsic->SetClimbRateFpsIC( get_Climb_Rate() ); + lsic->SetVNEDAirmassFpsIC( get_V_north_airmass(), + get_V_east_airmass(), + get_V_down_airmass() ); +} + +//Positions +void FGLaRCsim::set_Latitude(double lat) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Latitude: " << lat ); + snap_shot(); + lsic->SetLatitudeGDRadIC(lat); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_Longitude(double lon) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Longitude: " << lon ); + snap_shot(); + lsic->SetLongitudeRadIC(lon); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_Altitude(double alt) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Altitude: " << alt ); + snap_shot(); + lsic->SetAltitudeFtIC(alt); + set_ls(); + copy_from_LaRCsim(); //update the bus } +void FGLaRCsim::set_V_calibrated_kts(double vc) { + FG_LOG( FG_FLIGHT, FG_INFO, + "FGLaRCsim::set_V_calibrated_kts: " << vc ); + snap_shot(); + lsic->SetVcalibratedKtsIC(vc); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_Mach_number(double mach) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Mach_number: " << mach ); + snap_shot(); + lsic->SetMachIC(mach); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_Velocities_Local( double north, double east, double down ){ + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Velocities_local: " + << north << " " << east << " " << down ); + snap_shot(); + lsic->SetVNEDFpsIC(north, east, down); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_Velocities_Wind_Body( double u, double v, double w){ + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Velocities_Wind_Body: " + << u << " " << v << " " << w ); + snap_shot(); + lsic->SetUVWFpsIC(u,v,w); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +//Euler angles +void FGLaRCsim::set_Euler_Angles( double phi, double theta, double psi ) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Euler_angles: " + << phi << " " << theta << " " << psi ); + + snap_shot(); + lsic->SetPitchAngleRadIC(theta); + lsic->SetRollAngleRadIC(phi); + lsic->SetHeadingRadIC(psi); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +//Flight Path +void FGLaRCsim::set_Climb_Rate( double roc) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Climb_rate: " << roc ); + snap_shot(); + lsic->SetClimbRateFpsIC(roc); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_Gamma_vert_rad( double gamma) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Gamma_vert_rad: " << gamma ); + snap_shot(); + lsic->SetFlightPathAngleRadIC(gamma); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_Runway_altitude(double ralt) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Runway_altitude: " << ralt ); + snap_shot(); + lsic->SetRunwayAltitudeFtIC(ralt); + set_ls(); + copy_from_LaRCsim(); //update the bus +} + +void FGLaRCsim::set_AltitudeAGL(double altagl) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_AltitudeAGL: " << altagl ); + snap_shot(); + lsic->SetAltitudeAGLFtIC(altagl); + set_ls(); + copy_from_LaRCsim(); +} + +void FGLaRCsim::set_Velocities_Local_Airmass (double wnorth, + double weast, + double wdown ) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Velocities_Local_Airmass: " + << wnorth << " " << weast << " " << wdown ); + snap_shot(); + lsic->SetVNEDAirmassFpsIC( wnorth, weast, wdown ); + set_ls(); + copy_from_LaRCsim(); +} + +void FGLaRCsim::set_Static_pressure(double p) { + FG_LOG( FG_FLIGHT, FG_INFO, + "FGLaRCsim::set_Static_pressure: " << p ); + FG_LOG( FG_FLIGHT, FG_INFO, + "LaRCsim does not support externally supplied atmospheric data" ); +} + +void FGLaRCsim::set_Static_temperature(double T) { + FG_LOG( FG_FLIGHT, FG_INFO, + "FGLaRCsim::set_Static_temperature: " << T ); + FG_LOG( FG_FLIGHT, FG_INFO, + "LaRCsim does not support externally supplied atmospheric data" ); + +} + +void FGLaRCsim::set_Density(double rho) { + FG_LOG( FG_FLIGHT, FG_INFO, "FGLaRCsim::set_Density: " << rho ); + FG_LOG( FG_FLIGHT, FG_INFO, + "LaRCsim does not support externally supplied atmospheric data" ); + +}