geospacelab 0.11.0__py3-none-any.whl → 0.11.2__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.11.0"
+__version__ = "0.11.2"
 __maintainer__ = "Lei Cai"
 __email__ = "lei.cai@oulu.fi"
 __status__ = "Developing"

geospacelab/cs/_cs_base.py

@@ -215,13 +215,14 @@ class SpaceCartesianCS(SpaceCSBase):
 
         # phi: longitude, theta: co-latitude, r: radial distance
         r = np.sqrt(self.coords.x ** 2 + self.coords.y ** 2 + self.coords.z ** 2)
+        r_unit = self.coords.x_unit
         theta = np.arccos(self.coords.z / r)
         phi = npmath.trig_arctan_to_sph_lon(self.coords.x, self.coords.y)
         vector = None
         if self.vector is not None:
             mylog.StreamLogger.warning("The tranformation for vectors have not been implemented!")
 
-        coords = {'theta': theta, 'phi': phi, 'r': r}
+        coords = {'theta': theta, 'phi': phi, 'r': r, 'r_unit': r_unit}
         cs_new = set_cs(name=self.name, coords=coords,  vector=vector, kind='sph', ut=self.ut, **kwargs)
 
         return cs_new

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

@@ -159,47 +159,51 @@ class Dataset(datahub.DatasetSourced):
 
     def fix_geo_lon(self):
         from geospacelab.observatory.orbit.sc_orbit import OrbitPosition_SSCWS
-        from scipy.interpolate import interp1d
+        from scipy.interpolate import CubicSpline
+
+        dts_1 = self['SC_DATETIME'].value.flatten()
+        dt0 = dttool.get_start_of_the_day(self.dt_fr)
+        sectime_1 = [(dt - dt0).total_seconds() for dt in dts_1]
+        glat_1 = self['SC_GEO_LAT'].flatten()
+        glon_1 = self['SC_GEO_LON'].flatten()
+        alt_1 = self['SC_GEO_ALT'].flatten()
+        
         # check outliers
         orbit_obj = OrbitPosition_SSCWS(
             dt_fr=self.dt_fr - datetime.timedelta(minutes=30),
             dt_to=self.dt_to + datetime.timedelta(minutes=30),
             sat_id='dmsp' + self.sat_id.lower()
         )
-
-        glat_1 = self['SC_GEO_LAT'].value.flatten()
-        glon_1 = self['SC_GEO_LON'].value.flatten()
-        if glat_1.size < 2:
-            return
-
-        dts_1 = self['SC_DATETIME'].value.flatten()
-        dt0 = dttool.get_start_of_the_day(self.dt_fr)
-        sectime_1 = [(dt - dt0).total_seconds() for dt in dts_1]
-
-        glat_2 = orbit_obj['SC_GEO_LAT'].value.flatten()
-        glon_2 = orbit_obj['SC_GEO_LON'].value.flatten()
-        dts_2 = orbit_obj['SC_DATETIME'].value.flatten()
-        sectime_2 = [(dt - dt0).total_seconds() for dt in dts_2]
-
-        factor = np.pi / 180.
-        sin_glon_1 = np.sin(glon_1 * factor)
-        sin_glon_2 = np.sin(glon_2 * factor)
-        cos_glon_2 = np.cos(glon_2 * factor)
-        itpf_sin = interp1d(sectime_2, sin_glon_2, kind='cubic', bounds_error=False, fill_value='extrapolate')
-        itpf_cos = interp1d(sectime_2, cos_glon_2, kind='cubic', bounds_error=False, fill_value='extrapolate')
-        sin_glon_2_i = itpf_sin(sectime_1)
-        cos_glon_2_i = itpf_cos(sectime_1)
-        rad = np.sign(sin_glon_2_i) * (np.pi / 2 - np.arcsin(cos_glon_2_i))
-        glon_new = rad / factor
-        # rad = np.where((rad >= 0), rad, rad + 2 * numpy.pi)
-
-        ind_outliers = np.where(np.abs(sin_glon_1 - sin_glon_2_i) > 0.03)[0]
-
+        dts = orbit_obj['SC_DATETIME'].flatten()
+        sectimes, _ = dttool.convert_datetime_to_sectime(dts, dt0=dt0)
+        x = orbit_obj['SC_GEO_X'].flatten()
+        y = orbit_obj['SC_GEO_Y'].flatten()
+        z = orbit_obj['SC_GEO_Z'].flatten()
+        f_x = CubicSpline(sectimes, x)
+        x_i = f_x(sectime_1)
+        f_y = CubicSpline(sectimes, y)
+        y_i = f_y(sectime_1)
+        f_z = CubicSpline(sectimes, z)
+        z_i = f_z(sectime_1)
+        
+        
+        cs = geo_cs.GEOCCartesian(
+            coords={'x': x_i, 'y': y_i, 'z': z_i, 'x_unit': 'km', 'y_unit': 'km', 'z_unit': 'km'}
+        )
+        cs_new = cs.to_spherical()
+        
+        # glat_2 = cs_new['lat']
+        glon_2 = cs_new['lon']
+        # alt_2 = cs_new['height']
+        
+        delta = np.sin(glon_2 * np.pi / 180.) - np.sin(glon_2 * np.pi / 180.)
+        
         if self.replace_orbit:
-            glon_1 = glon_new
+            glon_1 = glon_2
         else:
-            glon_1[ind_outliers] = glon_new[ind_outliers]
-        self['SC_GEO_LON'].value = glon_1.reshape((glon_1.size, 1))
+            glon_1[np.abs(delta)>0.05] = glon_2[np.abs(delta)>0.05]
+            
+        self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
 
     def search_data_files(self, **kwargs):
 

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

@@ -19,6 +19,7 @@ from geospacelab.config import prf
 import geospacelab.toolbox.utilities.pybasic as basic
 import geospacelab.toolbox.utilities.pylogging as mylog
 import geospacelab.toolbox.utilities.pydatetime as dttool
+import geospacelab.cs as geo_cs
 from geospacelab.datahub.sources.madrigal.satellites.dmsp.s1.loader import Loader as default_Loader
 from geospacelab.datahub.sources.madrigal.satellites.dmsp.downloader import Downloader as default_Downloader
 import geospacelab.datahub.sources.madrigal.satellites.dmsp.s1.variable_config as var_config
@@ -179,53 +180,51 @@ class Dataset(datahub.DatasetSourced):
 
     def fix_geo_lon(self):
         from geospacelab.observatory.orbit.sc_orbit import OrbitPosition_SSCWS
-        from scipy.interpolate import interp1d
+        from scipy.interpolate import CubicSpline
+
+        dts_1 = self['SC_DATETIME'].value.flatten()
+        dt0 = dttool.get_start_of_the_day(self.dt_fr)
+        sectime_1 = [(dt - dt0).total_seconds() for dt in dts_1]
+        glat_1 = self['SC_GEO_LAT'].flatten()
+        glon_1 = self['SC_GEO_LON'].flatten()
+        alt_1 = self['SC_GEO_ALT'].flatten()
+        
         # check outliers
         orbit_obj = OrbitPosition_SSCWS(
             dt_fr=self.dt_fr - datetime.timedelta(minutes=30),
             dt_to=self.dt_to + datetime.timedelta(minutes=30),
             sat_id='dmsp' + self.sat_id.lower()
         )
-
-        glat_1 = self['SC_GEO_LAT'].value.flatten()
-        glon_1 = self['SC_GEO_LON'].value.flatten()
-        if glat_1.size < 2:
-            return
-
-        dts_1 = self['SC_DATETIME'].value.flatten()
-        dt0 = dttool.get_start_of_the_day(self.dt_fr)
-        sectime_1 = [(dt - dt0).total_seconds() for dt in dts_1]
-
-        glat_2 = orbit_obj['SC_GEO_LAT'].value.flatten()
-        glon_2 = orbit_obj['SC_GEO_LON'].value.flatten()
-        dts_2 = orbit_obj['SC_DATETIME'].value.flatten()
-        sectime_2 = [(dt - dt0).total_seconds() for dt in dts_2]
-
-        factor = np.pi / 180.
-        sin_glon_1 = np.sin(glon_1 * factor)
-        sin_glon_2 = np.sin(glon_2 * factor)
-        cos_glon_2 = np.cos(glon_2 * factor)
-        itpf_sin = interp1d(sectime_2, sin_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
-        itpf_cos = interp1d(sectime_2, cos_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
-        sin_glon_2_i = itpf_sin(sectime_1)
-        sin_glon_2_i = np.where(sin_glon_2_i > 1., 1., sin_glon_2_i)
-        sin_glon_2_i = np.where(sin_glon_2_i < -1., -1., sin_glon_2_i)
+        dts = orbit_obj['SC_DATETIME'].flatten()
+        sectimes, _ = dttool.convert_datetime_to_sectime(dts, dt0=dt0)
+        x = orbit_obj['SC_GEO_X'].flatten()
+        y = orbit_obj['SC_GEO_Y'].flatten()
+        z = orbit_obj['SC_GEO_Z'].flatten()
+        f_x = CubicSpline(sectimes, x)
+        x_i = f_x(sectime_1)
+        f_y = CubicSpline(sectimes, y)
+        y_i = f_y(sectime_1)
+        f_z = CubicSpline(sectimes, z)
+        z_i = f_z(sectime_1)
         
-        cos_glon_2_i = itpf_cos(sectime_1)
-        cos_glon_2_i = np.where(cos_glon_2_i > 1., 1., cos_glon_2_i)
-        cos_glon_2_i = np.where(cos_glon_2_i < -1., -1., cos_glon_2_i)
         
-        rad = np.sign(sin_glon_2_i) * (np.pi / 2 - np.arcsin(cos_glon_2_i))
-        glon_new = rad / factor
-        # rad = np.where((rad >= 0), rad, rad + 2 * numpy.pi)
-
-        ind_outliers = np.where(np.abs(sin_glon_1 - sin_glon_2_i) > 0.03)[0]
-
+        cs = geo_cs.GEOCCartesian(
+            coords={'x': x_i, 'y': y_i, 'z': z_i, 'x_unit': 'km', 'y_unit': 'km', 'z_unit': 'km'}
+        )
+        cs_new = cs.to_spherical()
+        
+        # glat_2 = cs_new['lat']
+        glon_2 = cs_new['lon']
+        # alt_2 = cs_new['height']
+        
+        delta = np.sin(glon_2 * np.pi / 180.) - np.sin(glon_2 * np.pi / 180.)
+        
         if self.replace_orbit:
-            glon_1 = glon_new
+            glon_1 = glon_2
         else:
-            glon_1[ind_outliers] = glon_new[ind_outliers]
-        self['SC_GEO_LON'].value = glon_1.reshape((glon_1.size, 1))
+            glon_1[np.abs(delta)>0.05] = glon_2[np.abs(delta)>0.05]
+            
+        self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
 
     def search_data_files(self, **kwargs):
 

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

@@ -151,47 +151,94 @@ class Dataset(datahub.DatasetSourced):
 
     def fix_geo_lon(self):
         from geospacelab.observatory.orbit.sc_orbit import OrbitPosition_SSCWS
-        from scipy.interpolate import interp1d
+        from scipy.interpolate import CubicSpline
+
+        dts_1 = self['SC_DATETIME'].value.flatten()
+        dt0 = dttool.get_start_of_the_day(self.dt_fr)
+        sectime_1 = [(dt - dt0).total_seconds() for dt in dts_1]
+        glat_1 = self['SC_GEO_LAT'].flatten()
+        glon_1 = self['SC_GEO_LON'].flatten()
+        alt_1 = self['SC_GEO_ALT'].flatten()
+        
         # check outliers
         orbit_obj = OrbitPosition_SSCWS(
             dt_fr=self.dt_fr - datetime.timedelta(minutes=30),
             dt_to=self.dt_to + datetime.timedelta(minutes=30),
             sat_id='dmsp' + self.sat_id.lower()
         )
-
-        glat_1 = self['SC_GEO_LAT'].value.flatten()
-        glon_1 = self['SC_GEO_LON'].value.flatten()
-        if glat_1.size < 2:
-            return
-
-        dts_1 = self['SC_DATETIME'].value.flatten()
-        dt0 = dttool.get_start_of_the_day(self.dt_fr)
-        sectime_1 = [(dt - dt0).total_seconds() for dt in dts_1]
-
-        glat_2 = orbit_obj['SC_GEO_LAT'].value.flatten()
-        glon_2 = orbit_obj['SC_GEO_LON'].value.flatten()
-        dts_2 = orbit_obj['SC_DATETIME'].value.flatten()
-        sectime_2 = [(dt - dt0).total_seconds() for dt in dts_2]
-
-        factor = np.pi / 180.
-        sin_glon_1 = np.sin(glon_1 * factor)
-        sin_glon_2 = np.sin(glon_2 * factor)
-        cos_glon_2 = np.cos(glon_2 * factor)
-        itpf_sin = interp1d(sectime_2, sin_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
-        itpf_cos = interp1d(sectime_2, cos_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
-        sin_glon_2_i = itpf_sin(sectime_1)
-        cos_glon_2_i = itpf_cos(sectime_1)
-        rad = np.sign(sin_glon_2_i) * (np.pi / 2 - np.arcsin(cos_glon_2_i))
-        glon_new = rad / factor
-        # rad = np.where((rad >= 0), rad, rad + 2 * numpy.pi)
-
-        ind_outliers = np.where(np.abs(sin_glon_1 - sin_glon_2_i) > 0.03)[0]
-
+        dts = orbit_obj['SC_DATETIME'].flatten()
+        sectimes, _ = dttool.convert_datetime_to_sectime(dts, dt0=dt0)
+        x = orbit_obj['SC_GEO_X'].flatten()
+        y = orbit_obj['SC_GEO_Y'].flatten()
+        z = orbit_obj['SC_GEO_Z'].flatten()
+        f_x = CubicSpline(sectimes, x)
+        x_i = f_x(sectime_1)
+        f_y = CubicSpline(sectimes, y)
+        y_i = f_y(sectime_1)
+        f_z = CubicSpline(sectimes, z)
+        z_i = f_z(sectime_1)
+        
+        
+        cs = geo_cs.GEOCCartesian(
+            coords={'x': x_i, 'y': y_i, 'z': z_i, 'x_unit': 'km', 'y_unit': 'km', 'z_unit': 'km'}
+        )
+        cs_new = cs.to_spherical()
+        
+        # glat_2 = cs_new['lat']
+        glon_2 = cs_new['lon']
+        # alt_2 = cs_new['height']
+        
+        delta = np.sin(glon_2 * np.pi / 180.) - np.sin(glon_2 * np.pi / 180.)
+        
         if self.replace_orbit:
-            glon_1 = glon_new
+            glon_1 = glon_2
         else:
-            glon_1[ind_outliers] = glon_new[ind_outliers]
-        self['SC_GEO_LON'].value = glon_1.reshape((glon_1.size, 1))
+            glon_1[np.abs(delta)>0.05] = glon_2[np.abs(delta)>0.05]
+            
+        self['SC_GEO_LON'].value = glon_1[:, np.newaxis]
+        
+        # from geospacelab.observatory.orbit.sc_orbit import OrbitPosition_SSCWS
+        # from scipy.interpolate import interp1d
+        # # check outliers
+        # orbit_obj = OrbitPosition_SSCWS(
+        #     dt_fr=self.dt_fr - datetime.timedelta(minutes=30),
+        #     dt_to=self.dt_to + datetime.timedelta(minutes=30),
+        #     sat_id='dmsp' + self.sat_id.lower()
+        # )
+
+        # glat_1 = self['SC_GEO_LAT'].value.flatten()
+        # glon_1 = self['SC_GEO_LON'].value.flatten()
+        # if glat_1.size < 2:
+        #     return
+
+        # dts_1 = self['SC_DATETIME'].value.flatten()
+        # dt0 = dttool.get_start_of_the_day(self.dt_fr)
+        # sectime_1 = [(dt - dt0).total_seconds() for dt in dts_1]
+
+        # glat_2 = orbit_obj['SC_GEO_LAT'].value.flatten()
+        # glon_2 = orbit_obj['SC_GEO_LON'].value.flatten()
+        # dts_2 = orbit_obj['SC_DATETIME'].value.flatten()
+        # sectime_2 = [(dt - dt0).total_seconds() for dt in dts_2]
+
+        # factor = np.pi / 180.
+        # sin_glon_1 = np.sin(glon_1 * factor)
+        # sin_glon_2 = np.sin(glon_2 * factor)
+        # cos_glon_2 = np.cos(glon_2 * factor)
+        # itpf_sin = interp1d(sectime_2, sin_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
+        # itpf_cos = interp1d(sectime_2, cos_glon_2, kind='linear', bounds_error=False, fill_value='extrapolate')
+        # sin_glon_2_i = itpf_sin(sectime_1)
+        # cos_glon_2_i = itpf_cos(sectime_1)
+        # rad = np.sign(sin_glon_2_i) * (np.pi / 2 - np.arcsin(cos_glon_2_i))
+        # glon_new = rad / factor
+        # # rad = np.where((rad >= 0), rad, rad + 2 * numpy.pi)
+
+        # ind_outliers = np.where(np.abs(sin_glon_1 - sin_glon_2_i) > 0.03)[0]
+
+        # if self.replace_orbit:
+        #     glon_1 = glon_new
+        # else:
+        #     glon_1[ind_outliers] = glon_new[ind_outliers]
+        # self['SC_GEO_LON'].value = glon_1.reshape((glon_1.size, 1))
 
     def search_data_files(self, **kwargs):
 

geospacelab/toolbox/utilities/numpyarray.py

@@ -270,6 +270,8 @@ def data_resample(
 #             grid_z_1[i, :] = z_i
 
 
+    
+
 def regridding_2d_xgaps(
     x, y, z,
     xtype=None, xres=None):
@@ -281,13 +283,28 @@ def regridding_2d_xgaps(
         x1 = sectime
     
     diff_x = numpy.diff(x1)
+    x_res_md = numpy.median(diff_x)
     if xres is None:
-        xres = numpy.median(diff_x)
+        xres = x_res_md
         
     if xtype == 'datetime' and basic.isnumeric(xres):
         xres = datetime.timedelta(seconds=xres)
+        x_res_md = datetime.timedelta(seconds=x_res_md)
+
+    xx_1 = np.empty_like(x)
+    xx_2 = np.empty_like(x)
+    for i, xxx in enumerate(x):
+        if i == 0:
+            xx_1[0] = xxx - np.min([datetime.timedelta(diff_x[0]), xres, x_res_md]) / 2
+        else:
+            xx_1[i] = xx_2[i-1] if xres > datetime.timedelta(seconds=diff_x[i-1]) else xxx - xres / 2
         
-    xx = numpy.hstack(((x.flatten() - xres / 2)[:, numpy.newaxis], (x.flatten() + xres / 2)[:, numpy.newaxis]))
+        if i == len(x) - 1:
+            xx_2[i] = xxx + np.min([datetime.timedelta(diff_x[i-1]), xres, x_res_md]) / 2 
+        else:
+            xx_2[i] = xxx + np.min([datetime.timedelta(diff_x[i]), xres, x_res_md]) / 2  
+    
+    xx = numpy.hstack((xx_1[:, numpy.newaxis], xx_2[:, numpy.newaxis]))
     xnew = xx.flatten()
     
     if len(y.shape) == 1 or 1 in y.shape:

{geospacelab-0.11.0.dist-info → geospacelab-0.11.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: geospacelab
-Version: 0.11.0
+Version: 0.11.2
 Summary: Collect, manage, and visualize geospace data.
 Home-page: https://github.com/JouleCai/geospacelab
 Author: Lei Cai

{geospacelab-0.11.0.dist-info → geospacelab-0.11.2.dist-info}/RECORD

@@ -1,4 +1,4 @@
-geospacelab/__init__.py,sha256=o9UGHuHohTJI5hQ6S62CZ35dSnAnkD-M5-ayloaqsLU,801
+geospacelab/__init__.py,sha256=z281cx-p62AcwL0TBdP-2pheA_u9DVw85sGurChtXrc,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
@@ -7,7 +7,7 @@ geospacelab/coords/geo_utilities.py,sha256=7_C98Cd29tzNl86kdWYmVC8UBS8E07oJIR6Ma
 geospacelab/cs/__init__.py,sha256=lzPhfypd0EfsOJnKsrrK2dQo3sxj8lpblYEhpz2i9xA,1282
 geospacelab/cs/_aacgm.py,sha256=6NmY_CVnECs-xZo1y2q3eglBc-ZrK-0SWtj6247Xq1E,723
 geospacelab/cs/_apex.py,sha256=9S-I0MMV5aU_MBO6Pi0AnFCGumfwIravW_VnGONvGYE,756
-geospacelab/cs/_cs_base.py,sha256=iZSlChh_tn6bw2k71c_LFs3jIvWy91ClxpOUO9MYYfg,17124
+geospacelab/cs/_cs_base.py,sha256=bERqTF1uNLUFA3F1ySccyJt-t1Dhu1_B1E7-prrgV5o,17178
 geospacelab/cs/_geo.py,sha256=cLrb2r-xD5aVwKENUzaFVV71tPpOC21zCmUIQmf67eg,18242
 geospacelab/cs/geo_utilities.py,sha256=RWsioSB6_n3bT502ONho8eyJaFfFVmNBv1T0RwOPD64,604
 geospacelab/datahub/__dataset_base__.py,sha256=XPcwlfIXNFeoqQoe_PNNNfKW8pPI7rEz3bdrVqLPSCQ,20616
@@ -196,13 +196,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=jnA8XpiD-bhi5lM1p1er93D-kS7aE8XgYZVVGcE0Aho,4886
-geospacelab/datahub/sources/madrigal/satellites/dmsp/e/__init__.py,sha256=2LmuPOzOH9lhMqXkrb4gj9Wxsmdq2mXSA0phTlrK0eU,10309
+geospacelab/datahub/sources/madrigal/satellites/dmsp/e/__init__.py,sha256=Csf-c864DxJXDjq9TL5YDF-8A97zVIdT8K8jerU45I4,10212
 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=u7Z0SFUaLNzPV7aIqcaEUTyIIYPHHJF3rXMWCDYbSUY,11536
+geospacelab/datahub/sources/madrigal/satellites/dmsp/s1/__init__.py,sha256=iPtRE0QJ-Fq6Bit5xRRDXYigQr_DKl6DA7F_GRFofuU,11171
 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=0Qkn18pmVBmc0y0Kk6TRJBynoU_ysYMsXLzAAXB37M8,10198
+geospacelab/datahub/sources/madrigal/satellites/dmsp/s4/__init__.py,sha256=85NIwKVTIB8nCmwzV_qk6dKeC4glIddjha6KQZXba64,12043
 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
@@ -331,7 +331,7 @@ geospacelab/toolbox/io/__init__.py,sha256=YcFkmOGWXNiOXWsMFxFPhn32UnbxmeSrluTp6O
 geospacelab/toolbox/io/dialog.py,sha256=y87PKCXQS5_axbrNuVrcYN2uIfr-ZHICr7yMQ0yH1M4,1067
 geospacelab/toolbox/unit/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 geospacelab/toolbox/utilities/__init__.py,sha256=AbpHGcgLb-kRsJGnwFEktk7uzpZOCcBY74-YBdrKVGs,1
-geospacelab/toolbox/utilities/numpyarray.py,sha256=olclT06DRPFSYaEdiD1uqwHvM0vHNDFjkLstck6_clE,11862
+geospacelab/toolbox/utilities/numpyarray.py,sha256=Sr7j5Kfe_vtc185GUX99olkML1Ld81jN6YmrNPX1wcM,12448
 geospacelab/toolbox/utilities/numpymath.py,sha256=xr-u9vsN5UMZ6pjomJ2O3A1n00Eg3yH5fN4lcY7rR_Q,3810
 geospacelab/toolbox/utilities/pybasic.py,sha256=-wLWoQMTroJ49nn8H_-hsFxDSiHXMiea86U6u0W8_cM,3310
 geospacelab/toolbox/utilities/pyclass.py,sha256=DnnzLH1ovlxArtqCpywZUfDoNzYi7HmMW2s17Xd30LI,2813
@@ -382,7 +382,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.11.0.dist-info/licenses/LICENSE,sha256=2yRlwLt4o5Z6OZAGcyvBj-zfFX1Uw7E6CzqODg7khqs,1515
+geospacelab-0.11.2.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
@@ -392,7 +392,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.11.0.dist-info/METADATA,sha256=IM_FXqyGJexx3lJm9sbC7H_HSJxbkRinuHpiAa6NxBg,24238
-geospacelab-0.11.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-geospacelab-0.11.0.dist-info/top_level.txt,sha256=98eDwrSNgyQFAtSA06QMP71gw9BzgIj0uvkTudpGly4,12
-geospacelab-0.11.0.dist-info/RECORD,,
+geospacelab-0.11.2.dist-info/METADATA,sha256=nnc9SDt62JaF_T7Aq4gqQrFiXIhdB4_zy-tEw-DAexU,24238
+geospacelab-0.11.2.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+geospacelab-0.11.2.dist-info/top_level.txt,sha256=98eDwrSNgyQFAtSA06QMP71gw9BzgIj0uvkTudpGly4,12
+geospacelab-0.11.2.dist-info/RECORD,,