1 // radiostack.cxx -- class to manage an instance of the radio stack
3 // Written by Curtis Olson, started April 2000.
5 // Copyright (C) 2000 Curtis L. Olson - curt@flightgear.org
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., 675 Mass Ave, Cambridge, MA 02139, USA.
24 #include <Aircraft/aircraft.hxx>
25 #include <Main/bfi.hxx>
26 #include <Navaids/ilslist.hxx>
27 #include <Navaids/navlist.hxx>
29 #include "radiostack.hxx"
32 * Boy, this is ugly! Make the VOR range vary by altitude difference.
35 kludgeRange (double stationElev, double aircraftElev, double nominalRange)
37 // Assume that the nominal range (usually
38 // 50nm) applies at a 5,000 ft difference.
40 double factor = ((aircraftElev*METER_TO_FEET) - stationElev) / 1000.0;
41 double range = fabs(nominalRange * factor);
43 // Clamp the range to keep it sane; for
44 // now, never less than 25% or more than
45 // 500% of nominal range.
46 if (range < nominalRange/4.0) {
47 range = nominalRange/4.0;
48 } else if (range > nominalRange*5.0) {
49 range = nominalRange*5.0;
56 FGRadioStack *current_radiostack;
60 FGRadioStack::FGRadioStack() {
67 FGRadioStack::~FGRadioStack() {
71 // Search the database for the current frequencies given current location
72 void FGRadioStack::update( double lon, double lat, double elev ) {
75 Point3D aircraft = fgGeodToCart( Point3D( lon, lat, elev ) );
80 station = Point3D( nav1_x, nav1_y, nav1_z );
81 nav1_loc_dist = aircraft.distance3D( station );
84 // staightline distance
85 station = Point3D( nav1_dme_x, nav1_dme_y, nav1_dme_z );
86 nav1_dme_dist = aircraft.distance3D( station );
92 station = Point3D( nav1_gs_x, nav1_gs_y, nav1_gs_z );
93 nav1_gs_dist = aircraft.distance3D( station );
98 if ( nav1_loc_dist < nav1_effective_range * NM_TO_METER ) {
102 geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
103 nav1_loclat, nav1_loclon,
106 // Alex: nav1_heading = - (az1 - FGBFI::getMagVar() / RAD_TO_DEG);
108 // cout << " heading = " << nav1_heading
109 // << " dist = " << nav1_dist << endl;
111 nav1_inrange = false;
116 nav1_radial = nav1_sel_radial;
119 nav1_inrange = false;
120 // cout << "not picking up vor. :-(" << endl;
124 station = Point3D( nav2_x, nav2_y, nav2_z );
125 nav2_loc_dist = aircraft.distance3D( station );
127 if ( nav2_has_dme ) {
128 // staightline distance
129 station = Point3D( nav2_dme_x, nav2_dme_y, nav2_dme_z );
130 nav2_dme_dist = aircraft.distance3D( station );
136 station = Point3D( nav2_gs_x, nav2_gs_y, nav2_gs_z );
137 nav2_gs_dist = aircraft.distance3D( station );
142 if ( nav2_loc_dist < nav2_effective_range * NM_TO_METER ) {
146 geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
147 nav2_loclat, nav2_loclon,
150 // Alex: nav2_heading = - (az1 - FGBFI::getMagVar() / RAD_TO_DEG);
152 // cout << " heading = " << nav2_heading
153 // << " dist = " << nav2_dist << endl;
155 nav2_inrange = false;
159 nav2_radial = nav2_sel_radial;
162 nav2_inrange = false;
163 // cout << "not picking up vor. :-(" << endl;
168 // staightline distance
169 station = Point3D( adf_x, adf_y, adf_z );
170 adf_dist = aircraft.distance3D( station );
172 if ( adf_dist < adf_effective_range * NM_TO_METER ) {
176 geo_inverse_wgs_84( elev, lat * RAD_TO_DEG, lon * RAD_TO_DEG,
181 // cout << " heading = " << nav2_heading
182 // << " dist = " << nav2_dist << endl;
192 // Update current nav/adf radio stations based on current postition
193 void FGRadioStack::search( double lon, double lat, double elev ) {
198 if ( current_ilslist->query( lon, lat, elev, nav1_freq, &ils ) ) {
201 nav1_has_dme = ils.get_has_dme();
202 nav1_has_gs = ils.get_has_gs();
204 nav1_loclon = ils.get_loclon();
205 nav1_loclat = ils.get_loclat();
206 nav1_gslon = ils.get_gslon();
207 nav1_gslat = ils.get_gslat();
208 nav1_dmelon = ils.get_dmelon();
209 nav1_dmelat = ils.get_dmelat();
210 nav1_elev = ils.get_gselev();
211 nav1_effective_range = FG_ILS_DEFAULT_RANGE;
212 nav1_target_gs = ils.get_gsangle();
213 nav1_radial = ils.get_locheading();
214 while ( nav1_radial < 0.0 ) { nav1_radial += 360.0; }
215 while ( nav1_radial > 360.0 ) { nav1_radial -= 360.0; }
216 nav1_x = ils.get_x();
217 nav1_y = ils.get_y();
218 nav1_z = ils.get_z();
219 nav1_gs_x = ils.get_gs_x();
220 nav1_gs_y = ils.get_gs_y();
221 nav1_gs_z = ils.get_gs_z();
222 nav1_dme_x = ils.get_dme_x();
223 nav1_dme_y = ils.get_dme_y();
224 nav1_dme_z = ils.get_dme_z();
225 // cout << "Found an ils station in range" << endl;
226 // cout << " id = " << ils.get_locident() << endl;
227 } else if ( current_navlist->query( lon, lat, elev, nav1_freq, &nav ) ) {
230 nav1_has_dme = nav.get_has_dme();
232 nav1_loclon = nav.get_lon();
233 nav1_loclat = nav.get_lat();
234 nav1_elev = nav.get_elev();
235 nav1_effective_range = kludgeRange(nav1_elev, elev, nav.get_range());
236 nav1_target_gs = 0.0;
237 nav1_radial = nav1_sel_radial;
238 nav1_x = nav1_dme_x = nav.get_x();
239 nav1_y = nav1_dme_y = nav.get_y();
240 nav1_z = nav1_dme_z = nav.get_z();
241 // cout << "Found a vor station in range" << endl;
242 // cout << " id = " << nav.get_ident() << endl;
245 // cout << "not picking up vor1. :-(" << endl;
248 if ( current_ilslist->query( lon, lat, elev, nav2_freq, &ils ) ) {
251 nav2_has_dme = ils.get_has_dme();
252 nav2_has_gs = ils.get_has_gs();
254 nav2_loclon = ils.get_loclon();
255 nav2_loclat = ils.get_loclat();
256 nav2_elev = ils.get_gselev();
257 nav2_effective_range = FG_ILS_DEFAULT_RANGE;
258 nav2_target_gs = ils.get_gsangle();
259 nav2_radial = ils.get_locheading();
260 while ( nav2_radial < 0.0 ) { nav2_radial += 360.0; }
261 while ( nav2_radial > 360.0 ) { nav2_radial -= 360.0; }
262 nav2_x = ils.get_x();
263 nav2_y = ils.get_y();
264 nav2_z = ils.get_z();
265 nav2_gs_x = ils.get_gs_x();
266 nav2_gs_y = ils.get_gs_y();
267 nav2_gs_z = ils.get_gs_z();
268 nav2_dme_x = ils.get_dme_x();
269 nav2_dme_y = ils.get_dme_y();
270 nav2_dme_z = ils.get_dme_z();
271 // cout << "Found an ils station in range" << endl;
272 // cout << " id = " << ils.get_locident() << endl;
273 } else if ( current_navlist->query( lon, lat, elev, nav2_freq, &nav ) ) {
276 nav2_has_dme = nav.get_has_dme();
277 nav2_has_dme = false;
278 nav2_loclon = nav.get_lon();
279 nav2_loclat = nav.get_lat();
280 nav2_elev = nav.get_elev();
281 nav2_effective_range = kludgeRange(nav2_elev, elev, nav.get_range());
282 nav2_target_gs = 0.0;
283 nav2_radial = nav2_sel_radial;
284 nav2_x = nav2_dme_x = nav.get_x();
285 nav2_y = nav2_dme_y = nav.get_y();
286 nav2_z = nav2_dme_z = nav.get_z();
287 // cout << "Found a vor station in range" << endl;
288 // cout << " id = " << nav.get_ident() << endl;
291 // cout << "not picking up vor2. :-(" << endl;
295 if ( current_navlist->query( lon, lat, elev, adf_freq, &nav ) ) {
298 adf_lon = nav.get_lon();
299 adf_lat = nav.get_lat();
300 adf_elev = nav.get_elev();
301 adf_effective_range = kludgeRange(adf_elev, elev, nav.get_range());
305 // cout << "Found an adf station in range" << endl;
306 // cout << " id = " << nav.get_ident() << endl;
309 // cout << "not picking up adf. :-(" << endl;
314 // periodic radio station search wrapper
315 void fgRadioSearch( void ) {
316 current_radiostack->search( cur_fdm_state->get_Longitude(),
317 cur_fdm_state->get_Latitude(),
318 cur_fdm_state->get_Altitude() * FEET_TO_METER );