1 // radiostack.cxx -- class to manage an instance of the radio stack
3 // Written by Curtis Olson, started April 2000.
5 // Copyright (C) 2000 Curtis L. Olson - curt@flightgear.org
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.
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.
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.
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>
30 #include "radiostack.hxx"
33 * Boy, this is ugly! Make the VOR range vary by altitude difference.
36 kludgeRange (double stationElev, double aircraftElev, double nominalRange)
38 // Assume that the nominal range (usually
39 // 50nm) applies at a 5,000 ft difference.
41 double factor = ((aircraftElev*METER_TO_FEET) - stationElev) / 1000.0;
42 double range = fabs(nominalRange * factor);
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;
57 FGRadioStack *current_radiostack;
60 // periodic radio station search wrapper
61 static void fgRadioSearch( void ) {
62 current_radiostack->search();
66 FGRadioStack::FGRadioStack() {
72 longitudeVal = fgGetValue("/position/longitude");
73 latitudeVal = fgGetValue("/position/latitude");
74 altitudeVal = fgGetValue("/position/altitude");
79 FGRadioStack::~FGRadioStack()
81 unbind(); // FIXME: should be called externally
91 // Search radio database once per second
92 global_events.Register( "fgRadioSearch()", fgRadioSearch,
93 fgEVENT::FG_EVENT_READY, 1000);
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);
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 fgTie("/radios/nav1/heading-needle-deflection", this,
117 &FGRadioStack::get_nav1_heading_needle_deflection);
118 fgTie("/radios/nav1/gs-needle-deflection", this,
119 &FGRadioStack::get_nav1_gs_needle_deflection);
122 fgTie("/radios/nav2/frequencies/selected", this,
123 &FGRadioStack::get_nav2_freq, &FGRadioStack::set_nav2_freq);
124 fgTie("/radios/nav2/frequencies/standby", this,
125 &FGRadioStack::get_nav2_alt_freq, &FGRadioStack::set_nav2_alt_freq);
126 fgTie("/radios/nav2/radials/selected", this,
127 &FGRadioStack::get_nav2_sel_radial,
128 &FGRadioStack::set_nav2_sel_radial);
131 fgTie("/radios/nav2/radials/actual", this, &FGRadioStack::get_nav2_radial);
132 fgTie("/radios/nav2/to-flag", this, &FGRadioStack::get_nav2_to_flag);
133 fgTie("/radios/nav2/from-flag", this, &FGRadioStack::get_nav2_from_flag);
134 fgTie("/radios/nav2/in-range", this, &FGRadioStack::get_nav2_inrange);
135 fgTie("/radios/nav2/dme/distance", this, &FGRadioStack::get_nav2_dme_dist);
136 fgTie("/radios/nav2/dme/in-range", this,
137 &FGRadioStack::get_nav2_dme_inrange);
138 fgTie("/radios/nav2/heading-needle-deflection", this,
139 &FGRadioStack::get_nav2_heading_needle_deflection);
140 fgTie("/radios/nav2/gs-needle-deflection", this,
141 &FGRadioStack::get_nav2_gs_needle_deflection);
144 fgTie("/radios/adf/frequencies/selected", this,
145 &FGRadioStack::get_adf_freq, &FGRadioStack::set_adf_freq);
146 fgTie("/radios/adf/frequencies/standby", this,
147 &FGRadioStack::get_adf_alt_freq, &FGRadioStack::set_adf_alt_freq);
148 fgTie("/radios/adf/rotation", this,
149 &FGRadioStack::get_adf_rotation, &FGRadioStack::set_adf_rotation);
153 FGRadioStack::unbind ()
155 fgUntie("/radios/nav1/frequencies/selected");
156 fgUntie("/radios/nav1/frequencies/standby");
157 fgUntie("/radios/nav1/radials/actual");
158 fgUntie("/radios/nav1/radials/selected");
159 fgUntie("/radios/nav1/to-flag");
160 fgUntie("/radios/nav1/from-flag");
161 fgUntie("/radios/nav1/in-range");
162 fgUntie("/radios/nav1/dme/distance");
163 fgUntie("/radios/nav1/dme/in-range");
164 fgUntie("/radios/nav1/heading-needle-deflection");
165 fgUntie("/radios/nav1/gs-needle-deflection");
167 fgUntie("/radios/nav2/frequencies/selected");
168 fgUntie("/radios/nav2/frequencies/standby");
169 fgUntie("/radios/nav2/radials/actual");
170 fgUntie("/radios/nav2/radials/selected");
171 fgUntie("/radios/nav2/to-flag");
172 fgUntie("/radios/nav2/from-flag");
173 fgUntie("/radios/nav2/in-range");
174 fgUntie("/radios/nav2/dme/distance");
175 fgUntie("/radios/nav2/dme/in-range");
176 fgUntie("/radios/nav2/heading-needle-deflection");
177 fgUntie("/radios/nav2/gs-needle-deflection");
179 fgUntie("/radios/adf/frequencies/selected");
180 fgUntie("/radios/adf/frequencies/standby");
181 fgUntie("/radios/adf/rotation");
184 // Search the database for the current frequencies given current location
186 FGRadioStack::update()
188 double lon = longitudeVal->getDoubleValue() * DEG_TO_RAD;
189 double lat = latitudeVal->getDoubleValue() * DEG_TO_RAD;
190 double elev = altitudeVal->getDoubleValue() * FEET_TO_METER;
194 Point3D aircraft = sgGeodToCart( Point3D( lon, lat, elev ) );
199 station = Point3D( nav1_x, nav1_y, nav1_z );
200 nav1_loc_dist = aircraft.distance3D( station );
202 if ( nav1_has_dme ) {
203 // staightline distance
204 station = Point3D( nav1_dme_x, nav1_dme_y, nav1_dme_z );
205 nav1_dme_dist = aircraft.distance3D( station );
211 station = Point3D( nav1_gs_x, nav1_gs_y, nav1_gs_z );
212 nav1_gs_dist = aircraft.distance3D( station );
217 if ( nav1_loc_dist < nav1_effective_range * NM_TO_METER ) {
221 geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
222 nav1_loclat, nav1_loclon,
224 // cout << "az1 = " << az1 << " magvar = " << nav1_magvar << endl;
225 nav1_heading = az1 - nav1_magvar;
226 // Alex: nav1_heading = - (az1 - FGBFI::getMagVar() / RAD_TO_DEG);
228 // cout << " heading = " << nav1_heading
229 // << " dist = " << nav1_dist << endl;
231 nav1_inrange = false;
236 nav1_radial = nav1_sel_radial;
239 nav1_inrange = false;
241 // cout << "not picking up vor. :-(" << endl;
245 station = Point3D( nav2_x, nav2_y, nav2_z );
246 nav2_loc_dist = aircraft.distance3D( station );
248 if ( nav2_has_dme ) {
249 // staightline distance
250 station = Point3D( nav2_dme_x, nav2_dme_y, nav2_dme_z );
251 nav2_dme_dist = aircraft.distance3D( station );
257 station = Point3D( nav2_gs_x, nav2_gs_y, nav2_gs_z );
258 nav2_gs_dist = aircraft.distance3D( station );
263 if ( nav2_loc_dist < nav2_effective_range * NM_TO_METER ) {
267 geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
268 nav2_loclat, nav2_loclon,
270 nav2_heading = az1 - nav2_magvar;
271 // Alex: nav2_heading = - (az1 - FGBFI::getMagVar() / RAD_TO_DEG);
273 // cout << " heading = " << nav2_heading
274 // << " dist = " << nav2_dist << endl;
276 nav2_inrange = false;
280 nav2_radial = nav2_sel_radial;
283 nav2_inrange = false;
284 // cout << "not picking up vor. :-(" << endl;
289 // staightline distance
290 station = Point3D( adf_x, adf_y, adf_z );
291 adf_dist = aircraft.distance3D( station );
293 if ( adf_dist < adf_effective_range * NM_TO_METER ) {
297 geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
302 // cout << " heading = " << nav2_heading
303 // << " dist = " << nav2_dist << endl;
313 // Update current nav/adf radio stations based on current postition
314 void FGRadioStack::search ()
316 double lon = longitudeVal->getDoubleValue() * DEG_TO_RAD;
317 double lat = latitudeVal->getDoubleValue() * DEG_TO_RAD;
318 double elev = altitudeVal->getDoubleValue() * FEET_TO_METER;
324 if ( current_ilslist->query( lon, lat, elev, nav1_freq, &ils ) ) {
327 nav1_has_dme = ils.get_has_dme();
328 nav1_has_gs = ils.get_has_gs();
330 nav1_loclon = ils.get_loclon();
331 nav1_loclat = ils.get_loclat();
332 nav1_gslon = ils.get_gslon();
333 nav1_gslat = ils.get_gslat();
334 nav1_dmelon = ils.get_dmelon();
335 nav1_dmelat = ils.get_dmelat();
336 nav1_elev = ils.get_gselev();
338 nav1_effective_range = FG_ILS_DEFAULT_RANGE;
339 nav1_target_gs = ils.get_gsangle();
340 nav1_radial = ils.get_locheading();
341 while ( nav1_radial < 0.0 ) { nav1_radial += 360.0; }
342 while ( nav1_radial > 360.0 ) { nav1_radial -= 360.0; }
343 nav1_x = ils.get_x();
344 nav1_y = ils.get_y();
345 nav1_z = ils.get_z();
346 nav1_gs_x = ils.get_gs_x();
347 nav1_gs_y = ils.get_gs_y();
348 nav1_gs_z = ils.get_gs_z();
349 nav1_dme_x = ils.get_dme_x();
350 nav1_dme_y = ils.get_dme_y();
351 nav1_dme_z = ils.get_dme_z();
352 // cout << "Found an ils station in range" << endl;
353 // cout << " id = " << ils.get_locident() << endl;
354 } else if ( current_navlist->query( lon, lat, elev, nav1_freq, &nav ) ) {
357 nav1_has_dme = nav.get_has_dme();
359 nav1_loclon = nav.get_lon();
360 nav1_loclat = nav.get_lat();
361 nav1_elev = nav.get_elev();
362 nav1_magvar = nav.get_magvar();
363 nav1_effective_range = kludgeRange(nav1_elev, elev, nav.get_range());
364 nav1_target_gs = 0.0;
365 nav1_radial = nav1_sel_radial;
366 nav1_x = nav1_dme_x = nav.get_x();
367 nav1_y = nav1_dme_y = nav.get_y();
368 nav1_z = nav1_dme_z = nav.get_z();
369 // cout << "Found a vor station in range" << endl;
370 // cout << " id = " << nav.get_ident() << endl;
375 // cout << "not picking up vor1. :-(" << endl;
378 if ( current_ilslist->query( lon, lat, elev, nav2_freq, &ils ) ) {
381 nav2_has_dme = ils.get_has_dme();
382 nav2_has_gs = ils.get_has_gs();
384 nav2_loclon = ils.get_loclon();
385 nav2_loclat = ils.get_loclat();
386 nav2_elev = ils.get_gselev();
388 nav2_effective_range = FG_ILS_DEFAULT_RANGE;
389 nav2_target_gs = ils.get_gsangle();
390 nav2_radial = ils.get_locheading();
391 while ( nav2_radial < 0.0 ) { nav2_radial += 360.0; }
392 while ( nav2_radial > 360.0 ) { nav2_radial -= 360.0; }
393 nav2_x = ils.get_x();
394 nav2_y = ils.get_y();
395 nav2_z = ils.get_z();
396 nav2_gs_x = ils.get_gs_x();
397 nav2_gs_y = ils.get_gs_y();
398 nav2_gs_z = ils.get_gs_z();
399 nav2_dme_x = ils.get_dme_x();
400 nav2_dme_y = ils.get_dme_y();
401 nav2_dme_z = ils.get_dme_z();
402 // cout << "Found an ils station in range" << endl;
403 // cout << " id = " << ils.get_locident() << endl;
404 } else if ( current_navlist->query( lon, lat, elev, nav2_freq, &nav ) ) {
407 nav2_has_dme = nav.get_has_dme();
408 nav2_has_dme = false;
409 nav2_loclon = nav.get_lon();
410 nav2_loclat = nav.get_lat();
411 nav2_elev = nav.get_elev();
412 nav2_magvar = nav.get_magvar();
413 nav2_effective_range = kludgeRange(nav2_elev, elev, nav.get_range());
414 nav2_target_gs = 0.0;
415 nav2_radial = nav2_sel_radial;
416 nav2_x = nav2_dme_x = nav.get_x();
417 nav2_y = nav2_dme_y = nav.get_y();
418 nav2_z = nav2_dme_z = nav.get_z();
419 // cout << "Found a vor station in range" << endl;
420 // cout << " id = " << nav.get_ident() << endl;
425 // cout << "not picking up vor2. :-(" << endl;
429 if ( current_navlist->query( lon, lat, elev, adf_freq, &nav ) ) {
432 adf_lon = nav.get_lon();
433 adf_lat = nav.get_lat();
434 adf_elev = nav.get_elev();
435 adf_effective_range = kludgeRange(adf_elev, elev, nav.get_range());
439 // cout << "Found an adf station in range" << endl;
440 // cout << " id = " << nav.get_ident() << endl;
443 // cout << "not picking up adf. :-(" << endl;
448 // return the amount of heading needle deflection, returns a value
449 // clamped to the range of ( -10 , 10 )
450 double FGRadioStack::get_nav1_heading_needle_deflection() const {
453 if ( nav1_inrange ) {
454 r = nav1_heading - nav1_radial;
455 // cout << "Radial = " << nav1_radial
456 // << " Bearing = " << nav1_heading << endl;
458 while ( r > 180.0 ) { r -= 360.0;}
459 while ( r < -180.0 ) { r += 360.0;}
460 if ( fabs(r) > 90.0 ) {
461 r = ( r<0.0 ? -r-180.0 : -r+180.0 );
467 // According to Robin Peel, the ILS is 4x more sensitive than a vor
468 if ( nav1_loc ) { r *= 4.0; }
469 if ( r < -10.0 ) { r = -10.0; }
470 if ( r > 10.0 ) { r = 10.0; }
478 // return the amount of heading needle deflection, returns a value
479 // clamped to the range of ( -10 , 10 )
480 double FGRadioStack::get_nav2_heading_needle_deflection() const {
483 if ( nav2_inrange ) {
484 r = nav2_heading - nav2_radial;
485 // cout << "Radial = " << nav2_radial
486 // << " Bearing = " << nav2_heading << endl;
488 while (r> 180.0) r-=360.0;
489 while (r<-180.0) r+=360.0;
490 if ( fabs(r) > 90.0 )
491 r = ( r<0.0 ? -r-180.0 : -r+180.0 );
492 // According to Robin Peel, the ILS is 4x more sensitive than a vor
493 if ( nav2_loc ) r *= 4.0;
494 if ( r < -10.0 ) r = -10.0;
495 if ( r > 10.0 ) r = 10.0;
503 // return the amount of glide slope needle deflection (.i.e. the
504 // number of degrees we are off the glide slope * 5.0
505 double FGRadioStack::get_nav1_gs_needle_deflection() const {
506 if ( nav1_inrange && nav1_has_gs ) {
507 double x = nav1_gs_dist;
508 double y = (FGBFI::getAltitude() - nav1_elev) * FEET_TO_METER;
509 double angle = atan2( y, x ) * RAD_TO_DEG;
510 return (nav1_target_gs - angle) * 5.0;
517 // return the amount of glide slope needle deflection (.i.e. the
518 // number of degrees we are off the glide slope * 5.0
519 double FGRadioStack::get_nav2_gs_needle_deflection() const {
520 if ( nav2_inrange && nav2_has_gs ) {
521 double x = nav2_gs_dist;
522 double y = (FGBFI::getAltitude() - nav2_elev) * FEET_TO_METER;
523 double angle = atan2( y, x ) * RAD_TO_DEG;
524 return (nav2_target_gs - angle) * 5.0;
532 * Return true if the NAV1 TO flag should be active.
535 FGRadioStack::get_nav1_to_flag () const
538 double offset = fabs(nav1_heading - nav1_radial);
542 return (offset <= 90.0 || offset >= 270.0);
550 * Return true if the NAV1 FROM flag should be active.
553 FGRadioStack::get_nav1_from_flag () const
556 double offset = fabs(nav1_heading - nav1_radial);
560 return (offset > 90.0 && offset < 270.0);
568 * Return true if the NAV2 TO flag should be active.
571 FGRadioStack::get_nav2_to_flag () const
574 double offset = fabs(nav2_heading - nav2_radial);
578 return (offset <= 90.0 || offset >= 270.0);
586 * Return true if the NAV2 FROM flag should be active.
589 FGRadioStack::get_nav2_from_flag () const
592 double offset = fabs(nav2_heading - nav2_radial);
596 return (offset > 90.0 && offset < 270.0);