getAltitudeMeanSeaLevelMeters(),
reverse=True)
arrivalTurnLeg.buildTurnLeg()
threeDegreesDescentGlideSlope.addGraph(arrivalTurnLeg)
print(' ================== great circle =============== ')
t1 = time.clock()
print("time end= ", t1, "duration= {0:3.10f} seconds".format(t1 - t0))
climbSlopeLastVertex = climbSlope.getVertex(
climbSlope.getNumberOfVertices() - 1)
climbSlopeLastWayPoint = climbSlopeLastVertex.getWeight()
print(climbSlopeLastWayPoint)
arrivalTurnLegFirstVertex = arrivalTurnLeg.getVertex(0)
arrivalTurnLegFirstWayPoint = arrivalTurnLegFirstVertex.getWeight()
finalGreatCircle = GreatCircleRoute(climbSlopeLastWayPoint,
arrivalTurnLegFirstWayPoint)
finalGreatCircle.computeGreatCircle()
print('great circle length= ' + str(greatCircle.computeLengthMeters()) +
' meters')
print(' ============= buildClimbRamp the final graph =================')
climbSlope.addGraph(finalGreatCircle)
climbSlope.addGraph(threeDegreesDescentGlideSlope)
print('route length= ' + str(climbSlope.computeLengthMeters()) + ' meters')
climbSlope.createKmlOutputFile()
t1 = time.clock()
def buildSimulatedArrivalPhase(self):
print(
self.className +
'=========== add final turn, descent and ground run ==================='
)
arrivalGroundRun = GroundRunLeg(runway=self.arrivalRunway,
aircraft=self.aircraft,
airport=self.arrivalAirport)
self.touchDownWayPoint = arrivalGroundRun.computeTouchDownWayPoint()
# add touch down to constraint list
self.constraintsList.append(
ArrivalRunWayTouchDownConstraint(self.touchDownWayPoint))
print(self.touchDownWayPoint)
''' distance from last fix to touch down '''
distanceToLastFixNautics = self.touchDownWayPoint.getDistanceMetersTo(
self.getLastWayPoint()) * Meter2NauticalMiles
print(
self.className +
'===================== final 3 degrees descending glide slope ================'
)
descentGlideSlope = DescentGlideSlope(
runway=self.arrivalRunway,
aircraft=self.aircraft,
arrivalAirport=self.arrivalAirport,
descentGlideSlopeDegrees=3.0)
''' if there is a fix nearer to 5 nautics of the touch-down then limit size of simulated glide slope '''
descentGlideSlopeSizeNautics = min(distanceToLastFixNautics / 2.0, 5.0)
''' build simulated glide slope '''
descentGlideSlope.buildSimulatedGlideSlope(
descentGlideSlopeSizeNautics)
self.firstGlideSlopeWayPoint = descentGlideSlope.getVertex(
v=0).getWeight()
print(self.className + ': top of arrival glide slope= {0}'.format(
self.firstGlideSlopeWayPoint))
print(
self.className +
' ================= need a turn leg to find the junction point the last way-point in the fix list to the top of the final glide slope'
)
'''
initial heading is the orientation of the run-way
'''
lastFixListWayPoint = self.wayPointsDict[self.fixList[-1]]
initialHeadingDegrees = self.arrivalRunway.getTrueHeadingDegrees()
print("=====> arrival runway - true heading degrees = {0}".format(
initialHeadingDegrees))
lastTurnLeg = TurnLeg(initialWayPoint=self.firstGlideSlopeWayPoint,
finalWayPoint=lastFixListWayPoint,
initialHeadingDegrees=initialHeadingDegrees,
aircraft=self.aircraft,
reverse=True)
lastTurnLeg.buildNewSimulatedArrivalTurnLeg(
deltaTimeSeconds=self.deltaTimeSeconds,
elapsedTimeSeconds=0.0,
distanceStillToFlyMeters=0.0,
simulatedAltitudeSeaLevelMeters=self.firstGlideSlopeWayPoint.
getAltitudeMeanSeaLevelMeters(),
flightPathAngleDegrees=3.0,
bankAngleDegrees=5.0)
descentGlideSlope.addGraph(lastTurnLeg)
#print self.className + ': compute arrival phase length= {0:.2f} meters'.format(descentGlideSlope.getLengthMeters())
#descentGlideSlope.createXlsxOutputFile()
#descentGlideSlope.createKmlOutputFile()
''' prepare next step '''
beginOfLastTurnLeg = lastTurnLeg.getVertex(v=0).getWeight()
print(self.className +
': begin of last turn= {0}'.format(beginOfLastTurnLeg))
''' add to constraint list '''
self.constraintsList.append(
TargetApproachConstraint(beginOfLastTurnLeg))
''' add the three last way-points in the fix list '''
self.insert(position='end', wayPoint=beginOfLastTurnLeg)
''' update the length of the flight path '''
self.distanceFromApproachToTouchDownMeters = descentGlideSlope.getLengthMeters(
)
self.flightLengthMeters = self.computeLengthMeters(
) + descentGlideSlope.getLengthMeters()
print(self.className +
': updated flight path length= {0:.2f} nautics'.format(
self.flightLengthMeters * Meter2NauticalMiles))
''' target approach fix is equal to the begin of the SIMULATED last turn leg '''
self.aircraft.setTargetApproachWayPoint(beginOfLastTurnLeg)
self.aircraft.setArrivalRunwayTouchDownWayPoint(self.touchDownWayPoint)
print(self.className + ': fix list= {0}'.format(self.fixList))
def buildSimulatedArrivalPhase(self):
print self.className + "=========== add final turn, descent and ground run ==================="
arrivalGroundRun = GroundRunLeg(runway=self.arrivalRunway, aircraft=self.aircraft, airport=self.arrivalAirport)
self.touchDownWayPoint = arrivalGroundRun.computeTouchDownWayPoint()
print self.touchDownWayPoint
""" distance from last fix to touch down """
distanceToLastFixNautics = (
self.touchDownWayPoint.getDistanceMetersTo(self.getLastWayPoint()) * Meter2NauticalMiles
)
print self.className + "===================== final 3 degrees descending glide slope ================"
descentGlideSlope = DescentGlideSlope(
runway=self.arrivalRunway,
aircraft=self.aircraft,
arrivalAirport=self.arrivalAirport,
descentGlideSlopeDegrees=3.0,
)
""" if there is a fix nearer to 5 nautics of the touch-down then limit size of simulated glide slope """
descentGlideSlopeSizeNautics = min(distanceToLastFixNautics / 2.0, 5.0)
""" build simulated glide slope """
descentGlideSlope.buildSimulatedGlideSlope(descentGlideSlopeSizeNautics)
self.firstGlideSlopeWayPoint = descentGlideSlope.getVertex(v=0).getWeight()
print self.className + ": top of arrival glide slope= {0}".format(self.firstGlideSlopeWayPoint)
print self.className + " ================= need a turn leg to find the junction point the last way-point in the fix list to the top of the final glide slope"
"""
initial heading is the orientation of the run-way
"""
lastFixListWayPoint = self.wayPointsDict[self.fixList[-1]]
initialHeadingDegrees = self.arrivalRunway.getTrueHeadingDegrees()
lastTurnLeg = TurnLeg(
initialWayPoint=self.firstGlideSlopeWayPoint,
finalWayPoint=lastFixListWayPoint,
initialHeadingDegrees=initialHeadingDegrees,
aircraft=self.aircraft,
reverse=True,
)
lastTurnLeg.buildNewSimulatedArrivalTurnLeg(
deltaTimeSeconds=self.deltaTimeSeconds,
elapsedTimeSeconds=0.0,
distanceStillToFlyMeters=0.0,
simulatedAltitudeSeaLevelMeters=self.firstGlideSlopeWayPoint.getAltitudeMeanSeaLevelMeters(),
flightPathAngleDegrees=3.0,
bankAngleDegrees=5.0,
)
descentGlideSlope.addGraph(lastTurnLeg)
# print self.className + ': compute arrival phase length= {0:.2f} meters'.format(descentGlideSlope.getLengthMeters())
# descentGlideSlope.createXlsxOutputFile()
# descentGlideSlope.createKmlOutputFile()
""" prepare next step """
beginOfLastTurnLeg = lastTurnLeg.getVertex(v=0).getWeight()
print self.className + ": begin of last turn= {0}".format(beginOfLastTurnLeg)
""" add the three last way-points in the fix list """
self.insert(position="end", wayPoint=beginOfLastTurnLeg)
""" update the length of the flight path """
self.distanceFromApproachToTouchDownMeters = descentGlideSlope.getLengthMeters()
self.flightLengthMeters = self.computeLengthMeters() + descentGlideSlope.getLengthMeters()
print self.className + ": updated flight path length= {0:.2f} nautics".format(
self.flightLengthMeters * Meter2NauticalMiles
)
""" target approach fix is equal to the begin of the SIMULATED last turn leg """
self.aircraft.setTargetApproachWayPoint(beginOfLastTurnLeg)
self.aircraft.setArrivalRunwayTouchDownWayPoint(self.touchDownWayPoint)
print self.className + ": fix list= {0}".format(self.fixList)