#!/usr/bin/python
# north_pass -- compute elevations of planes passing due north of me by the stars
# Copyright       : http://www.fsf.org/copyleft/gpl.html
# Author          : Dan Jacobson -- http://jidanni.org/comm/air/m750/
# Created On      : Mar 2003
# Last Modified On: Thu Mar 23 09:11:07 2006
# Update Count    : 41

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(25000*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