# m.drawcountries() # m.drawstates() #m.fillcontinents(color='coral', lake_color='aqua') m.drawparallels(np.arange(-90., 120., 30.)) m.drawmeridians(np.arange(0., 420., 60.)) if dn: m.nightshade(dn) #m.drawmapboundary(fill_color='aqua') return m time = GPSTime.from_datetime(datetime(2019,5,31,11,0,0)) dog = AstroDog() ionex_map = dog.get_ionex(time) # print(dog.get_ionex(time).get_TEC((5,5), time)) # Setup map: fig = plt.figure(1, figsize=(18, 12)) plt.clf() ax = fig.add_subplot(111) m = setup_map(ax) lon_bins = np.linspace(-180, 180, 73) lat_bins = np.linspace(-90, 90, 73) # print(lon_bins, lat_bins) tec = []
for idx, timestamp in enumerate(daterange.to_pydatetime()): save_idx = idx % num_save print(timestamp) i_kp = [kp_data.index.get_loc(t, method='nearest') for t in [timestamp]] kp = kp_data.iloc[i_kp, :]['kp'][0] i_f107 = [ f107_data.index.get_loc(t, method='nearest') for t in [timestamp] ] f107 = f107_data.iloc[i_f107, :]['f107'][0] time = GPSTime.from_datetime(timestamp) dog = AstroDog() ionex_map = dog.get_ionex(time) doy = timestamp.timetuple().tm_yday fTime = timestamp.hour + timestamp.minute / 60. sub_idx = 0 for lat in lat_bins: for lon in lon_bins: tec = dog.get_ionex(time).get_TEC((lat, lon), time) data[subbins * save_idx + sub_idx] = lat, lon, f107, kp, fTime, doy, tec sub_idx = sub_idx + 1 if (idx + 1) % num_save == 0: print("Storing results checkpoint...") try: