]> git.mxchange.org Git - flightgear.git/blob - src/Airports/runways.cxx
Added a <solve-weight> subtag of the approach/cruise parameters that can
[flightgear.git] / src / Airports / runways.cxx
1 // runways.cxx -- a simple class to manage airport runway info
2 //
3 // Written by Curtis Olson, started August 2000.
4 //
5 // Copyright (C) 2000  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 <math.h>               // fabs()
29 #include <stdio.h>              // sprintf()
30
31 #include <simgear/compiler.h>
32
33 #include <simgear/debug/logstream.hxx>
34 #include <simgear/misc/sgstream.hxx>
35
36 #include STL_STRING
37 #include STL_IOSTREAM
38 #include <map>
39
40 #include "runways.hxx"
41
42 SG_USING_NAMESPACE(std);
43 SG_USING_STD(istream);
44 SG_USING_STD(multimap);
45
46 inline istream&
47 operator >> ( istream& in, FGRunway& a )
48 {
49     string type;
50     int tmp;
51
52     in >> a.type;
53     if ( a.type == "R" ) {
54         in >> a.id >> a.rwy_no >> a.lat >> a.lon >> a.heading
55            >> a.length >> a.width >> a.surface_flags >> a.end1_flags
56            >> tmp >> tmp >> a.end2_flags >> tmp >> tmp;
57     } else if ( a.type == "T" ) {
58         // in >> a.id >> a.rwy_no >> a.lat >> a.lon >> a.heading
59         //    >> a.length >> a.width >> a.surface_flags;
60         in >> skipeol;
61     } else {
62         in >> skipeol;
63     }
64
65     return in;
66 }
67
68
69 FGRunwayList::FGRunwayList( const string& file ) {
70     SG_LOG( SG_GENERAL, SG_DEBUG, "Reading runway list: " << file );
71
72     // open the specified file for reading
73     sg_gzifstream in( file );
74     if ( !in.is_open() ) {
75         SG_LOG( SG_GENERAL, SG_ALERT, "Cannot open file: " << file );
76         exit(-1);
77     }
78
79     // skip header line
80     in >> skipeol;
81
82     FGRunway rwy;
83     while ( in ) {
84         in >> rwy;
85         if(rwy.type == "R") {
86             runways.insert(pair<const string, FGRunway>(rwy.id, rwy));
87         }
88     }
89 }
90
91
92 // Return reverse rwy number
93 // eg 01 -> 19
94 // 03L -> 21R
95 static string GetReverseRunwayNo(string rwyno) {        
96     // cout << "Original rwyno = " << rwyNo << '\n';
97     
98     // standardize input number
99     string tmp = rwyno.substr(1, 1);
100     if (( tmp == "L" || tmp == "R" || tmp == "C" ) || (rwyno.size() == 1)) {
101         tmp = rwyno;
102         rwyno = "0" + tmp;
103         SG_LOG( SG_GENERAL, SG_INFO, "Standardising rwy number from " << tmp
104                                      << " to " << rwyno );
105     }
106     
107     char buf[4];
108     int rn = atoi(rwyno.substr(0,2).c_str());
109     rn += 18;
110     while(rn > 36) {
111         rn -= 36;
112     }
113     sprintf(buf, "%02i", rn);
114     if(rwyno.size() == 3) {
115         if(rwyno.substr(2,1) == "L") {
116             buf[2] = 'R';
117             buf[3] = '\0';
118         } else if (rwyno.substr(2,1) == "R") {
119             buf[2] = 'L';
120             buf[3] = '\0';
121         } else if (rwyno.substr(2,1) == "C") {
122             buf[2] = 'C';
123             buf[3] = '\0';
124         } else {
125             SG_LOG(SG_GENERAL, SG_ALERT, "Unknown runway code "
126             << rwyno << " passed to GetReverseRunwayNo(...)");
127         }
128     }
129     return(buf);
130 }
131
132
133 // search for the specified apt id (wierd!)
134 bool FGRunwayList::search( const string& aptid, FGRunway* r ) {
135     runway_map_iterator pos;
136
137     pos = runways.lower_bound(aptid);
138     if ( pos != runways.end() ) {
139         current = pos;
140         *r = pos->second;
141         return true;
142     } else {
143         return false;
144     }
145 }
146
147
148 // search for the specified apt id and runway no
149 bool FGRunwayList::search( const string& aptid, const string& rwyno,
150                            FGRunway *r )
151 {
152     // standardize input number
153     string runwayno = rwyno;
154     string tmp = runwayno.substr(1, 1);
155     if (( tmp == "L" || tmp == "R" || tmp == "C" ) || (runwayno.size() == 1)) {
156         tmp = runwayno;
157         runwayno = "0" + tmp;
158         SG_LOG( SG_GENERAL, SG_INFO,
159                 "Standardising rwy number from " << tmp << " to " << runwayno );
160     }
161     string revrwyno = GetReverseRunwayNo(runwayno);
162
163     runway_map_iterator pos;
164     for ( pos = runways.lower_bound( aptid );
165           pos != runways.upper_bound( aptid ); ++pos)
166     {
167         if ( pos->second.rwy_no == runwayno ) {
168             current = pos;
169             *r = pos->second;
170             return true;
171         } else if ( pos->second.rwy_no == revrwyno ) {
172             // Search again with the other-end runway number.
173             // Remember we have to munge the heading and rwy_no
174             // results if this one matches
175             current = pos;
176             *r = pos->second;
177             // NOTE - matching revrwyno implies that runwayno was
178             // actually correct.
179             r->rwy_no = runwayno;
180             r->heading += 180.0;
181             string tmp = r->end1_flags;
182             r->end1_flags = r->end2_flags;
183             r->end2_flags = tmp;
184             return true;
185         }
186     }
187
188     return false;
189 }
190
191
192 // (wierd!)
193 FGRunway FGRunwayList::search( const string& aptid ) {
194     FGRunway a;
195     search( aptid, &a );
196     return a;
197 }
198
199
200 // Return the runway closest to a given heading
201 bool FGRunwayList::search( const string& aptid, const int tgt_hdg,
202                            FGRunway *runway )
203 {
204     string rwyNo = search(aptid, tgt_hdg);
205     return(rwyNo == "NN" ? false : search(aptid, rwyNo, runway));
206 }
207
208
209 // Return the runway number of the runway closest to a given heading
210 string FGRunwayList::search( const string& aptid, const int tgt_hdg ) {
211     FGRunway r;
212     FGRunway tmp_r;     
213     string rn;
214     double found_dir = 0.0;  
215  
216     if ( !search( aptid, &tmp_r ) ) {
217         SG_LOG( SG_GENERAL, SG_ALERT,
218                 "Failed to find " << aptid << " in database." );
219         return "NN";
220     }
221     
222     double diff;
223     double min_diff = 360.0;
224     
225     while ( tmp_r.id == aptid ) {
226         r = tmp_r;
227         
228         // forward direction
229         diff = tgt_hdg - r.heading;
230         while ( diff < -180.0 ) { diff += 360.0; }
231         while ( diff >  180.0 ) { diff -= 360.0; }
232         diff = fabs(diff);
233         // SG_LOG( SG_GENERAL, SG_INFO,
234         //         "Runway " << r.rwy_no << " heading = " << r.heading <<
235         //         " diff = " << diff );
236         if ( diff < min_diff ) {
237             min_diff = diff;
238             rn = r.rwy_no;
239             found_dir = 0;
240         }
241         
242         // reverse direction
243         diff = tgt_hdg - r.heading - 180.0;
244         while ( diff < -180.0 ) { diff += 360.0; }
245         while ( diff >  180.0 ) { diff -= 360.0; }
246         diff = fabs(diff);
247         // SG_LOG( SG_GENERAL, SG_INFO,
248         //         "Runway -" << r.rwy_no << " heading = " <<
249         //         r.heading + 180.0 <<
250         //         " diff = " << diff );
251         if ( diff < min_diff ) {
252             min_diff = diff;
253             rn = r.rwy_no;
254             found_dir = 180.0;
255         }
256         
257         next( &tmp_r );
258     }
259     
260     // SG_LOG( SG_GENERAL, SG_INFO, "closest runway = " << r.rwy_no
261     //         << " + " << found_dir );
262     rn = r.rwy_no;
263     // cout << "In search, rn = " << rn << endl;
264     if ( found_dir == 180 ) {
265         rn = GetReverseRunwayNo(rn);
266         //cout << "New rn = " << rn << '\n';
267     }   
268     
269     return rn;
270 }
271
272
273 bool FGRunwayList::next( FGRunway* runway ) {
274     ++current;
275     if ( current != runways.end() ) {
276         *runway = current->second;
277         return true;
278     } else {
279         return false;
280     }
281 }
282
283
284 FGRunway FGRunwayList::next() {
285     FGRunway result;
286
287     ++current;
288     if ( current != runways.end() ) {
289         result = current->second;
290     }
291
292     return result;
293 }
294
295
296 // Destructor
297 FGRunwayList::~FGRunwayList( void ) {
298 }