+
+ void parseRunwayLine850(const vector<string>& token)
+ {
+ double width = atof( token[1].c_str() ) * SG_METER_TO_FEET;
+ int surface_code = atoi( token[2].c_str() );
+
+ double lat_1 = atof( token[9].c_str() );
+ double lon_1 = atof( token[10].c_str() );
+ SGGeod pos_1(SGGeod::fromDegFt(lon_1, lat_1, 0.0));
+ rwy_lat_accum += lat_1;
+ rwy_lon_accum += lon_1;
+ rwy_count++;
+
+ double lat_2 = atof( token[18].c_str() );
+ double lon_2 = atof( token[19].c_str() );
+ SGGeod pos_2(SGGeod::fromDegFt(lon_2, lat_2, 0.0));
+ rwy_lat_accum += lat_2;
+ rwy_lon_accum += lon_2;
+ rwy_count++;
+
+ double length, heading_1, heading_2, dummy;
+ SGGeodesy::inverse( pos_1, pos_2, heading_1, heading_2, length );
+ SGGeod pos;
+ SGGeodesy::direct( pos_1, heading_1, length / 2.0, pos, dummy );
+ length *= SG_METER_TO_FEET;
+
+ last_rwy_heading = heading_1;
+
+ const string& rwy_no_1(token[8]);
+ const string& rwy_no_2(token[17]);
+ if ( rwy_no_1.size() == 0 || rwy_no_2.size() == 0 )
+ return;
+
+ double displ_thresh1 = atof( token[11].c_str() );
+ double displ_thresh2 = atof( token[20].c_str() );
+
+ double stopway1 = atof( token[12].c_str() );
+ double stopway2 = atof( token[21].c_str() );
+
+ FGRunway* rwy = new FGRunway(NULL, rwy_no_1, pos, heading_1, length,
+ width, displ_thresh1, stopway1, surface_code, false);
+ runways.push_back(rwy);
+
+ FGRunway* reciprocal = new FGRunway(NULL, rwy_no_2,
+ pos, heading_2, length, width,
+ displ_thresh2, stopway2, surface_code, true);
+ runways.push_back(reciprocal);
+
+ rwy->setReciprocalRunway(reciprocal);
+ reciprocal->setReciprocalRunway(rwy);
+ }
+
+ void parseWaterRunwayLine850(const vector<string>& token)
+ {
+ double width = atof( token[1].c_str() ) * SG_METER_TO_FEET;
+
+ double lat_1 = atof( token[4].c_str() );
+ double lon_1 = atof( token[5].c_str() );
+ SGGeod pos_1(SGGeod::fromDegFt(lon_1, lat_1, 0.0));
+ rwy_lat_accum += lat_1;
+ rwy_lon_accum += lon_1;
+ rwy_count++;
+
+ double lat_2 = atof( token[7].c_str() );
+ double lon_2 = atof( token[8].c_str() );
+ SGGeod pos_2(SGGeod::fromDegFt(lon_2, lat_2, 0.0));
+ rwy_lat_accum += lat_2;
+ rwy_lon_accum += lon_2;
+ rwy_count++;
+
+ double length, heading_1, heading_2, dummy;
+ SGGeodesy::inverse( pos_1, pos_2, heading_1, heading_2, length );
+ SGGeod pos;
+ SGGeodesy::direct( pos_1, heading_1, length / 2.0, pos, dummy );
+
+ last_rwy_heading = heading_1;
+
+ const string& rwy_no_1(token[3]);
+ const string& rwy_no_2(token[6]);
+
+ FGRunway* rwy = new FGRunway(NULL, rwy_no_1, pos, heading_1, length,
+ width, 0.0, 0.0, 13, false);
+ runways.push_back(rwy);
+
+ FGRunway* reciprocal = new FGRunway(NULL, rwy_no_2,
+ pos, heading_2, length, width,
+ 0.0, 0.0, 13, true);
+ runways.push_back(reciprocal);
+
+ rwy->setReciprocalRunway(reciprocal);
+ reciprocal->setReciprocalRunway(rwy);
+ }
+
+ void parseHelipadLine850(const vector<string>& token)
+ {
+ double length = atof( token[5].c_str() ) * SG_METER_TO_FEET;
+ double width = atof( token[6].c_str() ) * SG_METER_TO_FEET;
+
+ double lat = atof( token[2].c_str() );
+ double lon = atof( token[3].c_str() );
+ SGGeod pos(SGGeod::fromDegFt(lon, lat, 0.0));
+ rwy_lat_accum += lat;
+ rwy_lon_accum += lon;
+ rwy_count++;
+
+ double heading = atof( token[4].c_str() );
+
+ last_rwy_heading = heading;
+
+ const string& rwy_no(token[1]);
+ int surface_code = atoi( token[7].c_str() );
+
+ FGRunway* rwy = new FGRunway(NULL, rwy_no, pos, heading, length,
+ width, 0.0, 0.0, surface_code, false);
+ runways.push_back(rwy);
+ }
+
+ void parsePavementLine850(const vector<string>& token)
+ {
+ if ( token.size() >= 5 ) {
+ pavement_ident = token[4];
+ if ( !pavement_ident.empty() && pavement_ident[pavement_ident.size()-1] == '\r' )
+ pavement_ident.erase( pavement_ident.size()-1 );
+ } else {
+ pavement_ident = "xx";
+ }
+ }
+
+ void parsePavementNodeLine850(int num, const vector<string>& token)
+ {
+ double lat = atof( token[1].c_str() );
+ double lon = atof( token[2].c_str() );
+ SGGeod pos(SGGeod::fromDegFt(lon, lat, 0.0));
+
+ FGPavement* pvt = 0;
+ if ( !pavement_ident.empty() ) {
+ pvt = new FGPavement( pavement_ident, pos );
+ pavements.push_back( pvt );
+ pavement_ident = "";
+ } else {
+ pvt = pavements.back();
+ }
+ if ( num == 112 || num == 114 ) {
+ double lat_b = atof( token[3].c_str() );
+ double lon_b = atof( token[4].c_str() );
+ SGGeod pos_b(SGGeod::fromDegFt(lon_b, lat_b, 0.0));
+ pvt->addBezierNode(pos, pos_b, num == 114);
+ } else {
+ pvt->addNode(pos, num == 113);
+ }
+ }
+
+ void parseATISLine(FGCommList *comm_list, const vector<string>& token)
+ {
+ if ( rwy_count <= 0 ) {
+ SG_LOG( SG_GENERAL, SG_ALERT,
+ "No runways; skipping AWOS for " + last_apt_id);
+ }
+
+// This assumes/requires that any code-50 line (ATIS or AWOS)
+ // applies to the preceding code-1 line (airport ID and name)
+ // and that a full set of code-10 lines (runway descriptors)
+ // has come between the code-1 and code-50 lines.
+ // typical code-50 lines:
+ // 50 11770 ATIS
+ // 50 11770 AWOS 3
+ // This code parallels code found in "operator>>" in ATC.hxx;
+ // FIXME: unify the code.
+#if ENABLE_ATCDCL
+ ATCData a;
+ a.geod = SGGeod::fromDegFt(rwy_lon_accum / (double)rwy_count,
+ rwy_lat_accum / (double)rwy_count, last_apt_elev);
+ a.range = 50; // give all ATISs small range
+ a.ident = last_apt_id;
+ a.name = last_apt_name;
+ // short int representing tens of kHz:
+ a.freq = atoi(token[1].c_str());
+ if (token[2] == "ATIS") a.type = ATIS;
+ else a.type = AWOS; // ASOS same as AWOS
+
+ // generate cartesian coordinates
+ a.cart = SGVec3d::fromGeod(a.geod);
+ comm_list->commlist_freq[a.freq].push_back(a);
+
+ SGBucket bucket(a.geod);
+ int bucknum = bucket.gen_index();
+ comm_list->commlist_bck[bucknum].push_back(a);
+#else
+#endif
+#if 0
+ SG_LOG( SG_GENERAL, SG_ALERT,
+ "Loaded ATIS/AWOS for airport: " << a.ident
+ << " lat: " << a.geod.getLatitudeDeg()
+ << " lon: " << a.geod.getLongitudeDeg()
+ << " freq: " << a.freq
+ << " type: " << a.type );
+#endif
+ }
+