1 // dme.cxx - distance-measuring equipment.
2 // Written by David Megginson, started 2003.
4 // This file is in the Public Domain and comes with no warranty.
6 #include <simgear/compiler.h>
7 #include <simgear/math/sg_geodesy.hxx>
8 #include <simgear/math/sg_random.h>
10 #include <Main/fg_props.hxx>
11 #include <Navaids/navlist.hxx>
19 * Start by calculating the radar horizon based on the elevation
20 * difference, then clamp to the maximum, then add a fudge for
21 * borderline reception.
24 adjust_range (double transmitter_elevation_ft, double aircraft_altitude_ft,
27 double delta_elevation_ft =
28 fabs(aircraft_altitude_ft - transmitter_elevation_ft);
29 double range_nm = 1.23 * sqrt(delta_elevation_ft);
30 if (range_nm > max_range_nm)
31 range_nm = max_range_nm;
32 else if (range_nm < 20.0)
34 double rand = sg_random();
35 return range_nm + (range_nm * rand * rand);
39 DME::DME ( SGPropertyNode *node )
40 : _last_distance_nm(0),
41 _last_frequency_mhz(-1),
42 _time_before_search_sec(0),
43 _transmitter_valid(false),
44 _transmitter_elevation_ft(0),
45 _transmitter_range_nm(0),
46 _transmitter_bias(0.0),
51 for ( i = 0; i < node->nChildren(); ++i ) {
52 SGPropertyNode *child = node->getChild(i);
53 string cname = child->getName();
54 string cval = child->getStringValue();
55 if ( cname == "name" ) {
57 } else if ( cname == "number" ) {
58 num = child->getIntValue();
60 SG_LOG( SG_AUTOPILOT, SG_WARN, "Error in dme config logic" );
61 if ( name.length() ) {
62 SG_LOG( SG_AUTOPILOT, SG_WARN, "Section = " << name );
69 : _last_distance_nm(0),
70 _last_frequency_mhz(-1),
71 _time_before_search_sec(0),
72 _transmitter_valid(false),
73 _transmitter_elevation_ft(0),
74 _transmitter_range_nm(0),
75 _transmitter_bias(0.0)
87 branch = "/instrumentation/" + name;
89 SGPropertyNode *node = fgGetNode(branch.c_str(), num, true );
91 _longitude_node = fgGetNode("/position/longitude-deg", true);
92 _latitude_node = fgGetNode("/position/latitude-deg", true);
93 _altitude_node = fgGetNode("/position/altitude-ft", true);
94 _serviceable_node = node->getChild("serviceable", 0, true);
95 _electrical_node = fgGetNode("/systems/electrical/outputs/dme", true);
96 SGPropertyNode *fnode = node->getChild("frequencies", 0, true);
97 _source_node = fnode->getChild("source", 0, true);
98 _frequency_node = fnode->getChild("selected-mhz", 0, true);
99 _in_range_node = node->getChild("in-range", 0, true);
100 _distance_node = node->getChild("indicated-distance-nm", 0, true);
101 _speed_node = node->getChild("indicated-ground-speed", 0, true);
102 _time_node = node->getChild("indicated-time-min", 0, true);
104 _serviceable_node->setBoolValue(true);
108 DME::update (double delta_time_sec)
110 // Figure out the source
111 const char * source = _source_node->getStringValue();
112 if (source[0] == '\0') {
114 branch = "/instrumentation/" + name + "/frequencies/selected-mhz";
115 _source_node->setStringValue(branch.c_str());
116 source = _source_node->getStringValue();
120 double frequency_mhz = fgGetDouble(source, 108.0);
121 if (frequency_mhz != _last_frequency_mhz) {
122 _time_before_search_sec = 0;
123 _last_frequency_mhz = frequency_mhz;
126 // Get the aircraft position
127 double longitude_rad =
128 _longitude_node->getDoubleValue() * SGD_DEGREES_TO_RADIANS;
129 double latitude_rad =
130 _latitude_node->getDoubleValue() * SGD_DEGREES_TO_RADIANS;
132 _altitude_node->getDoubleValue() * SG_FEET_TO_METER;
134 // On timeout, scan again
135 _time_before_search_sec -= delta_time_sec;
136 if (_time_before_search_sec < 0)
137 search(frequency_mhz, longitude_rad,
138 latitude_rad, altitude_m);
140 // If it's off, don't bother.
141 if (!_serviceable_node->getBoolValue() ||
142 !_electrical_node->getBoolValue() ||
143 !_transmitter_valid) {
144 _last_distance_nm = 0;
145 _in_range_node->setBoolValue(false);
146 _distance_node->setDoubleValue(0);
147 _speed_node->setDoubleValue(0);
148 _time_node->setDoubleValue(0);
152 // Calculate the distance to the transmitter
154 sgGeodToCart(Point3D(longitude_rad, latitude_rad, altitude_m));
155 double distance_nm = _transmitter.distance3D(location) * SG_METER_TO_NM;
157 double range_nm = adjust_range(_transmitter_elevation_ft,
158 altitude_m * SG_METER_TO_FEET,
159 _transmitter_range_nm);
161 if (distance_nm <= range_nm) {
162 double speed_kt = (fabs(distance_nm - _last_distance_nm) *
163 ((1 / delta_time_sec) * 3600.0));
164 _last_distance_nm = distance_nm;
166 _in_range_node->setBoolValue(true);
167 double tmp_dist = distance_nm - _transmitter_bias;
168 if ( tmp_dist < 0.0 ) {
171 _distance_node->setDoubleValue( tmp_dist );
172 _speed_node->setDoubleValue(speed_kt);
173 _time_node->setDoubleValue(distance_nm/speed_kt*60.0);
176 _last_distance_nm = 0;
177 _in_range_node->setBoolValue(false);
178 _distance_node->setDoubleValue(0);
179 _speed_node->setDoubleValue(0);
180 _time_node->setDoubleValue(0);
185 DME::search (double frequency_mhz, double longitude_rad,
186 double latitude_rad, double altitude_m)
189 _time_before_search_sec = 1.0;
191 // try the ILS list first
193 = globals->get_dmelist()->findByFreq( frequency_mhz, longitude_rad,
194 latitude_rad, altitude_m);
196 _transmitter_valid = (dme != NULL);
198 if ( _transmitter_valid ) {
199 _transmitter = Point3D(dme->get_x(), dme->get_y(), dme->get_z());
200 _transmitter_elevation_ft = dme->get_elev_ft();
201 _transmitter_range_nm = dme->get_range();
202 _transmitter_bias = dme->get_multiuse();