src/Cockpit/navcom.cxx

   1 // navcom.cxx -- class to manage a navcom instance
   2 //
   3 // Written by Curtis Olson, started April 2000.
   4 //
   5 // Copyright (C) 2000 - 2002  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 #ifdef HAVE_CONFIG_H
  25 #  include <config.h>
  26 #endif
  27
  28 #include <stdio.h>      // snprintf
  29
  30 #include <simgear/compiler.h>
  31 #include <simgear/math/sg_random.h>
  32
  33 #include <Aircraft/aircraft.hxx>
  34 #include <Navaids/ilslist.hxx>
  35 #include <Navaids/navlist.hxx>
  36 #include <Time/FGEventMgr.hxx>
  37
  38 #include "navcom.hxx"
  39
  40 #include <string>
  41 SG_USING_STD(string);
  42
  43
  44 // Constructor
  45 FGNavCom::FGNavCom() :
  46     lon_node(fgGetNode("/position/longitude-deg", true)),
  47     lat_node(fgGetNode("/position/latitude-deg", true)),
  48     alt_node(fgGetNode("/position/altitude-ft", true)),
  49     last_nav_id(""),
  50     last_nav_vor(false),
  51     nav_play_count(0),
  52     nav_last_time(0),
  53     need_update(true),
  54     power_btn(true),
  55     audio_btn(true),
  56     comm_freq(0.0),
  57     comm_alt_freq(0.0),
  58     comm_vol_btn(0.0),
  59     nav_freq(0.0),
  60     nav_alt_freq(0.0),
  61     nav_radial(0.0),
  62     nav_vol_btn(0.0),
  63     nav_ident_btn(true)
  64 {
  65     SGPath path( globals->get_fg_root() );
  66     SGPath term = path;
  67     term.append( "Navaids/range.term" );
  68     SGPath low = path;
  69     low.append( "Navaids/range.low" );
  70     SGPath high = path;
  71     high.append( "Navaids/range.high" );
  72
  73     term_tbl = new SGInterpTable( term.str() );
  74     low_tbl = new SGInterpTable( low.str() );
  75     high_tbl = new SGInterpTable( high.str() );
  76
  77 }
  78
  79
  80 // Destructor
  81 FGNavCom::~FGNavCom()
  82 {
  83     delete term_tbl;
  84     delete low_tbl;
  85     delete high_tbl;
  86 }
  87
  88
  89 void
  90 FGNavCom::init ()
  91 {
  92     morse.init();
  93
  94     // We assume that index is valid now (it must be set before init()
  95     // is called.)
  96     char propname[256];
  97
  98     sprintf( propname, "/systems/electrical/outputs/navcom[%d]", index );
  99     // default to true in case no electrical system defined.
 100     fgSetDouble( propname, 60.0 );
 101     bus_power = fgGetNode( propname, true );
 102
 103     sprintf( propname, "/instrumentation/comm[%d]/servicable", index );
 104     com_servicable = fgGetNode( propname, true );
 105     com_servicable->setBoolValue( true );
 106
 107     sprintf( propname, "/instrumentation/nav[%d]/servicable", index );
 108     nav_servicable = fgGetNode( propname, true );
 109     nav_servicable->setBoolValue( true );
 110
 111     sprintf( propname, "/instrumentation/vor[%d]/cdi/servicable", index );
 112     cdi_servicable = fgGetNode( propname, true );
 113     cdi_servicable->setBoolValue( true );
 114
 115     sprintf( propname, "/instrumentation/vor[%d]/gs/servicable", index );
 116     gs_servicable = fgGetNode( propname, true );
 117     gs_servicable->setBoolValue( true );
 118
 119     sprintf( propname, "/instrumentation/vor[%d]/to-from/servicable", index );
 120     tofrom_servicable = fgGetNode( propname, true );
 121     tofrom_servicable->setBoolValue( true );
 122 }
 123
 124 void
 125 FGNavCom::bind ()
 126 {
 127     char propname[256];
 128
 129                                 // User inputs
 130     sprintf( propname, "/radios/comm[%d]/inputs/power-btn", index );
 131     fgTie( propname, this,
 132            &FGNavCom::get_power_btn, &FGNavCom::set_power_btn );
 133     fgSetArchivable( propname );
 134
 135     sprintf( propname, "/radios/comm[%d]/frequencies/selected-mhz", index );
 136     fgTie( propname, this, &FGNavCom::get_comm_freq, &FGNavCom::set_comm_freq );
 137     fgSetArchivable( propname );
 138
 139     sprintf( propname, "/radios/comm[%d]/frequencies/standby-mhz", index );
 140     fgTie( propname, this,
 141            &FGNavCom::get_comm_alt_freq, &FGNavCom::set_comm_alt_freq );
 142     fgSetArchivable( propname );
 143
 144     sprintf( propname, "/radios/comm[%d]/volume", index );
 145     fgTie( propname, this,
 146            &FGNavCom::get_comm_vol_btn, &FGNavCom::set_comm_vol_btn );
 147     fgSetArchivable( propname );
 148
 149     sprintf( propname, "/radios/nav[%d]/frequencies/selected-mhz", index );
 150     fgTie( propname, this,
 151           &FGNavCom::get_nav_freq, &FGNavCom::set_nav_freq );
 152     fgSetArchivable( propname );
 153
 154     sprintf( propname, "/radios/nav[%d]/frequencies/standby-mhz", index );
 155     fgTie( propname , this,
 156            &FGNavCom::get_nav_alt_freq, &FGNavCom::set_nav_alt_freq);
 157     fgSetArchivable( propname );
 158
 159     sprintf( propname, "/radios/nav[%d]/radials/selected-deg", index );
 160     fgTie( propname, this,
 161            &FGNavCom::get_nav_sel_radial, &FGNavCom::set_nav_sel_radial );
 162     fgSetArchivable( propname );
 163
 164     sprintf( propname, "/radios/nav[%d]/volume", index );
 165     fgTie( propname, this,
 166            &FGNavCom::get_nav_vol_btn, &FGNavCom::set_nav_vol_btn );
 167     fgSetArchivable( propname );
 168
 169     sprintf( propname, "/radios/nav[%d]/ident", index );
 170     fgTie( propname, this,
 171            &FGNavCom::get_nav_ident_btn, &FGNavCom::set_nav_ident_btn );
 172     fgSetArchivable( propname );
 173
 174                                 // Radio outputs
 175     sprintf( propname, "/radios/nav[%d]/audio-btn", index );
 176     fgTie( propname, this,
 177            &FGNavCom::get_audio_btn, &FGNavCom::set_audio_btn );
 178     fgSetArchivable( propname );
 179
 180     sprintf( propname, "/radios/nav[%d]/radials/actual-deg", index );
 181     fgTie( propname,  this, &FGNavCom::get_nav_radial );
 182
 183     sprintf( propname, "/radios/nav[%d]/to-flag", index );
 184     fgTie( propname, this, &FGNavCom::get_nav_to_flag );
 185
 186     sprintf( propname, "/radios/nav[%d]/from-flag", index );
 187     fgTie( propname, this, &FGNavCom::get_nav_from_flag );
 188
 189     sprintf( propname, "/radios/nav[%d]/in-range", index );
 190     fgTie( propname, this, &FGNavCom::get_nav_inrange );
 191
 192     sprintf( propname, "/radios/nav[%d]/heading-needle-deflection", index );
 193     fgTie( propname, this, &FGNavCom::get_nav_heading_needle_deflection );
 194
 195     sprintf( propname, "/radios/nav[%d]/gs-needle-deflection", index );
 196     fgTie( propname, this, &FGNavCom::get_nav_gs_needle_deflection );
 197
 198     // end of binding
 199 }
 200
 201
 202 void
 203 FGNavCom::unbind ()
 204 {
 205     char propname[256];
 206
 207     sprintf( propname, "/radios/comm[%d]/inputs/power-btn", index );
 208     fgUntie( propname );
 209     sprintf( propname, "/radios/comm[%d]/frequencies/selected-mhz", index );
 210     fgUntie( propname );
 211     sprintf( propname, "/radios/comm[%d]/frequencies/standby-mhz", index );
 212     fgUntie( propname );
 213
 214     sprintf( propname, "/radios/nav[%d]/frequencies/selected-mhz", index );
 215     fgUntie( propname );
 216     sprintf( propname, "/radios/nav[%d]/frequencies/standby-mhz", index );
 217     fgUntie( propname );
 218     sprintf( propname, "/radios/nav[%d]/radials/actual-deg", index );
 219     fgUntie( propname );
 220     sprintf( propname, "/radios/nav[%d]/radials/selected-deg", index );
 221     fgUntie( propname );
 222     sprintf( propname, "/radios/nav[%d]/ident", index );
 223     fgUntie( propname );
 224     sprintf( propname, "/radios/nav[%d]/to-flag", index );
 225     fgUntie( propname );
 226     sprintf( propname, "/radios/nav[%d]/from-flag", index );
 227     fgUntie( propname );
 228     sprintf( propname, "/radios/nav[%d]/in-range", index );
 229     fgUntie( propname );
 230     sprintf( propname, "/radios/nav[%d]/heading-needle-deflection", index );
 231     fgUntie( propname );
 232     sprintf( propname, "/radios/nav[%d]/gs-needle-deflection", index );
 233     fgUntie( propname );
 234 }
 235
 236
 237 // model standard VOR/DME/TACAN service volumes as per AIM 1-1-8
 238 double FGNavCom::adjustNavRange( double stationElev, double aircraftElev,
 239                                  double nominalRange )
 240 {
 241     // extend out actual usable range to be 1.3x the published safe range
 242     const double usability_factor = 1.3;
 243
 244     // assumptions we model the standard service volume, plus
 245     // ... rather than specifying a cylinder, we model a cone that
 246     // contains the cylinder.  Then we put an upside down cone on top
 247     // to model diminishing returns at too-high altitudes.
 248
 249     // altitude difference
 250     double alt = ( aircraftElev * SG_METER_TO_FEET - stationElev );
 251     // cout << "aircraft elev = " << aircraftElev * SG_METER_TO_FEET
 252     //      << " station elev = " << stationElev << endl;
 253
 254     if ( nominalRange < 25.0 + SG_EPSILON ) {
 255         // Standard Terminal Service Volume
 256         return term_tbl->interpolate( alt ) * usability_factor;
 257     } else if ( nominalRange < 50.0 + SG_EPSILON ) {
 258         // Standard Low Altitude Service Volume
 259         // table is based on range of 40, scale to actual range
 260         return low_tbl->interpolate( alt ) * nominalRange / 40.0
 261             * usability_factor;
 262     } else {
 263         // Standard High Altitude Service Volume
 264         // table is based on range of 130, scale to actual range
 265         return high_tbl->interpolate( alt ) * nominalRange / 130.0
 266             * usability_factor;
 267     }
 268 }
 269
 270
 271 // model standard ILS service volumes as per AIM 1-1-9
 272 double FGNavCom::adjustILSRange( double stationElev, double aircraftElev,
 273                                      double offsetDegrees, double distance )
 274 {
 275     // assumptions we model the standard service volume, plus
 276
 277     // altitude difference
 278     // double alt = ( aircraftElev * SG_METER_TO_FEET - stationElev );
 279     double offset = fabs( offsetDegrees );
 280
 281     if ( offset < 10 ) {
 282         return FG_ILS_DEFAULT_RANGE;
 283     } else if ( offset < 35 ) {
 284         return 10 + (35 - offset) * (FG_ILS_DEFAULT_RANGE - 10) / 25;
 285     } else if ( offset < 45 ) {
 286         return (45 - offset);
 287     } else if ( offset > 170 ) {
 288         return FG_ILS_DEFAULT_RANGE;
 289     } else if ( offset > 145 ) {
 290         return 10 + (offset - 145) * (FG_ILS_DEFAULT_RANGE - 10) / 25;
 291     } else if ( offset > 135 ) {
 292         return (offset - 135);
 293     } else {
 294         return 0;
 295     }
 296 }
 297
 298
 299 // Update the various nav values based on position and valid tuned in navs
 300 void
 301 FGNavCom::update(double dt)
 302 {
 303     double lon = lon_node->getDoubleValue() * SGD_DEGREES_TO_RADIANS;
 304     double lat = lat_node->getDoubleValue() * SGD_DEGREES_TO_RADIANS;
 305     double elev = alt_node->getDoubleValue() * SG_FEET_TO_METER;
 306
 307     need_update = false;
 308
 309     Point3D aircraft = sgGeodToCart( Point3D( lon, lat, elev ) );
 310     Point3D station;
 311     double az1, az2, s;
 312
 313     ////////////////////////////////////////////////////////////////////////
 314     // Nav.
 315     ////////////////////////////////////////////////////////////////////////
 316
 317     if ( nav_valid && power_btn && (bus_power->getDoubleValue() > 1.0)
 318          && nav_servicable->getBoolValue() )
 319     {
 320         station = Point3D( nav_x, nav_y, nav_z );
 321         nav_loc_dist = aircraft.distance3D( station );
 322
 323         if ( nav_has_gs ) {
 324             station = Point3D( nav_gs_x, nav_gs_y, nav_gs_z );
 325             nav_gs_dist = aircraft.distance3D( station );
 326         } else {
 327             nav_gs_dist = 0.0;
 328         }
 329
 330         // wgs84 heading
 331         geo_inverse_wgs_84( elev, lat * SGD_RADIANS_TO_DEGREES, lon * SGD_RADIANS_TO_DEGREES,
 332                             nav_loclat, nav_loclon,
 333                             &az1, &az2, &s );
 334         // cout << "az1 = " << az1 << " magvar = " << nav_magvar << endl;
 335         nav_heading = az1 - nav_magvar;
 336         // cout << " heading = " << nav_heading
 337         //      << " dist = " << nav_dist << endl;
 338
 339         if ( nav_loc ) {
 340             double offset = nav_heading - nav_radial;
 341             while ( offset < -180.0 ) { offset += 360.0; }
 342             while ( offset > 180.0 ) { offset -= 360.0; }
 343             // cout << "ils offset = " << offset << endl;
 344             nav_effective_range = adjustILSRange(nav_elev, elev, offset,
 345                                                   nav_loc_dist * SG_METER_TO_NM );
 346         } else {
 347             nav_effective_range = adjustNavRange(nav_elev, elev, nav_range);
 348         }
 349         // cout << "nav range = " << nav_effective_range
 350         //      << " (" << nav_range << ")" << endl;
 351
 352         if ( nav_loc_dist < nav_effective_range * SG_NM_TO_METER ) {
 353             nav_inrange = true;
 354         } else if ( nav_loc_dist < 2 * nav_effective_range * SG_NM_TO_METER ) {
 355             nav_inrange = sg_random() <
 356                 ( 2 * nav_effective_range * SG_NM_TO_METER - nav_loc_dist ) /
 357                 (nav_effective_range * SG_NM_TO_METER);
 358         } else {
 359             nav_inrange = false;
 360         }
 361
 362         if ( !nav_loc ) {
 363             nav_radial = nav_sel_radial;
 364         }
 365     } else {
 366         nav_inrange = false;
 367         // cout << "not picking up vor. :-(" << endl;
 368     }
 369
 370     if ( nav_valid && nav_inrange && nav_servicable->getBoolValue() ) {
 371         // play station ident via audio system if on + ident,
 372         // otherwise turn it off
 373         if ( power_btn && (bus_power->getDoubleValue() > 1.0)
 374              && nav_ident_btn && audio_btn )
 375         {
 376             FGSimpleSound *sound;
 377             sound = globals->get_soundmgr()->find( nav_fx_name );
 378             if ( sound != NULL ) {
 379                 sound->set_volume( nav_vol_btn );
 380             } else {
 381                 SG_LOG( SG_COCKPIT, SG_ALERT,
 382                         "Can't find nav-vor-ident sound" );
 383             }
 384             sound = globals->get_soundmgr()->find( dme_fx_name );
 385             if ( sound != NULL ) {
 386                 sound->set_volume( nav_vol_btn );
 387             } else {
 388                 SG_LOG( SG_COCKPIT, SG_ALERT,
 389                         "Can't find nav-dme-ident sound" );
 390             }
 391             // cout << "nav_last_time = " << nav_last_time << " ";
 392             // cout << "cur_time = "
 393             //      << globals->get_time_params()->get_cur_time();
 394             if ( nav_last_time <
 395                  globals->get_time_params()->get_cur_time() - 30 ) {
 396                 nav_last_time = globals->get_time_params()->get_cur_time();
 397                 nav_play_count = 0;
 398             }
 399             // cout << " nav_play_count = " << nav_play_count << endl;
 400             // cout << "playing = "
 401             //      << globals->get_soundmgr()->is_playing(nav_fx_name)
 402             //      << endl;
 403             if ( nav_play_count < 4 ) {
 404                 // play VOR ident
 405                 if ( !globals->get_soundmgr()->is_playing(nav_fx_name) ) {
 406                     globals->get_soundmgr()->play_once( nav_fx_name );
 407                     ++nav_play_count;
 408                 }
 409             } else if ( nav_play_count < 5 && nav_has_dme ) {
 410                 // play DME ident
 411                 if ( !globals->get_soundmgr()->is_playing(nav_fx_name) &&
 412                      !globals->get_soundmgr()->is_playing(dme_fx_name) ) {
 413                     globals->get_soundmgr()->play_once( dme_fx_name );
 414                     ++nav_play_count;
 415                 }
 416             }
 417         } else {
 418             globals->get_soundmgr()->stop( nav_fx_name );
 419             globals->get_soundmgr()->stop( dme_fx_name );
 420         }
 421     }
 422 }
 423
 424
 425 // Update current nav/adf radio stations based on current postition
 426 void FGNavCom::search()
 427 {
 428     double lon = lon_node->getDoubleValue() * SGD_DEGREES_TO_RADIANS;
 429     double lat = lat_node->getDoubleValue() * SGD_DEGREES_TO_RADIANS;
 430     double elev = alt_node->getDoubleValue() * SG_FEET_TO_METER;
 431
 432     FGILS ils;
 433     FGNav nav;
 434
 435     ////////////////////////////////////////////////////////////////////////
 436     // Nav.
 437     ////////////////////////////////////////////////////////////////////////
 438
 439     if ( current_ilslist->query( lon, lat, elev, nav_freq, &ils ) ) {
 440         nav_id = ils.get_locident();
 441         nav_valid = true;
 442         if ( last_nav_id != nav_id || last_nav_vor ) {
 443             nav_trans_ident = ils.get_trans_ident();
 444             last_nav_id = nav_id;
 445             last_nav_vor = false;
 446             nav_loc = true;
 447             nav_has_dme = ils.get_has_dme();
 448             nav_has_gs = ils.get_has_gs();
 449
 450             nav_loclon = ils.get_loclon();
 451             nav_loclat = ils.get_loclat();
 452             nav_gslon = ils.get_gslon();
 453             nav_gslat = ils.get_gslat();
 454             nav_elev = ils.get_gselev();
 455             nav_magvar = 0;
 456             nav_range = FG_ILS_DEFAULT_RANGE;
 457             nav_effective_range = nav_range;
 458             nav_target_gs = ils.get_gsangle();
 459             nav_radial = ils.get_locheading();
 460             while ( nav_radial <   0.0 ) { nav_radial += 360.0; }
 461             while ( nav_radial > 360.0 ) { nav_radial -= 360.0; }
 462             nav_x = ils.get_x();
 463             nav_y = ils.get_y();
 464             nav_z = ils.get_z();
 465             nav_gs_x = ils.get_gs_x();
 466             nav_gs_y = ils.get_gs_y();
 467             nav_gs_z = ils.get_gs_z();
 468
 469             if ( globals->get_soundmgr()->exists( nav_fx_name ) ) {
 470                 globals->get_soundmgr()->remove( nav_fx_name );
 471             }
 472             FGSimpleSound *sound;
 473             sound = morse.make_ident( nav_trans_ident, LO_FREQUENCY );
 474             sound->set_volume( 0.3 );
 475             globals->get_soundmgr()->add( sound, nav_fx_name );
 476
 477             if ( globals->get_soundmgr()->exists( dme_fx_name ) ) {
 478                 globals->get_soundmgr()->remove( dme_fx_name );
 479             }
 480             sound = morse.make_ident( nav_trans_ident, HI_FREQUENCY );
 481             sound->set_volume( 0.3 );
 482             globals->get_soundmgr()->add( sound, dme_fx_name );
 483
 484             int offset = (int)(sg_random() * 30.0);
 485             nav_play_count = offset / 4;
 486             nav_last_time = globals->get_time_params()->get_cur_time() -
 487                 offset;
 488             // cout << "offset = " << offset << " play_count = "
 489             //      << nav_play_count
 490             //      << " nav_last_time = " << nav_last_time
 491             //      << " current time = "
 492             //      << globals->get_time_params()->get_cur_time() << endl;
 493
 494             // cout << "Found an ils station in range" << endl;
 495             // cout << " id = " << ils.get_locident() << endl;
 496         }
 497     } else if ( current_navlist->query( lon, lat, elev, nav_freq, &nav ) ) {
 498         nav_id = nav.get_ident();
 499         nav_valid = true;
 500         if ( last_nav_id != nav_id || !last_nav_vor ) {
 501             last_nav_id = nav_id;
 502             last_nav_vor = true;
 503             nav_trans_ident = nav.get_trans_ident();
 504             nav_loc = false;
 505             nav_has_dme = nav.get_has_dme();
 506             nav_has_gs = false;
 507             nav_loclon = nav.get_lon();
 508             nav_loclat = nav.get_lat();
 509             nav_elev = nav.get_elev();
 510             nav_magvar = nav.get_magvar();
 511             nav_range = nav.get_range();
 512             nav_effective_range = adjustNavRange(nav_elev, elev, nav_range);
 513             nav_target_gs = 0.0;
 514             nav_radial = nav_sel_radial;
 515             nav_x = nav.get_x();
 516             nav_y = nav.get_y();
 517             nav_z = nav.get_z();
 518
 519             if ( globals->get_soundmgr()->exists( nav_fx_name ) ) {
 520                 globals->get_soundmgr()->remove( nav_fx_name );
 521             }
 522             FGSimpleSound *sound;
 523             sound = morse.make_ident( nav_trans_ident, LO_FREQUENCY );
 524             sound->set_volume( 0.3 );
 525             if ( globals->get_soundmgr()->add( sound, nav_fx_name ) ) {
 526                 // cout << "Added nav-vor-ident sound" << endl;
 527             } else {
 528                 cout << "Failed to add v1-vor-ident sound" << endl;
 529             }
 530
 531             if ( globals->get_soundmgr()->exists( dme_fx_name ) ) {
 532                 globals->get_soundmgr()->remove( dme_fx_name );
 533             }
 534             sound = morse.make_ident( nav_trans_ident, HI_FREQUENCY );
 535             sound->set_volume( 0.3 );
 536             globals->get_soundmgr()->add( sound, dme_fx_name );
 537
 538             int offset = (int)(sg_random() * 30.0);
 539             nav_play_count = offset / 4;
 540             nav_last_time = globals->get_time_params()->get_cur_time() -
 541                 offset;
 542             // cout << "offset = " << offset << " play_count = "
 543             //      << nav_play_count << " nav_last_time = "
 544             //      << nav_last_time << " current time = "
 545             //      << globals->get_time_params()->get_cur_time() << endl;
 546
 547             // cout << "Found a vor station in range" << endl;
 548             // cout << " id = " << nav.get_ident() << endl;
 549         }
 550     } else {
 551         nav_valid = false;
 552         nav_id = "";
 553         nav_radial = 0;
 554         nav_trans_ident = "";
 555         last_nav_id = "";
 556         if ( ! globals->get_soundmgr()->remove( nav_fx_name ) ) {
 557             cout << "Failed to remove nav-vor-ident sound" << endl;
 558         }
 559         globals->get_soundmgr()->remove( dme_fx_name );
 560         // cout << "not picking up vor1. :-(" << endl;
 561     }
 562 }
 563
 564
 565 // return the amount of heading needle deflection, returns a value
 566 // clamped to the range of ( -10 , 10 )
 567 double FGNavCom::get_nav_heading_needle_deflection() const {
 568     double r;
 569
 570     if ( nav_inrange
 571          && nav_servicable->getBoolValue() && cdi_servicable->getBoolValue() )
 572     {
 573         r = nav_heading - nav_radial;
 574         // cout << "Radial = " << nav_radial
 575         //      << "  Bearing = " << nav_heading << endl;
 576
 577         while ( r >  180.0 ) { r -= 360.0;}
 578         while ( r < -180.0 ) { r += 360.0;}
 579         if ( fabs(r) > 90.0 )
 580             r = ( r<0.0 ? -r-180.0 : -r+180.0 );
 581
 582         // According to Robin Peel, the ILS is 4x more sensitive than a vor
 583         if ( nav_loc ) { r *= 4.0; }
 584         if ( r < -10.0 ) { r = -10.0; }
 585         if ( r >  10.0 ) { r = 10.0; }
 586     } else {
 587         r = 0.0;
 588     }
 589
 590     return r;
 591 }
 592
 593
 594 // return the amount of glide slope needle deflection (.i.e. the
 595 // number of degrees we are off the glide slope * 5.0
 596 double FGNavCom::get_nav_gs_needle_deflection() const {
 597     if ( nav_inrange && nav_has_gs
 598          && nav_servicable->getBoolValue() && gs_servicable->getBoolValue() )
 599     {
 600         double x = nav_gs_dist;
 601         double y = (fgGetDouble("/position/altitude-ft") - nav_elev)
 602             * SG_FEET_TO_METER;
 603         double angle = atan2( y, x ) * SGD_RADIANS_TO_DEGREES;
 604         return (nav_target_gs - angle) * 5.0;
 605     } else {
 606         return 0.0;
 607     }
 608 }
 609
 610
 611 /**
 612  * Return true if the NAV TO flag should be active.
 613  */
 614 bool
 615 FGNavCom::get_nav_to_flag () const
 616 {
 617     if ( nav_inrange
 618          && nav_servicable->getBoolValue()
 619          && tofrom_servicable->getBoolValue() )
 620     {
 621         double offset = fabs(nav_heading - nav_radial);
 622         if (nav_loc) {
 623             return true;
 624         } else {
 625             return (offset <= 90.0 || offset >= 270.0);
 626         }
 627     } else {
 628         return false;
 629     }
 630 }
 631
 632
 633 /**
 634  * Return true if the NAV FROM flag should be active.
 635  */
 636 bool
 637 FGNavCom::get_nav_from_flag () const
 638 {
 639     if ( nav_inrange
 640          && nav_servicable->getBoolValue()
 641          && tofrom_servicable->getBoolValue() ) {
 642         double offset = fabs(nav_heading - nav_radial);
 643         if (nav_loc) {
 644             return false;
 645         } else {
 646           return (offset > 90.0 && offset < 270.0);
 647         }
 648     } else {
 649         return false;
 650     }
 651 }