src/Cockpit/radiostack.cxx

   1 // radiostack.cxx -- class to manage an instance of the radio stack
   2 //
   3 // Written by Curtis Olson, started April 2000.
   4 //
   5 // Copyright (C) 2000  Curtis L. Olson - curt@flightgear.org
   6 //
   7 // This program is free software; you can redistribute it and/or
   8 // modify it under the terms of the GNU General Public License as
   9 // published by the Free Software Foundation; either version 2 of the
  10 // License, or (at your option) any later version.
  11 //
  12 // This program is distributed in the hope that it will be useful, but
  13 // WITHOUT ANY WARRANTY; without even the implied warranty of
  14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 // General Public License for more details.
  16 //
  17 // You should have received a copy of the GNU General Public License
  18 // along with this program; if not, write to the Free Software
  19 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  20 //
  21 // $Id$
  22
  23
  24 #include <Aircraft/aircraft.hxx>
  25 #include <Main/bfi.hxx>
  26 #include <Navaids/ilslist.hxx>
  27 #include <Navaids/navlist.hxx>
  28 #include <Time/event.hxx>
  29
  30 #include "radiostack.hxx"
  31
  32 /**
  33  * Boy, this is ugly!  Make the VOR range vary by altitude difference.
  34  */
  35 static inline double
  36 kludgeRange (double stationElev, double aircraftElev, double nominalRange)
  37 {
  38                                 // Assume that the nominal range (usually
  39                                 // 50nm) applies at a 5,000 ft difference.
  40                                 // Just a wild guess!
  41   double factor = ((aircraftElev*METER_TO_FEET) - stationElev) / 1000.0;
  42   double range = fabs(nominalRange * factor);
  43
  44                                 // Clamp the range to keep it sane; for
  45                                 // now, never less than 25% or more than
  46                                 // 500% of nominal range.
  47   if (range < nominalRange/4.0) {
  48     range = nominalRange/4.0;
  49   } else if (range > nominalRange*5.0) {
  50     range = nominalRange*5.0;
  51   }
  52
  53   return range;
  54 }
  55
  56
  57 FGRadioStack *current_radiostack;
  58
  59
  60 // periodic radio station search wrapper
  61 static void fgRadioSearch( void ) {
  62     current_radiostack->search();
  63 }
  64
  65 // Constructor
  66 FGRadioStack::FGRadioStack() {
  67     nav1_radial = 0.0;
  68     nav1_dme_dist = 0.0;
  69     nav2_radial = 0.0;
  70     nav2_dme_dist = 0.0;
  71     need_update = true;
  72     longitudeVal = fgGetValue("/position/longitude");
  73     latitudeVal = fgGetValue("/position/latitude");
  74     altitudeVal = fgGetValue("/position/altitude");
  75 }
  76
  77
  78 // Destructor
  79 FGRadioStack::~FGRadioStack()
  80 {
  81     unbind();                   // FIXME: should be called externally
  82 }
  83
  84
  85 void
  86 FGRadioStack::init ()
  87 {
  88     search();
  89     update();
  90
  91     // Search radio database once per second
  92     global_events.Register( "fgRadioSearch()", fgRadioSearch,
  93                             fgEVENT::FG_EVENT_READY, 1000);
  94 }
  95
  96 void
  97 FGRadioStack::bind ()
  98 {
  99                                 // User inputs
 100     fgTie("/radios/nav1/frequencies/selected", this,
 101           &FGRadioStack::get_nav1_freq, &FGRadioStack::set_nav1_freq);
 102     fgTie("/radios/nav1/frequencies/standby", this,
 103           &FGRadioStack::get_nav1_alt_freq, &FGRadioStack::set_nav1_alt_freq);
 104     fgTie("/radios/nav1/radials/selected", this,
 105           &FGRadioStack::get_nav1_sel_radial,
 106           &FGRadioStack::set_nav1_sel_radial);
 107
 108                                 // Radio outputs
 109     fgTie("/radios/nav1/radials/actual", this, &FGRadioStack::get_nav1_radial);
 110     fgTie("/radios/nav1/to-flag", this, &FGRadioStack::get_nav1_to_flag);
 111     fgTie("/radios/nav1/from-flag", this, &FGRadioStack::get_nav1_from_flag);
 112     fgTie("/radios/nav1/in-range", this, &FGRadioStack::get_nav1_inrange);
 113     fgTie("/radios/nav1/dme/distance", this, &FGRadioStack::get_nav1_dme_dist);
 114     fgTie("/radios/nav1/dme/in-range", this,
 115           &FGRadioStack::get_nav1_dme_inrange);
 116
 117                                 // User inputs
 118     fgTie("/radios/nav2/frequencies/selected", this,
 119           &FGRadioStack::get_nav2_freq, &FGRadioStack::set_nav2_freq);
 120     fgTie("/radios/nav2/frequencies/standby", this,
 121           &FGRadioStack::get_nav2_alt_freq, &FGRadioStack::set_nav2_alt_freq);
 122     fgTie("/radios/nav2/radials/selected", this,
 123           &FGRadioStack::get_nav2_sel_radial,
 124           &FGRadioStack::set_nav2_sel_radial);
 125
 126                                 // Radio outputs
 127     fgTie("/radios/nav2/radials/actual", this, &FGRadioStack::get_nav2_radial);
 128     fgTie("/radios/nav2/to-flag", this, &FGRadioStack::get_nav2_to_flag);
 129     fgTie("/radios/nav2/from-flag", this, &FGRadioStack::get_nav2_from_flag);
 130     fgTie("/radios/nav2/in-range", this, &FGRadioStack::get_nav2_inrange);
 131     fgTie("/radios/nav2/dme/distance", this, &FGRadioStack::get_nav2_dme_dist);
 132     fgTie("/radios/nav2/dme/in-range", this,
 133           &FGRadioStack::get_nav2_dme_inrange);
 134
 135                                 // User inputs
 136     fgTie("/radios/adf/frequencies/selected", this,
 137           &FGRadioStack::get_adf_freq, &FGRadioStack::set_adf_freq);
 138     fgTie("/radios/adf/frequencies/standby", this,
 139           &FGRadioStack::get_adf_alt_freq, &FGRadioStack::set_adf_alt_freq);
 140     fgTie("/radios/adf/rotation", this,
 141           &FGRadioStack::get_adf_rotation, &FGRadioStack::set_adf_rotation);
 142 }
 143
 144 void
 145 FGRadioStack::unbind ()
 146 {
 147     fgUntie("/radios/nav1/frequencies/selected");
 148     fgUntie("/radios/nav1/frequencies/standby");
 149     fgUntie("/radios/nav1/radials/actual");
 150     fgUntie("/radios/nav1/radials/selected");
 151     fgUntie("/radios/nav1/to-flag");
 152     fgUntie("/radios/nav1/from-flag");
 153     fgUntie("/radios/nav1/in-range");
 154     fgUntie("/radios/nav1/dme/distance");
 155     fgUntie("/radios/nav1/dme/in-range");
 156
 157     fgUntie("/radios/nav2/frequencies/selected");
 158     fgUntie("/radios/nav2/frequencies/standby");
 159     fgUntie("/radios/nav2/radials/actual");
 160     fgUntie("/radios/nav2/radials/selected");
 161     fgUntie("/radios/nav2/to-flag");
 162     fgUntie("/radios/nav2/from-flag");
 163     fgUntie("/radios/nav2/in-range");
 164     fgUntie("/radios/nav2/dme/distance");
 165     fgUntie("/radios/nav2/dme/in-range");
 166
 167     fgUntie("/radios/adf/frequencies/selected");
 168     fgUntie("/radios/adf/frequencies/standby");
 169     fgUntie("/radios/adf/rotation");
 170 }
 171
 172 // Search the database for the current frequencies given current location
 173 void
 174 FGRadioStack::update()
 175 {
 176     double lon = longitudeVal->getDoubleValue() * DEG_TO_RAD;
 177     double lat = latitudeVal->getDoubleValue() * DEG_TO_RAD;
 178     double elev = altitudeVal->getDoubleValue() * FEET_TO_METER;
 179
 180     need_update = false;
 181
 182     Point3D aircraft = sgGeodToCart( Point3D( lon, lat, elev ) );
 183     Point3D station;
 184     double az1, az2, s;
 185
 186     if ( nav1_valid ) {
 187         station = Point3D( nav1_x, nav1_y, nav1_z );
 188         nav1_loc_dist = aircraft.distance3D( station );
 189
 190         if ( nav1_has_dme ) {
 191             // staightline distance
 192             station = Point3D( nav1_dme_x, nav1_dme_y, nav1_dme_z );
 193             nav1_dme_dist = aircraft.distance3D( station );
 194         } else {
 195             nav1_dme_dist = 0.0;
 196         }
 197
 198         if ( nav1_has_gs ) {
 199             station = Point3D( nav1_gs_x, nav1_gs_y, nav1_gs_z );
 200             nav1_gs_dist = aircraft.distance3D( station );
 201         } else {
 202             nav1_gs_dist = 0.0;
 203         }
 204
 205         if ( nav1_loc_dist < nav1_effective_range * NM_TO_METER ) {
 206             nav1_inrange = true;
 207
 208             // wgs84 heading
 209             geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
 210                                 nav1_loclat, nav1_loclon,
 211                                 &az1, &az2, &s );
 212             cout << "az1 = " << az1 << " magvar = " << nav1_magvar << endl;
 213             nav1_heading = az1 - nav1_magvar;
 214             // Alex: nav1_heading = - (az1 - FGBFI::getMagVar() / RAD_TO_DEG);
 215
 216             // cout << " heading = " << nav1_heading
 217             //      << " dist = " << nav1_dist << endl;
 218         } else {
 219             nav1_inrange = false;
 220         }
 221
 222         if ( nav1_loc ) {
 223         } else {
 224             nav1_radial = nav1_sel_radial;
 225         }
 226     } else {
 227         nav1_inrange = false;
 228         nav1_dme_dist = 0.0;
 229         // cout << "not picking up vor. :-(" << endl;
 230     }
 231
 232     if ( nav2_valid ) {
 233         station = Point3D( nav2_x, nav2_y, nav2_z );
 234         nav2_loc_dist = aircraft.distance3D( station );
 235
 236         if ( nav2_has_dme ) {
 237             // staightline distance
 238             station = Point3D( nav2_dme_x, nav2_dme_y, nav2_dme_z );
 239             nav2_dme_dist = aircraft.distance3D( station );
 240         } else {
 241             nav2_dme_dist = 0.0;
 242         }
 243
 244         if ( nav2_has_gs ) {
 245             station = Point3D( nav2_gs_x, nav2_gs_y, nav2_gs_z );
 246             nav2_gs_dist = aircraft.distance3D( station );
 247         } else {
 248             nav2_gs_dist = 0.0;
 249         }
 250
 251         if ( nav2_loc_dist < nav2_effective_range * NM_TO_METER ) {
 252             nav2_inrange = true;
 253
 254             // wgs84 heading
 255             geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
 256                                 nav2_loclat, nav2_loclon,
 257                                 &az1, &az2, &s );
 258             nav2_heading = az1 - nav2_magvar;
 259             // Alex: nav2_heading = - (az1 - FGBFI::getMagVar() / RAD_TO_DEG);
 260
 261             // cout << " heading = " << nav2_heading
 262             //      << " dist = " << nav2_dist << endl;
 263         } else {
 264             nav2_inrange = false;
 265         }
 266
 267         if ( !nav2_loc ) {
 268             nav2_radial = nav2_sel_radial;
 269         }
 270     } else {
 271         nav2_inrange = false;
 272         // cout << "not picking up vor. :-(" << endl;
 273     }
 274
 275     // adf
 276     if ( adf_valid ) {
 277         // staightline distance
 278         station = Point3D( adf_x, adf_y, adf_z );
 279         adf_dist = aircraft.distance3D( station );
 280
 281         if ( adf_dist < adf_effective_range * NM_TO_METER ) {
 282             adf_inrange = true;
 283
 284             // wgs84 heading
 285             geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
 286                                 adf_lat, adf_lon,
 287                                 &az1, &az2, &s );
 288             adf_heading = az1;
 289
 290             // cout << " heading = " << nav2_heading
 291             //      << " dist = " << nav2_dist << endl;
 292         } else {
 293             adf_inrange = false;
 294         }
 295     } else {
 296         adf_inrange = false;
 297     }
 298 }
 299
 300
 301 // Update current nav/adf radio stations based on current postition
 302 void FGRadioStack::search ()
 303 {
 304     double lon = longitudeVal->getDoubleValue() * DEG_TO_RAD;
 305     double lat = latitudeVal->getDoubleValue() * DEG_TO_RAD;
 306     double elev = altitudeVal->getDoubleValue() * FEET_TO_METER;
 307
 308     // nav1
 309     FGILS ils;
 310     FGNav nav;
 311
 312     if ( current_ilslist->query( lon, lat, elev, nav1_freq, &ils ) ) {
 313         nav1_valid = true;
 314         nav1_loc = true;
 315         nav1_has_dme = ils.get_has_dme();
 316         nav1_has_gs = ils.get_has_gs();
 317
 318         nav1_loclon = ils.get_loclon();
 319         nav1_loclat = ils.get_loclat();
 320         nav1_gslon = ils.get_gslon();
 321         nav1_gslat = ils.get_gslat();
 322         nav1_dmelon = ils.get_dmelon();
 323         nav1_dmelat = ils.get_dmelat();
 324         nav1_elev = ils.get_gselev();
 325         nav1_magvar = 0;
 326         nav1_effective_range = FG_ILS_DEFAULT_RANGE;
 327         nav1_target_gs = ils.get_gsangle();
 328         nav1_radial = ils.get_locheading();
 329         while ( nav1_radial <   0.0 ) { nav1_radial += 360.0; }
 330         while ( nav1_radial > 360.0 ) { nav1_radial -= 360.0; }
 331         nav1_x = ils.get_x();
 332         nav1_y = ils.get_y();
 333         nav1_z = ils.get_z();
 334         nav1_gs_x = ils.get_gs_x();
 335         nav1_gs_y = ils.get_gs_y();
 336         nav1_gs_z = ils.get_gs_z();
 337         nav1_dme_x = ils.get_dme_x();
 338         nav1_dme_y = ils.get_dme_y();
 339         nav1_dme_z = ils.get_dme_z();
 340         // cout << "Found an ils station in range" << endl;
 341         // cout << " id = " << ils.get_locident() << endl;
 342     } else if ( current_navlist->query( lon, lat, elev, nav1_freq, &nav ) ) {
 343         nav1_valid = true;
 344         nav1_loc = false;
 345         nav1_has_dme = nav.get_has_dme();
 346         nav1_has_gs = false;
 347         nav1_loclon = nav.get_lon();
 348         nav1_loclat = nav.get_lat();
 349         nav1_elev = nav.get_elev();
 350         nav1_magvar = nav.get_magvar();
 351         nav1_effective_range = kludgeRange(nav1_elev, elev, nav.get_range());
 352         nav1_target_gs = 0.0;
 353         nav1_radial = nav1_sel_radial;
 354         nav1_x = nav1_dme_x = nav.get_x();
 355         nav1_y = nav1_dme_y = nav.get_y();
 356         nav1_z = nav1_dme_z = nav.get_z();
 357         // cout << "Found a vor station in range" << endl;
 358         // cout << " id = " << nav.get_ident() << endl;
 359     } else {
 360         nav1_valid = false;
 361         nav1_radial = 0;
 362         nav1_dme_dist = 0;
 363         // cout << "not picking up vor1. :-(" << endl;
 364     }
 365
 366     if ( current_ilslist->query( lon, lat, elev, nav2_freq, &ils ) ) {
 367         nav2_valid = true;
 368         nav2_loc = true;
 369         nav2_has_dme = ils.get_has_dme();
 370         nav2_has_gs = ils.get_has_gs();
 371
 372         nav2_loclon = ils.get_loclon();
 373         nav2_loclat = ils.get_loclat();
 374         nav2_elev = ils.get_gselev();
 375         nav2_magvar = 0;
 376         nav2_effective_range = FG_ILS_DEFAULT_RANGE;
 377         nav2_target_gs = ils.get_gsangle();
 378         nav2_radial = ils.get_locheading();
 379         while ( nav2_radial <   0.0 ) { nav2_radial += 360.0; }
 380         while ( nav2_radial > 360.0 ) { nav2_radial -= 360.0; }
 381         nav2_x = ils.get_x();
 382         nav2_y = ils.get_y();
 383         nav2_z = ils.get_z();
 384         nav2_gs_x = ils.get_gs_x();
 385         nav2_gs_y = ils.get_gs_y();
 386         nav2_gs_z = ils.get_gs_z();
 387         nav2_dme_x = ils.get_dme_x();
 388         nav2_dme_y = ils.get_dme_y();
 389         nav2_dme_z = ils.get_dme_z();
 390         // cout << "Found an ils station in range" << endl;
 391         // cout << " id = " << ils.get_locident() << endl;
 392     } else if ( current_navlist->query( lon, lat, elev, nav2_freq, &nav ) ) {
 393         nav2_valid = true;
 394         nav2_loc = false;
 395         nav2_has_dme = nav.get_has_dme();
 396         nav2_has_dme = false;
 397         nav2_loclon = nav.get_lon();
 398         nav2_loclat = nav.get_lat();
 399         nav2_elev = nav.get_elev();
 400         nav2_magvar = nav.get_magvar();
 401         nav2_effective_range = kludgeRange(nav2_elev, elev, nav.get_range());
 402         nav2_target_gs = 0.0;
 403         nav2_radial = nav2_sel_radial;
 404         nav2_x = nav2_dme_x = nav.get_x();
 405         nav2_y = nav2_dme_y = nav.get_y();
 406         nav2_z = nav2_dme_z = nav.get_z();
 407         // cout << "Found a vor station in range" << endl;
 408         // cout << " id = " << nav.get_ident() << endl;
 409     } else {
 410         nav2_valid = false;
 411         nav2_radial = 0;
 412         nav2_dme_dist = 0;
 413         // cout << "not picking up vor2. :-(" << endl;
 414     }
 415
 416     // adf
 417     if ( current_navlist->query( lon, lat, elev, adf_freq, &nav ) ) {
 418         adf_valid = true;
 419
 420         adf_lon = nav.get_lon();
 421         adf_lat = nav.get_lat();
 422         adf_elev = nav.get_elev();
 423         adf_effective_range = kludgeRange(adf_elev, elev, nav.get_range());
 424         adf_x = nav.get_x();
 425         adf_y = nav.get_y();
 426         adf_z = nav.get_z();
 427         // cout << "Found an adf station in range" << endl;
 428         // cout << " id = " << nav.get_ident() << endl;
 429     } else {
 430         adf_valid = false;
 431         // cout << "not picking up adf. :-(" << endl;
 432     }
 433 }
 434
 435
 436 /**
 437  * Return true if the NAV1 TO flag should be active.
 438  */
 439 bool
 440 FGRadioStack::get_nav1_to_flag () const
 441 {
 442   if (nav1_inrange) {
 443     double offset = fabs(nav1_heading - nav1_radial);
 444     if (nav1_loc)
 445       return (offset <= 8.0 || offset >= 352.0);
 446     else
 447       return (offset <= 20.0 || offset >= 340.0);
 448   } else {
 449     return false;
 450   }
 451 }
 452
 453
 454 /**
 455  * Return true if the NAV1 FROM flag should be active.
 456  */
 457 bool
 458 FGRadioStack::get_nav1_from_flag () const
 459 {
 460   if (nav1_inrange) {
 461     double offset = fabs(nav1_heading - nav1_radial);
 462     if (nav1_loc)
 463       return (offset >= 172.0 && offset <= 188.0);
 464     else
 465       return (offset >= 160.0 && offset <= 200.0);
 466   } else {
 467     return false;
 468   }
 469 }
 470
 471
 472 /**
 473  * Return true if the NAV2 TO flag should be active.
 474  */
 475 bool
 476 FGRadioStack::get_nav2_to_flag () const
 477 {
 478   if (nav2_inrange) {
 479     double offset = fabs(nav2_heading - nav2_radial);
 480     if (nav2_loc)
 481       return (offset <= 8.0 || offset >= 352.0);
 482     else
 483       return (offset <= 20.0 || offset >= 340.0);
 484   } else {
 485     return false;
 486   }
 487 }
 488
 489
 490 /**
 491  * Return true if the NAV2 FROM flag should be active.
 492  */
 493 bool
 494 FGRadioStack::get_nav2_from_flag () const
 495 {
 496   if (nav2_inrange) {
 497     double offset = fabs(nav2_heading - nav2_radial);
 498     if (nav2_loc)
 499       return (offset >= 172.0 && offset <= 188.0);
 500     else
 501       return (offset >= 160.0 && offset <= 200.0);
 502   } else {
 503     return false;
 504   }
 505 }
 506
 507