]> git.mxchange.org Git - flightgear.git/blob - src/ATCDCL/tower.cxx
52d50f9a83c6ddc30656c0955e11cd260ae8b0f6
[flightgear.git] / src / ATCDCL / tower.cxx
1 // FGTower - a class to provide tower control at towered airports.
2 //
3 // Written by David Luff, started March 2002.
4 //
5 // Copyright (C) 2002  David C. Luff - david.luff@nottingham.ac.uk
6 // Copyright (C) 2008  Daniyar Atadjanov (ground clearance, gear check, weather, etc.)
7 //
8 // This program is free software; you can redistribute it and/or
9 // modify it under the terms of the GNU General Public License as
10 // published by the Free Software Foundation; either version 2 of the
11 // License, or (at your option) any later version.
12 //
13 // This program is distributed in the hope that it will be useful, but
14 // WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
16 // General Public License for more details.
17 //
18 // You should have received a copy of the GNU General Public License
19 // along with this program; if not, write to the Free Software
20 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
21
22 #ifdef HAVE_CONFIG_H
23 #  include <config.h>
24 #endif
25
26 #ifdef HAVE_STRINGS_H
27 #  include <strings.h>   // bcopy()
28 #else
29 #  include <string.h>    // MSVC doesn't have strings.h
30 #endif
31
32 #include <sstream>
33 #include <iomanip>
34 #include <iostream>
35
36 #include <simgear/debug/logstream.hxx>
37 #include <simgear/math/sg_geodesy.hxx>
38 #include <simgear/math/sg_random.h>
39 #include <simgear/misc/sg_path.hxx>
40
41 #include <Main/globals.hxx>
42 #include <Airports/runways.hxx>
43
44 #include "tower.hxx"
45 #include "ATCmgr.hxx"
46 #include "ATCutils.hxx"
47 #include "ATCDialog.hxx"
48 #include "commlist.hxx"
49 #include "AILocalTraffic.hxx"
50
51
52 using std::cout;
53
54 // TowerPlaneRec
55
56 TowerPlaneRec::TowerPlaneRec() :
57         planePtr(NULL),
58         clearedToLand(false),
59         clearedToLineUp(false),
60         clearedToTakeOff(false),
61         holdShortReported(false),
62         lineUpReported(false),
63         downwindReported(false),
64         longFinalReported(false),
65         longFinalAcknowledged(false),
66         finalReported(false),
67         finalAcknowledged(false),
68         rwyVacatedReported(false),
69         rwyVacatedAcknowledged(false),
70         goAroundReported(false),
71         instructedToGoAround(false),
72         onRwy(false),
73         nextOnRwy(false),
74         vfrArrivalReported(false),
75         vfrArrivalAcknowledged(false),
76         opType(TTT_UNKNOWN),
77         leg(LEG_UNKNOWN),
78         landingType(AIP_LT_UNKNOWN),
79         gearWasUp(false),
80         gearUpReported(false),
81         isUser(false)
82 {
83         plane.callsign = "UNKNOWN";
84 }
85
86 TowerPlaneRec::TowerPlaneRec(const PlaneRec& p) :
87         planePtr(NULL),
88         clearedToLand(false),
89         clearedToLineUp(false),
90         clearedToTakeOff(false),
91         holdShortReported(false),
92         lineUpReported(false),
93         downwindReported(false),
94         longFinalReported(false),
95         longFinalAcknowledged(false),
96         finalReported(false),
97         finalAcknowledged(false),
98         rwyVacatedReported(false),
99         rwyVacatedAcknowledged(false),
100         goAroundReported(false),
101         instructedToGoAround(false),
102         onRwy(false),
103         nextOnRwy(false),
104         vfrArrivalReported(false),
105         vfrArrivalAcknowledged(false),
106         opType(TTT_UNKNOWN),
107         leg(LEG_UNKNOWN),
108         landingType(AIP_LT_UNKNOWN),
109         gearWasUp(false),
110         gearUpReported(false),
111         isUser(false)
112 {
113         plane = p;
114 }
115
116 TowerPlaneRec::TowerPlaneRec(const Point3D& pt) :
117         planePtr(NULL),
118         clearedToLand(false),
119         clearedToLineUp(false),
120         clearedToTakeOff(false),
121         holdShortReported(false),
122         lineUpReported(false),
123         downwindReported(false),
124         longFinalReported(false),
125         longFinalAcknowledged(false),
126         finalReported(false),
127         finalAcknowledged(false),
128         rwyVacatedReported(false),
129         rwyVacatedAcknowledged(false),
130         goAroundReported(false),
131         instructedToGoAround(false),
132         onRwy(false),
133         nextOnRwy(false),
134         vfrArrivalReported(false),
135         vfrArrivalAcknowledged(false),
136         opType(TTT_UNKNOWN),
137         leg(LEG_UNKNOWN),
138         landingType(AIP_LT_UNKNOWN),
139         gearWasUp(false),
140         gearUpReported(false),
141         isUser(false)
142 {
143         plane.callsign = "UNKNOWN";
144         pos = pt;
145 }
146
147 TowerPlaneRec::TowerPlaneRec(const PlaneRec& p, const Point3D& pt) :
148         planePtr(NULL),
149         clearedToLand(false),
150         clearedToLineUp(false),
151         clearedToTakeOff(false),
152         holdShortReported(false),
153         lineUpReported(false),
154         downwindReported(false),
155         longFinalReported(false),
156         longFinalAcknowledged(false),
157         finalReported(false),
158         finalAcknowledged(false),
159         rwyVacatedReported(false),
160         rwyVacatedAcknowledged(false),
161         goAroundReported(false),
162         instructedToGoAround(false),
163         onRwy(false),
164         nextOnRwy(false),
165         vfrArrivalReported(false),
166         vfrArrivalAcknowledged(false),
167         opType(TTT_UNKNOWN),
168         leg(LEG_UNKNOWN),
169         landingType(AIP_LT_UNKNOWN),
170         gearWasUp(false),
171         gearUpReported(false),
172         isUser(false)
173 {
174         plane = p;
175         pos = pt;
176 }
177
178
179 // FGTower
180
181 /*******************************************
182                TODO List
183                            
184 Currently user is assumed to have taken off again when leaving the runway - check speed/elev for taxiing-in. (MAJOR)
185
186 Use track instead of heading to determine what leg of the circuit the user is flying. (MINOR)
187
188 Use altitude as well as position to try to determine if the user has left the circuit. (MEDIUM - other issues as well).
189
190 Currently HoldShortReported code assumes there will be only one plane holding for the runway at once and 
191 will break when planes start queueing. (CRITICAL)
192
193 Report-Runway-Vacated is left as only user ATC option following a go-around. (MAJOR)
194
195 Report-Downwind is not added as ATC option when user takes off to fly a circuit. (MAJOR)
196
197 eta of USER can be calculated very wrongly in circuit if flying straight out and turn4 etc are with +ve ortho y. 
198 This can then screw up circuit ordering for other planes (MEDIUM)
199
200 USER leaving circuit needs to be more robustly considered when intentions unknown
201 Currently only considered during climbout and breaks when user turns (MEDIUM).
202
203 GetPos() of the AI planes is called erratically - either too much or not enough. (MINOR)
204
205 GO-AROUND is instructed very late at < 12s to landing - possibly make more dependent on chance of rwy clearing before landing (FEATURE)
206
207 Need to make clear when TowerPlaneRecs do or don't get deleted in RemoveFromCircuitList etc. (MINOR until I misuse it - then CRITICAL!)
208
209 FGTower::RemoveAllUserDialogOptions() really ought to be replaced by an ATCDialog::clear_entries() function. (MINOR - efficiency).
210
211 At the moment planes in the lists are not guaranteed to always have a sensible ETA - it should be set as part of AddList functions, and lists should only be accessed this way. (FAIRLY MAJOR). 
212 *******************************************/
213
214 FGTower::FGTower() :
215         separateGround(true),
216         ground(0)
217 {
218         ATCmgr = globals->get_ATC_mgr();
219         
220         _type = TOWER;
221         
222         // Init the property nodes - TODO - need to make sure we're getting surface winds.
223         wind_from_hdg = fgGetNode("/environment/wind-from-heading-deg", true);
224         wind_speed_knots = fgGetNode("/environment/wind-speed-kt", true);
225         
226         update_count = 0;
227         update_count_max = 15;
228         
229         holdListItr = holdList.begin();
230         appListItr = appList.begin();
231         depListItr = depList.begin();
232         rwyListItr = rwyList.begin();
233         circuitListItr = circuitList.begin();
234         trafficListItr = trafficList.begin();
235         vacatedListItr = vacatedList.begin();
236         
237         freqClear = true;
238         
239         timeSinceLastDeparture = 9999;
240         departed = false;
241         
242         nominal_downwind_leg_pos = 1000.0;
243         nominal_base_leg_pos = -1000.0;
244         // TODO - set nominal crosswind leg pos based on minimum distance from takeoff end of rwy.
245         
246         _departureControlled = false;
247 }
248
249 FGTower::~FGTower() {
250         if(!separateGround) {
251                 delete ground;
252         }
253 }
254
255 void FGTower::Init() {
256         //cout << "Initialising tower " << ident << '\n';
257         
258         // Pointers to user's position
259         user_lon_node = fgGetNode("/position/longitude-deg", true);
260         user_lat_node = fgGetNode("/position/latitude-deg", true);
261         user_elev_node = fgGetNode("/position/altitude-ft", true);
262         user_hdg_node = fgGetNode("/orientation/heading-deg", true);
263         
264         // Need some way to initialise rwyOccupied flag correctly if the user is on the runway and to know its the user.
265         // I'll punt the startup issue for now though!!!
266         rwyOccupied = false;
267         
268         // Setup the ground control at this airport
269         AirportATC a;
270         //cout << "Tower ident = " << ident << '\n';
271         if(ATCmgr->GetAirportATCDetails(ident, &a)) {
272                 if(a.ground_freq) {             // Ground control
273                         ground = (FGGround*)ATCmgr->GetATCPointer(ident, GROUND);
274                         separateGround = true;
275                         if(ground == NULL) {
276                                 // Something has gone wrong :-(
277                                 SG_LOG(SG_ATC, SG_WARN, "ERROR - ground has frequency but can't get ground pointer :-(");
278                                 ground = new FGGround(ident);
279                                 separateGround = false;
280                                 ground->Init();
281                                 if(_display) {
282                                         ground->SetDisplay();
283                                 } else {
284                                         ground->SetNoDisplay();
285                                 }
286                         }
287                 } else {
288                         // Initialise ground anyway to do the shortest path stuff!
289                         // Note that we're now responsible for updating and deleting this - NOT the ATCMgr.
290                         ground = new FGGround(ident);
291                         separateGround = false;
292                         ground->Init();
293                         if(_display) {
294                                 ground->SetDisplay();
295                         } else {
296                                 ground->SetNoDisplay();
297                         }
298                 }
299         } else {
300                 SG_LOG(SG_ATC, SG_ALERT, "Unable to find airport details for " << ident << " in FGTower::Init()");
301                 // Initialise ground anyway to avoid segfault later
302                 ground = new FGGround(ident);
303                 separateGround = false;
304                 ground->Init();
305                 if(_display) {
306                         ground->SetDisplay();
307                 } else {
308                         ground->SetNoDisplay();
309                 }
310         }
311         
312         RemoveAllUserDialogOptions();
313         
314         // TODO - attempt to get a departure control pointer to see if we need to hand off departing traffic to departure.
315         
316         // Get the airport elevation
317         aptElev = fgGetAirportElev(ident.c_str());
318         
319         // TODO - this function only assumes one active rwy.
320         DoRwyDetails();
321         
322         // TODO - this currently assumes only one active runway.
323         rwyOccupied = OnActiveRunway(Point3D(user_lon_node->getDoubleValue(), user_lat_node->getDoubleValue(), 0.0));
324         
325         if(!OnAnyRunway(Point3D(user_lon_node->getDoubleValue(), user_lat_node->getDoubleValue(), 0.0), false)) {
326                 //cout << ident << "  ADD 0\n";
327                 current_atcdialog->add_entry(ident, "@AP Tower, @CS @MI miles @CD of the airport for full stop@AT",
328                                 "Contact tower for VFR arrival (full stop)", TOWER,
329                                 (int)USER_REQUEST_VFR_ARRIVAL_FULL_STOP);
330         } else {
331                 //cout << "User found on active runway\n";
332                 // Assume the user is started at the threshold ready to take-off
333                 TowerPlaneRec* t = new TowerPlaneRec;
334                 t->plane.callsign = fgGetString("/sim/user/callsign");
335                 t->plane.type = GA_SINGLE;      // FIXME - hardwired!!
336                 t->opType = TTT_UNKNOWN;        // We don't know if the user wants to do circuits or a departure...
337                 t->landingType = AIP_LT_UNKNOWN;
338                 t->leg = TAKEOFF_ROLL;
339                 t->isUser = true;
340                 t->planePtr = NULL;
341                 t->clearedToTakeOff = false;
342                 rwyList.push_back(t);
343                 rwyListItr = rwyList.begin();
344                 departed = false;
345                 current_atcdialog->add_entry(ident, "@CS @TO", "Request departure / take-off clearance",
346                                 TOWER, (int)USER_REQUEST_TAKE_OFF);
347         }
348 }
349
350 void FGTower::Update(double dt) {
351         //cout << "T" << endl;
352         // Each time step, what do we need to do?
353         // We need to go through the list of outstanding requests and acknowedgements
354         // and process at least one of them.
355         // We need to go through the list of planes under our control and check if
356         // any need to be addressed.
357         // We need to check for planes not under our control coming within our 
358         // control area and address if necessary.
359         
360         // TODO - a lot of the below probably doesn't need to be called every frame and should be staggered.
361         
362         // Sort the arriving planes
363         
364         /*
365         if(ident == "KEMT") {
366                 cout << update_count << "\ttL: " << trafficList.size() << "  cL: " << circuitList.size() << "  hL: " << holdList.size() << "  aL: " << appList.size() << '\n';
367         }
368         */
369         //if(ident == "EGNX") cout << display << '\n';
370         
371         if(departed != false) {
372                 timeSinceLastDeparture += dt;
373                 //if(ident == "KEMT") 
374                 //      cout << "  dt = " << dt << "  timeSinceLastDeparture = " << timeSinceLastDeparture << '\n';
375         }
376         
377         //cout << ident << " respond = " << respond << " responseReqd = " << responseReqd << '\n'; 
378         if(respond) {
379                 if(!responseReqd) SG_LOG(SG_ATC, SG_ALERT, "ERROR - respond is true and responseReqd is false in FGTower::Update(...)");
380                 Respond();
381                 respond = false;
382                 responseReqd = false;
383         }
384         
385         // Calculate the eta of each plane to the threshold.
386         // For ground traffic this is the fastest they can get there.
387         // For air traffic this is the middle approximation.
388         if(update_count == 1) {
389                 doThresholdETACalc();
390         }
391         
392         // Order the list of traffic as per expected threshold use and flag any conflicts
393         if(update_count == 2) {
394                 //bool conflicts = doThresholdUseOrder();
395                 doThresholdUseOrder();
396         }
397         
398         // sortConficts() !!!
399         
400         if(update_count == 4) {
401                 CheckHoldList(dt);
402         }
403         
404         // Uggh - HACK - why have we got rwyOccupied - wouldn't simply testing rwyList.size() do?
405         if(rwyList.size()) {
406                 rwyOccupied = true;
407         } else {
408                 rwyOccupied = false;
409         }
410         
411         if(update_count == 5 && rwyOccupied) {
412                 CheckRunwayList(dt);
413         }
414                 
415         if(update_count == 6) {
416                 CheckCircuitList(dt);
417         }
418         
419         if(update_count == 7) {
420                 CheckApproachList(dt);
421         }
422         
423         if(update_count == 8) {
424                 CheckDepartureList(dt);
425         }
426         
427         // TODO - do one plane from the departure list and set departed = false when out of consideration
428         
429         //doCommunication();
430         
431         if(!separateGround) {
432                 // The display stuff might have to get more clever than this when not separate 
433                 // since the tower and ground might try communicating simultaneously even though
434                 // they're mean't to be the same contoller/frequency!!
435                 // We could also get rid of this by overloading FGATC's Set(No)Display() functions.
436                 if(_display) {
437                         ground->SetDisplay();
438                 } else {
439                         ground->SetNoDisplay();
440                 }
441                 ground->Update(dt);
442         }
443         
444         ++update_count;
445         // How big should ii get - ie how long should the update cycle interval stretch?
446         if(update_count >= update_count_max) {
447                 update_count = 0;
448         }
449         
450         // Call the base class update for the response time handling.
451         FGATC::Update(dt);
452
453         /*
454         if(ident == "KEMT") {   
455                 // For AI debugging convienience - may be removed
456                 Point3D user_pos;
457                 user_pos.setlon(user_lon_node->getDoubleValue());
458                 user_pos.setlat(user_lat_node->getDoubleValue());
459                 user_pos.setelev(user_elev_node->getDoubleValue());
460                 Point3D user_ortho_pos = ortho.ConvertToLocal(user_pos);
461                 fgSetDouble("/AI/user/ortho-x", user_ortho_pos.x());
462                 fgSetDouble("/AI/user/ortho-y", user_ortho_pos.y());
463                 fgSetDouble("/AI/user/elev", user_elev_node->getDoubleValue());
464         }
465         */
466         
467         //cout << "Done T" << endl;
468 }
469
470 void FGTower::ReceiveUserCallback(int code) {
471         if(code == (int)USER_REQUEST_VFR_DEPARTURE) {
472                 RequestDepartureClearance("USER");
473         } else if(code == (int)USER_REQUEST_VFR_ARRIVAL) {
474                 VFRArrivalContact("USER");
475         } else if(code == (int)USER_REQUEST_VFR_ARRIVAL_FULL_STOP) {
476                 VFRArrivalContact("USER", FULL_STOP);
477         } else if(code == (int)USER_REQUEST_VFR_ARRIVAL_TOUCH_AND_GO) {
478                 VFRArrivalContact("USER", TOUCH_AND_GO);
479         } else if(code == (int)USER_REPORT_DOWNWIND) {
480                 ReportDownwind("USER");
481         } else if(code == (int)USER_REPORT_3_MILE_FINAL) {
482                 // For now we'll just call report final instead of report long final to avoid having to alter the response code
483                 ReportFinal("USER");
484         } else if(code == (int)USER_REPORT_RWY_VACATED) {
485                 ReportRunwayVacated("USER");
486         } else if(code == (int)USER_REPORT_GOING_AROUND) {
487                 ReportGoingAround("USER");
488         } else if(code == (int)USER_REQUEST_TAKE_OFF) {
489                 RequestTakeOffClearance("USER");
490         }
491 }
492
493 // **************** RESPONSE FUNCTIONS ****************
494
495 void FGTower::Respond() {
496         //cout << "\nEntering Respond, responseID = " << responseID << endl;
497         TowerPlaneRec* t = FindPlane(responseID);
498         if(t) {
499                 // This will grow!!!
500                 if(t->vfrArrivalReported && !t->vfrArrivalAcknowledged) {
501                         //cout << "Tower " << ident << " is responding to VFR arrival reported...\n";
502                         // Testing - hardwire straight in for now
503                         string trns = t->plane.callsign;
504                         trns += " ";
505                         trns += name;
506                         trns += " Tower,";
507                         // Should we clear staight in or for downwind entry?
508                         // For now we'll clear straight in if greater than 1km from a line drawn through the threshold perpendicular to the rwy.
509                         // Later on we might check the actual heading and direct some of those to enter on downwind or base.
510                         Point3D op = ortho.ConvertToLocal(t->pos);
511                         float gp = fgGetFloat("/gear/gear/position-norm");
512                         if(gp < 1)
513                                 t->gearWasUp = true; // This will be needed on final to tell "Gear down, ready to land."
514                         if(op.y() < -1000) {
515                                 trns += " Report three mile straight-in runway ";
516                                 t->opType = STRAIGHT_IN;
517                                 if(t->isUser) {
518                                         current_atcdialog->add_entry(ident, "@CS @MI mile final runway @RW@GR", "Report Final", TOWER, (int)USER_REPORT_3_MILE_FINAL);
519                                 } else {
520                                         t->planePtr->RegisterTransmission(14);
521                                 }
522                         } else {
523                                 // For now we'll just request reporting downwind.
524                                 // TODO - In real life something like 'report 2 miles southwest right downwind rwy 19R' might be used
525                                 // but I'm not sure how to handle all permutations of which direction to tell to report from yet.
526                                 trns += " Report ";
527                                 //cout << "Responding, rwy.patterDirection is " << rwy.patternDirection << '\n';
528                                 trns += ((rwy.patternDirection == 1) ? "right " : "left ");
529                                 trns += "downwind runway ";
530                                 t->opType = CIRCUIT;
531                                 // leave it in the app list until it gets into pattern though.
532                                 if(t->isUser) {
533                                         current_atcdialog->add_entry(ident, "@AP Tower, @CS Downwind @RW", "Report Downwind", TOWER, (int)USER_REPORT_DOWNWIND);
534                                 } else {
535                                         t->planePtr->RegisterTransmission(15);
536                                 }
537                         }
538                         trns += ConvertRwyNumToSpokenString(activeRwy);
539                         if(_display) {
540                                 pending_transmission = trns;
541                                 Transmit();
542                         } else {
543                                 //cout << "Not displaying, trns was " << trns << '\n';
544                         }
545                         t->vfrArrivalAcknowledged = true;
546                 } else if(t->downwindReported) {
547                         //cout << "Tower " << ident << " is responding to downwind reported...\n";
548                         ProcessDownwindReport(t);
549                         t->downwindReported = false;
550                 } else if(t->lineUpReported) {
551                         string trns = t->plane.callsign;
552                         if(rwyOccupied) {
553                                 double f = globals->get_ATC_mgr()->GetFrequency(ident, ATIS) / 100.0;
554                                 string wtr;
555                                 if(!f) {
556                                         wtr = ", " + GetWeather();
557                                 }
558                                 trns += " Cleared for take-off" + wtr;
559                                 t->clearedToTakeOff = true;
560                         } else {
561                                 if(!OnAnyRunway(Point3D(user_lon_node->getDoubleValue(), user_lat_node->getDoubleValue(), 0.0), true)) {
562                                         // TODO: Check if any AI Planes on final and tell something like: "After the landing CALLSIGN line up runway two eight right"
563                                         trns += " Line up runway " + ConvertRwyNumToSpokenString(activeRwy);
564                                         t->clearedToTakeOff = false;
565                                         current_atcdialog->add_entry(ident, "@CS @TO", "Report ready for take-off", TOWER, (int)USER_REQUEST_TAKE_OFF);
566
567                                 } else {
568                                         sg_srandom_time();
569                                         if((int(sg_random() * 10) + 1) != 3) {
570                                                 t->clearedToTakeOff = true;
571                                                 trns += " Cleared immediate take-off ";
572                                         } else {
573                                                 t->clearedToTakeOff = false;
574                                                 trns += " Negative, departure runway " + ConvertRwyNumToSpokenString(activeRwy);
575                                         }
576                                 }
577                         }
578                         if(_display) {
579                                 pending_transmission = trns;
580                                 Transmit();
581                         } else {
582                                 //cout << "Not displaying, trns was " << trns << '\n';
583                         }
584                         t->lineUpReported = false;
585                 } else if(t->holdShortReported) {
586                         //cout << "Tower " << ident << " is reponding to holdShortReported...\n";
587                         if(t->nextOnRwy) {
588                                 if(rwyOccupied) {       // TODO - ought to add a sanity check that it isn't this plane only on the runway (even though it shouldn't be!!)
589                                         // Do nothing for now - consider acknowloging hold short eventually
590                                 } else {
591                                         ClearHoldingPlane(t);
592                                         t->leg = TAKEOFF_ROLL;
593                                         rwyList.push_back(t);
594                                         rwyListItr = rwyList.begin();
595                                         rwyOccupied = true;
596                                         // WARNING - WE ARE ASSUMING ONLY ONE PLANE REPORTING HOLD AT A TIME BELOW
597                                         // FIXME TODO - FIX THIS!!!
598                                         if(!holdList.empty()) {
599                                                 if(holdListItr == holdList.end()) {
600                                                         holdListItr = holdList.begin();
601                                                 }
602                                                 holdList.erase(holdListItr);
603                                                 holdListItr = holdList.begin();
604                                         }
605                                 }
606                         } else {
607                                 // Tell him to hold and what position he is.
608                                 // Not currently sure under which circumstances we do or don't bother transmitting this.
609                                 string trns = t->plane.callsign;
610                                 trns += " hold position";
611                                 if(_display) {
612                                         pending_transmission = trns;
613                                         Transmit();
614                                 }
615                                 // TODO - add some idea of what traffic is blocking him.
616                         }
617                         t->holdShortReported = false;
618                 } else if(t->finalReported && !(t->finalAcknowledged)) {
619                         //cout << "Tower " << ident << " is responding to finalReported...\n";
620                         bool disp = true;
621                         string trns = t->plane.callsign;
622                         //cout << (t->nextOnRwy ? "Next on rwy " : "Not next!! ");
623                         //cout << (rwyOccupied ? "RWY OCCUPIED!!\n" : "Rwy not occupied\n");
624                         if(t->nextOnRwy && !rwyOccupied && !(t->instructedToGoAround)) {
625                                 if(t->landingType == FULL_STOP) {
626                                         trns += " cleared to land ";
627                                 } else {
628                                         double f = globals->get_ATC_mgr()->GetFrequency(ident, ATIS) / 100.0;
629                                         string wtr;
630                                         if(!f) {
631                                                 wtr = ", " + GetWeather();
632                                         } else {
633                                                 wtr = ", runway " + ConvertRwyNumToSpokenString(activeRwy);
634                                         }
635                                         trns += " cleared to land" + wtr;
636                                 }
637                                 // TODO - add winds
638                                 t->clearedToLand = true;
639                                 // Maybe remove report downwind from menu here as well incase user didn't bother to?
640                                 if(t->isUser) {
641                                         //cout << "ADD VACATED B\n";
642                                         // Put going around at the top (and hence default) since that'll be more desperate,
643                                         // or put rwy vacated at the top since that'll be more common?
644                                         current_atcdialog->add_entry(ident, "@CS Going Around", "Report going around", TOWER, USER_REPORT_GOING_AROUND);
645                                         current_atcdialog->add_entry(ident, "@CS Clear of the runway", "Report runway vacated", TOWER, USER_REPORT_RWY_VACATED);
646                                 } else {
647                                         t->planePtr->RegisterTransmission(7);
648                                 }
649                         } else if(t->eta < 20) {
650                                 // Do nothing - we'll be telling it to go around in less than 10 seconds if the
651                                 // runway doesn't clear so no point in calling "continue approach".
652                                 disp = false;
653                         } else {
654                                 trns += " continue approach";
655                                 trns += " and report ";
656                                 trns += ((rwy.patternDirection == 1) ? "right " : "left ");
657                                 trns += "downwind runway " + ConvertRwyNumToSpokenString(activeRwy);
658                                 t->opType = CIRCUIT;
659                                 if(t->isUser) {
660                                         current_atcdialog->add_entry(ident, "@AP Tower, @CS Downwind @RW", "Report Downwind", TOWER, (int)USER_REPORT_DOWNWIND);
661                                 } else {
662                                         t->planePtr->RegisterTransmission(15);
663                                 }
664                                 t->clearedToLand = false;
665                         }
666                         if(_display && disp) {
667                                 pending_transmission = trns;
668                                 Transmit();
669                         }
670                         t->finalAcknowledged = true;
671                 } else if(t->rwyVacatedReported && !(t->rwyVacatedAcknowledged)) {
672                         ProcessRunwayVacatedReport(t);
673                         t->rwyVacatedAcknowledged = true;
674                 }
675         }
676         //freqClear = true;     // FIXME - set this to come true after enough time to render the message
677         _releaseCounter = 0.0;
678         _releaseTime = 5.5;
679         _runReleaseCounter = true;
680         //cout << "Done Respond\n" << endl;
681 }
682
683 void FGTower::ProcessDownwindReport(TowerPlaneRec* t) {
684         int i = 1;
685         int a = 0;      // Count of preceding planes on approach
686         bool cf = false;        // conflicting traffic on final
687         bool cc = false;        // preceding traffic in circuit
688         TowerPlaneRec* tc = NULL;
689         for(tower_plane_rec_list_iterator twrItr = circuitList.begin(); twrItr != circuitList.end(); twrItr++) {
690                 if((*twrItr)->plane.callsign == responseID) break;
691                 tc = *twrItr;
692                 ++i;
693         }
694         if(i > 1) { cc = true; }
695         doThresholdETACalc();
696         TowerPlaneRec* tf = NULL;
697         for(tower_plane_rec_list_iterator twrItr = appList.begin(); twrItr != appList.end(); twrItr++) {
698                 if((*twrItr)->eta < (t->eta + 45) && strcmp((*twrItr)->plane.callsign.c_str(), t->plane.callsign.c_str()) != 0) { // don't let ATC ask you to follow yourself
699                         a++;
700                         tf = *twrItr;
701                         cf = true;
702                         // This should set the flagged plane to be the last conflicting one, and hence the one to follow.
703                         // It ignores the fact that we might have problems slotting into the approach traffic behind it - 
704                         // eventually we'll need some fancy algorithms for that!
705                 }
706         }
707         string trns = t->plane.callsign;
708         trns += " Number ";
709         trns += ConvertNumToSpokenDigits(i + a);
710         // This assumes that the number spoken is landing position, not circuit position, since some of the traffic might be on straight-in final.
711         trns += " ";
712         TowerPlaneRec* tt = NULL;
713         if((i == 1) && rwyList.empty() && (t->nextOnRwy) && (!cf)) {    // Unfortunately nextOnRwy currently doesn't handle circuit/straight-in ordering properly at present, hence the cf check below.
714                 trns += "Cleared to land";      // TODO - clear for the option if appropriate
715                 t->clearedToLand = true;
716                 if(!t->isUser) t->planePtr->RegisterTransmission(7);
717         } else if((i+a) > 1) {
718                 //First set tt to point to the correct preceding plane - final or circuit
719                 if(tc && tf) {
720                         tt = (tf->eta < tc->eta ? tf : tc);
721                 } else if(tc) {
722                         tt = tc;
723                 } else if(tf) {
724                         tt = tf;
725                 } else {
726                         // We should never get here!
727                         SG_LOG(SG_ATC, SG_ALERT, "ALERT - Logic error in FGTower::ProcessDownwindReport");
728                         return;
729                 }
730                 trns += "Follow the ";
731                 string s = tt->plane.callsign;
732                 int p = s.find('-');
733                 s = s.substr(0,p);
734                 trns += s;
735                 if((tt->opType) == CIRCUIT) {
736                         PatternLeg leg;
737                         if(tt->isUser) {
738                                 leg = tt->leg;
739                         } else {
740                                 leg = tt->planePtr->GetLeg();
741                         }
742                         if(leg == FINAL) {
743                                 trns += " on final";
744                         } else if(leg == TURN4) {
745                                 trns += " turning final";
746                         } else if(leg == BASE) {
747                                 trns += " on base";
748                         } else if(leg == TURN3) {
749                                 trns += " turning base";
750                         }
751                 } else {
752                         double miles_out = CalcDistOutMiles(tt);
753                         if(miles_out < 2) {
754                                 trns += " on short final";
755                         } else {
756                                 trns += " on ";
757                                 trns += ConvertNumToSpokenDigits((int)miles_out);
758                                 trns += " mile final";
759                         }
760                 }
761         }
762         if(_display) {
763                 pending_transmission = trns;
764                 Transmit();
765         }
766         if(t->isUser) {
767                 if(t->opType == TTT_UNKNOWN) t->opType = CIRCUIT;
768                 //cout << "ADD VACATED A\n";
769                 // Put going around at the top (and hence default) since that'll be more desperate,
770                 // or put rwy vacated at the top since that'll be more common?
771                 //cout << "ident = " << ident << ", adding go-around option\n";
772                 current_atcdialog->add_entry(ident, "@CS Going Around", "Report going around", TOWER, USER_REPORT_GOING_AROUND);
773                 current_atcdialog->add_entry(ident, "@CS Clear of the runway", "Report runway vacated", TOWER, USER_REPORT_RWY_VACATED);
774         }
775 }
776
777 void FGTower::ProcessRunwayVacatedReport(TowerPlaneRec* t) {
778         //cout << "Processing rwy vacated...\n";
779         if(t->isUser) current_atcdialog->remove_entry(ident, USER_REPORT_GOING_AROUND, TOWER);
780         string trns = t->plane.callsign;
781         if(separateGround) {
782                 trns += " Contact ground on ";
783                 double f = globals->get_ATC_mgr()->GetFrequency(ident, GROUND) / 100.0; 
784                 char buf[10];
785                 sprintf(buf, "%.2f", f);
786                 trns += buf;
787                 trns += " Good Day";
788                 if(!t->isUser) t->planePtr->RegisterTransmission(5);
789         } else {
790                 // Cop-out!!
791                 trns += " cleared for taxi to general aviation parking";
792                 if(!t->isUser) t->planePtr->RegisterTransmission(6);    // TODO - this is a mega-hack!!
793         }
794         //cout << "trns = " << trns << '\n';
795         if(_display) {
796                 pending_transmission = trns;
797                 Transmit();
798         }
799         RemoveFromRwyList(t->plane.callsign);
800         AddToVacatedList(t);
801         // Maybe we should check that the plane really *has* vacated the runway!
802 }
803
804 // *********** END RESPONSE FUNCTIONS *****************
805
806 // Currently this assumes we *are* next on the runway and doesn't check for planes about to land - 
807 // this should be done prior to calling this function.
808 void FGTower::ClearHoldingPlane(TowerPlaneRec* t) {
809         //cout << "Entering ClearHoldingPlane..." << endl;
810         // Lets Roll !!!!
811         string trns = t->plane.callsign;
812         //if(departed plane < some threshold in time away) {
813         if(0) {         // FIXME
814         //if(timeSinceLastDeparture <= 60.0 && departed == true) {
815                 trns += " line up runway " + ConvertRwyNumToSpokenString(activeRwy);
816                 t->clearedToLineUp = true;
817                 t->planePtr->RegisterTransmission(3);   // cleared to line-up
818         //} else if(arriving plane < some threshold away) {
819         } else if(GetTrafficETA(2) < 150.0 && (timeSinceLastDeparture > 60.0 || departed == false)) {   // Hack - hardwired time
820                 trns += " cleared immediate take-off";
821                 if(trafficList.size()) {
822                         tower_plane_rec_list_iterator trfcItr = trafficList.begin();
823                         trfcItr++;      // At the moment the holding plane should be first in trafficList.
824                         // Note though that this will break if holding planes aren't put in trafficList in the future.
825                         TowerPlaneRec* trfc = *trfcItr;
826                         trns += "... traffic is";
827                         switch(trfc->plane.type) {
828                         case UNKNOWN:
829                                 break;
830                         case GA_SINGLE:
831                                 trns += " a Cessna";    // TODO - add ability to specify actual plane type somewhere
832                                 break;
833                         case GA_HP_SINGLE:
834                                 trns += " a Piper";
835                                 break;
836                         case GA_TWIN:
837                                 trns += " a King-air";
838                                 break;
839                         case GA_JET:
840                                 trns += " a Learjet";
841                                 break;
842                         case MEDIUM:
843                                 trns += " a Regional";
844                                 break;
845                         case HEAVY:
846                                 trns += " a Heavy";
847                                 break;
848                         case MIL_JET:
849                                 trns += " Military";
850                                 break;
851                         }
852                         //if(trfc->opType == STRAIGHT_IN || trfc->opType == TTT_UNKNOWN) {
853                         if(trfc->opType == STRAIGHT_IN) {
854                                 double miles_out = CalcDistOutMiles(trfc);
855                                 if(miles_out < 2) {
856                                         trns += " on final";
857                                 } else {
858                                         trns += " on ";
859                                         trns += ConvertNumToSpokenDigits((int)miles_out);
860                                         trns += " mile final";
861                                 }
862                         } else if(trfc->opType == CIRCUIT) {
863                                 //cout << "Getting leg of " << trfc->plane.callsign << '\n';
864                                 switch(trfc->leg) {
865                                 case FINAL:
866                                         trns += " on final";
867                                         break;
868                                 case TURN4:
869                                         trns += " turning final";
870                                         break;
871                                 case BASE:
872                                         trns += " on base";
873                                         break;
874                                 case TURN3:
875                                         trns += " turning base";
876                                         break;
877                                 case DOWNWIND:
878                                         trns += " in circuit";  // At the moment the user plane is generally flagged as unknown opType when downwind incase its a downwind departure which means we won't get here.
879                                         break;
880                                 // And to eliminate compiler warnings...
881                                 case TAKEOFF_ROLL: break;
882                                 case CLIMBOUT:     break;
883                                 case TURN1:        break;
884                                 case CROSSWIND:    break;
885                                 case TURN2:        break;
886                                 case LANDING_ROLL: break;
887                                 case LEG_UNKNOWN:  break;
888                                 }
889                         }
890                 } else {
891                         // By definition there should be some arriving traffic if we're cleared for immediate takeoff
892                         SG_LOG(SG_ATC, SG_WARN, "Warning: Departing traffic cleared for *immediate* take-off despite no arriving traffic in FGTower");
893                 }
894                 t->clearedToTakeOff = true;
895                 t->planePtr->RegisterTransmission(4);   // cleared to take-off - TODO differentiate between immediate and normal take-off
896                 departed = false;
897                 timeSinceLastDeparture = 0.0;
898         } else {
899         //} else if(timeSinceLastDeparture > 60.0 || departed == false) {       // Hack - test for timeSinceLastDeparture should be in lineup block eventually
900                 trns += " cleared for take-off";
901                 // TODO - add traffic is... ?
902                 t->clearedToTakeOff = true;
903                 t->planePtr->RegisterTransmission(4);   // cleared to take-off
904                 departed = false;
905                 timeSinceLastDeparture = 0.0;
906         }
907         if(_display) {
908                 pending_transmission = trns;
909                 Transmit();
910         }
911         //cout << "Done ClearHoldingPlane " << endl;
912 }
913
914
915 // ***************************************************************************************
916 // ********** Functions to periodically check what the various traffic is doing **********
917
918 // Do one plane from the hold list
919 void FGTower::CheckHoldList(double dt) {
920         //cout << "Entering CheckHoldList..." << endl;
921         if(!holdList.empty()) {
922                 //cout << "*holdListItr = " << *holdListItr << endl;
923                 if(holdListItr == holdList.end()) {
924                         holdListItr = holdList.begin();
925                 }
926                 //cout << "*holdListItr = " << *holdListItr << endl;
927                 //Process(*holdListItr);
928                 TowerPlaneRec* t = *holdListItr;
929                 //cout << "t = " << t << endl;
930                 if(t->holdShortReported) {
931                         // NO-OP - leave it to the response handler.
932                 } else {        // not responding to report, but still need to clear if clear
933                         if(t->nextOnRwy) {
934                                 //cout << "departed = " << departed << '\n';
935                                 //cout << "timeSinceLastDeparture = " << timeSinceLastDeparture << '\n';
936                                 if(rwyOccupied) {
937                                         RemoveAllUserDialogOptions();
938                                         current_atcdialog->add_entry(ident, "@CS Ready for take-off", "Request take-off clearance", TOWER, (int)USER_REQUEST_TAKE_OFF);
939                                 } else if(timeSinceLastDeparture <= 60.0 && departed == true) {
940                                         // Do nothing - this is a bit of a hack - should maybe do line up be ready here
941                                 } else {
942                                         ClearHoldingPlane(t);
943                                         t->leg = TAKEOFF_ROLL;
944                                         rwyList.push_back(t);
945                                         rwyListItr = rwyList.begin();
946                                         rwyOccupied = true;
947                                         holdList.erase(holdListItr);
948                                         holdListItr = holdList.begin();
949                                         if (holdList.empty())
950                                           return;
951                                 }
952                         }
953                         // TODO - rationalise the considerable code duplication above!
954                 }
955                 ++holdListItr;
956         }
957         //cout << "Done CheckHoldList" << endl;
958 }
959
960 // do the ciruit list
961 void FGTower::CheckCircuitList(double dt) {
962         //cout << "Entering CheckCircuitList..." << endl;
963         // Clear the constraints - we recalculate here.
964         base_leg_pos = 0.0;
965         downwind_leg_pos = 0.0;
966         crosswind_leg_pos = 0.0;
967         
968         if(!circuitList.empty()) {      // Do one plane from the circuit
969                 if(circuitListItr == circuitList.end()) {
970                         circuitListItr = circuitList.begin();
971                 }
972                 TowerPlaneRec* t = *circuitListItr;
973                 //cout << ident <<  ' ' << circuitList.size() << ' ' << t->plane.callsign << " " << t->leg << " eta " << t->eta << '\n';
974                 if(t->isUser) {
975                         t->pos.setlon(user_lon_node->getDoubleValue());
976                         t->pos.setlat(user_lat_node->getDoubleValue());
977                         t->pos.setelev(user_elev_node->getDoubleValue());
978                         //cout << ident <<  ' ' << circuitList.size() << ' ' << t->plane.callsign << " " << t->leg << " eta " << t->eta << '\n';
979                 } else {
980                         t->pos = t->planePtr->GetPos();         // We should probably only set the pos's on one walk through the traffic list in the update function, to save a few CPU should we end up duplicating this.
981                         t->landingType = t->planePtr->GetLandingOption();
982                         //cout << "AI plane landing option is " << t->landingType << '\n';
983                 }
984                 Point3D tortho = ortho.ConvertToLocal(t->pos);
985                 if(t->isUser) {
986                         // Need to figure out which leg he's on
987                         //cout << "rwy.hdg = " << rwy.hdg << " user hdg = " << user_hdg_node->getDoubleValue();
988                         double ho = GetAngleDiff_deg(user_hdg_node->getDoubleValue(), rwy.hdg);
989                         //cout << " ho = " << ho << " abs(ho = " << abs(ho) << '\n';
990                         // TODO FIXME - get the wind and convert this to track, or otherwise use track somehow!!!
991                         // If it's gusty might need to filter the value, although we are leaving 30 degrees each way leeway!
992                         if(fabs(ho) < 30) {
993                                 // could be either takeoff, climbout or landing - check orthopos.y
994                                 //cout << "tortho.y = " << tortho.y() << '\n';
995                                 if((tortho.y() < 0) || (t->leg == TURN4) || (t->leg == FINAL)) {
996                                         t->leg = FINAL;
997                                         //cout << "Final\n";
998                                 } else {
999                                         t->leg = CLIMBOUT;      // TODO - check elev wrt. apt elev to differentiate takeoff roll and climbout
1000                                         //cout << "Climbout\n";
1001                                         // If it's the user we may be unsure of his/her intentions.
1002                                         // (Hopefully the AI planes won't try confusing the sim!!!)
1003                                         //cout << "tortho.y = " << tortho.y() << '\n';
1004                                         if(t->opType == TTT_UNKNOWN) {
1005                                                 if(tortho.y() > 5000) {
1006                                                         // 5 km out from threshold - assume it's a departure
1007                                                         t->opType = OUTBOUND;   // TODO - could check if the user has climbed significantly above circuit altitude as well.
1008                                                         // Since we are unknown operation we should be in depList already.
1009                                                         //cout << ident << " Removing user from circuitList (TTT_UNKNOWN)\n";
1010                                                         circuitListItr = circuitList.erase(circuitListItr);
1011                                                         RemoveFromTrafficList(t->plane.callsign);
1012                                                         if (circuitList.empty())
1013                                                             return;
1014                                                 }
1015                                         } else if(t->opType == CIRCUIT) {
1016                                                 if(tortho.y() > 10000) {
1017                                                         // 10 km out - assume the user has abandoned the circuit!!
1018                                                         t->opType = OUTBOUND;
1019                                                         depList.push_back(t);
1020                                                         depListItr = depList.begin();
1021                                                         //cout << ident << " removing user from circuitList (CIRCUIT)\n";
1022                                                         circuitListItr = circuitList.erase(circuitListItr);
1023                                                         if (circuitList.empty())
1024                                                           return;
1025                                                 }
1026                                         }
1027                                 }
1028                         } else if(fabs(ho) < 60) {
1029                                 // turn1 or turn 4
1030                                 // TODO - either fix or doublecheck this hack by looking at heading and pattern direction
1031                                 if((t->leg == CLIMBOUT) || (t->leg == TURN1)) {
1032                                         t->leg = TURN1;
1033                                         //cout << "Turn1\n";
1034                                 } else {
1035                                         t->leg = TURN4;
1036                                         //cout << "Turn4\n";
1037                                 }
1038                         } else if(fabs(ho) < 120) {
1039                                 // crosswind or base
1040                                 // TODO - either fix or doublecheck this hack by looking at heading and pattern direction
1041                                 if((t->leg == TURN1) || (t->leg == CROSSWIND)) {
1042                                         t->leg = CROSSWIND;
1043                                         //cout << "Crosswind\n";
1044                                 } else {
1045                                         t->leg = BASE;
1046                                         //cout << "Base\n";
1047                                 }
1048                         } else if(fabs(ho) < 150) {
1049                                 // turn2 or turn 3
1050                                 // TODO - either fix or doublecheck this hack by looking at heading and pattern direction
1051                                 if((t->leg == CROSSWIND) || (t->leg == TURN2)) {
1052                                         t->leg = TURN2;
1053                                         //cout << "Turn2\n";
1054                                 } else {
1055                                         t->leg = TURN3;
1056                                         // Probably safe now to assume the user is flying a circuit
1057                                         t->opType = CIRCUIT;
1058                                         //cout << "Turn3\n";
1059                                 }
1060                         } else {
1061                                 // downwind
1062                                 t->leg = DOWNWIND;
1063                                 //cout << "Downwind\n";
1064                         }
1065                         if(t->leg == FINAL) {
1066                                 if(OnActiveRunway(t->pos)) {
1067                                         t->leg = LANDING_ROLL;
1068                                 }
1069                         }
1070                 } else {
1071                         t->leg = t->planePtr->GetLeg();
1072                 }
1073                 
1074                 // Set the constraints IF this is the first plane in the circuit
1075                 // TODO - at the moment we're constraining plane 2 based on plane 1 - this won't (or might not) work for 3 planes in the circuit!!
1076                 if(circuitListItr == circuitList.begin()) {
1077                         switch(t->leg) {
1078                                 case FINAL:
1079                                 // Base leg must be at least as far out as the plane is - actually possibly not necessary for separation, but we'll use that for now.
1080                                 base_leg_pos = tortho.y();
1081                                 //cout << "base_leg_pos = " << base_leg_pos << '\n';
1082                                 break;
1083                                 case TURN4:
1084                                 // Fall through to base
1085                                 case BASE:
1086                                 base_leg_pos = tortho.y();
1087                                 //cout << "base_leg_pos = " << base_leg_pos << '\n';
1088                                 break;
1089                                 case TURN3:
1090                                 // Fall through to downwind
1091                                 case DOWNWIND:
1092                                 // Only have the downwind leg pos as turn-to-base constraint if more negative than we already have.
1093                                 base_leg_pos = (tortho.y() < base_leg_pos ? tortho.y() : base_leg_pos);
1094                                 //cout << "base_leg_pos = " << base_leg_pos;
1095                                 downwind_leg_pos = tortho.x();          // Assume that a following plane can simply be constrained by the immediately in front downwind plane
1096                                 //cout << " downwind_leg_pos = " << downwind_leg_pos << '\n';
1097                                 break;
1098                                 case TURN2:
1099                                 // Fall through to crosswind
1100                                 case CROSSWIND:
1101                                 crosswind_leg_pos = tortho.y();
1102                                 //cout << "crosswind_leg_pos = " << crosswind_leg_pos << '\n';
1103                                 t->instructedToGoAround = false;
1104                                 break;
1105                                 case TURN1:
1106                                 // Fall through to climbout
1107                                 case CLIMBOUT:
1108                                 // Only use current by constraint as largest
1109                                 crosswind_leg_pos = (tortho.y() > crosswind_leg_pos ? tortho.y() : crosswind_leg_pos);
1110                                 //cout << "crosswind_leg_pos = " << crosswind_leg_pos << '\n';
1111                                 break;
1112                                 case TAKEOFF_ROLL:
1113                                 break;
1114                                 case LEG_UNKNOWN:
1115                                 break;
1116                                 case LANDING_ROLL:
1117                                 break;
1118                                 default:
1119                                 break;
1120                         }
1121                 }
1122                 
1123                 if(t->leg == FINAL && !(t->instructedToGoAround)) {
1124                         doThresholdETACalc();
1125                         doThresholdUseOrder();
1126                         /*
1127                         if(t->isUser) {
1128                                 cout << "Checking USER on final... ";
1129                                 cout << "eta " << t->eta;
1130                                 if(t->clearedToLand) cout << " cleared to land\n";
1131                         }
1132                         */
1133                         //cout << "YES FINAL, t->eta = " << t->eta << ", rwyList.size() = " << rwyList.size() << '\n';
1134                         if(t->landingType == FULL_STOP) {
1135                                 t->opType = INBOUND;
1136                                 //cout << "\n******** SWITCHING TO INBOUND AT POINT AAA *********\n\n";
1137                         }
1138                         if(t->eta < 12 && rwyList.size()) {
1139                                 // TODO - need to make this more sophisticated 
1140                                 // eg. is the plane accelerating down the runway taking off [OK],
1141                                 // or stationary near the start [V. BAD!!].
1142                                 // For now this should stop the AI plane landing on top of the user.
1143                                 tower_plane_rec_list_iterator twrItr;
1144                                 twrItr = rwyList.begin();
1145                                 TowerPlaneRec* tpr = *twrItr;
1146                                 if(strcmp(tpr->plane.callsign.c_str(), t->plane.callsign.c_str()) == 0
1147                                                 && rwyList.size() == 1) {
1148                                         // Fixing bug when ATC says that we must go around because of traffic on rwy
1149                                         // but that traffic is our plane! In future we can use this expression
1150                                         // for other ATC-messages like "On ground at 46, vacate left."
1151
1152                                 } else {
1153                                         string trns = t->plane.callsign;
1154                                         trns += " GO AROUND TRAFFIC ON RUNWAY I REPEAT GO AROUND";
1155                                         pending_transmission = trns;
1156                                         ImmediateTransmit();
1157                                         t->instructedToGoAround = true;
1158                                         t->clearedToLand = false;
1159                                         // Assume it complies!!!
1160                                         t->opType = CIRCUIT;
1161                                         t->leg = CLIMBOUT;
1162                                         if(t->planePtr) {
1163                                                 //cout << "Registering Go-around transmission with AI plane\n";
1164                                                 t->planePtr->RegisterTransmission(13);
1165                                         }
1166                                 }
1167                         } else if(!t->clearedToLand) {
1168                                 // The whip through the appList is a hack since currently t->nextOnRwy doesn't always work
1169                                 // TODO - fix this!
1170                                 bool cf = false;
1171                                 for(tower_plane_rec_list_iterator twrItr = appList.begin(); twrItr != appList.end(); twrItr++) {
1172                                         if((*twrItr)->eta < t->eta) {
1173                                                 cf = true;
1174                                         }
1175                                 }
1176                                 if(t->nextOnRwy && !cf) {
1177                                         if(!rwyList.size()) {
1178                                                 string trns = t->plane.callsign;
1179                                                 trns += " Cleared to land";
1180                                                 pending_transmission = trns;
1181                                                 Transmit();
1182                                                 //if(t->isUser) cout << "Transmitting cleared to Land!!!\n";
1183                                                 t->clearedToLand = true;
1184                                                 if(!t->isUser) {
1185                                                         t->planePtr->RegisterTransmission(7);
1186                                                 }
1187                                         }
1188                                 } else {
1189                                         //if(t->isUser) cout << "Not next\n";
1190                                 }
1191                         }
1192                 } else if(t->leg == LANDING_ROLL) {
1193                         //cout << t->plane.callsign << " has landed - adding to rwyList\n";
1194                         rwyList.push_front(t);
1195                         // TODO - if(!clearedToLand) shout something!!
1196                         t->clearedToLand = false;
1197                         RemoveFromTrafficList(t->plane.callsign);
1198                         if(t->isUser) {
1199                                 t->opType = TTT_UNKNOWN;
1200                         }       // TODO - allow the user to specify opType via ATC menu
1201                         //cout << ident << " Removing " << t->plane.callsign << " from circuitList..." << endl;
1202                         circuitListItr = circuitList.erase(circuitListItr);
1203                         if(circuitListItr == circuitList.end() ) {
1204                                 circuitListItr = circuitList.begin();
1205                                 // avoid increment of circuitListItr (would increment to second element, or crash if no element left)
1206                                 return;
1207                         }
1208                 }
1209                 ++circuitListItr;
1210         }
1211         //cout << "Done CheckCircuitList" << endl;
1212 }
1213
1214 // Do the runway list - we'll do the whole runway list since it's important and there'll never be many planes on the rwy at once!!
1215 // FIXME - at the moment it looks like we're only doing the first plane from the rwy list.
1216 // (However, at the moment there should only be one airplane on the rwy at once, until we
1217 // start allowing planes to line up whilst previous arrival clears the rwy.)
1218 void FGTower::CheckRunwayList(double dt) {
1219         //cout << "Entering CheckRunwayList..." << endl;
1220         if(rwyOccupied) {
1221                 if(!rwyList.size()) {
1222                         rwyOccupied = false;
1223                 } else {
1224                         rwyListItr = rwyList.begin();
1225                         TowerPlaneRec* t = *rwyListItr;
1226                         if(t->isUser) {
1227                                 t->pos.setlon(user_lon_node->getDoubleValue());
1228                                 t->pos.setlat(user_lat_node->getDoubleValue());
1229                                 t->pos.setelev(user_elev_node->getDoubleValue());
1230                         } else {
1231                                 t->pos = t->planePtr->GetPos();         // We should probably only set the pos's on one walk through the traffic list in the update function, to save a few CPU should we end up duplicating this.
1232                         }
1233                         bool on_rwy = OnActiveRunway(t->pos);
1234                         if(!on_rwy) {
1235                                 // TODO - for all of these we need to check what the user is *actually* doing!
1236                                 if((t->opType == INBOUND) || (t->opType == STRAIGHT_IN)) {
1237                                         //cout << "Tower " << ident << " is removing plane " << t->plane.callsign << " from rwy list (vacated)\n";
1238                                         //cout << "Size of rwylist was " << rwyList.size() << '\n';
1239                                         //cout << "Size of vacatedList was " << vacatedList.size() << '\n';
1240                                         RemoveFromRwyList(t->plane.callsign);
1241                                         AddToVacatedList(t);
1242                                         //cout << "Size of rwylist is " << rwyList.size() << '\n';
1243                                         //cout << "Size of vacatedList is " << vacatedList.size() << '\n';
1244                                         // At the moment we wait until Runway Vacated is reported by the plane before telling to contact ground etc.
1245                                         // It's possible we could be a bit more proactive about this.
1246                                 } else if(t->opType == OUTBOUND) {
1247                                         depList.push_back(t);
1248                                         depListItr = depList.begin();
1249                                         rwyList.pop_front();
1250                                         departed = true;
1251                                         timeSinceLastDeparture = 0.0;
1252                                 } else if(t->opType == CIRCUIT) {
1253                                         //cout << ident << " adding " << t->plane.callsign << " to circuitList" << endl;
1254                                         circuitList.push_back(t);
1255                                         circuitListItr = circuitList.begin();
1256                                         AddToTrafficList(t);
1257                                         rwyList.pop_front();
1258                                         departed = true;
1259                                         timeSinceLastDeparture = 0.0;
1260                                 } else if(t->opType == TTT_UNKNOWN) {
1261                                         depList.push_back(t);
1262                                         depListItr = depList.begin();
1263                                         //cout << ident << " adding " << t->plane.callsign << " to circuitList" << endl;
1264                                         circuitList.push_back(t);
1265                                         circuitListItr = circuitList.begin();
1266                                         AddToTrafficList(t);
1267                                         rwyList.pop_front();
1268                                         departed = true;
1269                                         timeSinceLastDeparture = 0.0;   // TODO - we need to take into account that the user might taxi-in when flagged opType UNKNOWN - check speed/altitude etc to make decision as to what user is up to.
1270                                 } else {
1271                                         // HELP - we shouldn't ever get here!!!
1272                                 }
1273                         }
1274                 }
1275         }
1276         //cout << "Done CheckRunwayList" << endl;
1277 }
1278
1279 // Do one plane from the approach list
1280 void FGTower::CheckApproachList(double dt) {
1281         //cout << "CheckApproachList called for " << ident << endl;
1282         //cout << "AppList.size is " << appList.size() << endl;
1283         if(!appList.empty()) {
1284                 if(appListItr == appList.end()) {
1285                         appListItr = appList.begin();
1286                 }
1287                 TowerPlaneRec* t = *appListItr;
1288                 //cout << "t = " << t << endl;
1289                 //cout << "Checking " << t->plane.callsign << endl;
1290                 if(t->isUser) {
1291                         t->pos.setlon(user_lon_node->getDoubleValue());
1292                         t->pos.setlat(user_lat_node->getDoubleValue());
1293                         t->pos.setelev(user_elev_node->getDoubleValue());
1294                 } else {
1295                         // TODO - set/update the position if it's an AI plane
1296                 }
1297                 doThresholdETACalc();   // We need this here because planes in the lists are not guaranteed to *always* have the correct ETA
1298                 //cout << "eta is " << t->eta << ", rwy is " << (rwyList.size() ? "occupied " : "clear ") << '\n';
1299                 Point3D tortho = ortho.ConvertToLocal(t->pos);
1300                 if(t->eta < 12 && rwyList.size() && !(t->instructedToGoAround)) {
1301                         // TODO - need to make this more sophisticated 
1302                         // eg. is the plane accelerating down the runway taking off [OK],
1303                         // or stationary near the start [V. BAD!!].
1304                         // For now this should stop the AI plane landing on top of the user.
1305                         tower_plane_rec_list_iterator twrItr;
1306                         twrItr = rwyList.begin();
1307                         TowerPlaneRec* tpr = *twrItr;
1308                         if(strcmp ( tpr->plane.callsign.c_str(), t->plane.callsign.c_str() ) == 0 && rwyList.size() == 1) {
1309                                         // Fixing bug when ATC says that we must go around because of traffic on rwy
1310                                         // but that traffic is we! In future we can use this expression
1311                                         // for other ATC-messages like "On ground at 46, vacate left."
1312
1313                         } else {
1314                                 string trns = t->plane.callsign;
1315                                 trns += " GO AROUND TRAFFIC ON RUNWAY I REPEAT GO AROUND";
1316                                 pending_transmission = trns;
1317                                 ImmediateTransmit();
1318                                 t->instructedToGoAround = true;
1319                                 t->clearedToLand = false;
1320                                 t->nextOnRwy = false;   // But note this is recalculated so don't rely on it
1321                                 // Assume it complies!!!
1322                                 t->opType = CIRCUIT;
1323                                 t->leg = CLIMBOUT;
1324                                 if(!t->isUser) {
1325                                         if(t->planePtr) {
1326                                                 //cout << "Registering Go-around transmission with AI plane\n";
1327                                                 t->planePtr->RegisterTransmission(13);
1328                                         }
1329                                 } else {
1330                                         // TODO - add Go-around ack to comm options,
1331                                         // remove report rwy vacated. (possibly).
1332                                 }
1333                         }
1334                 } else if(t->isUser && t->eta < 90 && tortho.y() > -2500 && t->clearedToLand && t->gearUpReported == false) {
1335                         // Check if gear up or down
1336                         double gp = fgGetFloat("/gear/gear/position-norm");
1337                         if(gp < 1) {
1338                                 string trnsm = t->plane.callsign;
1339                                 sg_srandom_time();
1340                                 int rnd = int(sg_random() * 2) + 1;
1341                                 if(rnd == 2) {                          // Random message for more realistic ATC ;)
1342                                         trnsm += ", LANDING GEAR APPEARS UP!";
1343                                 } else {
1344                                         trnsm += ", Check wheels down and locked.";
1345                                 }
1346                                 pending_transmission = trnsm;
1347                                 ImmediateTransmit();
1348                                 t->gearUpReported = true;
1349                         }
1350                 } else if(t->eta < 90 && !t->clearedToLand) {
1351                         //doThresholdETACalc();
1352                         doThresholdUseOrder();
1353                         // The whip through the appList is a hack since currently t->nextOnRwy doesn't always work
1354                         // TODO - fix this!
1355                         bool cf = false;
1356                         for(tower_plane_rec_list_iterator twrItr = appList.begin(); twrItr != appList.end(); twrItr++) {
1357                                 if((*twrItr)->eta < t->eta) {
1358                                         cf = true;
1359                                 }
1360                         }
1361                         if(t->nextOnRwy && !cf) {
1362                                 if(!rwyList.size()) {
1363                                         string trns = t->plane.callsign;
1364                                         trns += " Cleared to land";
1365                                         pending_transmission = trns;
1366                                         Transmit();
1367                                         //if(t->isUser) cout << "Transmitting cleared to Land!!!\n";
1368                                         t->clearedToLand = true;
1369                                         if(!t->isUser) {
1370                                                 t->planePtr->RegisterTransmission(7);
1371                                         }
1372                                 }
1373                         } else {
1374                                 //if(t->isUser) cout << "Not next\n";
1375                         }
1376                 }
1377                 
1378                 // Check for landing...
1379                 bool landed = false;
1380                 if(!t->isUser) {
1381                         if(t->planePtr) {
1382                                 if(t->planePtr->GetLeg() == LANDING_ROLL) {
1383                                         landed = true;
1384                                 }
1385                         } else {
1386                                 SG_LOG(SG_ATC, SG_ALERT, "WARNING - not user and null planePtr in CheckApproachList!");
1387                         }
1388                 } else {        // user
1389                         if(OnActiveRunway(t->pos)) {
1390                                 landed = true;
1391                         }
1392                 }
1393                 
1394                 if(landed) {
1395                         // Duplicated in CheckCircuitList - must be able to rationalise this somehow!
1396                         //cout << "A " << t->plane.callsign << " has landed, adding to rwyList...\n";
1397                         rwyList.push_front(t);
1398                         // TODO - if(!clearedToLand) shout something!!
1399                         t->clearedToLand = false;
1400                         RemoveFromTrafficList(t->plane.callsign);
1401                         //if(t->isUser) {
1402                                 //      t->opType = TTT_UNKNOWN;
1403                         //}     // TODO - allow the user to specify opType via ATC menu
1404                         appListItr = appList.erase(appListItr);
1405                         if(appListItr == appList.end() ) {
1406                                 appListItr = appList.begin();
1407                         }
1408                         if (appList.empty())
1409                           return;
1410
1411                 }
1412                 
1413                 ++appListItr;
1414         }
1415         //cout << "Done" << endl;
1416 }
1417
1418 // Do one plane from the departure list
1419 void FGTower::CheckDepartureList(double dt) {
1420         if(!depList.empty()) {
1421                 if(depListItr == depList.end()) {
1422                         depListItr = depList.begin();
1423                 }
1424                 TowerPlaneRec* t = *depListItr;
1425                 //cout << "Dep list, checking " << t->plane.callsign;
1426                 
1427                 double distout; // meters
1428                 if(t->isUser) distout = dclGetHorizontalSeparation(Point3D(lon, lat, elev), Point3D(user_lon_node->getDoubleValue(), user_lat_node->getDoubleValue(), user_elev_node->getDoubleValue()));
1429                 else distout = dclGetHorizontalSeparation(Point3D(lon, lat, elev), t->planePtr->GetPos());
1430                 //cout << " distout = " << distout << '\n';
1431                 if(t->isUser && !(t->clearedToTakeOff)) {       // HACK - we use clearedToTakeOff to check if ATC already contacted with plane (and cleared take-off) or not
1432                         if(!OnAnyRunway(Point3D(user_lon_node->getDoubleValue(), user_lat_node->getDoubleValue(), 0.0), false)) {
1433                                 current_atcdialog->remove_entry(ident, USER_REQUEST_TAKE_OFF, TOWER);
1434                                 t->clearedToTakeOff = true;     // FIXME
1435                         }
1436                 }
1437                 if(distout > 10000) {
1438                         string trns = t->plane.callsign;
1439                         trns += " You are now clear of my airspace, good day";
1440                         pending_transmission = trns;
1441                         Transmit();
1442                         if(t->isUser) {
1443                                 // Change the communication options
1444                                 RemoveAllUserDialogOptions();
1445                                 //cout << "ADD A\n";
1446                                 current_atcdialog->add_entry(ident, "@AP Tower, @CS @MI miles @CD of the airport for full stop@AT", "Contact tower for VFR arrival (full stop)", TOWER, (int)USER_REQUEST_VFR_ARRIVAL_FULL_STOP);
1447                         } else {
1448                                 // Send a clear-of-airspace signal
1449                                 // TODO - implement this once we actually have departing AI traffic (currently all circuits or arrivals).
1450                         }
1451                         RemovePlane(t->plane.callsign);
1452                 } else {
1453                         ++depListItr;
1454                 }
1455         }
1456 }
1457
1458 // ********** End periodic check functions ***********************************************
1459 // ***************************************************************************************
1460
1461
1462 // Remove all dialog options for this tower.
1463 void FGTower::RemoveAllUserDialogOptions() {
1464         current_atcdialog->remove_entry(ident, USER_REQUEST_VFR_DEPARTURE, TOWER);
1465         current_atcdialog->remove_entry(ident, USER_REQUEST_VFR_ARRIVAL, TOWER);
1466         current_atcdialog->remove_entry(ident, USER_REQUEST_VFR_ARRIVAL_FULL_STOP, TOWER);
1467         current_atcdialog->remove_entry(ident, USER_REQUEST_VFR_ARRIVAL_TOUCH_AND_GO, TOWER);
1468         current_atcdialog->remove_entry(ident, USER_REPORT_3_MILE_FINAL, TOWER);
1469         current_atcdialog->remove_entry(ident, USER_REPORT_DOWNWIND, TOWER);
1470         current_atcdialog->remove_entry(ident, USER_REPORT_RWY_VACATED, TOWER);
1471         current_atcdialog->remove_entry(ident, USER_REPORT_GOING_AROUND, TOWER);        
1472         current_atcdialog->remove_entry(ident, USER_REQUEST_TAKE_OFF, TOWER);
1473 }
1474
1475 // Returns true if positions of crosswind/downwind/base leg turns should be constrained by previous traffic
1476 // plus the constraint position as a rwy orientated orthopos (meters)
1477 bool FGTower::GetCrosswindConstraint(double& cpos) {
1478         if(crosswind_leg_pos != 0.0) {
1479                 cpos = crosswind_leg_pos;
1480                 return(true);
1481         } else {
1482                 cpos = 0.0;
1483                 return(false);
1484         }
1485 }
1486 bool FGTower::GetDownwindConstraint(double& dpos) {
1487         if(fabs(downwind_leg_pos) > nominal_downwind_leg_pos) {
1488                 dpos = downwind_leg_pos;
1489                 return(true);
1490         } else {
1491                 dpos = 0.0;
1492                 return(false);
1493         }
1494 }
1495 bool FGTower::GetBaseConstraint(double& bpos) {
1496         if(base_leg_pos < nominal_base_leg_pos) {
1497                 bpos = base_leg_pos;
1498                 return(true);
1499         } else {
1500                 bpos = nominal_base_leg_pos;
1501                 return(false);
1502         }
1503 }
1504
1505
1506 // Figure out which runways are active.
1507 // For now we'll just be simple and do one active runway - eventually this will get much more complex
1508 // This is a private function - public interface to the results of this is through GetActiveRunway
1509 void FGTower::DoRwyDetails() {
1510         //cout << "GetRwyDetails called" << endl;
1511         
1512         // Based on the airport-id and wind get the active runway
1513         
1514   const FGAirport* apt = fgFindAirportID(ident);
1515   assert(apt);
1516         FGRunway* runway = apt->getActiveRunwayForUsage();
1517
1518   activeRwy = runway->ident();
1519   rwy.rwyID = runway->ident();
1520   SG_LOG(SG_ATC, SG_INFO, "In FGGround, active runway for airport " << ident << " is " << activeRwy);
1521   
1522   // Get the threshold position
1523   double other_way = runway->headingDeg() - 180.0;
1524   while(other_way <= 0.0) {
1525     other_way += 360.0;
1526   }
1527     // move to the +l end/center of the runway
1528   //cout << "Runway center is at " << runway._lon << ", " << runway._lat << '\n';
1529     Point3D origin = Point3D(runway->longitude(), runway->latitude(), aptElev);
1530   Point3D ref = origin;
1531     double tshlon, tshlat, tshr;
1532   double tolon, tolat, tor;
1533   rwy.length = runway->lengthM();
1534     geo_direct_wgs_84 ( aptElev, ref.lat(), ref.lon(), other_way, 
1535                         rwy.length / 2.0 - 25.0, &tshlat, &tshlon, &tshr );
1536     geo_direct_wgs_84 ( aptElev, ref.lat(), ref.lon(), runway->headingDeg(), 
1537                         rwy.length / 2.0 - 25.0, &tolat, &tolon, &tor );
1538
1539   // Note - 25 meters in from the runway end is a bit of a hack to put the plane ahead of the user.
1540   // now copy what we need out of runway into rwy
1541     rwy.threshold_pos = Point3D(tshlon, tshlat, aptElev);
1542   Point3D takeoff_end = Point3D(tolon, tolat, aptElev);
1543   //cout << "Threshold position = " << tshlon << ", " << tshlat << ", " << aptElev << '\n';
1544   //cout << "Takeoff position = " << tolon << ", " << tolat << ", " << aptElev << '\n';
1545   rwy.hdg = runway->headingDeg();
1546   // Set the projection for the local area based on this active runway
1547   ortho.Init(rwy.threshold_pos, rwy.hdg);       
1548   rwy.end1ortho = ortho.ConvertToLocal(rwy.threshold_pos);      // should come out as zero
1549   rwy.end2ortho = ortho.ConvertToLocal(takeoff_end);
1550   
1551   // Set the pattern direction
1552   // TODO - we'll check for a facilities file with this in eventually - for now assume left traffic except
1553   // for certain circumstances (RH parallel rwy).
1554   rwy.patternDirection = -1;            // Left
1555   if(rwy.rwyID.size() == 3) {
1556     rwy.patternDirection = (rwy.rwyID.substr(2,1) == "R" ? 1 : -1);
1557   }
1558   //cout << "Doing details, rwy.patterDirection is " << rwy.patternDirection << '\n';
1559 }
1560
1561
1562 // Figure out if a given position lies on the active runway
1563 // Might have to change when we consider more than one active rwy.
1564 bool FGTower::OnActiveRunway(const Point3D& pt) {
1565         // TODO - check that the centre calculation below isn't confused by displaced thesholds etc.
1566         Point3D xyc((rwy.end1ortho.x() + rwy.end2ortho.x())/2.0, (rwy.end1ortho.y() + rwy.end2ortho.y())/2.0, 0.0);
1567         Point3D xyp = ortho.ConvertToLocal(pt);
1568         
1569         //cout << "Length offset = " << fabs(xyp.y() - xyc.y()) << '\n';
1570         //cout << "Width offset = " << fabs(xyp.x() - xyc.x()) << '\n';
1571
1572         double rlen = rwy.length/2.0 + 5.0;
1573         double rwidth = rwy.width/2.0;
1574         double ldiff = fabs(xyp.y() - xyc.y());
1575         double wdiff = fabs(xyp.x() - xyc.x());
1576
1577         return((ldiff < rlen) && (wdiff < rwidth));
1578 }
1579
1580 // Figure out if a given position lies on any runway or not
1581 // Only call this at startup - reading the runways database is expensive and needs to be fixed!
1582 bool FGTower::OnAnyRunway(const Point3D& pt, bool onGround) {
1583         ATCData ad;
1584         double dist = current_commlist->FindClosest(lon, lat, elev, ad, TOWER, 7.0);
1585         if(dist < 0.0) {
1586                 return(false);
1587         }
1588         
1589         // Based on the airport-id, go through all the runways and check for a point in them
1590
1591   const FGAirport* apt = fgFindAirportID(ad.ident);
1592   assert(apt);
1593   
1594   for (unsigned int i=0; i<apt->numRunways(); ++i) {
1595     if (OnRunway(pt, apt->getRunwayByIndex(i))) {
1596       return true;
1597     }
1598   }
1599
1600   // if onGround is true, we only match real runways, so we're done
1601   if (onGround) return false;
1602
1603   // try taxiways as well
1604   for (unsigned int i=0; i<apt->numTaxiways(); ++i) {
1605     if (OnRunway(pt, apt->getTaxiwayByIndex(i))) {
1606       return true;
1607     }
1608   }
1609   
1610         return false;
1611 }
1612
1613
1614 // Returns true if successful
1615 bool FGTower::RemoveFromTrafficList(const string& id) {
1616         tower_plane_rec_list_iterator twrItr;
1617         for(twrItr = trafficList.begin(); twrItr != trafficList.end(); twrItr++) {
1618                 TowerPlaneRec* tpr = *twrItr;
1619                 if(tpr->plane.callsign == id) {
1620                         trafficList.erase(twrItr);
1621                         trafficListItr = trafficList.begin();
1622                         return(true);
1623                 }
1624         }       
1625         SG_LOG(SG_ATC, SG_WARN, "Warning - unable to remove aircraft " << id << " from trafficList in FGTower");
1626         return(false);
1627 }
1628
1629
1630 // Returns true if successful
1631 bool FGTower::RemoveFromAppList(const string& id) {
1632         tower_plane_rec_list_iterator twrItr;
1633         for(twrItr = appList.begin(); twrItr != appList.end(); twrItr++) {
1634                 TowerPlaneRec* tpr = *twrItr;
1635                 if(tpr->plane.callsign == id) {
1636                         appList.erase(twrItr);
1637                         appListItr = appList.begin();
1638                         return(true);
1639                 }
1640         }       
1641         //SG_LOG(SG_ATC, SG_WARN, "Warning - unable to remove aircraft " << id << " from appList in FGTower");
1642         return(false);
1643 }
1644
1645 // Returns true if successful
1646 bool FGTower::RemoveFromRwyList(const string& id) {
1647         tower_plane_rec_list_iterator twrItr;
1648         for(twrItr = rwyList.begin(); twrItr != rwyList.end(); twrItr++) {
1649                 TowerPlaneRec* tpr = *twrItr;
1650                 if(tpr->plane.callsign == id) {
1651                         rwyList.erase(twrItr);
1652                         rwyListItr = rwyList.begin();
1653                         return(true);
1654                 }
1655         }       
1656         //SG_LOG(SG_ATC, SG_WARN, "Warning - unable to remove aircraft " << id << " from rwyList in FGTower");
1657         return(false);
1658 }
1659
1660
1661 // Add a tower plane rec with ETA to the traffic list in the correct position ETA-wise
1662 // and set nextOnRwy if so.
1663 // Returns true if this could cause a threshold ETA conflict with other traffic, false otherwise.
1664 // For planes holding they are put in the first position with time to go, and the return value is
1665 // true if in the first position (nextOnRwy) and false otherwise.
1666 // See the comments in FGTower::doThresholdUseOrder for notes on the ordering
1667 bool FGTower::AddToTrafficList(TowerPlaneRec* t, bool holding) {
1668         //cout << "ADD: " << trafficList.size();
1669         //cout << "AddToTrafficList called, currently size = " << trafficList.size() << ", holding = " << holding << endl;
1670         double separation_time = 90.0;  // seconds - this is currently a guess for light plane separation, and includes a few seconds for a holding plane to taxi onto the rwy.
1671         double departure_sep_time = 60.0;       // Separation time behind departing airplanes.  Comments above also apply.
1672         bool conflict = false;
1673         double lastETA = 0.0;
1674         bool firstTime = true;
1675         // FIXME - make this more robust for different plane types eg. light following heavy.
1676         tower_plane_rec_list_iterator twrItr;
1677         //twrItr = trafficList.begin();
1678         //while(1) {
1679         for(twrItr = trafficList.begin(); twrItr != trafficList.end(); twrItr++) {
1680                 //if(twrItr == trafficList.end()) {
1681                 //      cout << "  END  ";
1682                 //      trafficList.push_back(t);
1683                 //      return(holding ? firstTime : conflict);
1684                 //} else {
1685                         TowerPlaneRec* tpr = *twrItr;
1686                         if(holding) {
1687                                 //cout << (tpr->isUser ? "USER!\n" : "NOT user\n");
1688                                 //cout << "tpr->eta - lastETA = " << tpr->eta - lastETA << '\n';
1689                                 double dep_allowance = (timeSinceLastDeparture < departure_sep_time ? departure_sep_time - timeSinceLastDeparture : 0.0); 
1690                                 double slot_time = (firstTime ? separation_time + dep_allowance : separation_time + departure_sep_time);
1691                                 // separation_time + departure_sep_time in the above accounts for the fact that the arrival could be touch and go,
1692                                 // and if not needs time to clear the rwy anyway.
1693                                 if(tpr->eta  - lastETA > slot_time) {
1694                                         t->nextOnRwy = firstTime;
1695                                         trafficList.insert(twrItr, t);
1696                                         //cout << "\tH\t" << trafficList.size() << '\n';
1697                                         return(firstTime);
1698                                 }
1699                                 firstTime = false;
1700                         } else {
1701                                 if(t->eta < tpr->eta) {
1702                                         // Ugg - this one's tricky.
1703                                         // It depends on what the two planes are doing and whether there's a conflict what we do.
1704                                         if(tpr->eta - t->eta > separation_time) {       // No probs, plane 2 can squeeze in before plane 1 with no apparent conflict
1705                                                 if(tpr->nextOnRwy) {
1706                                                         tpr->nextOnRwy = false;
1707                                                         t->nextOnRwy = true;
1708                                                 }
1709                                                 trafficList.insert(twrItr, t);
1710                                         } else {        // Ooops - this ones tricky - we have a potential conflict!
1711                                                 conflict = true;
1712                                                 // HACK - just add anyway for now and flag conflict - TODO - FIX THIS using CIRCUIT/STRAIGHT_IN and VFR/IFR precedence rules.
1713                                                 if(tpr->nextOnRwy) {
1714                                                         tpr->nextOnRwy = false;
1715                                                         t->nextOnRwy = true;
1716                                                 }
1717                                                 trafficList.insert(twrItr, t);
1718                                         }
1719                                         //cout << "\tC\t" << trafficList.size() << '\n';
1720                                         return(conflict);
1721                                 }
1722                         }
1723                 //}
1724                 //++twrItr;
1725         }
1726         // If we get here we must be at the end of the list, or maybe the list is empty.
1727         if(!trafficList.size()) {
1728                 t->nextOnRwy = true;
1729                 // conflict and firstTime should be false and true respectively in this case anyway.
1730         } else {
1731                 t->nextOnRwy = false;
1732         }
1733         trafficList.push_back(t);
1734         //cout << "\tE\t" << trafficList.size() << endl;
1735         return(holding ? firstTime : conflict);
1736 }
1737
1738 // Add a tower plane rec with ETA to the circuit list in the correct position ETA-wise
1739 // Returns true if this might cause a separation conflict (based on ETA) with other traffic, false otherwise.
1740 // Safe to add a plane that is already in - planes with the same callsign are not added.
1741 bool FGTower::AddToCircuitList(TowerPlaneRec* t) {
1742         if(!t) {
1743                 //cout << "**********************************************\n";
1744                 //cout << "AddToCircuitList called with NULL pointer!!!!!\n";
1745                 //cout << "**********************************************\n";
1746                 return false;
1747         }
1748         //cout << "ADD: " << circuitList.size();
1749         //cout << ident << " AddToCircuitList called for " << t->plane.callsign << ", currently size = " << circuitList.size() << endl;
1750         double separation_time = 60.0;  // seconds - this is currently a guess for light plane separation, and includes a few seconds for a holding plane to taxi onto the rwy.
1751         bool conflict = false;
1752         tower_plane_rec_list_iterator twrItr;
1753         // First check if the plane is already in the list
1754         //cout << "A" << endl;
1755         //cout << "Checking whether " << t->plane.callsign << " is already in circuit list..." << endl;
1756         //cout << "B" << endl;
1757         for(twrItr = circuitList.begin(); twrItr != circuitList.end(); twrItr++) {
1758                 if((*twrItr)->plane.callsign == t->plane.callsign) {
1759                         //cout << "In list - returning...\n";
1760                         return false;
1761                 }
1762         }
1763         //cout << "Not in list - adding..." << endl;
1764         
1765         for(twrItr = circuitList.begin(); twrItr != circuitList.end(); twrItr++) {
1766                 TowerPlaneRec* tpr = *twrItr;
1767                 //cout << tpr->plane.callsign << " eta is " << tpr->eta << '\n';
1768                 //cout << "New eta is " << t->eta << '\n';              
1769                 if(t->eta < tpr->eta) {
1770                         // Ugg - this one's tricky.
1771                         // It depends on what the two planes are doing and whether there's a conflict what we do.
1772                         if(tpr->eta - t->eta > separation_time) {       // No probs, plane 2 can squeeze in before plane 1 with no apparent conflict
1773                                 circuitList.insert(twrItr, t);
1774                                 circuitListItr = circuitList.begin();
1775                         } else {        // Ooops - this ones tricky - we have a potential conflict!
1776                                 conflict = true;
1777                                 // HACK - just add anyway for now and flag conflict.
1778                                 circuitList.insert(twrItr, t);
1779                                 circuitListItr = circuitList.begin();
1780                         }
1781                         //cout << "\tC\t" << circuitList.size() << '\n';
1782                         return(conflict);
1783                 }
1784         }
1785         // If we get here we must be at the end of the list, or maybe the list is empty.
1786         //cout << ident << " adding " << t->plane.callsign << " to circuitList" << endl;
1787         circuitList.push_back(t);       // TODO - check the separation with the preceding plane for the conflict flag.
1788         circuitListItr = circuitList.begin();
1789         //cout << "\tE\t" << circuitList.size() << endl;
1790         return(conflict);
1791 }
1792
1793 // Add to vacated list only if not already present
1794 void FGTower::AddToVacatedList(TowerPlaneRec* t) {
1795         tower_plane_rec_list_iterator twrItr;
1796         bool found = false;
1797         for(twrItr = vacatedList.begin(); twrItr != vacatedList.end(); twrItr++) {
1798                 if((*twrItr)->plane.callsign == t->plane.callsign) {
1799                         found = true;
1800                 }
1801         }
1802         if(found) return;
1803         vacatedList.push_back(t);
1804 }
1805
1806 void FGTower::AddToHoldingList(TowerPlaneRec* t) {
1807         tower_plane_rec_list_iterator it, end = holdList.end();
1808         for (it = holdList.begin(); it != end; ++it) {
1809                 if ((*it)->plane.callsign == t->plane.callsign)
1810                         return;
1811         
1812                 holdList.push_back(t);
1813         }
1814 }
1815
1816 // Calculate the eta of a plane to the threshold.
1817 // For ground traffic this is the fastest they can get there.
1818 // For air traffic this is the middle approximation.
1819 void FGTower::CalcETA(TowerPlaneRec* tpr, bool printout) {
1820         // For now we'll be very crude and hardwire expected speeds to C172-like values
1821         // The speeds below are specified in knots IAS and then converted to m/s
1822         double app_ias = 100.0 * 0.514444;                      // Speed during straight-in approach
1823         double circuit_ias = 80.0 * 0.514444;           // Speed around circuit
1824         double final_ias = 70.0 * 0.514444;             // Speed during final approach
1825         
1826         //if(printout) {
1827         //cout << "In CalcETA, airplane ident = " << tpr->plane.callsign << '\n';
1828         //cout << (tpr->isUser ? "USER\n" : "AI\n");
1829         //cout << flush;
1830         //}
1831         
1832         // Sign convention - dist_out is -ve for approaching planes and +ve for departing planes
1833         // dist_across is +ve in the pattern direction - ie a plane correctly on downwind will have a +ve dist_across
1834         
1835         Point3D op = ortho.ConvertToLocal(tpr->pos);
1836         //if(printout) {
1837         //if(!tpr->isUser) cout << "Orthopos is " << op.x() << ", " << op.y() << ' ';
1838         //cout << "opType is " << tpr->opType << '\n';
1839         //}
1840         double dist_out_m = op.y();
1841         double dist_across_m = fabs(op.x());    // The fabs is a hack to cope with the fact that we don't know the circuit direction yet
1842         //cout << "Doing ETA calc for " << tpr->plane.callsign << '\n';
1843         
1844         if(tpr->opType == STRAIGHT_IN || tpr->opType == INBOUND) {
1845                 //cout << "CASE 1\n";
1846                 double dist_to_go_m = sqrt((dist_out_m * dist_out_m) + (dist_across_m * dist_across_m));
1847                 if(dist_to_go_m < 1000) {
1848                         tpr->eta = dist_to_go_m / final_ias;
1849                 } else {
1850                         tpr->eta = (1000.0 / final_ias) + ((dist_to_go_m - 1000.0) / app_ias);
1851                 }
1852         } else if(tpr->opType == CIRCUIT || tpr->opType == TTT_UNKNOWN) {       // Hack alert - UNKNOWN has sort of been added here as a temporary hack.
1853                 //cout << "CASE 2\n";
1854                 // It's complicated - depends on if base leg is delayed or not
1855                 //if(printout) {
1856                 //cout << "Leg = " << tpr->leg << '\n';
1857                 //}
1858                 if(tpr->leg == LANDING_ROLL) {
1859                         tpr->eta = 0;
1860                 } else if((tpr->leg == FINAL) || (tpr->leg == TURN4)) {
1861                         //cout << "dist_out_m = " << dist_out_m << '\n';
1862                         tpr->eta = fabs(dist_out_m) / final_ias;
1863                 } else if((tpr->leg == BASE) || (tpr->leg == TURN3)) {
1864                         tpr->eta = (fabs(dist_out_m) / final_ias) + (dist_across_m / circuit_ias);
1865                 } else {
1866                         // Need to calculate where base leg is likely to be
1867                         // FIXME - for now I'll hardwire it to 1000m which is what AILocalTraffic uses!!!
1868                         // TODO - as a matter of design - AILocalTraffic should get the nominal no-traffic base turn distance from Tower, since in real life the published pattern might differ from airport to airport
1869                         double nominal_base_dist_out_m = -1000;
1870                         double current_base_dist_out_m;
1871                         if(!GetBaseConstraint(current_base_dist_out_m)) {
1872                                 current_base_dist_out_m = nominal_base_dist_out_m;
1873                         }
1874                         //cout << "current_base_dist_out_m = " << current_base_dist_out_m << '\n';
1875                         double nominal_dist_across_m = 1000;    // Hardwired value from AILocalTraffic
1876                         double current_dist_across_m;
1877                         if(!GetDownwindConstraint(current_dist_across_m)) {
1878                                 current_dist_across_m = nominal_dist_across_m;
1879                         }
1880                         double nominal_cross_dist_out_m = 2000; // Bit of a guess - AI plane turns to crosswind at 700ft agl.
1881                         tpr->eta = fabs(current_base_dist_out_m) / final_ias;   // final
1882                         //cout << "a = " << tpr->eta << '\n';
1883                         if((tpr->leg == DOWNWIND) || (tpr->leg == TURN2)) {
1884                                 tpr->eta += dist_across_m / circuit_ias;
1885                                 //cout << "b = " << tpr->eta << '\n';
1886                                 tpr->eta += fabs(current_base_dist_out_m - dist_out_m) / circuit_ias;
1887                                 //cout << "c = " << tpr->eta << '\n';
1888                         } else if((tpr->leg == CROSSWIND) || (tpr->leg == TURN1)) {
1889                                 //cout << "CROSSWIND calc: ";
1890                                 //cout << tpr->eta << ' ';
1891                                 if(dist_across_m > nominal_dist_across_m) {
1892                                         tpr->eta += dist_across_m / circuit_ias;
1893                                         //cout << "a ";
1894                                 } else {
1895                                         tpr->eta += nominal_dist_across_m / circuit_ias;
1896                                         //cout << "b ";
1897                                 }
1898                                 //cout << tpr->eta << ' ';
1899                                 // should we use the dist across of the previous plane if there is previous still on downwind?
1900                                 //if(printout) cout << "bb = " << tpr->eta << '\n';
1901                                 if(dist_out_m > nominal_cross_dist_out_m) {
1902                                         tpr->eta += fabs(current_base_dist_out_m - dist_out_m) / circuit_ias;
1903                                         //cout << "c ";
1904                                 } else {
1905                                         tpr->eta += fabs(current_base_dist_out_m - nominal_cross_dist_out_m) / circuit_ias;
1906                                         //cout << "d ";
1907                                 }
1908                                 //cout << tpr->eta << ' ';
1909                                 //if(printout) cout << "cc = " << tpr->eta << '\n';
1910                                 if(nominal_dist_across_m > dist_across_m) {
1911                                         tpr->eta += (nominal_dist_across_m - dist_across_m) / circuit_ias;
1912                                         //cout << "e ";
1913                                 } else {
1914                                         // Nothing to add
1915                                         //cout << "f ";
1916                                 }
1917                                 //cout << tpr->eta << '\n';
1918                                 //if(printout) cout << "dd = " << tpr->eta << '\n';
1919                         } else {
1920                                 // We've only just started - why not use a generic estimate?
1921                                 tpr->eta = 240.0;
1922                         }
1923                 }
1924                 //if(printout) {
1925                 //      cout << "ETA = " << tpr->eta << '\n';
1926                 //}
1927                 //if(!tpr->isUser) cout << tpr->plane.callsign << '\t' << tpr->eta << '\n';
1928         } else {
1929                 tpr->eta = 99999;
1930         }       
1931 }
1932
1933
1934 // Calculate the distance of a plane to the threshold in meters
1935 // TODO - Modify to calculate flying distance of a plane in the circuit
1936 double FGTower::CalcDistOutM(TowerPlaneRec* tpr) {
1937         return(dclGetHorizontalSeparation(rwy.threshold_pos, tpr->pos));
1938 }
1939
1940
1941 // Calculate the distance of a plane to the threshold in miles
1942 // TODO - Modify to calculate flying distance of a plane in the circuit
1943 double FGTower::CalcDistOutMiles(TowerPlaneRec* tpr) {
1944         return(CalcDistOutM(tpr) / 1600.0);             // FIXME - use a proper constant if possible.
1945 }
1946
1947
1948 // Iterate through all the lists, update the position of, and call CalcETA for all the planes.
1949 void FGTower::doThresholdETACalc() {
1950         //cout << "Entering doThresholdETACalc..." << endl;
1951         tower_plane_rec_list_iterator twrItr;
1952         // Do the approach list first
1953         for(twrItr = appList.begin(); twrItr != appList.end(); twrItr++) {
1954                 TowerPlaneRec* tpr = *twrItr;
1955                 if(!(tpr->isUser)) tpr->pos = tpr->planePtr->GetPos();
1956                 //cout << "APP: ";
1957                 CalcETA(tpr);
1958         }       
1959         // Then the circuit list
1960         //cout << "Circuit list size is " << circuitList.size() << '\n';
1961         for(twrItr = circuitList.begin(); twrItr != circuitList.end(); twrItr++) {
1962                 TowerPlaneRec* tpr = *twrItr;
1963                 if(!(tpr->isUser)) tpr->pos = tpr->planePtr->GetPos();
1964                 //cout << "CIRC: ";
1965                 CalcETA(tpr);
1966         }
1967         //cout << "Done doThresholdETCCalc" << endl;
1968 }
1969                 
1970
1971 // Check that the planes in traffic list are correctly ordered,
1972 // that the nearest (timewise) is flagged next on rwy, and return
1973 // true if any threshold use conflicts are detected, false otherwise.
1974 bool FGTower::doThresholdUseOrder() {
1975         //cout << "Entering doThresholdUseOrder..." << endl;
1976         bool conflict = false;
1977         
1978         // Wipe out traffic list, go through circuit, app and hold list, and reorder them in traffic list.
1979         // Here's the rather simplistic assumptions we're using:
1980         // Currently all planes are assumed to be GA light singles with corresponding speeds and separation times.
1981         // In order of priority for runway use:
1982         // STRAIGHT_IN > CIRCUIT > HOLDING_FOR_DEPARTURE
1983         // No modification of planes speeds occurs - conflicts are resolved by delaying turn for base,
1984         // and holding planes until a space.
1985         // When calculating if a holding plane can use the runway, time clearance from last departure
1986         // as well as time clearance to next arrival must be considered.
1987         
1988         trafficList.clear();
1989         
1990         tower_plane_rec_list_iterator twrItr;
1991         // Do the approach list first
1992         //if(ident == "KRHV") cout << "A" << flush;
1993         for(twrItr = appList.begin(); twrItr != appList.end(); twrItr++) {
1994                 TowerPlaneRec* tpr = *twrItr;
1995                 //if(ident == "KRHV") cout << tpr->plane.callsign << '\n';
1996                 conflict = AddToTrafficList(tpr);
1997         }       
1998         // Then the circuit list
1999         //if(ident == "KRHV") cout << "C" << flush;
2000         for(twrItr = circuitList.begin(); twrItr != circuitList.end(); twrItr++) {
2001                 TowerPlaneRec* tpr = *twrItr;
2002                 //if(ident == "KRHV") cout << tpr->plane.callsign << '\n';
2003                 conflict = AddToTrafficList(tpr);
2004         }
2005         // And finally the hold list
2006         //cout << "H" << endl;
2007         for(twrItr = holdList.begin(); twrItr != holdList.end(); twrItr++) {
2008                 TowerPlaneRec* tpr = *twrItr;
2009                 AddToTrafficList(tpr, true);
2010         }
2011         
2012         
2013         if(0) {
2014         //if(ident == "KRHV") {
2015                 cout << "T\n";
2016                 for(twrItr = trafficList.begin(); twrItr != trafficList.end(); twrItr++) {
2017                         TowerPlaneRec* tpr = *twrItr;
2018                         cout << tpr->plane.callsign << '\t' << tpr->eta << '\t';
2019                 }
2020                 cout << endl;
2021         }
2022         
2023         //cout << "Done doThresholdUseOrder" << endl;
2024         return(conflict);
2025 }
2026
2027
2028 // Return the ETA of plane no. list_pos (1-based) in the traffic list.
2029 // i.e. list_pos = 1 implies next to use runway.
2030 double FGTower::GetTrafficETA(unsigned int list_pos, bool printout) {
2031         if(trafficList.size() < list_pos) {
2032                 return(99999);
2033         }
2034
2035         tower_plane_rec_list_iterator twrItr;
2036         twrItr = trafficList.begin();
2037         for(unsigned int i = 1; i < list_pos; i++, twrItr++);
2038         TowerPlaneRec* tpr = *twrItr;
2039         CalcETA(tpr, printout);
2040         //cout << "ETA returned = " << tpr->eta << '\n';
2041         return(tpr->eta);
2042 }
2043         
2044
2045 void FGTower::ContactAtHoldShort(const PlaneRec& plane, FGAIPlane* requestee, tower_traffic_type operation) {
2046         // HACK - assume that anything contacting at hold short is new for now - FIXME LATER
2047         TowerPlaneRec* t = new TowerPlaneRec;
2048         t->plane = plane;
2049         t->planePtr = requestee;
2050         t->holdShortReported = true;
2051         t->clearedToLineUp = false;
2052         t->clearedToTakeOff = false;
2053         t->opType = operation;
2054         t->pos = requestee->GetPos();
2055         
2056         //cout << "Hold Short reported by " << plane.callsign << '\n';
2057         SG_LOG(SG_ATC, SG_BULK, "Hold Short reported by " << plane.callsign);
2058
2059 /*      
2060         bool next = AddToTrafficList(t, true);
2061         if(next) {
2062                 double teta = GetTrafficETA(2);
2063                 if(teta < 150.0) {
2064                         t->clearanceCounter = 7.0;      // This reduces the delay before response to 3 secs if an immediate takeoff is reqd
2065                         //cout << "Reducing response time to request due imminent traffic\n";
2066                 }
2067         } else {
2068         }
2069 */
2070         // TODO - possibly add the reduced interval to clearance when immediate back in under the new scheme
2071
2072         holdList.push_back(t);
2073         
2074         responseReqd = true;
2075 }
2076
2077 // Register the presence of an AI plane at a point where contact would already have been made in real life
2078 // CAUTION - currently it is assumed that this plane's callsign is unique - it is up to AIMgr to generate unique callsigns.
2079 void FGTower::RegisterAIPlane(const PlaneRec& plane, FGAIPlane* ai, const tower_traffic_type& op, const PatternLeg& lg) {
2080         // At the moment this is only going to be tested with inserting an AI plane on downwind
2081         TowerPlaneRec* t = new TowerPlaneRec;
2082         t->plane = plane;
2083         t->planePtr = ai;
2084         t->opType = op;
2085         t->leg = lg;
2086         t->pos = ai->GetPos();
2087         
2088         CalcETA(t);
2089         
2090         if(op == CIRCUIT && lg != LEG_UNKNOWN) {
2091                 AddToCircuitList(t);
2092         } else {
2093                 // FLAG A WARNING
2094         }
2095         
2096         doThresholdUseOrder();
2097 }
2098
2099 void FGTower::DeregisterAIPlane(const string& id) {
2100         RemovePlane(id);
2101 }
2102
2103 // Contact tower for VFR approach
2104 // eg "Cessna Charlie Foxtrot Golf Foxtrot Sierra eight miles South of the airport for full stop with Bravo"
2105 // This function probably only called via user interaction - AI planes will have an overloaded function taking a planerec.
2106 // opt defaults to AIP_LT_UNKNOWN
2107 void FGTower::VFRArrivalContact(const string& ID, const LandingType& opt) {
2108         //cout << "USER Request Landing Clearance called for ID " << ID << '\n';
2109         
2110         // For now we'll assume that the user is a light plane and can get him/her to join the circuit if necessary.
2111
2112         TowerPlaneRec* t;       
2113         string usercall = fgGetString("/sim/user/callsign");
2114         if(ID == "USER" || ID == usercall) {
2115                 t = FindPlane(usercall);
2116                 if(!t) {
2117                         //cout << "NOT t\n";
2118                         t = new TowerPlaneRec;
2119                         t->isUser = true;
2120                         t->pos.setlon(user_lon_node->getDoubleValue());
2121                         t->pos.setlat(user_lat_node->getDoubleValue());
2122                         t->pos.setelev(user_elev_node->getDoubleValue());
2123                 } else {
2124                         //cout << "IS t\n";
2125                         // Oops - the plane is already registered with this tower - maybe we took off and flew a giant circuit without
2126                         // quite getting out of tower airspace - just ignore for now and treat as new arrival.
2127                         // TODO - Maybe should remove from departure and circuit list if in there though!!
2128                 }
2129         } else {
2130                 // Oops - something has gone wrong - put out a warning
2131                 cout << "WARNING - FGTower::VFRContact(string ID, LandingType lt) called with ID " << ID << " which does not appear to be the user.\n";
2132                 return;
2133         }
2134                 
2135         
2136         // TODO
2137         // Calculate where the plane is in relation to the active runway and it's circuit
2138         // and set the op-type as appropriate.
2139         
2140         // HACK - to get up and running I'm going to assume that the user contacts tower on a staight-in final for now.
2141         t->opType = STRAIGHT_IN;
2142         
2143         t->plane.type = GA_SINGLE;      // FIXME - Another assumption!
2144         t->plane.callsign = usercall;
2145         
2146         t->vfrArrivalReported = true;
2147         responseReqd = true;
2148         
2149         appList.push_back(t);   // Not necessarily permanent
2150         appListItr = appList.begin();
2151         AddToTrafficList(t);
2152         
2153         current_atcdialog->remove_entry(ident, USER_REQUEST_VFR_ARRIVAL, TOWER);
2154         current_atcdialog->remove_entry(ident, USER_REQUEST_VFR_ARRIVAL_FULL_STOP, TOWER);
2155         current_atcdialog->remove_entry(ident, USER_REQUEST_VFR_ARRIVAL_TOUCH_AND_GO, TOWER);
2156 }
2157
2158 // landingType defaults to AIP_LT_UNKNOWN
2159 void FGTower::VFRArrivalContact(const PlaneRec& plane, FGAIPlane* requestee, const LandingType& lt) {
2160         //cout << "VFRArrivalContact called for plane " << plane.callsign << " at " << ident << '\n';
2161         // Possible hack - assume this plane is new for now - TODO - should check really
2162         TowerPlaneRec* t = new TowerPlaneRec;
2163         t->plane = plane;
2164         t->planePtr = requestee;
2165         t->landingType = lt;
2166         t->pos = requestee->GetPos();
2167         
2168         //cout << "Hold Short reported by " << plane.callsign << '\n';
2169         SG_LOG(SG_ATC, SG_BULK, "VFR arrival contact made by " << plane.callsign);
2170         //cout << "VFR arrival contact made by " << plane.callsign << '\n';
2171
2172         // HACK - to get up and running I'm going to assume a staight-in final for now.
2173         t->opType = STRAIGHT_IN;
2174         
2175         t->vfrArrivalReported = true;
2176         responseReqd = true;
2177         
2178         //cout << "Before adding, appList.size = " << appList.size() << " at " << ident << '\n';
2179         appList.push_back(t);   // Not necessarily permanent
2180         appListItr = appList.begin();
2181         //cout << "After adding, appList.size = " << appList.size() << " at " << ident << '\n';
2182         AddToTrafficList(t);
2183 }
2184
2185 void FGTower::RequestDepartureClearance(const string& ID) {
2186         //cout << "Request Departure Clearance called...\n";
2187 }
2188
2189 void FGTower::RequestTakeOffClearance(const string& ID) {
2190         string uid=ID;
2191         if(ID == "USER") {
2192                 uid = fgGetString("/sim/user/callsign");
2193                 current_atcdialog->remove_entry(ident, USER_REQUEST_TAKE_OFF, TOWER);
2194         }
2195         TowerPlaneRec* t = FindPlane(uid);
2196         if(t) {
2197                 if(!(t->clearedToTakeOff)) {
2198                         departed = false;
2199                         t->lineUpReported=true;
2200                         responseReqd = true;
2201                 }
2202         }
2203         else {
2204                 SG_LOG(SG_ATC, SG_WARN, "WARNING: Unable to find plane " << ID << " in FGTower::RequestTakeOffClearance(...)");
2205         }
2206 }
2207         
2208 void FGTower::ReportFinal(const string& ID) {
2209         //cout << "Report Final Called at tower " << ident << " by plane " << ID << '\n';
2210         string uid=ID;
2211         if(ID == "USER") {
2212                 uid = fgGetString("/sim/user/callsign");
2213                 current_atcdialog->remove_entry(ident, USER_REPORT_3_MILE_FINAL, TOWER);
2214         }
2215         TowerPlaneRec* t = FindPlane(uid);
2216         if(t) {
2217                 t->finalReported = true;
2218                 t->finalAcknowledged = false;
2219                 if(!(t->clearedToLand)) {
2220                         responseReqd = true;
2221                 } else {
2222                         // possibly respond with wind even if already cleared to land?
2223                         t->finalReported = false;
2224                         t->finalAcknowledged = true;
2225                         // HACK!! - prevents next reporting being misinterpreted as this one.
2226                 }
2227         } else {
2228                 SG_LOG(SG_ATC, SG_WARN, "WARNING: Unable to find plane " << ID << " in FGTower::ReportFinal(...)");
2229         }
2230 }
2231
2232 void FGTower::ReportLongFinal(const string& ID) {
2233         string uid=ID;
2234         if(ID == "USER") {
2235                 uid = fgGetString("/sim/user/callsign");
2236                 current_atcdialog->remove_entry(ident, USER_REPORT_3_MILE_FINAL, TOWER);
2237         }
2238         TowerPlaneRec* t = FindPlane(uid);
2239         if(t) {
2240                 t->longFinalReported = true;
2241                 t->longFinalAcknowledged = false;
2242                 if(!(t->clearedToLand)) {
2243                         responseReqd = true;
2244                 } // possibly respond with wind even if already cleared to land?
2245         } else {
2246                 SG_LOG(SG_ATC, SG_WARN, "WARNING: Unable to find plane " << ID << " in FGTower::ReportLongFinal(...)");
2247         }
2248 }
2249
2250 //void FGTower::ReportOuterMarker(string ID);
2251 //void FGTower::ReportMiddleMarker(string ID);
2252 //void FGTower::ReportInnerMarker(string ID);
2253
2254 void FGTower::ReportRunwayVacated(const string& ID) {
2255         //cout << "Report Runway Vacated Called at tower " << ident << " by plane " << ID << '\n';
2256         string uid=ID;
2257         if(ID == "USER") {
2258                 uid = fgGetString("/sim/user/callsign");
2259                 current_atcdialog->remove_entry(ident, USER_REPORT_RWY_VACATED, TOWER);
2260         }
2261         TowerPlaneRec* t = FindPlane(uid);
2262         if(t) {
2263                 //cout << "Found it...\n";
2264                 t->rwyVacatedReported = true;
2265                 responseReqd = true;
2266         } else {
2267                 SG_LOG(SG_ATC, SG_WARN, "WARNING: Unable to find plane " << ID << " in FGTower::ReportRunwayVacated(...)");
2268                 SG_LOG(SG_ATC, SG_ALERT, "A WARNING: Unable to find plane " << ID << " in FGTower::ReportRunwayVacated(...)");
2269                 //cout << "WARNING: Unable to find plane " << ID << " in FGTower::ReportRunwayVacated(...)\n";
2270         }
2271 }
2272
2273 TowerPlaneRec* FGTower::FindPlane(const string& ID) {
2274         //cout << "FindPlane called for " << ID << "...\n";
2275         tower_plane_rec_list_iterator twrItr;
2276         // Do the approach list first
2277         for(twrItr = appList.begin(); twrItr != appList.end(); twrItr++) {
2278                 //cout << "appList callsign is " << (*twrItr)->plane.callsign << '\n';
2279                 if((*twrItr)->plane.callsign == ID) return(*twrItr);
2280         }       
2281         // Then the circuit list
2282         for(twrItr = circuitList.begin(); twrItr != circuitList.end(); twrItr++) {
2283                 //cout << "circuitList callsign is " << (*twrItr)->plane.callsign << '\n';
2284                 if((*twrItr)->plane.callsign == ID) return(*twrItr);
2285         }
2286         // Then the runway list
2287         //cout << "rwyList.size() is " << rwyList.size() << '\n';
2288         for(twrItr = rwyList.begin(); twrItr != rwyList.end(); twrItr++) {
2289                 //cout << "rwyList callsign is " << (*twrItr)->plane.callsign << '\n';
2290                 if((*twrItr)->plane.callsign == ID) return(*twrItr);
2291         }
2292         // The hold list
2293         for(twrItr = holdList.begin(); twrItr != holdList.end(); twrItr++) {
2294                 if((*twrItr)->plane.callsign == ID) return(*twrItr);
2295         }
2296         // And finally the vacated list
2297         for(twrItr = vacatedList.begin(); twrItr != vacatedList.end(); twrItr++) {
2298                 //cout << "vacatedList callsign is " << (*twrItr)->plane.callsign << '\n';
2299                 if((*twrItr)->plane.callsign == ID) return(*twrItr);
2300         }
2301         SG_LOG(SG_ATC, SG_WARN, "Unable to find " << ID << " in FGTower::FindPlane(...)");
2302         //exit(-1);
2303         return(NULL);
2304 }
2305
2306 void FGTower::RemovePlane(const string& ID) {
2307         //cout << ident << " RemovePlane called for " << ID << '\n';
2308         // We have to be careful here - we want to erase the plane from all lists it is in,
2309         // but we can only delete it once, AT THE END.
2310         TowerPlaneRec* t = NULL;
2311         tower_plane_rec_list_iterator twrItr;
2312         for(twrItr = appList.begin(); twrItr != appList.end();) {
2313                 if((*twrItr)->plane.callsign == ID) {
2314                         t = *twrItr;
2315                         twrItr = appList.erase(twrItr);
2316                         appListItr = appList.begin();
2317                         // HACK: aircraft are sometimes more than once in a list, so we need to
2318                         // remove them all before we can delete the TowerPlaneRec class
2319                         //break;
2320                 } else
2321                         ++twrItr;
2322         }
2323         for(twrItr = depList.begin(); twrItr != depList.end();) {
2324                 if((*twrItr)->plane.callsign == ID) {
2325                         t = *twrItr;
2326                         twrItr = depList.erase(twrItr);
2327                         depListItr = depList.begin();
2328                 } else
2329                         ++twrItr;
2330         }
2331         for(twrItr = circuitList.begin(); twrItr != circuitList.end();) {
2332                 if((*twrItr)->plane.callsign == ID) {
2333                         t = *twrItr;
2334                         twrItr = circuitList.erase(twrItr);
2335                         circuitListItr = circuitList.begin();
2336                 } else
2337                         ++twrItr;
2338         }
2339         for(twrItr = holdList.begin(); twrItr != holdList.end();) {
2340                 if((*twrItr)->plane.callsign == ID) {
2341                         t = *twrItr;
2342                         twrItr = holdList.erase(twrItr);
2343                         holdListItr = holdList.begin();
2344                 } else
2345                         ++twrItr;
2346         }
2347         for(twrItr = rwyList.begin(); twrItr != rwyList.end();) {
2348                 if((*twrItr)->plane.callsign == ID) {
2349                         t = *twrItr;
2350                         twrItr = rwyList.erase(twrItr);
2351                         rwyListItr = rwyList.begin();
2352                 } else
2353                         ++twrItr;
2354         }
2355         for(twrItr = vacatedList.begin(); twrItr != vacatedList.end();) {
2356                 if((*twrItr)->plane.callsign == ID) {
2357                         t = *twrItr;
2358                         twrItr = vacatedList.erase(twrItr);
2359                         vacatedListItr = vacatedList.begin();
2360                 } else
2361                         ++twrItr;
2362         }
2363         for(twrItr = trafficList.begin(); twrItr != trafficList.end();) {
2364                 if((*twrItr)->plane.callsign == ID) {
2365                         t = *twrItr;
2366                         twrItr = trafficList.erase(twrItr);
2367                         trafficListItr = trafficList.begin();
2368                 } else
2369                         ++twrItr;
2370         }
2371         // And finally, delete the record.
2372         delete t;
2373 }
2374
2375 void FGTower::ReportDownwind(const string& ID) {
2376         //cout << "ReportDownwind(...) called\n";
2377         string uid=ID;
2378         if(ID == "USER") {
2379                 uid = fgGetString("/sim/user/callsign");
2380                 current_atcdialog->remove_entry(ident, USER_REPORT_DOWNWIND, TOWER);
2381         }
2382         TowerPlaneRec* t = FindPlane(uid);
2383         if(t) {
2384                 t->downwindReported = true;
2385                 responseReqd = true;
2386                 // If the plane is in the app list, remove it and put it in the circuit list instead.
2387                 // Ideally we might want to do this at the 2 mile report prior to 45 deg entry, but at
2388                 // the moment that would b&gg?r up the constraint position calculations.
2389                 RemoveFromAppList(ID);
2390                 t->leg = DOWNWIND;
2391                 if(t->isUser) {
2392                         t->pos.setlon(user_lon_node->getDoubleValue());
2393                         t->pos.setlat(user_lat_node->getDoubleValue());
2394                         t->pos.setelev(user_elev_node->getDoubleValue());
2395                 } else {
2396                         // ASSERT(t->planePtr != NULL);
2397                         t->pos = t->planePtr->GetPos();
2398                 }
2399                 CalcETA(t);
2400                 AddToCircuitList(t);
2401         } else {
2402                 SG_LOG(SG_ATC, SG_WARN, "WARNING: Unable to find plane " << ID << " in FGTower::ReportDownwind(...)");
2403         }
2404 }
2405
2406 void FGTower::ReportGoingAround(const string& ID) {
2407         string uid=ID;
2408         if(ID == "USER") {
2409                 uid = fgGetString("/sim/user/callsign");
2410                 RemoveAllUserDialogOptions();   // TODO - it would be much more efficient if ATCDialog simply had a clear() function!!!
2411                 current_atcdialog->add_entry(ident, "@AP Tower, @CS Downwind @RW", "Report Downwind", TOWER, (int)USER_REPORT_DOWNWIND);
2412         }
2413         TowerPlaneRec* t = FindPlane(uid);
2414         if(t) {
2415                 //t->goAroundReported = true;  // No need to set this until we start responding to it.
2416                 responseReqd = false;   // might change in the future but for now we'll not distract them during the go-around.
2417                 // If the plane is in the app list, remove it and put it in the circuit list instead.
2418                 RemoveFromAppList(ID);
2419                 t->leg = CLIMBOUT;
2420                 if(t->isUser) {
2421                         t->pos.setlon(user_lon_node->getDoubleValue());
2422                         t->pos.setlat(user_lat_node->getDoubleValue());
2423                         t->pos.setelev(user_elev_node->getDoubleValue());
2424                 } else {
2425                         // ASSERT(t->planePtr != NULL);
2426                         t->pos = t->planePtr->GetPos();
2427                 }
2428                 CalcETA(t);
2429                 AddToCircuitList(t);
2430         } else {
2431                 SG_LOG(SG_ATC, SG_WARN, "WARNING: Unable to find plane " << ID << " in FGTower::ReportDownwind(...)");
2432         }
2433 }
2434
2435 string FGTower::GenText(const string& m, int c) {
2436         const int cmax = 300;
2437         //string message;
2438         char tag[4];
2439         char crej = '@';
2440         char mes[cmax];
2441         char dum[cmax];
2442         //char buf[10];
2443         char *pos;
2444         int len;
2445         //FGTransmission t;
2446         string usercall = fgGetString("/sim/user/callsign");
2447         TowerPlaneRec* t = FindPlane(responseID);
2448         
2449         //transmission_list_type     tmissions = transmissionlist_station[station];
2450         //transmission_list_iterator current   = tmissions.begin();
2451         //transmission_list_iterator last      = tmissions.end();
2452         
2453         //for ( ; current != last ; ++current ) {
2454         //      if ( current->get_code().c1 == code.c1 &&  
2455         //              current->get_code().c2 == code.c2 &&
2456         //          current->get_code().c3 == code.c3 ) {
2457                         
2458                         //if ( ttext ) message = current->get_transtext();
2459                         //else message = current->get_menutext();
2460                         strcpy( &mes[0], m.c_str() ); 
2461                         
2462                         // Replace all the '@' parameters with the actual text.
2463                         int check = 0;  // If mes gets overflowed the while loop can go infinite
2464                         double gp = fgGetFloat("/gear/gear/position-norm");
2465                         while ( strchr(&mes[0], crej) != NULL  ) {      // ie. loop until no more occurances of crej ('@') found
2466                                 pos = strchr( &mes[0], crej );
2467                                 memmove(&tag[0], pos, 3);
2468                                 tag[3] = '\0';
2469                                 int i;
2470                                 len = 0;
2471                                 for ( i=0; i<cmax; i++ ) {
2472                                         if ( mes[i] == crej ) {
2473                                                 len = i; 
2474                                                 break;
2475                                         }
2476                                 }
2477                                 strncpy( &dum[0], &mes[0], len );
2478                                 dum[len] = '\0';
2479                                 
2480                                 if ( strcmp ( tag, "@ST" ) == 0 )
2481                                         //strcat( &dum[0], tpars.station.c_str() );
2482                                         strcat(&dum[0], ident.c_str());
2483                                 else if ( strcmp ( tag, "@AP" ) == 0 )
2484                                         //strcat( &dum[0], tpars.airport.c_str() );
2485                                         strcat(&dum[0], name.c_str());
2486                                 else if ( strcmp ( tag, "@CS" ) == 0 ) 
2487                                         //strcat( &dum[0], tpars.callsign.c_str() );
2488                                         strcat(&dum[0], usercall.c_str());
2489                                 else if ( strcmp ( tag, "@TD" ) == 0 ) {
2490                                         /*
2491                                         if ( tpars.tdir == 1 ) {
2492                                                 char buf[] = "left";
2493                                                 strcat( &dum[0], &buf[0] );
2494                                         }
2495                                         else {
2496                                                 char buf[] = "right";
2497                                                 strcat( &dum[0], &buf[0] );
2498                                         }
2499                                         */
2500                                 }
2501                                 else if ( strcmp ( tag, "@HE" ) == 0 ) {
2502                                         /*
2503                                         char buf[10];
2504                                         sprintf( buf, "%i", (int)(tpars.heading) );
2505                                         strcat( &dum[0], &buf[0] );
2506                                         */
2507                                 }
2508                                 else if ( strcmp ( tag, "@AT" ) == 0 ) {        // ATIS ID
2509                                         /*
2510                                         char buf[10];
2511                                         sprintf( buf, "%i", (int)(tpars.heading) );
2512                                         strcat( &dum[0], &buf[0] );
2513                                         */
2514                                         double f = globals->get_ATC_mgr()->GetFrequency(ident, ATIS) / 100.0;
2515                                         if(f) {
2516                                                 string atis_id;
2517                                                 atis_id = ", information " + GetATISID();
2518                                                 strcat( &dum[0], atis_id.c_str() );
2519                                         }
2520                                 }
2521                                 else if ( strcmp ( tag, "@VD" ) == 0 ) {
2522                                         /*
2523                                         if ( tpars.VDir == 1 ) {
2524                                                 char buf[] = "Descend and maintain";
2525                                                 strcat( &dum[0], &buf[0] );
2526                                         }
2527                                         else if ( tpars.VDir == 2 ) {
2528                                                 char buf[] = "Maintain";
2529                                                 strcat( &dum[0], &buf[0] );
2530                                         }
2531                                         else if ( tpars.VDir == 3 ) {
2532                                                 char buf[] = "Climb and maintain";
2533                                                 strcat( &dum[0], &buf[0] );
2534                                         } 
2535                                         */
2536                                 }
2537                                 else if ( strcmp ( tag, "@AL" ) == 0 ) {
2538                                         /*
2539                                         char buf[10];
2540                                         sprintf( buf, "%i", (int)(tpars.alt) );
2541                                         strcat( &dum[0], &buf[0] );
2542                                         */
2543                                 }
2544                                 else if ( strcmp ( tag, "@TO" ) == 0 ) {      // Requesting take-off or departure clearance
2545                                         string tmp;
2546                                         if (rwyOccupied) {
2547                                                 tmp = "Ready for take-off";
2548                                         } else {
2549                                                 if (OnAnyRunway(Point3D(user_lon_node->getDoubleValue(),
2550                                                                 user_lat_node->getDoubleValue(), 0.0),true)) {
2551                                                         tmp = "Request take-off clearance";
2552                                                 } else {
2553                                                         tmp = "Request departure clearance";
2554                                                 }
2555                                         }
2556                                         strcat(&dum[0], tmp.c_str());
2557                                 }
2558                                 else if ( strcmp ( tag, "@MI" ) == 0 ) {
2559                                         char buf[10];
2560                                         //sprintf( buf, "%3.1f", tpars.miles );
2561                                         int dist_miles = (int)dclGetHorizontalSeparation(Point3D(lon, lat, elev), Point3D(user_lon_node->getDoubleValue(), user_lat_node->getDoubleValue(), user_elev_node->getDoubleValue())) / 1600;
2562                                         sprintf(buf, "%i", dist_miles);
2563                                         strcat( &dum[0], &buf[0] );
2564                                 }
2565                                 else if ( strcmp ( tag, "@FR" ) == 0 ) {
2566                                         /*
2567                                         char buf[10];
2568                                         sprintf( buf, "%6.2f", tpars.freq );
2569                                         strcat( &dum[0], &buf[0] );
2570                                         */
2571                                 }
2572                                 else if ( strcmp ( tag, "@RW" ) == 0 ) {
2573                                         strcat(&dum[0], ConvertRwyNumToSpokenString(activeRwy).c_str());
2574                                 }
2575                                 else if ( strcmp ( tag, "@GR" ) == 0 ) {        // Gear position (on final)
2576                                         if(t->gearWasUp && gp > 0.99) {
2577                                                 strcat(&dum[0], ", gear down, ready to land.");
2578                                         }
2579                                 }
2580                                 else if(strcmp(tag, "@CD") == 0) {      // @CD = compass direction
2581                                         double h = GetHeadingFromTo(Point3D(lon, lat, elev), Point3D(user_lon_node->getDoubleValue(), user_lat_node->getDoubleValue(), user_elev_node->getDoubleValue()));
2582                                         while(h < 0.0) h += 360.0;
2583                                         while(h > 360.0) h -= 360.0;
2584                                         if(h < 22.5 || h > 337.5) {
2585                                                 strcat(&dum[0], "North");
2586                                         } else if(h < 67.5) {
2587                                                 strcat(&dum[0], "North-East");
2588                                         } else if(h < 112.5) {
2589                                                 strcat(&dum[0], "East");
2590                                         } else if(h < 157.5) {
2591                                                 strcat(&dum[0], "South-East");
2592                                         } else if(h < 202.5) {
2593                                                 strcat(&dum[0], "South");
2594                                         } else if(h < 247.5) {
2595                                                 strcat(&dum[0], "South-West");
2596                                         } else if(h < 292.5) {
2597                                                 strcat(&dum[0], "West");
2598                                         } else {
2599                                                 strcat(&dum[0], "North-West");
2600                                         }
2601                                 } else {
2602                                         cout << "Tag " << tag << " not found" << endl;
2603                                         break;
2604                                 }
2605                                 strcat( &dum[0], &mes[len+3] );
2606                                 strcpy( &mes[0], &dum[0] );
2607                                 
2608                                 ++check;
2609                                 if(check > 10) {
2610                                         SG_LOG(SG_ATC, SG_WARN, "WARNING: Possibly endless loop terminated in FGTransmissionlist::gen_text(...)"); 
2611                                         break;
2612                                 }
2613                         }
2614                         
2615                         //cout << mes  << endl;  
2616                         //break;
2617                 //}
2618         //}
2619         return mes[0] ? mes : "No transmission found";
2620 }
2621
2622 string FGTower::GetWeather() {
2623         std::ostringstream msg;
2624
2625         // wind
2626         double hdg = wind_from_hdg->getDoubleValue();
2627         double speed = wind_speed_knots->getDoubleValue();
2628         if (speed < 1)
2629                 msg << "wind calm";
2630         else
2631                 msg << "wind " << int(hdg) << " degrees at " << int(speed) << " knots";
2632
2633         // visibility
2634         double visibility = fgGetDouble("/environment/visibility-m");
2635         if (visibility < 10000)
2636                 msg << ", visibility " << int(visibility / 1609) << " miles";
2637
2638         // pressure / altimeter
2639         double pressure = fgGetDouble("/environment/pressure-sea-level-inhg");
2640         msg << ", QFE " << fixed << setprecision(2) << pressure << ".";
2641
2642         return msg.str();
2643 }
2644
2645 string FGTower::GetATISID() {
2646         int hours = fgGetInt("/sim/time/utc/hour");
2647         int phonetic_id = current_commlist->GetCallSign(ident, hours, 0);
2648         return GetPhoneticIdent(phonetic_id);
2649 }
2650
2651 ostream& operator << (ostream& os, tower_traffic_type ttt) {
2652         switch(ttt) {
2653         case(CIRCUIT):      return(os << "CIRCUIT");
2654         case(INBOUND):      return(os << "INBOUND");
2655         case(OUTBOUND):     return(os << "OUTBOUND");
2656         case(TTT_UNKNOWN):  return(os << "UNKNOWN");
2657         case(STRAIGHT_IN):  return(os << "STRAIGHT_IN");
2658         }
2659         return(os << "ERROR - Unknown switch in tower_traffic_type operator << ");
2660 }
2661