imap-processing 0.9.0__py3-none-any.whl → 0.11.0__py3-none-any.whl

imap_processing/cdf/config/imap_swe_l2_variable_attrs.yaml

@@ -0,0 +1,238 @@
+# <=== Coordinates ===>
+esa_step:
+  CATDESC: Energy step id in lookup table
+  FIELDNAM: energy step
+  FILLVAL: -9223370000000000000
+  FORMAT: I2
+  LABLAXIS: Energy Step
+  SCALE_TYP: linear
+  UNITS: " "
+  VALIDMAX: 23
+  VALIDMIN: 0
+  VAR_TYPE: support_data
+
+energy:
+  CATDESC: Energy of the particle in eV
+  FIELDNAM: Energy (eV)
+  FILLVAL: -9223370000000000000
+  FORMAT: E5.4
+  LABLAXIS: Energy
+  SCALE_TYP: linear
+  UNITS: eV
+  VALIDMAX: 2000.0
+  VALIDMIN: 0.0
+  VAR_TYPE: support_data
+
+spin_sector:
+  CATDESC: Spin sector. 30 measurements nominally
+  FIELDNAM: Spin Sector
+  FILLVAL: -9223370000000000000
+  FORMAT: I2
+  LABLAXIS: Spin Sector
+  SCALE_TYP: linear
+  UNITS: " "
+  VALIDMAX: 29
+  VALIDMIN: 0
+  VAR_TYPE: support_data
+
+inst_az:
+  CATDESC: Spin angle bins in instrument frame. Angle resolution is 12 degree nominally
+  FIELDNAM: spin angle
+  FILLVAL: -9223370000000000000
+  FORMAT: E5.4
+  LABLAXIS: Angle
+  SCALE_TYP: linear
+  UNITS: "Degree"
+  VALIDMAX: 359.999
+  VALIDMIN: 0.00
+  VAR_TYPE: support_data
+  VAR_NOTES: >
+    Bin CENTERS  - Assuming 0.5 sec measurements and a spin starting
+    at Angle=0, these are the angles at the middle of each measurement
+    bin.
+
+cem_id:
+  CATDESC: Data of each CEM detector
+  FIELDNAM: CEM data
+  FILLVAL: -9223370000000000000
+  FORMAT: E14.7
+  LABLAXIS: Rates
+  SCALE_TYP: linear
+  UNITS: " "
+  VALIDMAX: 0.000015514
+  VALIDMIN: 0.0
+  VAR_NOTES: " "
+  VAR_TYPE: support_data
+
+inst_el:
+  CATDESC: Angle of each CEM detector
+  FIELDNAM: CEM Angle
+  FILLVAL: -9223370000000000000
+  FORMAT: I2
+  LABLAXIS: Angle
+  SCALE_TYP: linear
+  UNITS: Degree
+  VALIDMAX: 63
+  VALIDMIN: -63
+  VAR_TYPE: support_data
+
+cycle:
+  CATDESC: Full cycle data takes 4 spins' data
+  FIELDNAM: quarter cycle
+  FILLVAL: -9223370000000000000
+  FORMAT: I2
+  LABLAXIS: Cycle
+  SCALE_TYP: linear
+  UNITS: " "
+  VALIDMAX: 3
+  VALIDMIN: 0
+  VAR_TYPE: support_data
+
+# <=== Label Attributes ===>
+# LABL_PTR_i expects VAR_TYPE of metadata with char data type
+esa_step_label:
+  CATDESC: Energy step id in lookup table
+  FIELDNAM: Energy Step ID
+  FORMAT: A2
+  VAR_TYPE: metadata
+
+spin_sector_label:
+  CATDESC: Spin sector. 30 measurements nominally
+  FIELDNAM: Measurement
+  FORMAT: A2
+  VAR_TYPE: metadata
+
+cem_id_label:
+  CATDESC: Data rates of each CEM detector
+  FIELDNAM: CEM Rates
+  FORMAT: A1
+  VAR_TYPE: metadata
+
+energy_label:
+  CATDESC: Energy of the particle in eV
+  FIELDNAM: Energy
+  FORMAT: A5
+  VAR_TYPE: metadata
+
+inst_az_label:
+  CATDESC: Spin angle bins in instrument frame. Angle resolution is 12 degree nominally
+  FIELDNAM: Spin Angle
+  FORMAT: A3
+  VAR_TYPE: metadata
+
+inst_el_label:
+  CATDESC: Angle of each CEM detector
+  FIELDNAM: CEM Angle
+  FORMAT: A3
+  VAR_TYPE: metadata
+
+# <=== Data Variables ===>
+
+inst_az_spin_sector:
+  CATDESC: Spin angle organized by voltage step and spin sector
+  DEPEND_0: epoch
+  DEPEND_1: esa_step
+  DEPEND_2: spin_sector
+  LABL_PTR_1: esa_step_label
+  LABL_PTR_2: spin_sector_label
+  DISPLAY_TYPE: spectrogram
+  FIELDNAM: Spin Angle
+  FILLVAL: -1.0000000E+31
+  FORMAT: E4.5
+  UNITS: Degree
+  VALIDMAX: 359.99
+  VALIDMIN: 0.0
+  VAR_TYPE: data
+
+phase_space_density_spin_sector:
+  CATDESC: Phase space density organized by voltage step and spin sector and CEM
+  DEPEND_0: epoch
+  DEPEND_1: esa_step
+  DEPEND_2: spin_sector
+  DEPEND_3: cem_id
+  LABL_PTR_1: esa_step_label
+  LABL_PTR_2: spin_sector_label
+  LABL_PTR_3: cem_id_label
+  DISPLAY_TYPE: spectrogram
+  FIELDNAM: Phase Space Density
+  FILLVAL: -1.0000000E+31
+  UNITS: s^3 / (cm^6 * ster)
+  FORMAT: E14.7
+  VALIDMAX: 0.000015514
+  VALIDMIN: 0
+  VAR_TYPE: data
+
+phase_space_density:
+  CATDESC: Phase space density organized by energy (eV), spin angle, CEM angle
+  DEPEND_0: epoch
+  DEPEND_1: energy
+  DEPEND_2: inst_az
+  DEPEND_3: inst_el
+  LABL_PTR_1: energy_label
+  LABL_PTR_2: inst_az_label
+  LABL_PTR_3: inst_el_label
+  DISPLAY_TYPE: spectrogram
+  FIELDNAM: Phase Space Density
+  FILLVAL: -1.0000000E+31
+  UNITS: s^3 / (cm^6 * ster)
+  FORMAT: E14.7
+  VALIDMAX: 0.000015514
+  VALIDMIN: 0
+  VAR_TYPE: data
+
+flux_spin_sector:
+  CATDESC: Flux organized by voltage step and spin sector and CEM
+  DEPEND_0: epoch
+  DEPEND_1: esa_step
+  DEPEND_2: spin_sector
+  DEPEND_3: cem_id
+  LABL_PTR_1: esa_step_label
+  LABL_PTR_2: spin_sector_label
+  LABL_PTR_3: cem_id_label
+  DISPLAY_TYPE: spectrogram
+  FIELDNAM: Flux
+  FILLVAL: -1.0000000E+31
+  UNITS: 1 / (2 * eV * cm^2 * s * ster)
+  FORMAT: E14.7
+  VALIDMAX: 0.000015514
+  VALIDMIN: 0
+  VAR_TYPE: data
+
+flux:
+  CATDESC: Flux organized by energy (eV), spin angle, CEM angle
+  DEPEND_0: epoch
+  DEPEND_1: energy
+  DEPEND_2: inst_az
+  DEPEND_3: inst_el
+  LABL_PTR_1: energy_label
+  LABL_PTR_2: inst_az_label
+  LABL_PTR_3: inst_el_label
+  DISPLAY_TYPE: spectrogram
+  FIELDNAM: Flux
+  FILLVAL: -1.0000000E+31
+  UNITS: 1 / (2 * eV * cm^2 * s * ster)
+  FORMAT: E14.7
+  VALIDMAX: 0.000015514
+  VALIDMIN: 0
+  VAR_TYPE: data
+
+acquisition_time:
+  CATDESC: Acquisition time organized by voltage step and spin sector
+  DEPEND_0: epoch
+  DEPEND_1: esa_step
+  DEPEND_2: spin_sector
+  LABL_PTR_1: esa_step_label
+  LABL_PTR_2: spin_sector_label
+  DISPLAY_TYPE: spectrogram
+  FIELDNAM: Acquisition time by volt step and spin sector
+  FILLVAL: -1.0000000E+31
+  FORMAT: I20
+  VALIDMAX: 9223372036854775807
+  VALIDMIN: 0
+  UNITS: sec
+  VAR_TYPE: support_data
+  VAR_NOTES: >
+    This stores the acquisition times of each measurement. This time is
+    calculated by combining ACQ_START_COARSE, ACQ_START_FINE,
+    acquisition duration and settle duration. It is time of each energy
+    step and each science measurement.

imap_processing/cdf/config/imap_ultra_l1b_variable_attrs.yaml

@@ -9,6 +9,13 @@ default_attrs: &default
   VAR_TYPE: data
   UNITS: " "
 
+component:
+  CATDESC: Velocity components (vx, vy, vz)
+  FIELDNAM: component
+  LABLAXIS: component
+  FORMAT: A3
+  VAR_TYPE: metadata
+
 default_uint8_attrs: &default_uint8
   <<: *default
   FILLVAL: 255
@@ -92,6 +99,13 @@ x_coin:
   LABLAXIS: x coincidence position
   UNITS: mm / 100
 
+event_times:
+  <<: *default_float32
+  CATDESC: time of event as calculated using aux and de l1a packets
+  FIELDNAM: event_times
+  LABLAXIS: event times
+  UNITS: seconds
+
 tof_start_stop:
   <<: *default_float32
   CATDESC: Particle time of flight from start to stop
@@ -134,26 +148,20 @@ start_type:
   LABLAXIS: start type
   UNITS: " "
 
-vx_ultra:
+direct_event_velocity:
   <<: *default_float32
-  CATDESC: Normalized component of the velocity vector in x direction in the instrument frame.
-  FIELDNAM: vx_ultra
-  LABLAXIS: vx ultra
-  UNITS: " "
+  CATDESC: Direct event velocity vector in x, y, and z directions in the instrument frame.
+  DEPEND_1: component
+  FIELDNAM: v_ultra
+  LABLAXIS: v ultra
+  UNITS: "km/s"
 
-vy_ultra:
+velocity_magnitude:
   <<: *default_float32
-  CATDESC: Normalized component of the velocity vector in y direction in the instrument frame.
-  FIELDNAM: vy_ultra
-  LABLAXIS: vy ultra
-  UNITS: " "
-
-vz_ultra:
-  <<: *default_float32
-  CATDESC: Normalized component of the velocity vector in z direction in the instrument frame.
-  FIELDNAM: vz_ultra
-  LABLAXIS: vz ultra
-  UNITS: " "
+  CATDESC: Magnitude of the particle velocity.
+  FIELDNAM: vmag_ultra
+  LABLAXIS: vmag ultra
+  UNITS: "km/s"
 
 energy:
   <<: *default_float32
@@ -162,6 +170,13 @@ energy:
   LABLAXIS: energy
   UNITS: keV
 
+tof_energy:
+  <<: *default_float32
+  CATDESC: Energy measured using kinetic energy equation.
+  FIELDNAM: energy
+  LABLAXIS: energy
+  UNITS: keV
+
 species:
   <<: *default
   DISPLAY_TYPE: no_plot
@@ -186,74 +201,49 @@ path_length:
   LABLAXIS: path_length
   UNITS: mm / 100
 
-event_efficiency:
-  <<: *default_float32
-  CATDESC: Estimated event efficiency for this path through the instrument.
-  FIELDNAM: event_efficiency
-  LABLAXIS: event efficiency
-  UNITS: " "
-
-vx_sc:
-  <<: *default_float32
-  CATDESC: x-component of the velocity vector of the ena in the spacecraft frame.
-  FIELDNAM: vx_sc
-  LABLAXIS: vx sc
-  UNITS: " "
-
-vy_sc:
-  <<: *default_float32
-  CATDESC: y-component of the velocity vector of the ena in the spacecraft frame.
-  FIELDNAM: vy_sc
-  LABLAXIS: vy sc
-  UNITS: " "
-
-vz_sc:
+azimuth:
   <<: *default_float32
-  CATDESC: z-component of the velocity vector of the ena in the spacecraft frame.
-  FIELDNAM: vz_sc
-  LABLAXIS: vz sc
-  UNITS: " "
+  CATDESC: Azimuth.
+  FIELDNAM: azimuth
+  LABLAXIS: azimuth
+  UNITS: radians
 
-vx_dps_sc:
+elevation:
   <<: *default_float32
-  CATDESC: x-component of the velocity vector of the ena in the DPS frame at rest WRT spacecraft.
-  FIELDNAM: vx_dps_sc
-  LABLAXIS: vx dps sc
-  UNITS: " "
+  CATDESC: Elevation.
+  FIELDNAM: elevation
+  LABLAXIS: elevation
+  UNITS: radians
 
-vy_dps_sc:
-  <<: *default_float32
-  CATDESC: y-component of the velocity vector of the ena in the DPS frame at rest WRT spacecraft.
-  FIELDNAM: vy_dps_sc
-  LABLAXIS: vy dps sc
-  UNITS: " "
-
-vz_dps_sc:
+event_efficiency:
   <<: *default_float32
-  CATDESC: z-component of the velocity vector of the ena in the DPS frame at rest WRT spacecraft.
-  FIELDNAM: vz_dps_sc
-  LABLAXIS: vz dps sc
+  CATDESC: Estimated event efficiency for this path through the instrument.
+  FIELDNAM: event_efficiency
+  LABLAXIS: event efficiency
   UNITS: " "
 
-vx_dps_helio:
+velocity_sc:
   <<: *default_float32
-  CATDESC: x-component of the velocity vector of the ena in the DPS frame at rest WRT heliosphere.
-  FIELDNAM: vx_dps_helio
-  LABLAXIS: vx dps helio
+  CATDESC: x,y,z-components of the velocity vector of the ena in the spacecraft frame.
+  DEPEND_1: component
+  FIELDNAM: v_sc
+  LABLAXIS: v sc
   UNITS: " "
 
-vy_dps_helio:
+velocity_dps_sc:
   <<: *default_float32
-  CATDESC: y-component of the velocity vector of the ena in the DPS frame at rest WRT heliosphere.
-  FIELDNAM: vy_dps_helio
-  LABLAXIS: vy dps helio
+  CATDESC: x,y,z-components of the velocity vector of the ena in the DPS frame at rest WRT spacecraft.
+  DEPEND_1: component
+  FIELDNAM: v_dps_sc
+  LABLAXIS: v dps sc
   UNITS: " "
 
-vz_dps_helio:
+velocity_dps_helio:
   <<: *default_float32
-  CATDESC: z-component of the velocity vector of the ena in the DPS frame at rest WRT heliosphere.
-  FIELDNAM: vz_dps_helio
-  LABLAXIS: vz dps helio
+  CATDESC: x,y,z-components of the velocity vector of the ena in the DPS frame at rest WRT heliosphere.
+  DEPEND_1: component
+  FIELDNAM: v_dps_helio
+  LABLAXIS: v dps helio
   UNITS: " "
 
 de_event_met:
@@ -266,14 +256,6 @@ de_event_met:
   VAR_TYPE: support_data
   SCALE_TYP: linear
 
-spin_number:
-  <<: *default
-  CATDESC: Spin number from Universal Spin Table.
-  FIELDNAM: spin_number
-  LABLAXIS: spin_number
-  # TODO: come back to format
-  UNITS: " "
-
 spin_start_time:
   <<: *default
   CATDESC: Spin start time from Universal Spin Table.
@@ -282,14 +264,20 @@ spin_start_time:
   # TODO: come back to format
   UNITS: s
 
-avg_spin_period:
+spin_period:
   <<: *default
-  CATDESC: Average spin period from Universal Spin Table.
-  FIELDNAM: avg_spin_period
-  LABLAXIS: avg_spin_period
-  # TODO: come back to format
+  CATDESC: Spin period from Universal Spin Table.
+  FIELDNAM: spin_period
+  LABLAXIS: spin_period
   UNITS: s
 
+spin_rate:
+  <<: *default
+  CATDESC: Spin rate from Universal Spin Table.
+  FIELDNAM: spin_rate
+  LABLAXIS: spin_rate
+  UNITS: rpm
+
 rate_start_pulses:
   <<: *default
   CATDESC: Rate of start pulses (/s).
@@ -330,16 +318,36 @@ rate_rejected_events:
   # TODO: come back to format
   UNITS: 1/s
 
-quality_hk:
-  <<: *default
-  CATDESC: Spin filter derived from Ultra housekeeping. Bitwise flagging used to filter the spin.
-  FIELDNAM: quality_hk
-  LABLAXIS: quality hk
-  # TODO: come back to format
+# Extendedspin attributes
+ena_rates:
+  <<: *default_float32
+  CATDESC: Rates calculated from de packet.
+  FIELDNAM: ena_rates
+  LABLAXIS: ena_rates
+  DEPEND_0: spin_number
+  DEPEND_1: energy_bin_geometric_mean
+  UNITS: " "
+
+ena_rates_threshold:
+  <<: *default_float32
+  CATDESC: Rates threshold used for flagging data.
+  FIELDNAM: ena_rates_threshold
+  LABLAXIS: ena_rates_threshold
+  DEPEND_0: spin_number
+  DEPEND_1: energy_bin_geometric_mean
+  UNITS: " "
+
+quality_ena_rates:
+  <<: *default_uint16
+  CATDESC: Spin filter derived from Ultra ena rates. Bitwise flagging used to filter the spin.
+  FIELDNAM: quality_ena_rates
+  LABLAXIS: quality_ena_rates
+  DEPEND_0: spin_number
+  DEPEND_1: energy_bin_geometric_mean
   UNITS: " "
 
 quality_attitude:
-  <<: *default
+  <<: *default_uint16
   CATDESC: Spin filter derived from IMAP attitude. Bitwise flagging used to filter the spin.
   FIELDNAM: quality_attitude
   LABLAXIS: quality attitude
@@ -347,7 +355,7 @@ quality_attitude:
   UNITS: " "
 
 quality_instruments:
-  <<: *default
+  <<: *default_uint16
   CATDESC: Spin filter derived from instruments other than Ultra. Bitwise flagging used to filter the spin.
   FIELDNAM: quality_instruments
   LABLAXIS: quality instruments

imap_processing/cdf/utils.py

@@ -13,7 +13,7 @@ from cdflib.xarray.cdf_to_xarray import ISTP_TO_XARRAY_ATTRS
 
 import imap_processing
 from imap_processing._version import __version__, __version_tuple__  # noqa: F401
-from imap_processing.spice.time import J2000_EPOCH
+from imap_processing.spice.time import TTJ2000_EPOCH
 
 logger = logging.getLogger(__name__)
 
@@ -94,7 +94,7 @@ def write_cdf(
     # Convert J2000 epoch referenced data to datetime64
     # TODO: This implementation of epoch to time string results in an error of
     #       5 seconds due to 5 leap-second occurrences since the J2000 epoch.
-    dt64 = J2000_EPOCH + dataset["epoch"].values[0].astype("timedelta64[ns]")
+    dt64 = TTJ2000_EPOCH + dataset["epoch"].values[0].astype("timedelta64[ns]")
     start_time = np.datetime_as_string(dt64, unit="D").replace("-", "")
 
     # Will now accept vXXX or XXX formats, as batch starter sends versions as vXXX.

imap_processing/cli.py

@@ -38,6 +38,7 @@ from imap_processing.cdf.utils import load_cdf, write_cdf
 from imap_processing.codice import codice_l1a, codice_l1b
 from imap_processing.glows.l1a.glows_l1a import glows_l1a
 from imap_processing.glows.l1b.glows_l1b import glows_l1b
+from imap_processing.glows.l2.glows_l2 import glows_l2
 from imap_processing.hi.l1a import hi_l1a
 from imap_processing.hi.l1b import hi_l1b
 from imap_processing.hi.l1c import hi_l1c
@@ -226,6 +227,8 @@ class ProcessInstrument(ABC):
     ----------
     data_level : str
         The data level to process (e.g. ``l1a``).
+    data_descriptor : str
+        The descriptor of the data to process (e.g. ``sci``).
     dependency_str : str
         A string representation of the dependencies for the instrument in the
         format: "[{
@@ -248,6 +251,7 @@ class ProcessInstrument(ABC):
     def __init__(
         self,
         data_level: str,
+        data_descriptor: str,
         dependency_str: str,
         start_date: str,
         end_date: str,
@@ -255,6 +259,7 @@ class ProcessInstrument(ABC):
         upload_to_sdc: bool,
     ) -> None:
         self.data_level = data_level
+        self.descriptor = data_descriptor
 
         # Convert string into a dictionary
         self.dependencies = loads(dependency_str.replace("'", '"'))
@@ -283,8 +288,8 @@ class ProcessInstrument(ABC):
                 # TODO: Validate dep dict
                 # TODO: determine what dependency information is optional
                 return_query = imap_data_access.query(
-                    start_date=self.start_date,
-                    end_date=self.end_date,
+                    start_date=dependency["start_date"],
+                    end_date=dependency.get("end_date", None),
                     instrument=dependency["instrument"],
                     data_level=dependency["data_level"],
                     version=dependency["version"],
@@ -395,7 +400,13 @@ class ProcessInstrument(ABC):
         datasets : list[xarray.Dataset]
             A list of datasets (products) produced by do_processing method.
         """
+        if len(datasets) == 0:
+            logger.info("No products to write to CDF file.")
+            return
+
         logger.info("Writing products to local storage")
+        logger.info("Parent files: %s", self._dependency_list)
+
         products = [
             write_cdf(dataset, parent_files=self._dependency_list)
             for dataset in datasets
@@ -482,6 +493,15 @@ class Glows(ProcessInstrument):
             input_dataset = load_cdf(dependencies[0])
             datasets = [glows_l1b(input_dataset, self.version)]
 
+        if self.data_level == "l2":
+            if len(dependencies) > 1:
+                raise ValueError(
+                    f"Unexpected dependencies found for GLOWS L2:"
+                    f"{dependencies}. Expected only one input dependency."
+                )
+            input_dataset = load_cdf(dependencies[0])
+            datasets = glows_l2(input_dataset, self.version)
+
         return datasets
 
 
@@ -517,7 +537,14 @@ class Hi(ProcessInstrument):
             dependencies = [load_cdf(dependency) for dependency in dependencies]
             datasets = [hi_l1b.hi_l1b(dependencies[0], self.version)]
         elif self.data_level == "l1c":
-            dependencies = [load_cdf(dependency) for dependency in dependencies]
+            # TODO: Add PSET calibration product config file dependency and remove
+            #    below injected dependency
+            dependencies.append(
+                Path(__file__).parent
+                / "tests/hi/test_data/l1"
+                / "imap_his_pset-calibration-prod-config_20240101_v001.csv"
+            )
+            dependencies[0] = load_cdf(dependencies[0])
             datasets = [hi_l1c.hi_l1c(dependencies, self.version)]
         else:
             raise NotImplementedError(
@@ -734,13 +761,20 @@ class Swapi(ProcessInstrument):
         datasets: list[xr.Dataset] = []
 
         if self.data_level == "l1":
-            if len(dependencies) > 1:
+            # For science, we expect l0 raw file and L1 housekeeping file
+            if self.descriptor == "sci" and len(dependencies) != 2:
                 raise ValueError(
-                    f"Unexpected dependencies found for SWAPI L1:"
+                    f"Unexpected dependencies found for SWAPI L1 science:"
+                    f"{dependencies}. Expected only two dependencies."
+                )
+            # For housekeeping, we expect only L0 raw file
+            if self.descriptor == "hk" and len(dependencies) != 1:
+                raise ValueError(
+                    f"Unexpected dependencies found for SWAPI L1 housekeeping:"
                     f"{dependencies}. Expected only one dependency."
                 )
-            # process data
-            datasets = [swapi_l1(dependencies[0], self.version)]
+            # process science or housekeeping data
+            datasets = swapi_l1(dependencies, self.version)
         elif self.data_level == "l2":
             if len(dependencies) > 1:
                 raise ValueError(
@@ -780,7 +814,7 @@ class Swe(ProcessInstrument):
                     f"Unexpected dependencies found for SWE L1A:"
                     f"{dependencies}. Expected only one dependency."
                 )
-            datasets = [swe_l1a(str(dependencies[0]), data_version=self.version)]
+            datasets = swe_l1a(str(dependencies[0]), data_version=self.version)
             # Right now, we only process science data. Therefore,
             # we expect only one dataset to be returned.
 
@@ -792,7 +826,8 @@ class Swe(ProcessInstrument):
                 )
             # read CDF file
             l1a_dataset = load_cdf(dependencies[0])
-            datasets = [swe_l1b(l1a_dataset, data_version=self.version)]
+            # TODO: read lookup table and in-flight calibration data here.
+            datasets = swe_l1b(l1a_dataset, data_version=self.version)
         else:
             print("Did not recognize data level. No processing done.")
 
@@ -861,6 +896,7 @@ def main() -> None:
     cls = getattr(sys.modules[__name__], args.instrument.capitalize())
     instrument = cls(
         args.data_level,
+        args.descriptor,
         args.dependency,
         args.start_date,
         args.end_date,