geospacelab 0.12.1__py3-none-any.whl → 0.12.3__py3-none-any.whl

geospacelab/__init__.py

@@ -6,7 +6,7 @@ __author__ = "Lei Cai"
 __copyright__ = "Copyright 2021, GeospaceLAB"
 __credits__ = ["Lei Cai"]
 __license__ = "BSD-3-Clause License"
-__version__ = "0.12.1"
+__version__ = "0.12.3"
 __maintainer__ = "Lei Cai"
 __email__ = "lei.cai@oulu.fi"
 __status__ = "Developing"

geospacelab/datahub/sources/cdaweb/downloader.py

@@ -213,8 +213,11 @@ class DownloaderFromFTPBase(DownloaderBase):
         if file_paths_local is None:
             self.file_paths_local = []
             for fp_remote in self.file_paths_remote:
+                fp_local = self.root_dir_local
+                for p in fp_remote.replace(root_dir_remote, '').split('/'):
+                    fp_local = fp_local / p
                 self.file_paths_local.append(
-                    self.root_dir_local / fp_remote.replace(root_dir_remote, '')
+                    fp_local
                 )
         self.done = [False] * len(self.file_paths_remote)
         
@@ -347,8 +350,11 @@ class DownloaderFromHTTPBase(DownloaderBase):
         if file_paths_local is None:
             self.file_paths_local = []
             for fp_remote in self.file_paths_remote:
+                fp_local = self.root_dir_local
+                for p in fp_remote.replace(root_url, '').split('/'):
+                    fp_local = fp_local / p
                 self.file_paths_local.append(
-                    self.root_dir_local / fp_remote.replace(root_url, '')
+                    fp_local
                 )
         else:
             self.file_paths_local = file_paths_local

geospacelab/datahub/sources/madrigal/satellites/dmsp/e/__init__.py

@@ -200,11 +200,12 @@ class Dataset(datahub.DatasetSourced):
         delta = np.sin(glon_2 * np.pi / 180.) - np.sin(glon_2 * np.pi / 180.)
         
         if self.replace_orbit:
-            glon_1 = glon_2
+            self['SC_GEO_LON'].value = glon_2[:, np.newaxis]
+            self['SC_GEO_LAT'].value = cs_new['lat'][:, np.newaxis]
         else:
             glon_1[np.abs(delta)>0.001] = glon_2[np.abs(delta)>0.001]
-            
-        self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
+            self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
+        
 
     def search_data_files(self, **kwargs):
 

geospacelab/datahub/sources/madrigal/satellites/dmsp/s1/__init__.py

@@ -134,7 +134,7 @@ class Dataset(datahub.DatasetSourced):
             self.time_filter_by_range(var_datetime_name='SC_DATETIME')
 
         if self.calib_orbit:
-            self.fix_geo_lon()
+           self.fix_geo_lon()
 
         if self.add_AACGM:
             self.convert_to_AACGM()
@@ -220,11 +220,11 @@ class Dataset(datahub.DatasetSourced):
         delta = np.sin(glon_2 * np.pi / 180.) - np.sin(glon_1 * np.pi / 180.)
         
         if self.replace_orbit:
-            glon_1 = glon_2
+            self['SC_GEO_LON'].value = glon_2[:, np.newaxis]
+            self['SC_GEO_LAT'].value = cs_new['lat'][:, np.newaxis]
         else:
             glon_1[np.abs(delta)>0.001] = glon_2[np.abs(delta)>0.001]
-            
-        self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
+            self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
 
     def search_data_files(self, **kwargs):
 

geospacelab/datahub/sources/madrigal/satellites/dmsp/s4/__init__.py

@@ -192,11 +192,11 @@ class Dataset(datahub.DatasetSourced):
         delta = np.sin(glon_2 * np.pi / 180.) - np.sin(glon_2 * np.pi / 180.)
         
         if self.replace_orbit:
-            glon_1 = glon_2
+            self['SC_GEO_LON'].value = glon_2[:, np.newaxis]
+            self['SC_GEO_LAT'].value = cs_new['lat'][:, np.newaxis]
         else:
             glon_1[np.abs(delta)>0.001] = glon_2[np.abs(delta)>0.001]
-            
-        self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
+            self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
         
         # from geospacelab.observatory.orbit.sc_orbit import OrbitPosition_SSCWS
         # from scipy.interpolate import interp1d

geospacelab/datahub/sources/tud/downloader.py

@@ -106,11 +106,10 @@ class TUDownloader(DownloaderBase):
                 version_dict[version_in_fn[:3]], 
                 self.mission + '_data'
             ]]) + '/'
-
-            fp_local = fp_remote.replace(replacement, '')
-            file_paths_local.append(
-                self.root_dir_local / self.mission.upper() /
-                self.product.upper() / version_in_fn / fp_local)
+            fp_local = self.root_dir_local / self.mission.upper() / self.product.upper() / version_in_fn
+            for p in fp_remote.replace(replacement, '').split('/'):
+                fp_local = fp_local / p
+            file_paths_local.append(fp_local)
         super().save_files_from_http(file_paths_local, root_dir_remote)  
         
         for i, (done, file_path) in enumerate(zip(self.done, self.file_paths_local)):

geospacelab/observatory/orbit/utilities.py

@@ -8,7 +8,7 @@ import scipy.signal as sig
 from geospacelab.datahub import DatasetUser
 from geospacelab.toolbox.utilities import pydatetime as dttool
 import geospacelab.toolbox.utilities.numpymath as npmath
-from geospacelab.cs import GEOCSpherical
+import geospacelab.cs as geo_cs
 
 
 class LEOToolbox(DatasetUser):
@@ -34,6 +34,27 @@ class LEOToolbox(DatasetUser):
         self.sector_cs = 'GEO'
         self.sectors = {}
 
+    def add_GEO_cartesian(self, cs_in='GEO', ):
+        
+        if cs_in == 'GEO':
+            coords={
+                'lat': self['SC_GEO_LAT'].value.flatten(),
+                'lon': self['SC_GEO_LON'].value.flatten(),
+                'height': self['SC_GEO_ALT'].value.flatten(),
+                'lat_unit': 'deg', 'lon_unit': 'deg', 'height_unit': 'km'
+            }
+        else:
+            raise NotImplementedError
+        
+        cs = geo_cs.GEOSpherical(
+            coords=coords
+        )
+        cs_new = cs.to_cartesian()
+        self.add_variable(var_name='SC_GEO_X', value=cs_new['x'].reshape((cs_new['x'].size, 1)))
+        self.add_variable(var_name='SC_GEO_Y', value=cs_new['y'].reshape((cs_new['y'].size, 1)))
+        self.add_variable(var_name='SC_GEO_Z', value=cs_new['z'].reshape((cs_new['z'].size, 1)))
+        return 
+    
     def search_orbit_nodes(self, data_interval=1, t_res=None):
         def check_nodes(nodes):
             dts_nodes = nodes['DATETIME']
@@ -330,7 +351,266 @@ class LEOToolbox(DatasetUser):
             }
         else:
             raise NotImplementedError
+    
+    def griddata_by_sector_new(
+        self, 
+        variable_names=None,
+        sector_name=None,
+        x_grid_res=20*60,
+        y_grid_res=0.5,
+        x_data_res=None,
+        y_data_res=None,
+        grid_method='linear',
+        along_x_interp_method='linear',
+        along_track_interp=True,
+        along_track_interp_method='linear',
+        along_track_binning=False,
+        along_track_binning_method='mean',
+        along_track_binning_size=None,
+        visual='on',
+    ):      
+        # Initialize parameters
+        if along_track_binning and along_track_interp:
+            raise ValueError('Only one of along_track_binning and along_track_interp can be True.') 
+        if along_track_binning:
+            if along_track_binning_size is None:
+                along_track_binning_size = y_grid_res
+        
+        self.visual = visual
+        dts_c = self['_'.join(('SECTOR', sector_name, 'DATETIME'))].value.flatten()
+        dts = self['SC_DATETIME'].value.flatten()
+        sector = self['SECTOR_' + sector_name].value.flatten()
+        pseudo_lat = self['_'.join(('SECTOR', sector_name, 'PSEUDO_LAT'))].value.flatten()
+        lat_range = self.sectors[sector_name]['PSEUDO_LAT_RANGE']
+        n_tracks = int(np.max(sector))
+        
+        if along_track_interp:
+            x_data = dts[sector > 0]
+        else:
+            x_data = dts_c[sector > 0]    
+        x_data_c = dts_c[sector>0]
+        x = np.array([(dt - dttool.get_start_of_the_day(self.dt_fr)).total_seconds() for dt in x_data])
+        x_c = np.array([(dt - dttool.get_start_of_the_day(self.dt_fr)).total_seconds() for dt in x_data_c])
+        x_unique_c = np.unique(x_c)
+        if x_data_res is None:
+            x_data_res = np.nanmedian(np.diff(x_unique_c))
+        
+        y_data = pseudo_lat[sector > 0]
+        y= y_data
+        if y_data_res is None:
+            diff_ys = []
+            for i in range(n_tracks):
+                sector_i = sector[sector>0]
+                y_track = y_data[sector_i == i+1]
+                if len(y_track) > 3:
+                    diff_ys.append(np.abs(y_track[len(y_track) // 2] - y_track[len(y_track) // 2 - 1]))
+            y_data_res = np.nanmedian(diff_ys)
+            if y_grid_res > y_data_res:
+                y_data_res = y_grid_res
+                
+        # Set the grids
+        dt0 = dttool.get_start_of_the_day(self.dt_fr)
+        ny = int(np.ceil(np.diff(lat_range) / y_grid_res) + 1)
+        ys = np.linspace(lat_range[0], lat_range[1], ny)
+        if x_grid_res is None:
+            xs = x_unique_c
+            # ys = ys[:-1] + y_grid_res / 2 
+            grid_x, grid_y = np.meshgrid(xs, ys)
+            x_interp = False
+            x_grid_res=x_data_res
+        else:
+            min_x = np.floor((self.dt_fr - dt0).total_seconds() / x_grid_res) * x_grid_res
+            max_x = np.ceil((self.dt_to - dt0).total_seconds() / x_grid_res) * x_grid_res
+            grid_x, grid_y = np.meshgrid(
+                np.arange(min_x, max_x, x_grid_res),
+                ys 
+            )
+            x_interp = True
+        if along_track_binning:
+            y_bins = [
+                ys - along_track_binning_size / 2,
+                ys + along_track_binning_size / 2
+            ]
+
+        grid_lat = griddata((x_c, y), y, (grid_x, grid_y), method='nearest')
+
+        if self.sector_cs == 'AACGM':
+            mask_y = np.abs(grid_y - grid_lat) > y_data_res * 2
+            which_lat = 'AACGM_LAT'
+        if self.sector_cs == 'APEX':
+            mask_y = np.abs(grid_y - grid_lat) > y_data_res * 2
+            which_lat = 'APEX_LAT' 
+        else:
+            mask_y = np.abs(grid_y - grid_lat) > y_data_res * 1.2
+            which_lat = 'GEO_LAT'
+
+        grid_sectime = griddata((x_c, y), x_c, (grid_x, grid_y), method='nearest')
+        mask_x = np.abs(grid_x - grid_sectime) > x_data_res * 1.2
+
+        # remove gaps
+        # i_mask_x = np.where(mask_x)
+        # rec_bounds_ind = {}
+        # if list(i_mask_x[0]):
+        #     for ii, jj in zip(*i_mask_x):
+        #         xx = grid_x[ii, jj]
+        #         x_ii = grid_x[ii, :].flatten()
+        #         if xx in rec_bounds_ind.keys():
+        #             bound_x_1 = rec_bounds_ind[xx][0]
+        #             bound_x_2 = rec_bounds_ind[xx][1]
+        #         else:
+        #             iii = np.where(x_unique_c<xx)[0]
+        #             if not list(iii):
+        #                 bound_x_1 = x_ii[0] - x_data_res
+        #             else:
+        #                 bound_x_1 = x_unique_c[iii[-1]] - x_grid_res
+        #             iii = np.where(x_unique_c > xx)[0]
+        #             if not list(iii):
+        #                 bound_x_2 = x_ii[-1] + x_grid_res
+        #             else:
+        #                 bound_x_2 = x_unique_c[iii[0]] + x_data_res
+        #             rec_bounds_ind[xx] = [bound_x_1, bound_x_2]
+        #         mask_x[ii, (x_ii >= bound_x_1) & (x_ii <= bound_x_2)] = True
+        #     pass
+        
+        # 
+        
+        for var_name in variable_names:
+            # print(var_name)
+            vrb = self[var_name].value.flatten()[sector>0]
+            if along_track_binning:
+                sector_1 = sector[sector>0]
+                n_lats = ny
+                xx_1 = np.full((n_tracks, n_lats), np.nan)
+                yy_1 = grid_y[:, 0].flatten()
+                zz_1 = np.full((n_tracks, n_lats), np.nan)
+                grid_z = np.ones_like(grid_x)*np.nan
                 
+                for ii in np.arange(1, n_tracks+1):
+                    xd = x[sector_1 == ii]
+                    yd = y[sector_1 == ii]
+                    zd = vrb[sector_1 == ii]
+                    if (not list(xd)) or (not list(yd)):
+                        continue
+                    inds_sorted = np.argsort(yd)
+                    yd_sorted = yd[inds_sorted]
+                    inds_y_bin_0 = np.searchsorted(yd_sorted, y_bins[0])
+                    inds_y_bin_1 = np.searchsorted(yd_sorted, y_bins[1])
+                    for ib0, ib1 in zip(inds_y_bin_0, inds_y_bin_1):
+                        if along_track_binning_method == 'mean':
+                            z_i = np.nanmean(zd[inds_sorted[ib0:ib1]])
+                        elif along_track_binning_method == 'median':
+                            z_i = np.nanmedian(zd[inds_sorted[ib0:ib1]])
+                        elif along_track_binning_method == 'max':
+                            z_i = np.nanmax(zd[inds_sorted[ib0:ib1]])
+                        elif along_track_binning_method == 'min':
+                            z_i = np.nanmin(zd[inds_sorted[ib0:ib1]])
+                        
+                            if list(z_i):
+                                z_bin = np.nanmean(z_i)
+                            else:
+                                z_bin = np.nan
+                    f = interp1d(yd, xd, bounds_error=False, fill_value= 'extrapolate')
+                    x_i = f(yy_1)
+                    xx_1[ii-1, :] = x_i
+
+                    if 'LON' in var_name:
+                        z_i = npmath.interp_period_data(yd, zd, yy_1,  period=360., method='linear', bounds_error=False)
+                    elif 'LST' in var_name:
+                        z_i = npmath.interp_period_data(yd, zd, yy_1,  period=24., method='linear', bounds_error=False)
+                    else:
+                        # f = interp1d(yd, zd, bounds_error=False, fill_value='extrapolate')
+                        # z_i = f(yy_1)
+                        z_i = np.ones_like(yy_1) * np.nan       
+                        for j, yyy in enumerate(yy_1):
+                            inds_y = np.where((yd >= yyy-y_grid_res/2) & (yd < yyy+y_grid_res/2))[0]
+                            if list(inds_y):
+                                z_i[j] = np.nanmean(zd[inds_y])
+                    zz_1[ii-1, :] = z_i 
+                for ii in range(grid_x.shape[0]):
+                    xd = xx_1[:, ii].flatten()
+                    zd = zz_1[:, ii].flatten()
+                    if x_interp:
+                        if 'LON' in var_name:
+                            z_i = npmath.interp_period_data(xd, zd, grid_x[0],  period=360., method='linear', bounds_error=False)
+                        elif 'LST' in var_name:
+                            z_i = npmath.interp_period_data(xd, zd, grid_x[0],  period=24., method='linear', bounds_error=False)
+                        else:
+                            f = interp1d(xd, zd, bounds_error=False, fill_value='extrapolate')
+                            z_i = f(grid_x[0])
+                        grid_z[ii, :] = z_i
+                    else:
+                        grid_z[ii, :] = zd 
+            elif along_track_interp:
+                sector_1 = sector[sector>0]
+                n_tracks = int(np.max(sector))
+                n_lats = ny
+                xx_1 = np.empty((n_tracks, n_lats))
+                zz_1 = np.empty((n_tracks, n_lats))
+                yy_1 = grid_y[:, 0].flatten()
+                grid_z = np.ones_like(grid_x) * np.nan
+                for ii in np.arange(1, n_tracks+1):
+                    xd = x[sector_1 == ii]
+                    yd = y[sector_1 == ii]
+                    zd = vrb[sector_1 == ii]
+                    
+                    if (not list(xd)) or (not list(yd)):
+                        continue
+
+                    f = interp1d(yd, xd, bounds_error=False, fill_value= 'extrapolate')
+                    x_i = f(yy_1)
+                    xx_1[ii-1, :] = x_i
+
+                    if 'LON' in var_name:
+                        z_i = npmath.interp_period_data(yd, zd, yy_1,  period=360., method='linear', bounds_error=False)
+                    elif 'LST' in var_name:
+                        z_i = npmath.interp_period_data(yd, zd, yy_1,  period=24., method='linear', bounds_error=False)
+                    else:
+                        f = interp1d(yd, zd, bounds_error=False, fill_value='extrapolate')
+                        z_i = f(yy_1)
+                    zz_1[ii-1, :] = z_i 
+                for ii in range(grid_x.shape[0]):
+                    xd = xx_1[:, ii].flatten()
+                    zd = zz_1[:, ii].flatten()
+                    if x_interp:
+                        if 'LON' in var_name:
+                            z_i = npmath.interp_period_data(xd, zd, grid_x[0],  period=360., method='linear', bounds_error=False)
+                        elif 'LST' in var_name:
+                            z_i = npmath.interp_period_data(xd, zd, grid_x[0],  period=24., method='linear', bounds_error=False)
+                        else:
+                            f = interp1d(xd, zd, bounds_error=False, fill_value='extrapolate')
+                            z_i = f(grid_x[0])
+                        grid_z[ii, :] = z_i
+                    else:
+                        grid_z[ii, :] = zd
+            else:
+                grid_z = griddata((x, y), vrb, (grid_x, grid_y), method=method)
+            grid_z[mask_y] = np.nan
+            grid_z[mask_x] = np.nan
+            var_name_in = '_'.join(('SECTOR', sector_name, 'GRID', var_name))
+            self.sectors[sector_name]['VARIABLE_NAMES'].append(var_name_in) 
+            
+            # self.add_variable(var_name=var_name_in, value=grid_z.T)
+            self[var_name_in] = self[var_name].clone(omit_attrs='visual')
+            self[var_name_in].visual = 'new'
+            var = self[var_name_in]
+            var.value = grid_z.T
+            var.set_depend(0, {'UT': '_'.join(('SECTOR', sector_name, 'GRID_DATETIME'))})
+            var.set_depend(1, {which_lat: '_'.join(('SECTOR', sector_name, 'GRID_LAT'))})
+            var.visual.axis[0].data = '@d.' + '_'.join(('SECTOR', sector_name, 'GRID_DATETIME'))
+            var.visual.axis[1].data = '@d.' + '_'.join(('SECTOR', sector_name, 'GRID_LAT'))
+            var.visual.axis[1].label = 'GLAT' if self.sector_cs == 'GEO' else 'MLAT'
+            var.visual.axis[1].unit = r'$^\circ$'
+            var.visual.axis[1].lim = lat_range
+            var.visual.axis[2].data = '@v.value'
+            var.visual.axis[2].data_scale = self[var_name].visual.axis[1].data_scale
+            var.visual.axis[2].label = '@v.label'
+            var.visual.axis[2].unit = '@v.unit'
+            var.visual.plot_config.style = '2P'
+            var.visual.plot_config.pcolormesh.update(cmap='jet')
+        
+           
+        return
+    
     def griddata_by_sector(
             self, sector_name=None, variable_names=None, x_grid_res=20*60, y_grid_res=0.5, along_track_interp=True,
             x_data_res=None, y_data_res=None, along_track_binning=False, 

{geospacelab-0.12.1.dist-info → geospacelab-0.12.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: geospacelab
-Version: 0.12.1
+Version: 0.12.3
 Summary: Collect, manage, and visualize geospace data.
 Home-page: https://github.com/JouleCai/geospacelab
 Author: Lei Cai
@@ -27,8 +27,8 @@ Requires-Dist: beautifulsoup4>=4.9.3
 Requires-Dist: natsort>=7.1.1
 Requires-Dist: numpy<2.4.0
 Requires-Dist: scipy>=1.6.0
-Requires-Dist: h5py>=3.2.1
 Requires-Dist: netcdf4>=1.5.7
+Requires-Dist: h5py>=3.2.1
 Requires-Dist: matplotlib>=3.5
 Requires-Dist: madrigalweb>=3.3
 Requires-Dist: aacgmv2>=2.6.2
@@ -376,6 +376,12 @@ Output:
 From then on, the data source has been switched from JHUAPL to CDAWeb.
  
 ```python
+
+import datetime
+import matplotlib.pyplot as plt
+import geospacelab.visualization.mpl.geomap.geodashboards as geomap
+
+
 dt_fr = datetime.datetime(2015, 9, 8, 8)
 dt_to = datetime.datetime(2015, 9, 8, 23, 59)
 time_c = datetime.datetime(2015, 9, 8, 20, 21)

{geospacelab-0.12.1.dist-info → geospacelab-0.12.3.dist-info}/RECORD

@@ -1,4 +1,4 @@
-geospacelab/__init__.py,sha256=wQvOWtd7JdoSlfT_9BWr10LK2nbINk2GnJSok-4mppU,801
+geospacelab/__init__.py,sha256=k0D5pyAhRJ3aOpvHg8OpbijUV_a_Fx4BJzB7U1iCUMA,801
 geospacelab/config/__init__.py,sha256=--F2bcKRCNIbPaFz4bySMkuxGg1ZxF1j1uyxMA0t3xA,660
 geospacelab/config/__mpl__.py,sha256=bO10-mtYDF1EhnRFuXX4H8TUU3esxOZe-qfTrzRhBk0,4705
 geospacelab/config/_preferences.py,sha256=VKXXWNvFzIMXGmvb7Udd900DoubMlqrcu-gLmzHmzMk,4754
@@ -23,7 +23,7 @@ geospacelab/datahub/_variable_base.py,sha256=dUm15Y012FQ0hAdeA48ky0tKEc9il59FHvO
 geospacelab/datahub/test.py,sha256=ZXtv9nhZd3RKr8BPsKqSB_boCpB2P4SYF6hFKdo8aAg,682
 geospacelab/datahub/sources/__init__.py,sha256=VZtk8qV_O3LQtUQOOwvBwbqmcI_labsWdml-iqpH8u4,327
 geospacelab/datahub/sources/cdaweb/__init__.py,sha256=fP1ziTZFxTOz-HF07az67F5AfvKJsXdJluEK-NZi6yU,680
-geospacelab/datahub/sources/cdaweb/downloader.py,sha256=-09HuoIa48zMVHn3-IDZVn_59WnzJoSx1w17jtCQVFk,15057
+geospacelab/datahub/sources/cdaweb/downloader.py,sha256=yR4clbgB69zzy0xcfUSKyYIzo_odD9Yu4aY9J8691ew,15287
 geospacelab/datahub/sources/cdaweb/downloader_backup.py,sha256=LCn0tupzAkir9Z_280tBulHEPlqoX6DWe8HKFk1h4IE,3191
 geospacelab/datahub/sources/cdaweb/dmsp/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 geospacelab/datahub/sources/cdaweb/dmsp/ssusi/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -202,13 +202,13 @@ geospacelab/datahub/sources/madrigal/isr/risr_n/vi/variable_config.py,sha256=dE-
 geospacelab/datahub/sources/madrigal/satellites/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 geospacelab/datahub/sources/madrigal/satellites/dmsp/__init__.py,sha256=YnuaB5J_ER8LrLUFghlty0YHYnmd9_2IHNW9TBQLaSI,487
 geospacelab/datahub/sources/madrigal/satellites/dmsp/downloader.py,sha256=WUtJN1OPF_CzH1ZxvM01TG1O71Vnr5sx1KcyfVWGPbE,4913
-geospacelab/datahub/sources/madrigal/satellites/dmsp/e/__init__.py,sha256=_l1rcUrGWUV8XJg0D1KzQEeA3bLEWg_y5w2be2wDhM4,10246
+geospacelab/datahub/sources/madrigal/satellites/dmsp/e/__init__.py,sha256=fVcmj9Jyym_OiVG-aT-lawJMzmIEdykdzQK4KTwrXbc,10347
 geospacelab/datahub/sources/madrigal/satellites/dmsp/e/loader.py,sha256=kZ9AeI0LsNWS26Kwfp6azvDl2IutcWKBPDd6o0tKASc,4902
 geospacelab/datahub/sources/madrigal/satellites/dmsp/e/variable_config.py,sha256=exp1bYVmDBVlIbraeNgAuVEnjxg5IC2ED1cZIEiPj1U,10280
-geospacelab/datahub/sources/madrigal/satellites/dmsp/s1/__init__.py,sha256=le8raC-lW_vUTXLvdsmf5-jhUGS7tlZ7ZeLH__zD9NI,11173
+geospacelab/datahub/sources/madrigal/satellites/dmsp/s1/__init__.py,sha256=otMfzaf9vTgultEXXMDnwVtXVO1HCPqNOrPYuFRLaPw,11264
 geospacelab/datahub/sources/madrigal/satellites/dmsp/s1/loader.py,sha256=o5zEkosq7LyNlBIjlefKUpYpePCDbyL_yd3HrrUrk5Q,3360
 geospacelab/datahub/sources/madrigal/satellites/dmsp/s1/variable_config.py,sha256=Amp90wHogRQOESLH2b4D0jfeBozIKvvOnwx-_waU0us,6745
-geospacelab/datahub/sources/madrigal/satellites/dmsp/s4/__init__.py,sha256=6lA2b-pE_h9ofuWMgTJqpEdDRGY2rGlVAKvIqoGvXe4,12077
+geospacelab/datahub/sources/madrigal/satellites/dmsp/s4/__init__.py,sha256=eqSbH10r6e650ShTwOLGcqXSYNGYDmiak2eZFKkSFeA,12169
 geospacelab/datahub/sources/madrigal/satellites/dmsp/s4/loader.py,sha256=LL8IAMwmgQX06lgVECh6Rx6TGXEoHmAkFDmouJ6kQOo,3229
 geospacelab/datahub/sources/madrigal/satellites/dmsp/s4/variable_config.py,sha256=4DhvQxFonoYiU_3f9ObE2OkJOXDkHfw0RUMHa3FLTv4,4054
 geospacelab/datahub/sources/ncei/__init__.py,sha256=MHpNJ-PJBjywheG5CZDgj89fsxVG5qSxi_vkTsxibd8,658
@@ -246,7 +246,7 @@ geospacelab/datahub/sources/supermag/magnetometer/downloader.py,sha256=P9o888vds
 geospacelab/datahub/sources/supermag/magnetometer/loader.py,sha256=Jef9eydQDPNcCVHpji6nUWGdtBMR5Phbmd0i7PPEd3M,2002
 geospacelab/datahub/sources/supermag/magnetometer/variable_config.py,sha256=DgXzDgrz8cLQj8tlxk_Fqmb5vvmHQdUfNrmkBaytgM8,5106
 geospacelab/datahub/sources/tud/__init__.py,sha256=83o4oiR4K8ZYQdTZ8ZwvZyyYkEhEzXjIo44861XTMSg,2706
-geospacelab/datahub/sources/tud/downloader.py,sha256=eQSuPpJ2fYFtL3dMONg2mcKIfPI5rR1AfI3xksD0ZBw,12837
+geospacelab/datahub/sources/tud/downloader.py,sha256=FtX4ae3v-VEmMcIVwouW-xt6mPDbFPU3UOQ5p3QtdpU,12874
 geospacelab/datahub/sources/tud/champ/__init__.py,sha256=iGqgkK3TM-vGgySrlQR_6ry1EGRqcTDXQTAewf-cgNU,559
 geospacelab/datahub/sources/tud/champ/dns_acc/__init__.py,sha256=hEVXH7H7XtjgzBfu3HAvbkqRqrYMr0nUTCpAM6ll3V4,10276
 geospacelab/datahub/sources/tud/champ/dns_acc/downloader.py,sha256=vjpKCZqa5-zTMAo7LHNUMqWCqRutnu-loRGZ-SfG1AU,1121
@@ -335,7 +335,7 @@ geospacelab/observatory/earth/geodesy.py,sha256=6ouOEgiqlTrDBHdPjhZ5p8BScfG8-NYr
 geospacelab/observatory/earth/sun_position.py,sha256=gy-zPPoX_92r6Cz7BT7JQyiNO_tte3osCQFkNq9lWe0,1434
 geospacelab/observatory/orbit/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 geospacelab/observatory/orbit/sc_orbit.py,sha256=FgTAFnoL6HRNlP8C9o2IDE0TpfGwaAP_3_2_8x4Euxo,15588
-geospacelab/observatory/orbit/utilities.py,sha256=ludGybXzPaeVJOS_XfO7M0TEnf-7MiRf_b6CYgGlLLI,34083
+geospacelab/observatory/orbit/utilities.py,sha256=zhSBcSWAbfx_vD3Gd4zHS0AuZzfwPp-WYyjMEuajRZc,47245
 geospacelab/observatory/site/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 geospacelab/observatory/site/mixins.py,sha256=rhMJngBVgv0wR1m8e5cvnPxHZQnntbcJIgp4bEyTMFE,457
 geospacelab/quantity/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -400,7 +400,7 @@ geospacelab/wrapper/geopack/geopack/t89.py,sha256=zDVNPrmtK1NnNHgohQEPqOOJDsm2Z-
 geospacelab/wrapper/geopack/geopack/t96.py,sha256=ktcoo1R7Z3NtkWHENuseu48ub4-JfQGqFV0ZOtd0zH8,65292
 geospacelab/wrapper/geopack/geopack/test_geopack1.md,sha256=dMUY0O1BgZsKpmJ6BLSQ80B6p6DZcB7OceFeyPOlFK0,15324
 geospacelab/wrapper/geopack/geopack/test_geopack1.py,sha256=qjLz6O3BAk3H58IpmxXyftwZTkh3vPGp49C-al4hjf0,6669
-geospacelab-0.12.1.dist-info/licenses/LICENSE,sha256=2yRlwLt4o5Z6OZAGcyvBj-zfFX1Uw7E6CzqODg7khqs,1515
+geospacelab-0.12.3.dist-info/licenses/LICENSE,sha256=2yRlwLt4o5Z6OZAGcyvBj-zfFX1Uw7E6CzqODg7khqs,1515
 test/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 test/test_ampere.py,sha256=0-HZURubpv1mBK3bJ_qTqx39L1jezgRoU5neXMPYgZQ,2968
 test/test_dmsp_s1.py,sha256=5m_7mjdDGja8ovshNPV3LKW_6q6mIwT9XKqoyRiH79A,3588
@@ -410,7 +410,7 @@ test/test_omni.py,sha256=Zk1LZozPiY5V0aSRmK6GTQuB01hHn_j2j3Brm6Ea_po,1632
 test/test_superdarn.py,sha256=uP55muvXryPzNGHinWkiGv2PxvRs4f9M9h1WIBEBW7k,2846
 test/test_swarm.py,sha256=PDDE9nUshhQpXZbV_ZwcsjbMhI73fRaojTZv9rtRzZE,15568
 test/test_swarm_new.py,sha256=mzhMAx-M9W3Ue5noTyfBx4c3Vtc3b_ZUEvGgL9v8UE4,853
-geospacelab-0.12.1.dist-info/METADATA,sha256=LalBLB6dH4RnIFpRHkp7QF66g4EnNwj6X_J8FPecPRo,24190
-geospacelab-0.12.1.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
-geospacelab-0.12.1.dist-info/top_level.txt,sha256=98eDwrSNgyQFAtSA06QMP71gw9BzgIj0uvkTudpGly4,12
-geospacelab-0.12.1.dist-info/RECORD,,
+geospacelab-0.12.3.dist-info/METADATA,sha256=woYG0JznX74bBp-IBp_dGuKltHt5cdAzfEn0fYK9Fys,24309
+geospacelab-0.12.3.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
+geospacelab-0.12.3.dist-info/top_level.txt,sha256=98eDwrSNgyQFAtSA06QMP71gw9BzgIj0uvkTudpGly4,12
+geospacelab-0.12.3.dist-info/RECORD,,