2 // simple.cxx -- a really simplistic class to manage airport ID,
3 // lat, lon of the center of one of it's runways, and
6 // Written by Curtis Olson, started April 1998.
7 // Updated by Durk Talsma, started December, 2004.
9 // Copyright (C) 1998 Curtis L. Olson - http://www.flightgear.org/~curt
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.
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.
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.
34 #include <simgear/compiler.h>
36 #include <Environment/environment_mgr.hxx>
37 #include <Environment/environment.hxx>
38 #include <simgear/misc/sg_path.hxx>
39 #include <simgear/props/props.hxx>
40 #include <simgear/structure/subsystem_mgr.hxx>
41 #include <simgear/debug/logstream.hxx>
42 #include <Main/globals.hxx>
43 #include <Main/fg_props.hxx>
44 #include <Airports/runways.hxx>
45 #include <Airports/dynamics.hxx>
50 #include "xmlloader.hxx"
53 using std::random_shuffle;
58 /***************************************************************************
60 ***************************************************************************/
62 FGAirport::FGAirport(const string &id, const SGGeod& location, const SGGeod& tower_location,
63 const string &name, bool has_metar, Type aType) :
64 FGPositioned(aType, id, location),
65 _tower_location(tower_location),
67 _has_metar(has_metar),
73 FGAirport::~FGAirport()
78 bool FGAirport::isAirport() const
80 return type() == AIRPORT;
83 bool FGAirport::isSeaport() const
85 return type() == SEAPORT;
88 bool FGAirport::isHeliport() const
90 return type() == HELIPORT;
93 FGAirportDynamics * FGAirport::getDynamics()
98 //cerr << "Trying to load dynamics for " << _id << endl;
99 _dynamics = new FGAirportDynamics(this);
100 XMLLoader::load(_dynamics);
102 FGRunwayPreference rwyPrefs(this);
103 XMLLoader::load(&rwyPrefs);
104 _dynamics->setRwyUse(rwyPrefs);
109 unsigned int FGAirport::numRunways() const
111 return mRunways.size();
114 FGRunway* FGAirport::getRunwayByIndex(unsigned int aIndex) const
116 assert(aIndex >= 0 && aIndex < mRunways.size());
117 return mRunways[aIndex];
120 bool FGAirport::hasRunwayWithIdent(const string& aIdent) const
122 return (getIteratorForRunwayIdent(aIdent) != mRunways.end());
125 FGRunway* FGAirport::getRunwayByIdent(const string& aIdent) const
127 Runway_iterator it = getIteratorForRunwayIdent(aIdent);
128 if (it == mRunways.end()) {
129 SG_LOG(SG_GENERAL, SG_ALERT, "no such runway '" << aIdent << "' at airport " << ident());
130 throw sg_range_exception("unknown runway " + aIdent + " at airport:" + ident(), "FGAirport::getRunwayByIdent");
136 FGAirport::Runway_iterator
137 FGAirport::getIteratorForRunwayIdent(const string& aIdent) const
139 string ident(aIdent);
140 if ((aIdent.size() == 1) || !isdigit(aIdent[1])) {
141 ident = "0" + aIdent;
144 Runway_iterator it = mRunways.begin();
145 for (; it != mRunways.end(); ++it) {
146 if ((*it)->ident() == ident) {
154 static double normaliseBearing(double aBearing)
156 while (aBearing < -180) {
160 while (aBearing > 180.0) {
167 FGRunway* FGAirport::findBestRunwayForHeading(double aHeading) const
169 Runway_iterator it = mRunways.begin();
170 FGRunway* result = NULL;
171 double currentBestQuality = 0.0;
173 SGPropertyNode *param = fgGetNode("/sim/airport/runways/search", true);
174 double lengthWeight = param->getDoubleValue("length-weight", 0.01);
175 double widthWeight = param->getDoubleValue("width-weight", 0.01);
176 double surfaceWeight = param->getDoubleValue("surface-weight", 10);
177 double deviationWeight = param->getDoubleValue("deviation-weight", 1);
179 for (; it != mRunways.end(); ++it) {
180 double good = (*it)->score(lengthWeight, widthWeight, surfaceWeight);
182 double dev = normaliseBearing(aHeading - (*it)->headingDeg());
183 double bad = fabs(deviationWeight * dev) + 1e-20;
184 double quality = good / bad;
186 if (quality > currentBestQuality) {
187 currentBestQuality = quality;
195 bool FGAirport::hasHardRunwayOfLengthFt(double aLengthFt) const
197 unsigned int numRunways(mRunways.size());
198 for (unsigned int r=0; r<numRunways; ++r) {
199 FGRunway* rwy = mRunways[r];
200 if (rwy->isReciprocal()) {
201 continue; // we only care about lengths, so don't do work twice
204 if (rwy->isHardSurface() && (rwy->lengthFt() >= aLengthFt)) {
205 return true; // we're done!
207 } // of runways iteration
212 unsigned int FGAirport::numTaxiways() const
214 return mTaxiways.size();
217 FGRunway* FGAirport::getTaxiwayByIndex(unsigned int aIndex) const
219 assert(aIndex >= 0 && aIndex < mTaxiways.size());
220 return mTaxiways[aIndex];
223 void FGAirport::addRunway(FGRunway* aRunway)
225 aRunway->setAirport(this);
227 if (aRunway->isTaxiway()) {
228 mTaxiways.push_back(aRunway);
230 mRunways.push_back(aRunway);
234 FGRunway* FGAirport::getActiveRunwayForUsage() const
236 static FGEnvironmentMgr* envMgr = NULL;
238 envMgr = (FGEnvironmentMgr *) globals->get_subsystem("environment");
241 FGEnvironment stationWeather(envMgr->getEnvironment(mPosition));
243 double windSpeed = stationWeather.get_wind_speed_kt();
244 double hdg = stationWeather.get_wind_from_heading_deg();
245 if (windSpeed <= 0.0) {
246 hdg = 270; // This forces West-facing rwys to be used in no-wind situations
247 // which is consistent with Flightgear's initial setup.
250 return findBestRunwayForHeading(hdg);
253 FGAirport* FGAirport::findClosest(const SGGeod& aPos, double aCuttofNm, Filter* filter)
255 AirportFilter aptFilter;
256 if (filter == NULL) {
260 FGPositionedRef r = FGPositioned::findClosest(aPos, aCuttofNm, filter);
265 return static_cast<FGAirport*>(r.ptr());
268 FGAirport::HardSurfaceFilter::HardSurfaceFilter(double minLengthFt) :
269 mMinLengthFt(minLengthFt)
273 bool FGAirport::HardSurfaceFilter::pass(FGPositioned* aPos) const
275 if (aPos->type() != AIRPORT) {
276 return false; // exclude seaports and heliports as well, we need a runways
279 return static_cast<FGAirport*>(aPos)->hasHardRunwayOfLengthFt(mMinLengthFt);
282 /******************************************************************************
284 *****************************************************************************/
286 // Populates a list of subdirectories of $FG_ROOT/Airports/AI so that
287 // the add() method doesn't have to try opening 2 XML files in each of
288 // thousands of non-existent directories. FIXME: should probably add
289 // code to free this list after parsing of apt.dat is finished;
290 // non-issue at the moment, however, as there are no AI subdirectories
291 // in the base package.
293 // Note: 2005/12/23: This is probably not necessary anymore, because I'm
294 // Switching to runtime airport dynamics loading (DT).
295 FGAirportList::FGAirportList()
299 // SGPath aid( globals->get_fg_root() );
300 // aid.append( "/Airports/AI" );
301 // if((d = ulOpenDir(aid.c_str())) == NULL)
303 // while((dent = ulReadDir(d)) != NULL) {
304 // SG_LOG( SG_GENERAL, SG_DEBUG, "Dent: " << dent->d_name );
305 // ai_dirs.insert(dent->d_name);
311 FGAirportList::~FGAirportList( void )
313 for (unsigned int i = 0; i < airports_array.size(); ++i) {
314 delete airports_array[i];
319 // add an entry to the list
320 FGAirport* FGAirportList::add( const string &id, const SGGeod& location, const SGGeod& tower_location,
321 const string &name, bool has_metar, FGPositioned::Type aType)
323 FGAirport* a = new FGAirport(id, location, tower_location, name, has_metar, aType);
324 airports_by_id[a->getId()] = a;
325 // try and read in an auxilary file
327 airports_array.push_back( a );
331 // search for the specified id
332 FGAirport* FGAirportList::search( const string& id)
334 airport_map_iterator itr = airports_by_id.find(id);
335 return (itr == airports_by_id.end() ? NULL : itr->second);
339 FGAirportList::size () const
341 return airports_array.size();
345 const FGAirport *FGAirportList::getAirport( unsigned int index ) const
347 if (index < airports_array.size()) {
348 return(airports_array[index]);
354 // find basic airport location info from airport database
355 const FGAirport *fgFindAirportID( const string& id)
357 const FGAirport* result = NULL;
359 SG_LOG( SG_GENERAL, SG_BULK, "Searching for airport code = " << id );
361 result = globals->get_airports()->search( id );
363 if ( result == NULL ) {
364 SG_LOG( SG_GENERAL, SG_ALERT,
365 "Failed to find " << id << " in apt.dat.gz" );
371 SG_LOG( SG_GENERAL, SG_BULK,
372 "Position for " << id << " is ("
373 << result->getLongitude() << ", "
374 << result->getLatitude() << ")" );
380 // get airport elevation
381 double fgGetAirportElev( const string& id )
383 SG_LOG( SG_GENERAL, SG_BULK,
384 "Finding elevation for airport: " << id );
386 const FGAirport *a=fgFindAirportID( id);
388 return a->getElevation();
395 // get airport position
396 Point3D fgGetAirportPos( const string& id )
398 SG_LOG( SG_ATC, SG_BULK,
399 "Finding position for airport: " << id );
401 const FGAirport *a = fgFindAirportID( id);
404 return Point3D(a->getLongitude(), a->getLatitude(), a->getElevation());
406 return Point3D(0.0, 0.0, -9999.0);