src/Airports/airport.cxx

   1 //
   2 // simple.cxx -- a really simplistic class to manage airport ID,
   3 //               lat, lon of the center of one of it's runways, and
   4 //               elevation in feet.
   5 //
   6 // Written by Curtis Olson, started April 1998.
   7 // Updated by Durk Talsma, started December, 2004.
   8 //
   9 // Copyright (C) 1998  Curtis L. Olson  - http://www.flightgear.org/~curt
  10 //
  11 // This program is free software; you can redistribute it and/or
  12 // modify it under the terms of the GNU General Public License as
  13 // published by the Free Software Foundation; either version 2 of the
  14 // License, or (at your option) any later version.
  15 //
  16 // This program is distributed in the hope that it will be useful, but
  17 // WITHOUT ANY WARRANTY; without even the implied warranty of
  18 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19 // General Public License for more details.
  20 //
  21 // You should have received a copy of the GNU General Public License
  22 // along with this program; if not, write to the Free Software
  23 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  24 //
  25 // $Id$
  26
  27 #ifdef HAVE_CONFIG_H
  28 #  include <config.h>
  29 #endif
  30
  31 #include "airport.hxx"
  32
  33 #include <algorithm>
  34 #include <cassert>
  35 #include <boost/foreach.hpp>
  36
  37 #include <simgear/misc/sg_path.hxx>
  38 #include <simgear/props/props.hxx>
  39 #include <simgear/props/props_io.hxx>
  40 #include <simgear/debug/logstream.hxx>
  41 #include <simgear/sg_inlines.h>
  42 #include <simgear/structure/exception.hxx>
  43
  44 #include <Environment/environment_mgr.hxx>
  45 #include <Environment/environment.hxx>
  46 #include <Main/fg_props.hxx>
  47 #include <Airports/runways.hxx>
  48 #include <Airports/pavement.hxx>
  49 #include <Airports/dynamics.hxx>
  50 #include <Airports/xmlloader.hxx>
  51 #include <Navaids/procedure.hxx>
  52 #include <Navaids/waypoint.hxx>
  53 #include <ATC/CommStation.hxx>
  54 #include <Navaids/NavDataCache.hxx>
  55 #include <Navaids/navrecord.hxx>
  56
  57 using std::vector;
  58 using std::pair;
  59
  60 using namespace flightgear;
  61
  62 /***************************************************************************
  63  * FGAirport
  64  ***************************************************************************/
  65
  66 AirportCache FGAirport::airportCache;
  67
  68 FGAirport::FGAirport( PositionedID aGuid,
  69                       const std::string &id,
  70                       const SGGeod& location,
  71                       const std::string &name,
  72                       bool has_metar,
  73                       Type aType ):
  74     FGPositioned(aGuid, aType, id, location),
  75     _name(name),
  76     _has_metar(has_metar),
  77     _dynamics(0),
  78     mTowerDataLoaded(false),
  79     mRunwaysLoaded(false),
  80     mHelipadsLoaded(false),
  81     mTaxiwaysLoaded(false),
  82     mProceduresLoaded(false),
  83     mThresholdDataLoaded(false),
  84     mILSDataLoaded(false)
  85 {
  86 }
  87
  88
  89 FGAirport::~FGAirport()
  90 {
  91     delete _dynamics;
  92 }
  93
  94 bool FGAirport::isAirport() const
  95 {
  96   return type() == AIRPORT;
  97 }
  98
  99 bool FGAirport::isSeaport() const
 100 {
 101   return type() == SEAPORT;
 102 }
 103
 104 bool FGAirport::isHeliport() const
 105 {
 106   return type() == HELIPORT;
 107 }
 108
 109 bool FGAirport::isAirportType(FGPositioned* pos)
 110 {
 111     if (!pos) {
 112         return false;
 113     }
 114
 115     return (pos->type() >= AIRPORT) && (pos->type() <= SEAPORT);
 116 }
 117
 118 FGAirportDynamics * FGAirport::getDynamics()
 119 {
 120     if (_dynamics) {
 121         return _dynamics;
 122     }
 123
 124     _dynamics = new FGAirportDynamics(this);
 125     XMLLoader::load(_dynamics);
 126     _dynamics->init();
 127
 128     FGRunwayPreference rwyPrefs(this);
 129     XMLLoader::load(&rwyPrefs);
 130     _dynamics->setRwyUse(rwyPrefs);
 131
 132     return _dynamics;
 133 }
 134
 135 //------------------------------------------------------------------------------
 136 unsigned int FGAirport::numRunways() const
 137 {
 138   loadRunways();
 139   return mRunways.size();
 140 }
 141
 142 //------------------------------------------------------------------------------
 143 unsigned int FGAirport::numHelipads() const
 144 {
 145   loadHelipads();
 146   return mHelipads.size();
 147 }
 148
 149 //------------------------------------------------------------------------------
 150 FGRunwayRef FGAirport::getRunwayByIndex(unsigned int aIndex) const
 151 {
 152   loadRunways();
 153   return mRunways.at(aIndex);
 154 }
 155
 156 //------------------------------------------------------------------------------
 157 FGHelipadRef FGAirport::getHelipadByIndex(unsigned int aIndex) const
 158 {
 159   loadHelipads();
 160   return loadById<FGHelipad>(mHelipads, aIndex);
 161 }
 162
 163 //------------------------------------------------------------------------------
 164 FGRunwayMap FGAirport::getRunwayMap() const
 165 {
 166   loadRunways();
 167   FGRunwayMap map;
 168
 169   double minLengthFt = fgGetDouble("/sim/navdb/min-runway-length-ft");
 170
 171   BOOST_FOREACH(FGRunwayRef rwy, mRunways)
 172   {
 173     // ignore unusably short runways
 174     // TODO other methods don't check this...
 175     if( rwy->lengthFt() >= minLengthFt )
 176       map[ rwy->ident() ] = rwy;
 177   }
 178
 179   return map;
 180 }
 181
 182 //------------------------------------------------------------------------------
 183 FGHelipadMap FGAirport::getHelipadMap() const
 184 {
 185   loadHelipads();
 186   FGHelipadMap map;
 187
 188   BOOST_FOREACH(PositionedID id, mHelipads)
 189   {
 190     FGHelipad* rwy = loadById<FGHelipad>(id);
 191     map[ rwy->ident() ] = rwy;
 192   }
 193
 194   return map;
 195 }
 196
 197 //------------------------------------------------------------------------------
 198 bool FGAirport::hasRunwayWithIdent(const std::string& aIdent) const
 199 {
 200   loadRunways();
 201   BOOST_FOREACH(FGRunwayRef rwy, mRunways) {
 202     if (rwy->ident() == aIdent) {
 203       return true;
 204     }
 205   }
 206
 207   return false;
 208 }
 209
 210 //------------------------------------------------------------------------------
 211 bool FGAirport::hasHelipadWithIdent(const std::string& aIdent) const
 212 {
 213   return flightgear::NavDataCache::instance()
 214     ->airportItemWithIdent(guid(), FGPositioned::HELIPAD, aIdent) != 0;
 215 }
 216
 217 //------------------------------------------------------------------------------
 218 FGRunwayRef FGAirport::getRunwayByIdent(const std::string& aIdent) const
 219 {
 220   loadRunways();
 221   BOOST_FOREACH(FGRunwayRef rwy, mRunways) {
 222     if (rwy->ident() == aIdent) {
 223       return rwy;
 224     }
 225   }
 226
 227   SG_LOG(SG_GENERAL, SG_ALERT, "no such runway '" << aIdent << "' at airport " << ident());
 228   throw sg_range_exception("unknown runway " + aIdent + " at airport:" + ident(), "FGAirport::getRunwayByIdent");
 229 }
 230
 231 //------------------------------------------------------------------------------
 232 FGHelipadRef FGAirport::getHelipadByIdent(const std::string& aIdent) const
 233 {
 234   PositionedID id = flightgear::NavDataCache::instance()->airportItemWithIdent(guid(), FGPositioned::HELIPAD, aIdent);
 235   if (id == 0) {
 236     SG_LOG(SG_GENERAL, SG_ALERT, "no such helipad '" << aIdent << "' at airport " << ident());
 237     throw sg_range_exception("unknown helipad " + aIdent + " at airport:" + ident(), "FGAirport::getRunwayByIdent");
 238   }
 239
 240   return loadById<FGHelipad>(id);
 241 }
 242
 243 //------------------------------------------------------------------------------
 244 FGRunwayRef FGAirport::findBestRunwayForHeading(double aHeading, struct FindBestRunwayForHeadingParams * parms ) const
 245 {
 246   loadRunways();
 247
 248   FGRunway* result = NULL;
 249   double currentBestQuality = 0.0;
 250
 251   struct FindBestRunwayForHeadingParams fbrfhp;
 252   if( NULL != parms ) fbrfhp = *parms;
 253
 254   SGPropertyNode_ptr searchNode = fgGetNode("/sim/airport/runways/search");
 255   if( searchNode.valid() ) {
 256     fbrfhp.lengthWeight = searchNode->getDoubleValue("length-weight", fbrfhp.lengthWeight );
 257     fbrfhp.widthWeight = searchNode->getDoubleValue("width-weight", fbrfhp.widthWeight );
 258     fbrfhp.surfaceWeight = searchNode->getDoubleValue("surface-weight", fbrfhp.surfaceWeight );
 259     fbrfhp.deviationWeight = searchNode->getDoubleValue("deviation-weight", fbrfhp.deviationWeight );
 260     fbrfhp.ilsWeight = searchNode->getDoubleValue("ils-weight", fbrfhp.ilsWeight );
 261   }
 262
 263   BOOST_FOREACH(FGRunwayRef rwy, mRunways) {
 264     double good = rwy->score( fbrfhp.lengthWeight,  fbrfhp.widthWeight,  fbrfhp.surfaceWeight,  fbrfhp.ilsWeight );
 265     double dev = aHeading - rwy->headingDeg();
 266     SG_NORMALIZE_RANGE(dev, -180.0, 180.0);
 267     double bad = fabs( fbrfhp.deviationWeight * dev) + 1e-20;
 268     double quality = good / bad;
 269
 270     if (quality > currentBestQuality) {
 271       currentBestQuality = quality;
 272       result = rwy;
 273     }
 274   }
 275
 276   return result;
 277 }
 278
 279 //------------------------------------------------------------------------------
 280 FGRunwayRef FGAirport::findBestRunwayForPos(const SGGeod& aPos) const
 281 {
 282   loadRunways();
 283
 284   FGRunway* result = NULL;
 285   double currentLowestDev = 180.0;
 286
 287   BOOST_FOREACH(FGRunwayRef rwy, mRunways) {
 288     double inboundCourse = SGGeodesy::courseDeg(aPos, rwy->end());
 289     double dev = inboundCourse - rwy->headingDeg();
 290     SG_NORMALIZE_RANGE(dev, -180.0, 180.0);
 291
 292     dev = fabs(dev);
 293     if (dev < currentLowestDev) { // new best match
 294       currentLowestDev = dev;
 295       result = rwy;
 296     }
 297   } // of runway iteration
 298
 299   return result;
 300
 301 }
 302
 303 //------------------------------------------------------------------------------
 304 bool FGAirport::hasHardRunwayOfLengthFt(double aLengthFt) const
 305 {
 306   loadRunways();
 307
 308   BOOST_FOREACH(FGRunwayRef rwy, mRunways) {
 309     if (rwy->isHardSurface() && (rwy->lengthFt() >= aLengthFt)) {
 310       return true; // we're done!
 311     }
 312   } // of runways iteration
 313
 314   return false;
 315 }
 316
 317 //------------------------------------------------------------------------------
 318 FGRunwayList FGAirport::getRunways() const
 319 {
 320   loadRunways();
 321
 322   return mRunways;
 323 }
 324
 325 //------------------------------------------------------------------------------
 326 FGRunwayList FGAirport::getRunwaysWithoutReciprocals() const
 327 {
 328   loadRunways();
 329
 330   FGRunwayList r;
 331
 332   BOOST_FOREACH(FGRunwayRef rwy, mRunways) {
 333     FGRunway* recip = rwy->reciprocalRunway();
 334     if (recip) {
 335       FGRunwayList::iterator it = std::find(r.begin(), r.end(), recip);
 336       if (it != r.end()) {
 337         continue; // reciprocal already in result set, don't include us
 338       }
 339     }
 340
 341     r.push_back(rwy);
 342   }
 343
 344   return r;
 345 }
 346
 347 //------------------------------------------------------------------------------
 348 unsigned int FGAirport::numTaxiways() const
 349 {
 350   loadTaxiways();
 351   return mTaxiways.size();
 352 }
 353
 354 //------------------------------------------------------------------------------
 355 FGTaxiwayRef FGAirport::getTaxiwayByIndex(unsigned int aIndex) const
 356 {
 357   loadTaxiways();
 358   return loadById<FGTaxiway>(mTaxiways, aIndex);
 359 }
 360
 361 //------------------------------------------------------------------------------
 362 FGTaxiwayList FGAirport::getTaxiways() const
 363 {
 364   loadTaxiways();
 365   return loadAllById<FGTaxiway>(mTaxiways);
 366 }
 367
 368 //------------------------------------------------------------------------------
 369 unsigned int FGAirport::numPavements() const
 370 {
 371   loadTaxiways();
 372   return mPavements.size();
 373 }
 374
 375 //------------------------------------------------------------------------------
 376 FGPavementRef FGAirport::getPavementByIndex(unsigned int aIndex) const
 377 {
 378   loadTaxiways();
 379   return loadById<FGPavement>(mPavements, aIndex);
 380 }
 381
 382 //------------------------------------------------------------------------------
 383 FGPavementList FGAirport::getPavements() const
 384 {
 385   loadTaxiways();
 386   return loadAllById<FGPavement>(mPavements);
 387 }
 388
 389 //------------------------------------------------------------------------------
 390 FGRunwayRef FGAirport::getActiveRunwayForUsage() const
 391 {
 392   FGEnvironmentMgr* envMgr = (FGEnvironmentMgr *) globals->get_subsystem("environment");
 393
 394   // This forces West-facing rwys to be used in no-wind situations
 395   // which is consistent with Flightgear's initial setup.
 396   double hdg = 270;
 397
 398   if (envMgr) {
 399     FGEnvironment stationWeather(envMgr->getEnvironment(geod()));
 400
 401     double windSpeed = stationWeather.get_wind_speed_kt();
 402     if (windSpeed > 0.0) {
 403       hdg = stationWeather.get_wind_from_heading_deg();
 404     }
 405   }
 406
 407   return findBestRunwayForHeading(hdg);
 408 }
 409
 410 //------------------------------------------------------------------------------
 411 FGAirportRef FGAirport::findClosest( const SGGeod& aPos,
 412                                      double aCuttofNm,
 413                                      Filter* filter )
 414 {
 415   AirportFilter aptFilter;
 416   if( !filter )
 417     filter = &aptFilter;
 418
 419   return static_pointer_cast<FGAirport>
 420   (
 421     FGPositioned::findClosest(aPos, aCuttofNm, filter)
 422   );
 423 }
 424
 425 FGAirport::HardSurfaceFilter::HardSurfaceFilter(double minLengthFt) :
 426   mMinLengthFt(minLengthFt)
 427 {
 428   if (minLengthFt < 0.0) {
 429     mMinLengthFt = fgGetDouble("/sim/navdb/min-runway-length-ft", 0.0);
 430   }
 431 }
 432
 433 bool FGAirport::HardSurfaceFilter::passAirport(FGAirport* aApt) const
 434 {
 435   return aApt->hasHardRunwayOfLengthFt(mMinLengthFt);
 436 }
 437
 438 //------------------------------------------------------------------------------
 439 FGAirport::TypeRunwayFilter::TypeRunwayFilter():
 440   _type(FGPositioned::AIRPORT),
 441   _min_runway_length_ft( fgGetDouble("/sim/navdb/min-runway-length-ft", 0.0) )
 442 {
 443
 444 }
 445
 446 //------------------------------------------------------------------------------
 447 bool FGAirport::TypeRunwayFilter::fromTypeString(const std::string& type)
 448 {
 449   if(      type == "heliport" ) _type = FGPositioned::HELIPORT;
 450   else if( type == "seaport"  ) _type = FGPositioned::SEAPORT;
 451   else if( type == "airport"  ) _type = FGPositioned::AIRPORT;
 452   else                          return false;
 453
 454   return true;
 455 }
 456
 457 //------------------------------------------------------------------------------
 458 bool FGAirport::TypeRunwayFilter::pass(FGPositioned* pos) const
 459 {
 460   FGAirport* apt = static_cast<FGAirport*>(pos);
 461   if(  (apt->type() == FGPositioned::AIRPORT)
 462     && !apt->hasHardRunwayOfLengthFt(_min_runway_length_ft)
 463     )
 464     return false;
 465
 466   return true;
 467 }
 468
 469 //------------------------------------------------------------------------------
 470 FGAirportRef FGAirport::findByIdent(const std::string& aIdent)
 471 {
 472   AirportCache::iterator it = airportCache.find(aIdent);
 473   if (it != airportCache.end())
 474    return it->second;
 475
 476   PortsFilter filter;
 477   FGAirportRef r = static_pointer_cast<FGAirport>
 478   (
 479     FGPositioned::findFirstWithIdent(aIdent, &filter)
 480   );
 481
 482   // add airport to the cache (even when it's NULL, so we don't need to search in vain again)
 483   airportCache[aIdent] = r;
 484
 485   // we don't warn here when r==NULL, let the caller do that
 486   return r;
 487 }
 488
 489 //------------------------------------------------------------------------------
 490 FGAirportRef FGAirport::getByIdent(const std::string& aIdent)
 491 {
 492   FGAirportRef r = findByIdent(aIdent);
 493   if (!r)
 494     throw sg_range_exception("No such airport with ident: " + aIdent);
 495   return r;
 496 }
 497
 498 char** FGAirport::searchNamesAndIdents(const std::string& aFilter)
 499 {
 500   return NavDataCache::instance()->searchAirportNamesAndIdents(aFilter);
 501 }
 502
 503 // find basic airport location info from airport database
 504 const FGAirport *fgFindAirportID( const std::string& id)
 505 {
 506     if ( id.empty() ) {
 507         return NULL;
 508     }
 509
 510     return FGAirport::findByIdent(id);
 511 }
 512
 513 PositionedIDVec FGAirport::itemsOfType(FGPositioned::Type ty) const
 514 {
 515   flightgear::NavDataCache* cache = flightgear::NavDataCache::instance();
 516   return cache->airportItemsOfType(guid(), ty);
 517 }
 518
 519 void FGAirport::loadRunways() const
 520 {
 521   if (mRunwaysLoaded) {
 522     return; // already loaded, great
 523   }
 524
 525   loadSceneryDefinitions();
 526
 527   mRunwaysLoaded = true;
 528   PositionedIDVec rwys(itemsOfType(FGPositioned::RUNWAY));
 529   BOOST_FOREACH(PositionedID id, rwys) {
 530     mRunways.push_back(loadById<FGRunway>(id));
 531   }
 532 }
 533
 534 void FGAirport::loadHelipads() const
 535 {
 536   if (mHelipadsLoaded) {
 537     return; // already loaded, great
 538   }
 539
 540   mHelipadsLoaded = true;
 541   mHelipads = itemsOfType(FGPositioned::HELIPAD);
 542 }
 543
 544 void FGAirport::loadTaxiways() const
 545 {
 546   if (mTaxiwaysLoaded) {
 547     return; // already loaded, great
 548   }
 549
 550   mTaxiwaysLoaded =  true;
 551   mTaxiways = itemsOfType(FGPositioned::TAXIWAY);
 552 }
 553
 554 void FGAirport::loadProcedures() const
 555 {
 556   if (mProceduresLoaded) {
 557     return;
 558   }
 559
 560   mProceduresLoaded = true;
 561   SGPath path;
 562   if (!XMLLoader::findAirportData(ident(), "procedures", path)) {
 563     SG_LOG(SG_GENERAL, SG_INFO, "no procedures data available for " << ident());
 564     return;
 565   }
 566
 567   SG_LOG(SG_GENERAL, SG_INFO, ident() << ": loading procedures from " << path.str());
 568   RouteBase::loadAirportProcedures(path, const_cast<FGAirport*>(this));
 569 }
 570
 571 void FGAirport::loadSceneryDefinitions() const
 572 {
 573   if (mThresholdDataLoaded) {
 574     return;
 575   }
 576
 577   mThresholdDataLoaded = true;
 578
 579   SGPath path;
 580   if (!XMLLoader::findAirportData(ident(), "threshold", path)) {
 581     return; // no XML threshold data
 582   }
 583
 584   try {
 585     SGPropertyNode_ptr rootNode = new SGPropertyNode;
 586     readProperties(path.str(), rootNode);
 587     const_cast<FGAirport*>(this)->readThresholdData(rootNode);
 588   } catch (sg_exception& e) {
 589     SG_LOG(SG_NAVAID, SG_WARN, ident() << "loading threshold XML failed:" << e.getFormattedMessage());
 590   }
 591 }
 592
 593 void FGAirport::readThresholdData(SGPropertyNode* aRoot)
 594 {
 595   SGPropertyNode* runway;
 596   int runwayIndex = 0;
 597   for (; (runway = aRoot->getChild("runway", runwayIndex)) != NULL; ++runwayIndex) {
 598     SGPropertyNode* t0 = runway->getChild("threshold", 0),
 599       *t1 = runway->getChild("threshold", 1);
 600     assert(t0);
 601     assert(t1); // too strict? maybe we should finally allow single-ended runways
 602
 603     processThreshold(t0);
 604     processThreshold(t1);
 605   } // of runways iteration
 606 }
 607
 608 void FGAirport::processThreshold(SGPropertyNode* aThreshold)
 609 {
 610   // first, let's identify the current runway
 611   std::string rwyIdent(aThreshold->getStringValue("rwy"));
 612   NavDataCache* cache = NavDataCache::instance();
 613   PositionedID id = cache->airportItemWithIdent(guid(), FGPositioned::RUNWAY, rwyIdent);
 614
 615   double lon = aThreshold->getDoubleValue("lon"),
 616   lat = aThreshold->getDoubleValue("lat");
 617   SGGeod newThreshold(SGGeod::fromDegM(lon, lat, elevationM()));
 618
 619   double newHeading = aThreshold->getDoubleValue("hdg-deg");
 620   double newDisplacedThreshold = aThreshold->getDoubleValue("displ-m");
 621   double newStopway = aThreshold->getDoubleValue("stopw-m");
 622
 623   if (id == 0) {
 624     SG_LOG(SG_GENERAL, SG_DEBUG, "FGAirport::processThreshold: "
 625            "found runway not defined in the global data:" << ident() << "/" << rwyIdent);
 626     // enable this code when threshold.xml contains sufficient data to
 627     // fully specify a new runway, *and* we figure out how to assign runtime
 628     // Positioned IDs and insert temporary items into the spatial map.
 629 #if 0
 630     double newLength = 0.0, newWidth = 0.0;
 631     int surfaceCode = 0;
 632     FGRunway* rwy = new FGRunway(id, guid(), rwyIdent, newThreshold,
 633                        newHeading,
 634                        newLength, newWidth,
 635                        newDisplacedThreshold, newStopway,
 636                        surfaceCode);
 637     // insert into the spatial map too
 638     mRunways.push_back(rwy);
 639 #endif
 640   } else {
 641     FGRunway* rwy = loadById<FGRunway>(id);
 642     rwy->updateThreshold(newThreshold, newHeading,
 643                          newDisplacedThreshold, newStopway);
 644
 645   }
 646 }
 647
 648 SGGeod FGAirport::getTowerLocation() const
 649 {
 650   validateTowerData();
 651   return mTowerPosition;
 652 }
 653
 654 void FGAirport::validateTowerData() const
 655 {
 656   if (mTowerDataLoaded) {
 657     return;
 658   }
 659
 660   mTowerDataLoaded = true;
 661
 662 // first, load data from the cache (apt.dat)
 663   NavDataCache* cache = NavDataCache::instance();
 664   PositionedIDVec towers = cache->airportItemsOfType(guid(), FGPositioned::TOWER);
 665   if (towers.empty()) {
 666     SG_LOG(SG_GENERAL, SG_ALERT, "No towers defined for:" <<ident());
 667     mTowerPosition = geod(); // use airport position
 668     // increase tower elevation by 20 metres above the field elevation
 669     mTowerPosition.setElevationM(geod().getElevationM() + 20.0);
 670   } else {
 671     FGPositionedRef tower = cache->loadById(towers.front());
 672     mTowerPosition = tower->geod();
 673   }
 674
 675   SGPath path;
 676   if (!XMLLoader::findAirportData(ident(), "twr", path)) {
 677     return; // no XML tower data, base position is fine
 678   }
 679
 680   try {
 681     SGPropertyNode_ptr rootNode = new SGPropertyNode;
 682     readProperties(path.str(), rootNode);
 683     const_cast<FGAirport*>(this)->readTowerData(rootNode);
 684   } catch (sg_exception& e){
 685     SG_LOG(SG_NAVAID, SG_WARN, ident() << "loading twr XML failed:" << e.getFormattedMessage());
 686   }
 687 }
 688
 689 void FGAirport::readTowerData(SGPropertyNode* aRoot)
 690 {
 691   SGPropertyNode* twrNode = aRoot->getChild("tower")->getChild("twr");
 692   double lat = twrNode->getDoubleValue("lat"),
 693     lon = twrNode->getDoubleValue("lon"),
 694     elevM = twrNode->getDoubleValue("elev-m");
 695 // tower elevation is AGL, not AMSL. Since we don't want to depend on the
 696 // scenery for a precise terrain elevation, we use the field elevation
 697 // (this is also what the apt.dat code does)
 698   double fieldElevationM = geod().getElevationM();
 699   mTowerPosition = SGGeod::fromDegM(lon, lat, fieldElevationM + elevM);
 700 }
 701
 702 void FGAirport::validateILSData()
 703 {
 704   if (mILSDataLoaded) {
 705     return;
 706   }
 707
 708   // to avoid re-entrancy on this code-path, ensure we set loaded
 709   // immediately.
 710   mILSDataLoaded = true;
 711
 712   SGPath path;
 713   if (!XMLLoader::findAirportData(ident(), "ils", path)) {
 714     return; // no XML tower data
 715   }
 716
 717   try {
 718       SGPropertyNode_ptr rootNode = new SGPropertyNode;
 719       readProperties(path.str(), rootNode);
 720       readILSData(rootNode);
 721   } catch (sg_exception& e){
 722       SG_LOG(SG_NAVAID, SG_WARN, ident() << "loading ils XML failed:" << e.getFormattedMessage());
 723   }
 724 }
 725
 726 void FGAirport::readILSData(SGPropertyNode* aRoot)
 727 {
 728   NavDataCache* cache = NavDataCache::instance();
 729   // find the entry matching the runway
 730   SGPropertyNode* runwayNode, *ilsNode;
 731   for (int i=0; (runwayNode = aRoot->getChild("runway", i)) != NULL; ++i) {
 732     for (int j=0; (ilsNode = runwayNode->getChild("ils", j)) != NULL; ++j) {
 733       // must match on both nav-ident and runway ident, to support the following:
 734       // - runways with multiple distinct ILS installations (KEWD, for example)
 735       // - runways where both ends share the same nav ident (LFAT, for example)
 736       PositionedID ils = cache->findILS(guid(), ilsNode->getStringValue("rwy"),
 737                                         ilsNode->getStringValue("nav-id"));
 738       if (ils == 0) {
 739         SG_LOG(SG_GENERAL, SG_INFO, "reading ILS data for " << ident() <<
 740                ", couldn't find runway/navaid for:" <<
 741                ilsNode->getStringValue("rwy") << "/" <<
 742                ilsNode->getStringValue("nav-id"));
 743         continue;
 744       }
 745
 746       double hdgDeg = ilsNode->getDoubleValue("hdg-deg"),
 747         lon = ilsNode->getDoubleValue("lon"),
 748         lat = ilsNode->getDoubleValue("lat"),
 749         elevM = ilsNode->getDoubleValue("elev-m");
 750
 751       FGNavRecordRef nav(FGPositioned::loadById<FGNavRecord>(ils));
 752       assert(nav.valid());
 753       nav->updateFromXML(SGGeod::fromDegM(lon, lat, elevM), hdgDeg);
 754     } // of ILS iteration
 755   } // of runway iteration
 756 }
 757
 758 void FGAirport::addSID(flightgear::SID* aSid)
 759 {
 760   mSIDs.push_back(aSid);
 761 }
 762
 763 void FGAirport::addSTAR(STAR* aStar)
 764 {
 765   mSTARs.push_back(aStar);
 766 }
 767
 768 void FGAirport::addApproach(Approach* aApp)
 769 {
 770   mApproaches.push_back(aApp);
 771 }
 772
 773 //------------------------------------------------------------------------------
 774 unsigned int FGAirport::numSIDs() const
 775 {
 776   loadProcedures();
 777   return mSIDs.size();
 778 }
 779
 780 //------------------------------------------------------------------------------
 781 flightgear::SID* FGAirport::getSIDByIndex(unsigned int aIndex) const
 782 {
 783   loadProcedures();
 784   return mSIDs[aIndex];
 785 }
 786
 787 //------------------------------------------------------------------------------
 788 flightgear::SID* FGAirport::findSIDWithIdent(const std::string& aIdent) const
 789 {
 790   loadProcedures();
 791   for (unsigned int i=0; i<mSIDs.size(); ++i) {
 792     if (mSIDs[i]->ident() == aIdent) {
 793       return mSIDs[i];
 794     }
 795   }
 796
 797   return NULL;
 798 }
 799
 800 //------------------------------------------------------------------------------
 801 flightgear::SIDList FGAirport::getSIDs() const
 802 {
 803   loadProcedures();
 804   return flightgear::SIDList(mSIDs.begin(), mSIDs.end());
 805 }
 806
 807 //------------------------------------------------------------------------------
 808 unsigned int FGAirport::numSTARs() const
 809 {
 810   loadProcedures();
 811   return mSTARs.size();
 812 }
 813
 814 //------------------------------------------------------------------------------
 815 STAR* FGAirport::getSTARByIndex(unsigned int aIndex) const
 816 {
 817   loadProcedures();
 818   return mSTARs[aIndex];
 819 }
 820
 821 //------------------------------------------------------------------------------
 822 STAR* FGAirport::findSTARWithIdent(const std::string& aIdent) const
 823 {
 824   loadProcedures();
 825   for (unsigned int i=0; i<mSTARs.size(); ++i) {
 826     if (mSTARs[i]->ident() == aIdent) {
 827       return mSTARs[i];
 828     }
 829   }
 830
 831   return NULL;
 832 }
 833
 834 //------------------------------------------------------------------------------
 835 STARList FGAirport::getSTARs() const
 836 {
 837   loadProcedures();
 838   return STARList(mSTARs.begin(), mSTARs.end());
 839 }
 840
 841 unsigned int FGAirport::numApproaches() const
 842 {
 843   loadProcedures();
 844   return mApproaches.size();
 845 }
 846
 847 //------------------------------------------------------------------------------
 848 Approach* FGAirport::getApproachByIndex(unsigned int aIndex) const
 849 {
 850   loadProcedures();
 851   return mApproaches[aIndex];
 852 }
 853
 854 //------------------------------------------------------------------------------
 855 Approach* FGAirport::findApproachWithIdent(const std::string& aIdent) const
 856 {
 857   loadProcedures();
 858   for (unsigned int i=0; i<mApproaches.size(); ++i) {
 859     if (mApproaches[i]->ident() == aIdent) {
 860       return mApproaches[i];
 861     }
 862   }
 863
 864   return NULL;
 865 }
 866
 867 //------------------------------------------------------------------------------
 868 ApproachList FGAirport::getApproaches(ProcedureType type) const
 869 {
 870   loadProcedures();
 871   if( type == PROCEDURE_INVALID )
 872     return ApproachList(mApproaches.begin(), mApproaches.end());
 873
 874   ApproachList ret;
 875   for(size_t i = 0; i < mApproaches.size(); ++i)
 876   {
 877     if( mApproaches[i]->type() == type )
 878       ret.push_back(mApproaches[i]);
 879   }
 880   return ret;
 881 }
 882
 883 CommStationList
 884 FGAirport::commStations() const
 885 {
 886   NavDataCache* cache = NavDataCache::instance();
 887   CommStationList result;
 888   BOOST_FOREACH(PositionedID pos, cache->airportItemsOfType(guid(),
 889                                                             FGPositioned::FREQ_GROUND,
 890                                                             FGPositioned::FREQ_UNICOM))
 891   {
 892     result.push_back( loadById<CommStation>(pos) );
 893   }
 894
 895   return result;
 896 }
 897
 898 CommStationList
 899 FGAirport::commStationsOfType(FGPositioned::Type aTy) const
 900 {
 901   NavDataCache* cache = NavDataCache::instance();
 902   CommStationList result;
 903   BOOST_FOREACH(PositionedID pos, cache->airportItemsOfType(guid(), aTy)) {
 904     result.push_back( loadById<CommStation>(pos) );
 905   }
 906
 907   return result;
 908 }
 909
 910 class AirportWithSize
 911 {
 912 public:
 913     AirportWithSize(FGPositionedRef pos) :
 914         _pos(pos),
 915         _sizeMetric(0)
 916     {
 917         assert(pos->type() == FGPositioned::AIRPORT);
 918         FGAirport* apt = static_cast<FGAirport*>(pos.get());
 919         BOOST_FOREACH(FGRunway* rwy, apt->getRunwaysWithoutReciprocals()) {
 920             _sizeMetric += static_cast<int>(rwy->lengthFt());
 921         }
 922     }
 923
 924     bool operator<(const AirportWithSize& other) const
 925     {
 926         return _sizeMetric < other._sizeMetric;
 927     }
 928
 929     FGPositionedRef pos() const
 930     { return _pos; }
 931 private:
 932     FGPositionedRef _pos;
 933     unsigned int _sizeMetric;
 934
 935 };
 936
 937 void FGAirport::sortBySize(FGPositionedList& airportList)
 938 {
 939     std::vector<AirportWithSize> annotated;
 940     BOOST_FOREACH(FGPositionedRef p, airportList) {
 941         annotated.push_back(AirportWithSize(p));
 942     }
 943     std::sort(annotated.begin(), annotated.end());
 944
 945     for (unsigned int i=0; i<annotated.size(); ++i) {
 946         airportList[i] = annotated[i].pos();
 947     }
 948 }
 949
 950 // get airport elevation
 951 double fgGetAirportElev( const std::string& id )
 952 {
 953     const FGAirport *a=fgFindAirportID( id);
 954     if (a) {
 955         return a->getElevation();
 956     } else {
 957         return -9999.0;
 958     }
 959 }
 960
 961
 962 // get airport position
 963 SGGeod fgGetAirportPos( const std::string& id )
 964 {
 965     const FGAirport *a = fgFindAirportID( id);
 966
 967     if (a) {
 968         return SGGeod::fromDegM(a->getLongitude(), a->getLatitude(), a->getElevation());
 969     } else {
 970         return SGGeod::fromDegM(0.0, 0.0, -9999.0);
 971     }
 972 }