#!/usr/bin/python3
# north_pass -- compute elevations of planes passing due north of me by the stars
# Copyright: https://www.gnu.org/licenses/gpl.html
# Author: Dan Jacobson https://www.jidanni.org/
# Created: Mar 2003
# Last-Updated: 2024-06-13T12:05:43+0000
#      Update #: 50

latitude=24.18170479 #mine, also the elevation of the pole (star)
north=9338.0 #horizontal distance M-750 air route passes north of me in meters
my_elevation=777
f2m=.3048 #feet to meters
import math
fl=range(1000*(int(16000*f2m/1000)),int(46000*f2m),1000)
for height in fl:
    elevation=math.atan(height/north)*180/math.pi
    declination=90-(elevation-latitude)
    print ("%5d %.1f %.1f %5d %5d" %(
        height, #Meters above surface
        elevation, #Angle above my horizon
        declination, #Declination on star chart at the north-south meridian
        height+my_elevation, #Meters above mean sea level
        (height+my_elevation)/f2m)) #Feet above mean sea level