imap-processing 0.15.0__py3-none-any.whl → 0.16.1__py3-none-any.whl

imap_processing/ialirt/l0/process_swapi.py

@@ -1,20 +1,144 @@
 """Functions to support I-ALiRT SWAPI processing."""
 
 import logging
+from decimal import Decimal
+from typing import Optional
 
 import numpy as np
+import pandas as pd
 import xarray as xr
-from xarray import DataArray
+from scipy.optimize import curve_fit
+from scipy.special import erf
 
+from imap_processing import imap_module_directory
+from imap_processing.ialirt.constants import IalirtSwapiConstants as Consts
 from imap_processing.ialirt.utils.grouping import find_groups
-
-# from imap_processing.swapi.l1.swapi_l1 import process_sweep_data
-# from imap_processing.swapi.l2.swapi_l2 import TIME_PER_BIN
+from imap_processing.ialirt.utils.time import calculate_time
+from imap_processing.spice.time import met_to_ttj2000ns, met_to_utc
+from imap_processing.swapi.l1.swapi_l1 import process_sweep_data
+from imap_processing.swapi.l2.swapi_l2 import TIME_PER_BIN
 
 logger = logging.getLogger(__name__)
 
 
-def process_swapi_ialirt(unpacked_data: xr.Dataset) -> dict[str, DataArray]:
+def count_rate(
+    energy_pass: float, speed: float, density: float, temp: float
+) -> float | np.ndarray:
+    """
+    Compute SWAPI count rate for provided energy passband, speed, density and temp.
+
+    This model for coincidence count rate was developed by the SWAPI instrument
+    science team, detailed on page 52 of the IMAP SWAPI Instrument Algorithms Document.
+
+    Parameters
+    ----------
+    energy_pass : float
+        Energy passband [eV].
+    speed : float
+        Bulk solar wind speed [km/s].
+    density : float
+        Proton density [cm^-3].
+    temp : float
+        Temperature [K].
+
+    Returns
+    -------
+    count_rate : float | np.ndarray
+        Particle coincidence count rate.
+    """
+    # thermal velocity of solar wind ions
+    thermal_velocity = np.sqrt(2 * Consts.boltz * temp / Consts.prot_mass)
+    beta = 1 / (thermal_velocity**2)
+    # convert energy to Joules
+    center_speed = np.sqrt(2 * energy_pass * 1.60218e-19 / Consts.prot_mass)
+    speed = speed * 1000  # convert km/s to m/s
+    density = density * 1e6  # convert 1/cm**3 to 1/m**3
+
+    return (
+        (density * Consts.eff_area * (beta / np.pi) ** (3 / 2))
+        * (np.exp(-beta * (center_speed**2 + speed**2 - 2 * center_speed * speed)))
+        * np.sqrt(np.pi / (beta * speed * center_speed))
+        * erf(np.sqrt(beta * speed * center_speed) * (Consts.az_fov / 2))
+        * (
+            center_speed**4
+            * Consts.speed_ew
+            * np.arcsin(thermal_velocity / center_speed)
+        )
+    )
+
+
+def optimize_pseudo_parameters(
+    count_rates: np.ndarray,
+    count_rate_error: np.ndarray,
+    energy_passbands: Optional[np.ndarray] = None,
+) -> (dict)[str, list[float]]:
+    """
+    Find the pseudo speed (u), density (n) and temperature (T) of solar wind particles.
+
+    Fit a curve to calculated count rate values as a function of energy passband.
+
+    Parameters
+    ----------
+    count_rates : np.ndarray
+        Particle coincidence count rates.
+    count_rate_error : np.ndarray
+        Standard deviation of the coincidence count rates parameter.
+    energy_passbands : np.ndarray, default None
+        Energy passbands, passed in only for testing purposes.
+
+    Returns
+    -------
+    solution_dict : dict
+        Dictionary containing the optimized speed, density, and temperature values for
+        each sweep included in the input count_rates array.
+    """
+    if not energy_passbands:
+        # Read in energy passbands
+        energy_data = pd.read_csv(
+            f"{imap_module_directory}/tests/swapi/lut/imap_swapi_esa-unit"
+            f"-conversion_20250211_v000.csv"
+        )
+        energy_passbands = (
+            energy_data["Energy"][0:63]
+            .replace(",", "", regex=True)
+            .to_numpy()
+            .astype(float)
+        )
+
+    # Initial guess pulled from page 52 of the IMAP SWAPI Instrument Algorithms Document
+    initial_param_guess = np.array([550, 5.27, 1e5])
+    solution_dict = {  # type: ignore
+        "pseudo_speed": [],
+        "pseudo_density": [],
+        "pseudo_temperature": [],
+    }
+
+    for sweep in np.arange(count_rates.shape[0]):
+        current_sweep_count_rates = count_rates[sweep, :]
+        current_sweep_count_rate_errors = count_rate_error[sweep, :]
+        # Find the max count rate, and use the 6 points surrounding it (inclusive)
+        max_index = np.argmax(current_sweep_count_rates)
+        sol = curve_fit(
+            f=count_rate,
+            xdata=energy_passbands.take(
+                range(max_index - 3, max_index + 3), mode="wrap"
+            ),
+            ydata=current_sweep_count_rates.take(
+                range(max_index - 3, max_index + 3), mode="wrap"
+            ),
+            sigma=current_sweep_count_rate_errors.take(
+                range(max_index - 3, max_index + 3), mode="wrap"
+            ),
+            p0=initial_param_guess,
+        )
+        solution_dict["pseudo_speed"].append(sol[0][0])
+        solution_dict["pseudo_density"].append(sol[0][1])
+        solution_dict["pseudo_temperature"].append(sol[0][2])
+
+    return solution_dict
+
+
+def process_swapi_ialirt(unpacked_data: xr.Dataset) -> list[dict]:
     """
     Extract I-ALiRT variables and calculate coincidence count rate.
 
@@ -32,7 +156,21 @@ def process_swapi_ialirt(unpacked_data: xr.Dataset) -> dict[str, DataArray]:
 
     sci_dataset = unpacked_data.sortby("epoch", ascending=True)
 
-    grouped_dataset = find_groups(sci_dataset, (0, 11), "swapi_seq_number", "swapi_acq")
+    met = calculate_time(
+        sci_dataset["sc_sclk_sec"], sci_dataset["sc_sclk_sub_sec"], 256
+    )
+
+    # Add required parameters.
+    sci_dataset["met"] = met
+    met_values = []
+
+    grouped_dataset = find_groups(sci_dataset, (0, 11), "swapi_seq_number", "met")
+
+    if grouped_dataset.group.size == 0:
+        logger.warning(
+            "There was an issue with the SWAPI grouping process, returning empty data."
+        )
+        return []
 
     for group in np.unique(grouped_dataset["group"]):
         # Sequence values for the group should be 0-11 with no duplicates.
@@ -40,6 +178,10 @@ def process_swapi_ialirt(unpacked_data: xr.Dataset) -> dict[str, DataArray]:
             (grouped_dataset["group"] == group)
         ]
 
+        met_values.append(
+            int(grouped_dataset["met"][(grouped_dataset["group"] == group).values][0])
+        )
+
         # Ensure no duplicates and all values from 0 to 11 are present
         if not np.array_equal(seq_values.astype(int), np.arange(12)):
             logger.info(
@@ -48,22 +190,29 @@ def process_swapi_ialirt(unpacked_data: xr.Dataset) -> dict[str, DataArray]:
             )
             continue
 
-    total_packets = len(grouped_dataset["swapi_seq_number"].data)
-
-    # It takes 12 sequence data to make one full SWAPI sweep
-    total_sequence = 12
-    total_full_sweeps = total_packets // total_sequence
-
-    met_values = grouped_dataset["swapi_shcoarse"].data.reshape(total_full_sweeps, 12)[
-        :, 0
-    ]
-
-    # raw_coin_count = process_sweep_data(grouped_dataset, "coin_cnt")
-    # raw_coin_rate = raw_coin_count / TIME_PER_BIN
-
-    swapi_data = {
-        "met": met_values
-        # more variables to go here
-    }
+    raw_coin_count = process_sweep_data(grouped_dataset, "swapi_coin_cnt")
+    raw_coin_rate = raw_coin_count / TIME_PER_BIN
+    count_rate_error = np.sqrt(raw_coin_count) / TIME_PER_BIN
+
+    solution = optimize_pseudo_parameters(raw_coin_rate, count_rate_error)
+
+    swapi_data = []
+
+    for entry in np.arange(0, len(solution["pseudo_speed"])):
+        swapi_data.append(
+            {
+                "apid": 478,
+                "met": int(met_values[entry]),
+                "met_in_utc": met_to_utc(met_values[entry]).split(".")[0],
+                "ttj2000ns": int(met_to_ttj2000ns(met_values[entry])),
+                "swapi_pseudo_proton_speed": Decimal(solution["pseudo_speed"][entry]),
+                "swapi_pseudo_proton_density": Decimal(
+                    solution["pseudo_density"][entry]
+                ),
+                "swapi_pseudo_proton_temperature": Decimal(
+                    solution["pseudo_temperature"][entry]
+                ),
+            }
+        )
 
     return swapi_data

imap_processing/ialirt/l0/process_swe.py

@@ -1,7 +1,6 @@
 """Functions to support I-ALiRT SWE processing."""
 
 import logging
-from decimal import Decimal
 
 import numpy as np
 import pandas as pd
@@ -547,16 +546,10 @@ def process_swe(accumulated_data: xr.Dataset, in_flight_cal_files: list) -> list
             {
                 "apid": 478,
                 "met": int(grouped["met"].min()),
-                "utc": met_to_utc(grouped["met"].min()).split(".")[0],
+                "met_in_utc": met_to_utc(grouped["met"].min()).split(".")[0],
                 "ttj2000ns": int(met_to_ttj2000ns(grouped["met"].min())),
-                **{
-                    f"swe_normalized_counts_half_1_esa_{i}": Decimal(str(val))
-                    for i, val in enumerate(summed_first)
-                },
-                **{
-                    f"swe_normalized_counts_half_2_esa_{i}": Decimal(str(val))
-                    for i, val in enumerate(summed_second)
-                },
+                "swe_normalized_counts_half_1_esa": [int(val) for val in summed_first],
+                "swe_normalized_counts_half_2_esa": [int(val) for val in summed_second],
                 "swe_counterstreaming_electrons": max(bde_first_half, bde_second_half),
             }
         )

imap_processing/ialirt/utils/constants.py

@@ -1,6 +1,20 @@
 """Keys for I-ALiRT data products."""
 
 IALIRT_KEYS = [
+    # H intensities in 15 energy ranges and binned into 4 azimuths and 4 spin angle bins
+    "codicehi_h",
+    # C/O abundance ratio
+    "codicelo_c_over_o_abundance",
+    # Mg/O abundance ratio
+    "codicelo_mg_over_o_abundance",
+    # Fe/O abundance ratio
+    "codicelo_fe_over_o_abundance",
+    # C+6/C+5 charge state ratio
+    "codicelo_c_plus_6_over_c_plus_5_ratio",
+    # O+7/O+6 charge state ratio
+    "codicelo_o_plus_7_over_o_plus_6_ratio",
+    # Fe low/Fe high charge state ratio
+    "codicelo_fe_low_over_fe_high_ratio",
     # Low energy (~300 keV) electrons (A-side)
     "hit_e_a_side_low_en",
     # Medium energy (~3 MeV) electrons (A-side)
@@ -40,9 +54,9 @@ IALIRT_KEYS = [
     # Pseudo temperature of solar wind protons in plasma frame
     "swapi_pseudo_proton_temperature",
     # SWE Normalized Counts - Half Cycle 1
-    *[f"swe_normalized_counts_half_1_esa_{i}" for i in range(8)],
+    "swe_normalized_counts_half_1",
     # SWE Normalized Counts - Half Cycle 2
-    *[f"swe_normalized_counts_half_2_esa_{i}" for i in range(8)],
+    "swe_normalized_counts_half_2",
     # SWE Counterstreaming flag
     "swe_counterstreaming_electrons",
 ]

imap_processing/ialirt/utils/create_xarray.py

@@ -7,7 +7,7 @@ from imap_processing.cdf.imap_cdf_manager import ImapCdfAttributes
 from imap_processing.ialirt.utils.constants import IALIRT_KEYS
 
 
-def create_xarray_from_records(records: list[dict]) -> xr.Dataset:
+def create_xarray_from_records(records: list[dict]) -> xr.Dataset:  # noqa: PLR0912
     """
     Create dataset from a list of records.
 
@@ -39,16 +39,51 @@ def create_xarray_from_records(records: list[dict]) -> xr.Dataset:
         data=ttj2000ns_values,
         name="epoch",
         dims=["epoch"],
-        attrs=cdf_manager.get_variable_attributes("epoch"),
+        attrs=cdf_manager.get_variable_attributes("epoch", check_schema=False),
     )
     component = xr.DataArray(
         ["x", "y", "z"],
         name="component",
         dims=["component"],
-        attrs=cdf_manager.get_variable_attributes("component"),
+        attrs=cdf_manager.get_variable_attributes("component", check_schema=False),
     )
 
-    coords = {"epoch": epoch, "component": component}
+    esa_step = xr.DataArray(
+        data=np.arange(8, dtype=np.uint8),
+        name="esa_step",
+        dims=["esa_step"],
+        attrs=cdf_manager.get_variable_attributes("esa_step", check_schema=False),
+    )
+
+    energy_ranges = xr.DataArray(
+        data=np.arange(15, dtype=np.uint8),
+        name="energy_ranges",
+        dims=["energy_ranges"],
+        attrs=cdf_manager.get_variable_attributes("energy_ranges", check_schema=False),
+    )
+
+    azimuth = xr.DataArray(
+        data=np.arange(4, dtype=np.uint8),
+        name="azimuth",
+        dims=["azimuth"],
+        attrs=cdf_manager.get_variable_attributes("azimuth", check_schema=False),
+    )
+
+    spin_angle_bin = xr.DataArray(
+        data=np.arange(4, dtype=np.uint8),
+        name="spin_angle_bin",
+        dims=["spin_angle_bin"],
+        attrs=cdf_manager.get_variable_attributes("spin_angle_bin", check_schema=False),
+    )
+
+    coords = {
+        "epoch": epoch,
+        "component": component,
+        "esa_step": esa_step,
+        "energy_ranges": energy_ranges,
+        "azimuth": azimuth,
+        "spin_angle_bin": spin_angle_bin,
+    }
     dataset = xr.Dataset(
         coords=coords,
         attrs=cdf_manager.get_global_attributes("imap_ialirt_l1_realtime"),
@@ -56,17 +91,25 @@ def create_xarray_from_records(records: list[dict]) -> xr.Dataset:
 
     # Create empty dataset for each key.
     for key in instrument_keys:
-        attrs = cdf_manager.get_variable_attributes(key)
+        attrs = cdf_manager.get_variable_attributes(key, check_schema=False)
         fillval = attrs.get("FILLVAL")
         if key.startswith("mag"):
             data = np.full((n, 3), fillval, dtype=np.float32)
             dims = ["epoch", "component"]
             dataset[key] = xr.DataArray(data, dims=dims, attrs=attrs)
+        elif key.startswith("codicehi"):
+            data = np.full((n, 15, 4, 4), fillval, dtype=np.float32)
+            dims = ["epoch", "energy", "azimuth", "spin_angle_bin"]
+            dataset[key] = xr.DataArray(data, dims=dims, attrs=attrs)
         elif key == "swe_counterstreaming_electrons":
             data = np.full(n, fillval, dtype=np.uint8)
             dims = ["epoch"]
             dataset[key] = xr.DataArray(data, dims=dims, attrs=attrs)
-        elif key.startswith(("hit", "swe")):
+        elif key.startswith("swe"):
+            data = np.full((n, 8), fillval, dtype=np.uint32)
+            dims = ["epoch", "esa_step"]
+            dataset[key] = xr.DataArray(data, dims=dims, attrs=attrs)
+        elif key.startswith("hit"):
             data = np.full(n, fillval, dtype=np.uint32)
             dims = ["epoch"]
             dataset[key] = xr.DataArray(data, dims=dims, attrs=attrs)
@@ -77,11 +120,16 @@ def create_xarray_from_records(records: list[dict]) -> xr.Dataset:
 
     # Populate the dataset variables
     for i, record in enumerate(records):
-        for key in record.keys():
-            val = record[key]
-            if key.startswith("mag"):
-                dataset[key].data[i] = [direction for direction in val]
-            elif key in instrument_keys:
+        for key, val in record.items():
+            if key in ["apid", "met", "met_in_utc", "ttj2000ns"]:
+                continue
+            elif key.startswith("mag"):
+                dataset[key].data[i, :] = val
+            elif key.startswith("swe_normalized_counts"):
+                dataset[key].data[i, :] = val
+            elif key.startswith("codicehi"):
+                dataset[key].data[i, :, :, :] = val
+            else:
                 dataset[key].data[i] = val
 
     return dataset

imap_processing/ultra/utils/ultra_l1_utils.py

@@ -57,14 +57,16 @@ def create_dataset(  # noqa: PLR0912
             data_dict["epoch"],
             name="epoch",
             dims=["epoch"],
-            attrs=cdf_manager.get_variable_attributes("epoch"),
+            attrs=cdf_manager.get_variable_attributes("epoch", check_schema=False),
         )
         if "sensor-de" in name:
             component = xr.DataArray(
                 ["vx", "vy", "vz"],
                 name="component",
                 dims=["component"],
-                attrs=cdf_manager.get_variable_attributes("component"),
+                attrs=cdf_manager.get_variable_attributes(
+                    "component", check_schema=False
+                ),
             )
             coords = {"epoch": epoch_time, "component": component}
         else:

{imap_processing-0.15.0.dist-info → imap_processing-0.16.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: imap-processing
-Version: 0.15.0
+Version: 0.16.1
 Summary: IMAP Science Operations Center Processing
 License: MIT
 Keywords: IMAP,SDC,SOC,Science Operations

{imap_processing-0.15.0.dist-info → imap_processing-0.16.1.dist-info}/RECORD

@@ -1,5 +1,5 @@
 imap_processing/__init__.py,sha256=b9xHlf8_0OtN_OyhPlrC6ayahYR0QTt_e713NjxZObw,1305
-imap_processing/_version.py,sha256=irZ9t9UpzoARk751rz8rjjmcOizvIGEr5-xJe0qW6KE,127
+imap_processing/_version.py,sha256=ab-yNmCvJwJgN5bAYEJTLlnA0k6zw7i42yU_-iPzIzg,127
 imap_processing/ancillary/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 imap_processing/ancillary/ancillary_dataset_combiner.py,sha256=Pxg1wQLjPKzEkgE3a4Tart5TfeuH9rqyCKTd5GbgiL0,9602
 imap_processing/ccsds/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -22,11 +22,11 @@ imap_processing/cdf/config/imap_glows_l2_variable_attrs.yaml,sha256=eHGzS5_ctJ3g
 imap_processing/cdf/config/imap_hi_global_cdf_attrs.yaml,sha256=28VVPcB_dRsk24Rba-nVsyivixnCOBEEHqL6oUu61Vc,2509
 imap_processing/cdf/config/imap_hi_variable_attrs.yaml,sha256=VsIpaOP7SrwsWqiFo-Dhcz1iBVnZE7PmzzVaGwKZ_YI,14678
 imap_processing/cdf/config/imap_hit_global_cdf_attrs.yaml,sha256=zjHAXLlXrrQPgzKh3aTuoY1qb41J0i8zQ81xLTf7Y7g,3121
-imap_processing/cdf/config/imap_hit_l1a_variable_attrs.yaml,sha256=5X8V6Hzg7UGP6HeqHGKy-ogJP3K-V4sLi5271N8fBNw,11281
+imap_processing/cdf/config/imap_hit_l1a_variable_attrs.yaml,sha256=nVphxZP1hc8EHgR9K0Qh7HZd-G8iM0C0gSXeEyRdx0M,48050
 imap_processing/cdf/config/imap_hit_l1b_variable_attrs.yaml,sha256=gxeC7-SUKkWiy1OE6MgfKxeV-0bdXuauodOIfFtQ3JQ,11287
 imap_processing/cdf/config/imap_hit_l2_variable_attrs.yaml,sha256=f8ZLlF9EIes8BcFrpKTO8Mb_VVnP2d9HoLtQVOgz23Q,61402
 imap_processing/cdf/config/imap_ialirt_global_cdf_attrs.yaml,sha256=AGQ9J0zEdWRZ_AR6_AYYHkRb57mYIe2jVflhYv8SEsQ,961
-imap_processing/cdf/config/imap_ialirt_l1_variable_attrs.yaml,sha256=Hl3VE_JRPPDWu9_KDE_zUu0CVGTSYOvIR31gFoyD0D0,9891
+imap_processing/cdf/config/imap_ialirt_l1_variable_attrs.yaml,sha256=IWJwYEeIIWiL1EN8EufKc4O8BThHNDwL2CcfQWRcSPA,8362
 imap_processing/cdf/config/imap_idex_global_cdf_attrs.yaml,sha256=b0dbQ_IJtgTHmwdD754ZHr2-kwo8KIoG3FUZdK5vMEA,2920
 imap_processing/cdf/config/imap_idex_l1a_variable_attrs.yaml,sha256=KaRu-HPqAdAJ2vrX7m54pKUJsJtPaCfofi2OBs906cs,22493
 imap_processing/cdf/config/imap_idex_l1b_variable_attrs.yaml,sha256=Cx-s8hGmUtM_7pon1EbJTP3Im4uFfNeN6gZlGW5DudE,11153
@@ -57,13 +57,13 @@ imap_processing/cdf/config/imap_ultra_l1c_variable_attrs.yaml,sha256=D4uKUqZ8G5h
 imap_processing/cdf/config/imap_variable_schema.yaml,sha256=oKfSQKJ7vYh7SC9dwFeE5PzVQNXYCpcp3AZSclKUVeg,728
 imap_processing/cdf/imap_cdf_manager.py,sha256=K2W_5SzRTxJij2D-ZoKA4UdUxcVC1H90soEloJehsUk,2086
 imap_processing/cdf/utils.py,sha256=Go43WgJbp-BOH_5cLNAru76cCb6jW0ZzwwBCEmeV19E,7789
-imap_processing/cli.py,sha256=CzuVj86GDcBYBMsTt3_dXowiNN8eu9JK9H-vQtSP8Ls,53026
+imap_processing/cli.py,sha256=6kvGorFaUWQeWt43INKbFCH6J7JUwZHSCMQdDTQs8-k,52921
 imap_processing/codice/__init__.py,sha256=Ga6Ioj68HnQ4moi7SxJdOyYmGHcRvKOh4uZUEPJZHl4,20
 imap_processing/codice/codice_l0.py,sha256=_StmC58uRLBuUKSpVaqv6Bo49CeexkcOqnt_askKtHk,1701
-imap_processing/codice/codice_l1a.py,sha256=fCCwAdqGU3RWflHTYl8QYq3J9UvpPaR8IWSQ4V-7L9w,56366
+imap_processing/codice/codice_l1a.py,sha256=3sOzARFDNpL2JhNH2PugPMuMFkOUNMcI0I3T6jt8QhQ,57489
 imap_processing/codice/codice_l1b.py,sha256=_x--ZXHxxSknEC546hHQGqhRTanYtJTfKDsB8d9Pz_Q,5408
-imap_processing/codice/codice_l2.py,sha256=2UsNIWpsd6DRHlaFYF4AykaaCqPcoRln9dLhzyYMQcw,2947
-imap_processing/codice/constants.py,sha256=QOo0YVbZuTlWopRgbx4SPLhiSqXXeohMBa6bLACZrE0,45328
+imap_processing/codice/codice_l2.py,sha256=6bWXaEZxE-LuLv9vlvTIDlCxRWhkYMva0Re7X7-MS_4,2603
+imap_processing/codice/constants.py,sha256=nND2C9KVVkQ1_ZmmL3T9edqIF_gdqQ8h4NyEDVr7_HY,45621
 imap_processing/codice/data/esa_sweep_values.csv,sha256=RuRTBG_TftVKGDEslxDo5SVwLIWq2uh3WrlVsXABpfg,5336
 imap_processing/codice/data/lo_stepping_values.csv,sha256=Z-hdv2nxTzl-RXBWTa7Y_NRXR1HyW1mJrYMmX6vDnHo,5318
 imap_processing/codice/decompress.py,sha256=-EA7vlGIKx6d7gltfDLr8XmU-D6Q9_7FZRtzIUgW54M,5232
@@ -96,41 +96,35 @@ imap_processing/glows/packet_definitions/P_GLX_TMSCHIST.xml,sha256=Q7X5gF9jM4TA5
 imap_processing/glows/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 imap_processing/glows/utils/constants.py,sha256=REjvd0NgnST_rsfFQNGDo0F3_KKDZz20Dyt-19T2txM,2837
 imap_processing/hi/__init__.py,sha256=Ga6Ioj68HnQ4moi7SxJdOyYmGHcRvKOh4uZUEPJZHl4,20
-imap_processing/hi/l1a/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/hi/l1a/hi_l1a.py,sha256=Nj0V9XJZXchdMwsZEyztS2IY_Y0ENMHEPqWUum2VZOQ,3247
-imap_processing/hi/l1a/histogram.py,sha256=ptwPrVDSVjVWJKww3W1XGUNh0Hw9A56_aRZpSANCAwI,4791
-imap_processing/hi/l1a/science_direct_event.py,sha256=dd0-Nf1kTd35GMqVUbidLSLP6b8JVlTzBrJ4yyCNUWA,7150
-imap_processing/hi/l1b/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/hi/l1b/hi_l1b.py,sha256=IvHoTqtAf4u8zemWQ41frARR_gwbQboskp-7hk7AltY,14777
-imap_processing/hi/l1c/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/hi/l1c/hi_l1c.py,sha256=iHOmVZkKEhxwwne1T5tFMxJZjD2n5wI0-WueumTwA50,28459
-imap_processing/hi/l2/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/hi/l2/hi_l2.py,sha256=9Cd1FqdKk3ya6_u9f7gal-9pIxrLgEQT9EBREJ_DhOA,9288
+imap_processing/hi/hi_l1a.py,sha256=m-KGuWw6NzRNuJJp8RVQI1zy636GS7HkXjJNSedPlJ0,14731
+imap_processing/hi/hi_l1b.py,sha256=PegKVakMz_M73eenywiNmna_LlHPLCqRzAKkoEA5uR8,14759
+imap_processing/hi/hi_l1c.py,sha256=lT_XFSRME1ol7Bo2y38RCHghEEXeW21oQ2Rcj2fmjDo,28529
+imap_processing/hi/hi_l2.py,sha256=veO2FdinuiRcONQ5JucnRj8lvs6wwLVrA_oveqsFzFE,9360
 imap_processing/hi/packet_definitions/TLM_HI_COMBINED_SCI.xml,sha256=Cwu_sE-b6lpXtuikQKuU5uQ_pMFTG9HoZeH9c-QDQgU,251743
 imap_processing/hi/packet_definitions/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/hi/utils.py,sha256=J_LpWqnhZnukhNk5EVyGaU6ZNZjcCosQ-RUSUTFjXFc,8098
+imap_processing/hi/utils.py,sha256=5_JowCx2jQsI58M0vCHOgtib_a6Jll9U2vOylpB26A4,8166
 imap_processing/hit/__init__.py,sha256=Ga6Ioj68HnQ4moi7SxJdOyYmGHcRvKOh4uZUEPJZHl4,20
 imap_processing/hit/hit_utils.py,sha256=-O8xt90w6sMPQ2t4CYmqrwkBw3rF1h1dMn58Q3Ioz3g,13371
 imap_processing/hit/l0/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/hit/l0/constants.py,sha256=MlLaYv73tTxVi4MNZL8oP-Ddnj92Xq_f1Ep0Jo2hkHU,5605
-imap_processing/hit/l0/decom_hit.py,sha256=ur0M3AVhuptWRvUVUvugxfAMf6yKI0Xee_7RNqJyF4E,15391
+imap_processing/hit/l0/constants.py,sha256=hJnOMvRVt3gm9mZY8iTVc1mBYq_MNuTz9twkuF4G1F0,5656
+imap_processing/hit/l0/decom_hit.py,sha256=o8rPVfGm67v1M4GUIjtSF7FpaQnYZbVRHJR79olDZOs,16436
 imap_processing/hit/l1a/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/hit/l1a/hit_l1a.py,sha256=BWxm73t6CeNXBt23NArnkXufeRH_OkhVOReLXEpiD3M,11241
+imap_processing/hit/l1a/hit_l1a.py,sha256=gy1wGEJZOxb-sZnlQACv7FVD966fek1EBvE8jzpZdWg,11345
 imap_processing/hit/l1b/constants.py,sha256=ZL5h2GgC5j_s11YgMxtY4EHpmkfo6IXnniCv-iXdG84,9692
-imap_processing/hit/l1b/hit_l1b.py,sha256=7CwzKpWD8HgsdQN0T72Dj9tP-jOE5_4kMViJ5yzFpYg,17495
+imap_processing/hit/l1b/hit_l1b.py,sha256=nXy8SUXJujFcWOrOBufmdrxq9xzdmB9dyZQOYz3-xiM,18489
 imap_processing/hit/l2/constants.py,sha256=jvs7Uic3Bl1wyRC01s1-F5LAsTzTwi089TpN19fvlFg,18813
-imap_processing/hit/l2/hit_l2.py,sha256=E_3hyn_XzmnUbAFqTSHeM_QQtCYqMnOzLm6Km_crmOk,27102
+imap_processing/hit/l2/hit_l2.py,sha256=qVGTHGO6tXNmz0nqKg3MEEqjlVGHhVsff5uS4jNajsk,26998
 imap_processing/hit/packet_definitions/hit_packet_definitions.xml,sha256=CyYmMnuObmk2y1nlrd7xOZnwZ_qZ3X__AhhKmRyNOvg,136147
 imap_processing/ialirt/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+imap_processing/ialirt/constants.py,sha256=mwVPgSVNuPpTO4zx3LeUW7L5a6aJG2ZrezObPvR9s1c,1241
 imap_processing/ialirt/l0/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 imap_processing/ialirt/l0/ialirt_spice.py,sha256=WNuMlxA5UDGiBd8ap2k_dHxFmi2UCFgnQ4_wEZUhVU0,5649
 imap_processing/ialirt/l0/mag_l0_ialirt_data.py,sha256=vgIerXkk4ZoOxxVaNEgvM1ESWVkGusSZ-3k73-Cl_tI,5276
-imap_processing/ialirt/l0/parse_mag.py,sha256=x51q2w_GPToSfKSnYzsscJkvsC2B8VSdL0jUuSo3FmQ,12780
-imap_processing/ialirt/l0/process_codicehi.py,sha256=cYkm_mKa502g4IeetTH-fVwxjFEZ2tRYFJ9WjrTEp40,4778
-imap_processing/ialirt/l0/process_codicelo.py,sha256=pF-1vBmIv2JhKzblSHgJMTXK23psJqCdFgWpL-ppl_Y,1037
-imap_processing/ialirt/l0/process_hit.py,sha256=Hws13pxQWRh7hYV2kl4bvg9R0sAtUo3DfnGAuaN5xSk,5680
-imap_processing/ialirt/l0/process_swapi.py,sha256=CLS-NVTddI1dGXEXdjSor2RNGnmdxMTe1_-WlU1PmyI,2113
-imap_processing/ialirt/l0/process_swe.py,sha256=5GGQ_f6uQlaEZB4vxIxvWCUUdN28coZHStxgqleyQWg,18150
+imap_processing/ialirt/l0/parse_mag.py,sha256=v5Efczml4s-TFdREhjg3A6PzEFTYAx4dbJyeACmkbPc,12787
+imap_processing/ialirt/l0/process_codice.py,sha256=guLrLsk1_lk7LJY_kJ4lMQocpTWIwrx50RcYnXL6Cig,3041
+imap_processing/ialirt/l0/process_hit.py,sha256=HJ_h8nd5Xs0NC-mWNS7Ce3D86H1408lT6ZCdXlNQv48,5408
+imap_processing/ialirt/l0/process_swapi.py,sha256=spcxacRDdHBpga1ZoAHN5EMUz-pDHdfazuYOrvzdywE,7493
+imap_processing/ialirt/l0/process_swe.py,sha256=QXp6rmlxw8K5M0bapKwAArCd_99qHFS8sxOnhHLVkhU,17953
 imap_processing/ialirt/packet_definitions/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 imap_processing/ialirt/packet_definitions/ialirt.xml,sha256=-hlT16Mr2Q6uCUfJ3jdpw0mPdeGCcZiWiorbkiEanwA,40784
 imap_processing/ialirt/packet_definitions/ialirt_codicehi.xml,sha256=PwzGTJj4dROXJBqbT32UZ6swaPbjvXdXXn7E_Lxch2A,11367
@@ -140,8 +134,8 @@ imap_processing/ialirt/packet_definitions/ialirt_mag.xml,sha256=Kt3t3S1gk7kqv006
 imap_processing/ialirt/packet_definitions/ialirt_swapi.xml,sha256=tGC4V-IV7WO6U3HtC2rDWA-500oOubcXgidMNyeGOE4,8400
 imap_processing/ialirt/packet_definitions/ialirt_swe.xml,sha256=FJNhfWjBjK1Ze3BCW86hvNJIiWoIK5dujt2Dw76pjF0,13442
 imap_processing/ialirt/process_ephemeris.py,sha256=RN6rzEY4KBc0V-ZKrsKffu8MEo7k0q3zJATBqaat9OY,8073
-imap_processing/ialirt/utils/constants.py,sha256=tdsY51p0qTOtTVrEdqPbHgBbXJMYQ80dpnVIK6MddQw,1852
-imap_processing/ialirt/utils/create_xarray.py,sha256=e4yLK725-BAkBjVAZWQPW9TLUNB_ZiqoByH4fezzaXg,2938
+imap_processing/ialirt/utils/constants.py,sha256=rF9KlnVig5Xk82lwo4OGbFPgf60lYNKHRP3FrlDQ4ls,2324
+imap_processing/ialirt/utils/create_xarray.py,sha256=mNMiFk06YTl7dUxqguNtX_qytISWfp9VqVaVm6vF-NI,4743
 imap_processing/ialirt/utils/grouping.py,sha256=Hk5u_TC0r72fN6s2CUEsD6opGig31EiwnRepAH9oOEE,3795
 imap_processing/ialirt/utils/time.py,sha256=fXBDKxNp2IXYJ1laSPPAIw_GTKvjFIgPwh54fnNQMEI,636
 imap_processing/idex/__init__.py,sha256=jO9AUm6ORgQhIHszOzLb3Vb2RVyz_SDUV64peGcr20A,52
@@ -262,10 +256,10 @@ imap_processing/ultra/packet_definitions/README.md,sha256=Iam1bd6UVkRLayApXlHG_y
 imap_processing/ultra/packet_definitions/ULTRA_SCI_COMBINED.xml,sha256=q1Bmz5vJTdqjBcKPExBZ_LCismlQFHOFyTUI1dJsjiY,939289
 imap_processing/ultra/packet_definitions/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 imap_processing/ultra/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-imap_processing/ultra/utils/ultra_l1_utils.py,sha256=QbDKjBo2m78AZMCorDdPpCL-ayQa7Fkm2PEbwmThCGU,5243
+imap_processing/ultra/utils/ultra_l1_utils.py,sha256=W1wABph96AJVkFNGrmq_0Vf6i__fIFxIVTvjCPbysY0,5321
 imap_processing/utils.py,sha256=zdcBXBQKL2NnElJcEbd-2QPeEDz0H3Yy0hVjRy9xcCE,14062
-imap_processing-0.15.0.dist-info/LICENSE,sha256=F2rxhvc6auEI0Dk9IGjglQSQQk60EvTe8M1dORMZPOg,1098
-imap_processing-0.15.0.dist-info/METADATA,sha256=o3cpLf1n619mCGzVt3RgsHIBuLiQo7z4eIjStubAg5U,9102
-imap_processing-0.15.0.dist-info/WHEEL,sha256=Nq82e9rUAnEjt98J6MlVmMCZb-t9cYE2Ir1kpBmnWfs,88
-imap_processing-0.15.0.dist-info/entry_points.txt,sha256=5r8ijLImHSNJxr-SGDC8kJy81BtXjmeUOmNfWSfLuRs,104
-imap_processing-0.15.0.dist-info/RECORD,,
+imap_processing-0.16.1.dist-info/LICENSE,sha256=F2rxhvc6auEI0Dk9IGjglQSQQk60EvTe8M1dORMZPOg,1098
+imap_processing-0.16.1.dist-info/METADATA,sha256=zjyVEsFl8ALuxnGSQ8rsqCWTvwsq3zUsPaDzTCSOmKQ,9102
+imap_processing-0.16.1.dist-info/WHEEL,sha256=Nq82e9rUAnEjt98J6MlVmMCZb-t9cYE2Ir1kpBmnWfs,88
+imap_processing-0.16.1.dist-info/entry_points.txt,sha256=5r8ijLImHSNJxr-SGDC8kJy81BtXjmeUOmNfWSfLuRs,104
+imap_processing-0.16.1.dist-info/RECORD,,