1 // runways.cxx -- a simple class to manage airport runway info
3 // Written by Curtis Olson, started August 2000.
5 // Copyright (C) 2000 Curtis L. Olson - http://www.flightgear.org/~curt
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., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
28 #include <cstdio> // sprintf()
29 #include <cstdlib> // atoi()
32 #include <simgear/compiler.h>
34 #include <simgear/props/props.hxx>
38 #include "runways.hxx"
42 static std::string cleanRunwayNo(const std::string& aRwyNo)
44 if (aRwyNo[0] == 'x') {
45 return std::string(); // no ident for taxiways
48 string result(aRwyNo);
49 // canonicalise runway ident
50 if ((aRwyNo.size() == 1) || !isdigit(aRwyNo[1])) {
51 result = "0" + aRwyNo;
54 // trim off trailing garbage
55 if (result.size() > 2) {
56 char suffix = toupper(result[2]);
58 result = result.substr(0, 2);
65 FGRunway::FGRunway(FGAirport* aAirport, const string& aIdent,
67 const double heading, const double length,
69 const double displ_thresh,
71 const int surface_code,
73 FGRunwayBase(RUNWAY, cleanRunwayNo(aIdent), aGeod, heading, length, width, surface_code, true),
75 _isReciprocal(reciprocal),
77 _displ_thresh(displ_thresh),
83 string FGRunway::reverseIdent(const string& aRunwayIdent)
85 // Helipads don't have a seperate number per end
86 if (aRunwayIdent.size() && (aRunwayIdent[0] == 'H' || aRunwayIdent[0] == 'h' || aRunwayIdent[0] == 'x')) {
90 std::string ident(aRunwayIdent);
93 int rn = atoi(ident.substr(0,2).c_str());
99 sprintf(buf, "%02i", rn);
100 if(ident.size() == 3) {
101 char suffix = toupper(ident[2]);
104 } else if (suffix == 'R') {
107 // something else, just copy
117 double FGRunway::score(double aLengthWt, double aWidthWt, double aSurfaceWt) const
120 if (_surface_code == 12 || _surface_code == 5) // dry lakebed & gravel
122 else if (_surface_code == 1 || _surface_code == 2) // asphalt & concrete
125 return _length * aLengthWt + _width * aWidthWt + surface * aSurfaceWt + 1e-20;
128 SGGeod FGRunway::begin() const
130 return pointOnCenterline(0.0);
133 SGGeod FGRunway::end() const
135 return pointOnCenterline(lengthM());
138 SGGeod FGRunway::threshold() const
140 return pointOnCenterline(_displ_thresh * SG_FEET_TO_METER);
143 void FGRunway::processThreshold(SGPropertyNode* aThreshold)
145 assert(ident() == aThreshold->getStringValue("rwy"));
147 double lon = aThreshold->getDoubleValue("lon"),
148 lat = aThreshold->getDoubleValue("lat");
149 SGGeod newThreshold(SGGeod::fromDegM(lon, lat, mPosition.getElevationM()));
151 _heading = aThreshold->getDoubleValue("hdg-deg");
152 _displ_thresh = aThreshold->getDoubleValue("displ-m") * SG_METER_TO_FEET;
153 _stopway = aThreshold->getDoubleValue("stopw-m") * SG_METER_TO_FEET;
155 // compute the new runway center, based on the threshold lat/lon and length,
156 double offsetFt = (0.5 * _length);
159 SGGeodesy::direct(newThreshold, _heading, offsetFt * SG_FEET_TO_METER, newCenter, dummy);
160 mPosition = newCenter;
163 void FGRunway::setReciprocalRunway(FGRunway* other)
165 assert(_reciprocal==NULL);
166 assert((other->_reciprocal == NULL) || (other->_reciprocal == this));