import numpy
import math
import requests
import datetime

deg2rad = math.pi/180
rad2deg = 180/math.pi

def cartesian_x(f,l):
    #f = latitude, l = longitude
    return (math.cos(f*deg2rad)*math.cos(l*deg2rad))

def cartesian_y(f,l):
    #f = latitude, l = longitude
    return (math.cos(f*deg2rad)*math.sin(l*deg2rad))

def cartesian_z(f,l):
    #f = latitude, l = longitude
    return (math.sin(f*deg2rad))

def cross_x(x, y, z, i,j,k):
    return ((y*k)-(z*j))
def cross_y(x, y, z, i,j,k):
    return ((z*i)-(x*k))
def cross_z(x, y, z, i,j,k):
    return ((x*j)-(y*i))

def spherical_lat(x,y,z):
    r = math.sqrt(x*x + y*y)
    #Omitting the special cases because points will always
    #be separated for this application
    return (math.atan2(z, r)*rad2deg) # return degrees

def spherical_lng(x,y,z):
    #Omitting the special cases because points will always
    #be separated for this application
    return (math.atan2(y, x)*rad2deg) # return degrees


#Let's make a sphere for the Earth

earth = sphere(color='pink', center=(0,0,0), size=1, opacity=0.4)
earth

#Now let's get the height of the nearby ionosphere so that we can make another sphere
#First, the radius of the Earth in km. We'll need this
earth_rad_km = 6371

#What was the height of the f2 when the spot happened?
fromDate = "2023-08-21T03:43:00"
endDate =  "2023-08-21T05:43:00"
start_date_win = fromDate#.strftime("%Y.%m.%d %H:%M:%S")
start_date_win = start_date_win.replace(" ", "%20")
    #print(start_date_win)
end_date_win = endDate#.strftime("%Y.%m.%d %H:%M:%S")
end_date_win = end_date_win.replace(" ", "%20")
    #print("start_date_win " + start_date_win)
print("end_date_win " + end_date_win)
#https://github.com/hcarter333/rm-rbn-history/blob/0e4ed3b1867917be88a9122839c91b781d2fb54b/ionodata.py#L16
iono_url = "https://lgdc.uml.edu/common/DIDBGetValues?ursiCode=PA836&charName=MUFD,hF2,hmF2&fromDate=" + start_date_win + "&toDate=" + end_date_win

#iono_url = "https://lgdc.uml.edu/common/DIDBGetValues?ursiCode=PA836&charName=MUFD,hmF2&fromDate=2024.01.27T11:00&toDate=2024.01.27T14:23"

iodata=requests.get(iono_url)
#print(iodata.text)
iono_data = iodata.text
iono_lines = iono_data.split("\n")
data_found = 0
for data_line in iono_lines:
    print(data_line)
#    if(data_found == 1):
#        #We've arrived at the data
#        #Split on space and print the max F2
#        data_fields = data_line.split()
#        hmF2 = data_fields[4]
#        print("hmF2 " + data_fields[4])
        #break
#we're going to use the closest time to the beginning of the range for now
#    if(data_line.find("#Time") != -1):
#        data_found = 1

#show(earth)



#This might not be the most efficient way, but let's try to
#get the 18 degree more accurately
#lat, lng for tx
t_lat = 37.7248952200944
t_lng = -122.422936174405
#lat,lng for landing point
r_lat = 55.24556451189041
r_lng = -110.29305571520665

tx_x = cartesian_x(t_lat,t_lng)
tx_y = cartesian_y(t_lat,t_lng)
tx_z = cartesian_z(t_lat,t_lng)
rx_x = cartesian_x(r_lat,r_lng)
rx_y = cartesian_y(r_lat,r_lng)
rx_z = cartesian_z(r_lat,r_lng)

#Take the cross product of the two positions
g_x = cross_x(tx_x, tx_y, tx_z, rx_x, rx_y, rx_z)
g_y = cross_y(tx_x, tx_y, tx_z, rx_x, rx_y, rx_z)
g_z = cross_z(tx_x, tx_y, tx_z, rx_x, rx_y, rx_z)

#Now, get the magnitude of the cross product
mag_g = math.sqrt(g_x**2+g_y**2+g_z**2)
mag_tx = math.sqrt(tx_x**2+tx_y**2+tx_z**2)
print("magnitude of cross product" + str(mag_g))
print("magnitude of cross product" + str(mag_tx))
theta = math.asin(mag_g)*rad2deg
print("swept_angle = " + str(theta))

tx_cos = math.cos(deg2rad*theta)
F2_min = earth_rad_km*((1-tx_cos)/tx_cos)
print("Calculated minimum height of F2", str(F2_min))
print("radius of earth in km ", earth_rad_km)

print("Now let's look at the F2 height at Idaho")
iono_url = "https://lgdc.uml.edu/common/DIDBGetValues?ursiCode=IF843&charName=MUFD,hF2,hmF2&fromDate=" + start_date_win + "&toDate=" + end_date_win

#iono_url = "https://lgdc.uml.edu/common/DIDBGetValues?ursiCode=PA836&charName=MUFD,hmF2&fromDate=2024.01.27T11:00&toDate=2024.01.27T14:23"

iodata=requests.get(iono_url)
#print(iodata.text)
iono_data = iodata.text
iono_lines = iono_data.split("\n")
data_found = 0
for data_line in iono_lines:
    print(data_line)

print("\n\n\nAnd then from Thule")
iono_url = "https://lgdc.uml.edu/common/DIDBGetValues?ursiCode=THJ76&charName=MUFD,hF2,hmF2&fromDate=" + start_date_win + "&toDate=" + end_date_win

#iono_url = "https://lgdc.uml.edu/common/DIDBGetValues?ursiCode=PA836&charName=MUFD,hmF2&fromDate=2024.01.27T11:00&toDate=2024.01.27T14:23"

iodata=requests.get(iono_url)
#print(iodata.text)
iono_data = iodata.text
iono_lines = iono_data.split("\n")
data_found = 0
for data_line in iono_lines:
    print(data_line)

#catenary(a,x) = a*cosh(x/a)
#a2 = plot(catenary(0.2,x), 0, 0.5, color='red', legend_label='a=0.2', legend_color='red')
#a3 = plot(catenary(0.3,x), 0, 0.5, color='green', legend_label='a=0.3', legend_color='green')
#a4 = plot(catenary(0.4,x), 0, 0.5, color='blue', le
exam(x)=-2*(x^4)-(5*(x^3))-(x^2)-(8*x)-6
hamie=plot(exam(x), -3, 3, color='green', ymin=-30,ymax=40)
hamie