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

imap_processing/codice/codice_l1a.py

@@ -27,6 +27,7 @@ from imap_processing.codice import constants
 from imap_processing.codice.codice_l0 import decom_packets
 from imap_processing.codice.decompress import decompress
 from imap_processing.codice.utils import CODICEAPID
+from imap_processing.spice.time import met_to_ttj2000ns
 
 logger = logging.getLogger(__name__)
 logger.setLevel(logging.INFO)
@@ -57,6 +58,8 @@ class CoDICEL1aPipeline:
 
     Methods
     -------
+    calculate_epoch_values()
+        Calculate and return the values to be used for `epoch`.
     decompress_data(science_values)
         Perform decompression on the data.
     define_coordinates()
@@ -84,6 +87,26 @@ class CoDICEL1aPipeline:
         self.plan_step = plan_step
         self.view_id = view_id
 
+    def calculate_epoch_values(self) -> NDArray[int]:
+        """
+        Calculate and return the values to be used for `epoch`.
+
+        On CoDICE, the epoch values are derived from the `acq_start_seconds` and
+        `acq_start_subseconds` fields in the packet. Note that the
+        `acq_start_subseconds` field needs to be converted from microseconds to
+        seconds
+
+        Returns
+        -------
+        epoch : NDArray[int]
+            List of epoch values.
+        """
+        epoch = met_to_ttj2000ns(
+            self.dataset.acq_start_seconds + self.dataset.acq_start_subseconds / 1e6
+        )
+
+        return epoch
+
     def decompress_data(self, science_values: list[str]) -> None:
         """
         Perform decompression on the data.
@@ -126,17 +149,34 @@ class CoDICEL1aPipeline:
         """
         self.coords = {}
 
-        for name in self.config["coords"]:
+        coord_names = ["epoch", *list(self.config["dims"].keys())]
+
+        # These are labels unique to lo-counters products coordinates
+        if self.config["dataset_name"] in [
+            "imap_codice_l1a_lo-counters-aggregated",
+            "imap_codice_l1a_lo-counters-singles",
+        ]:
+            coord_names.append("spin_sector_pairs_label")
+
+        for name in coord_names:
             if name == "epoch":
-                values = self.dataset.epoch.data
-            elif name == "esa_step":
-                values = np.arange(self.config["num_energy_steps"])
-            elif name == "inst_az":
-                values = np.arange(self.config["num_positions"])
-            elif name == "spin_sector":
-                values = np.arange(self.config["num_spin_sectors"])
+                values = self.calculate_epoch_values()
+            elif name in ["esa_step", "inst_az", "spin_sector", "spin_sector_pairs"]:
+                values = np.arange(self.config["dims"][name])
+            elif name == "spin_sector_pairs_label":
+                values = np.array(
+                    [
+                        "0-30 deg",
+                        "30-60 deg",
+                        "60-90 deg",
+                        "90-120 deg",
+                        "120-150 deg",
+                        "150-180 deg",
+                    ]
+                )
             else:
-                # TODO: Need to implement other types of coords
+                # TODO: May need to implement other types of coords for Hi
+                #       and/or event data products
                 continue
 
             coord = xr.DataArray(
@@ -170,26 +210,31 @@ class CoDICEL1aPipeline:
         # Stack the data so that it is easier to reshape and iterate over
         all_data = np.stack(self.data)
 
-        # The dimension of all data is (epoch, num_counters, num_energy_steps,
-        # num_positions, num_spin_sectors) (or may be slightly different
-        # depending on the data product). In any case, iterate over the
-        # num_counters dimension to isolate the data for each counter so
-        # that it can be placed in a CDF data variable.
+        # The dimension of all_data is something like (epoch, num_counters,
+        # num_energy_steps, num_positions, num_spin_sectors) (or may be slightly
+        # different depending on the data product). In any case, iterate over
+        # the num_counters dimension to isolate the data for each counter so
+        # each counter's data can be placed in a separate CDF data variable.
         for counter, variable_name in zip(
             range(all_data.shape[1]), self.config["variable_names"]
         ):
-            counter_data = all_data[:, counter, :, :, :]
+            # Extract the counter data
+            counter_data = all_data[:, counter, ...]
 
             # Get the CDF attributes
             descriptor = self.config["dataset_name"].split("imap_codice_l1a_")[-1]
             cdf_attrs_key = f"{descriptor}-{variable_name}"
             attrs = self.cdf_attrs.get_variable_attributes(cdf_attrs_key)
 
+            # The final CDF dimensions always has "epoch" as the first dimension,
+            # followed by the dimensions for the specific data product
+            dims = ["epoch", *list(self.config["dims"].keys())]
+
             # Create the CDF data variable
             dataset[variable_name] = xr.DataArray(
                 counter_data,
                 name=variable_name,
-                dims=self.config["dims"],
+                dims=dims,
                 attrs=attrs,
             )
 
@@ -216,11 +261,20 @@ class CoDICEL1aPipeline:
         dataset : xarray.Dataset
             ``xarray`` dataset for the data product, with added support variables.
         """
+        # These variables can be gathered from the packet data
+        packet_data_variables = [
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
+        ]
+
         for variable_name in self.config["support_variables"]:
+            # These variables require reading in external tables
             if variable_name == "energy_table":
                 variable_data = self.get_energy_table()
                 dims = ["esa_step"]
-                attrs = self.cdf_attrs.get_variable_attributes("esa_step")
+                attrs = self.cdf_attrs.get_variable_attributes("energy_table")
 
             elif variable_name == "acquisition_time_per_step":
                 variable_data = self.get_acquisition_times()
@@ -229,10 +283,24 @@ class CoDICEL1aPipeline:
                     "acquisition_time_per_step"
                 )
 
-            else:
-                # TODO: Need to implement methods to gather and set other
-                #       support attributes
-                continue
+            elif variable_name in packet_data_variables:
+                variable_data = self.dataset[variable_name].data
+                dims = ["epoch"]
+                attrs = self.cdf_attrs.get_variable_attributes(variable_name)
+
+            # Data quality is named differently in packet data and needs to be
+            # treated slightly differently
+            elif variable_name == "data_quality":
+                variable_data = self.dataset.suspect.data
+                dims = ["epoch"]
+                attrs = self.cdf_attrs.get_variable_attributes("data_quality")
+
+            # Spin period requires the application of a conversion factor
+            # See Table B.5 in the algorithm document
+            elif variable_name == "spin_period":
+                variable_data = self.dataset.spin_period.data * 0.00032
+                dims = ["epoch"]
+                attrs = self.cdf_attrs.get_variable_attributes("spin_period")
 
             # Add variable to the dataset
             dataset[variable_name] = xr.DataArray(
@@ -341,45 +409,23 @@ class CoDICEL1aPipeline:
         3D arrays representing dimensions such as spin sectors, positions, and
         energies (depending on the data product).
         """
+        # This will contain the reshaped data for all counters
         self.data = []
 
-        # For CoDICE-lo, data are a 3D arrays with a shape representing
-        # [<num_positions>,<num_spin_sectors>,<num_energy_steps>]
-        if self.config["instrument"] == "lo":
-            for packet_data in self.raw_data:
-                if packet_data:
-                    reshaped_packet_data = np.array(
-                        packet_data, dtype=np.uint32
-                    ).reshape(
-                        (
-                            self.config["num_counters"],
-                            self.config["num_energy_steps"],
-                            self.config["num_positions"],
-                            self.config["num_spin_sectors"],
-                        )
-                    )
-                    self.data.append(reshaped_packet_data)
-                else:
-                    self.data.append(None)
-
-        # For CoDICE-hi, data are a 3D array with a shape representing
-        # [<num_energy_steps>,<num_positions>,<num_spin_sectors>]
-        elif self.config["instrument"] == "hi":
-            for packet_data in self.raw_data:
-                if packet_data:
-                    reshaped_packet_data = np.array(
-                        packet_data, dtype=np.uint32
-                    ).reshape(
-                        (
-                            self.config["num_counters"],
-                            self.config["num_energy_steps"],
-                            self.config["num_positions"],
-                            self.config["num_spin_sectors"],
-                        )
-                    )
-                    self.data.append(reshaped_packet_data)
-                else:
-                    self.data.append(None)
+        # Typically, data are a 4D arrays with a shape representing some
+        # combination of <num_counters>, <num_positions>, <num_spin_sectors>,
+        # and <num_energy_steps>. However, the existence and order of these
+        # dimensions can vary depending on the specific data product, so we
+        # define this in the "dims" value configuration dictionary. The number
+        # of counters is the first dimension/axis.
+        reshape_dims = (self.config["num_counters"], *self.config["dims"].values())
+
+        # For each packet/epoch, reshape the data along these dimensions
+        for packet_data in self.raw_data:
+            reshaped_packet_data = np.array(packet_data, dtype=np.uint32).reshape(
+                reshape_dims
+            )
+            self.data.append(reshaped_packet_data)
 
         # No longer need to keep the raw data around
         del self.raw_data

imap_processing/codice/constants.py

@@ -12,9 +12,6 @@ PUI = PickUp Ion
 ESA = ElectroStatic Analyzer
 """
 
-# TODO: What to do in the case of a value of 255 in LOSSY_A and LOSSY_B
-#       compression? (Joey uses 0x100000000)
-
 from imap_processing.codice.utils import CODICEAPID, CoDICECompression
 
 APIDS_FOR_SCIENCE_PROCESSING = [
@@ -34,8 +31,7 @@ APIDS_FOR_SCIENCE_PROCESSING = [
 
 
 # CDF-friendly names for lo data products
-LO_INST_COUNTS_AGGREGATED_VARIABLE_NAMES = ["aggregated"]
-LO_INST_COUNTS_SINGLES_VARIABLE_NAMES = ["apd_singles"]
+LO_COUNTERS_SINGLES_VARIABLE_NAMES = ["apd_singles"]
 LO_SW_ANGULAR_VARIABLE_NAMES = ["hplus", "heplusplus", "oplus6", "fe_loq"]
 LO_NSW_ANGULAR_VARIABLE_NAMES = ["heplusplus"]
 LO_SW_PRIORITY_VARIABLE_NAMES = [
@@ -81,270 +77,230 @@ HI_INST_COUNTS_SINGLES_VARIABLE_NAMES = ["tcr", "ssdo", "stssd"]
 HI_OMNI_SPECIES_VARIABLE_NAMES = ["h", "he3", "he4", "c", "o", "ne_mg_si", "fe", "uh"]
 HI_SECT_SPECIES_VARIABLE_NAMES = ["h", "he3he4", "cno", "fe"]
 
-# TODO: Hi products shape should be energy x ssd index x spin sector (8*12*12)
+# lo-counters-aggregated data product variables are dynamically determined
+# based on the number of active counters
+# TODO: Try to convince Joey to move to lower case variable names with
+#       underscores?
+LO_COUNTERS_AGGREGATED_ACTIVE_VARIABLES = {
+    "TCR": True,
+    "DCR": True,
+    "TOFPlusAPD": False,
+    "TOFOnly": False,
+    "PositionPlusAPD": False,
+    "PositionOnly": False,
+    "STAorSTBPlusAPD": False,
+    "STAorSTBOnly": False,
+    "Reserved1": False,
+    "Reserved2": False,
+    "SPOnly": False,
+    "APDOnly": False,
+    "LowTOFCutoff": False,
+    "STA": True,
+    "STB": True,
+    "SP": True,
+    "TotalPositionCount": True,
+    "InvalidPositionCount": False,
+    "ASIC1FlagInvalid": False,
+    "ASIC2FlagInvalid": False,
+    "ASIC1ChannelInvalid": False,
+    "ASIC2ChannelInvalid": False,
+    "TEC4TimeoutTOFNoPos": False,
+    "TEC4TimeoutPosNoTOF": False,
+    "TEC4TimeoutNoPosTOF": False,
+    "TEC5TimeoutTOFNoPos": False,
+    "TEC5TimeoutPosNoTOF": False,
+    "TEC5TimeoutNoPosTOF": False,
+}
+LO_COUNTERS_AGGREGATED_VARIABLE_NAMES = [
+    name
+    for name, is_active in LO_COUNTERS_AGGREGATED_ACTIVE_VARIABLES.items()
+    if is_active
+]
+
+# TODO: Possibly move to consistent order of dimensions with other instruments
+#       TBD after discussion with Joey and at the Science Team Meeting in Feb
 DATA_PRODUCT_CONFIGURATIONS = {
     CODICEAPID.COD_HI_INST_COUNTS_AGGREGATED: {
-        "coords": [
-            "epoch",
-            "inst_az",
-            "spin_sector",
-            "esa_step",
-            "energy_label",
-        ],  # TODO: These will likely change
         "dataset_name": "imap_codice_l1a_hi-counters-aggregated",
-        "dims": [
-            "epoch",
-            "esa_step",
-            "inst_az",
-            "spin_sector",
-        ],  # TODO: These will likely change
+        "dims": {
+            "esa_step": 1,
+            "inst_az": 6,
+            "spin_sector": 1,
+        },  # TODO: Double check with Joey
         "instrument": "hi",
         "num_counters": 1,
-        "num_energy_steps": 1,  # TODO: Double check with Joey
-        "num_positions": 6,  # TODO: Double check with Joey
-        "num_spin_sectors": 1,
         "support_variables": [],  # TODO: Double check with Joey
         "variable_names": HI_INST_COUNTS_AGGREGATED_VARIABLE_NAMES,
     },
     CODICEAPID.COD_HI_INST_COUNTS_SINGLES: {
-        "coords": [
-            "epoch",
-            "inst_az",
-            "spin_sector",
-            "esa_step",
-            "energy_label",
-        ],  # TODO: These will likely change
         "dataset_name": "imap_codice_l1a_hi-counters-singles",
-        "dims": [
-            "epoch",
-            "esa_step",
-            "inst_az",
-            "spin_sector",
-        ],  # TODO: These will likely change
+        "dims": {
+            "esa_step": 1,
+            "inst_az": 12,
+            "spin_sector": 1,
+        },  # TODO: Double check with Joey
         "instrument": "hi",
         "num_counters": 3,
-        "num_energy_steps": 1,  # TODO: Double check with Joey
-        "num_positions": 12,  # TODO: Double check with Joey
-        "num_spin_sectors": 1,
         "support_variables": [],  # No support variables for this one
         "variable_names": HI_INST_COUNTS_SINGLES_VARIABLE_NAMES,
     },
     CODICEAPID.COD_HI_OMNI_SPECIES_COUNTS: {
-        "coords": [
-            "epoch",
-            "inst_az",
-            "spin_sector",
-            "esa_step",
-            "energy_label",
-        ],  # TODO: These will likely change
         "dataset_name": "imap_codice_l1a_hi-omni",
-        "dims": [
-            "epoch",
-            "esa_step",
-            "inst_az",
-            "spin_sector",
-        ],  # TODO: These will likely change
+        "dims": {
+            "esa_step": 15,
+            "inst_az": 4,
+            "spin_sector": 1,
+        },  # TODO: Double check with Joey
         "instrument": "hi",
         "num_counters": 8,
-        "num_energy_steps": 15,  # TODO: Double check with Joey
-        "num_positions": 4,  # TODO: Double check with Joey
-        "num_spin_sectors": 1,
         "support_variables": ["data_quality", "spin_period"],
         "variable_names": HI_OMNI_SPECIES_VARIABLE_NAMES,
     },
     CODICEAPID.COD_HI_SECT_SPECIES_COUNTS: {
-        "coords": [
-            "epoch",
-            "inst_az",
-            "spin_sector",
-            "esa_step",
-            "energy_label",
-        ],  # TODO: These will likely change
         "dataset_name": "imap_codice_l1a_hi-sectored",
-        "dims": [
-            "epoch",
-            "esa_step",
-            "inst_az",
-            "spin_sector",
-        ],  # TODO: These will likely change
+        "dims": {
+            "esa_step": 8,
+            "inst_az": 12,
+            "spin_sector": 12,
+        },  # TODO: Double check with Joey
         "instrument": "hi",
         "num_counters": 4,
-        "num_energy_steps": 8,
-        "num_positions": 12,
-        "num_spin_sectors": 12,
         "support_variables": ["data_quality", "spin_period"],
         "variable_names": HI_SECT_SPECIES_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_INST_COUNTS_AGGREGATED: {
-        "coords": [
-            "epoch",
-            "inst_az",
-            "spin_sector",
-            "esa_step",
-            "energy_label",
-        ],  # TODO: These will likely change
         "dataset_name": "imap_codice_l1a_lo-counters-aggregated",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"spin_sector_pairs": 6, "esa_step": 128},
         "instrument": "lo",
-        "num_counters": 1,
-        "num_energy_steps": 128,
-        "num_positions": 6,
-        "num_spin_sectors": 6,
+        "num_counters": len(
+            LO_COUNTERS_AGGREGATED_VARIABLE_NAMES
+        ),  # The number of counters depencds on the number of active counters
         "support_variables": [
             "energy_table",
             "acquisition_time_per_step",
-        ],  # TODO: Double check with Joey
-        "variable_names": LO_INST_COUNTS_AGGREGATED_VARIABLE_NAMES,
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
+            "data_quality",
+            "spin_period",
+        ],
+        "variable_names": LO_COUNTERS_AGGREGATED_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_INST_COUNTS_SINGLES: {
-        "coords": [
-            "epoch",
-            "inst_az",
-            "spin_sector",
-            "esa_step",
-            "energy_label",
-        ],  # TODO: These will likely change
         "dataset_name": "imap_codice_l1a_lo-counters-singles",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"inst_az": 24, "spin_sector_pairs": 6, "esa_step": 128},
         "instrument": "lo",
         "num_counters": 1,
-        "num_energy_steps": 128,
-        "num_positions": 24,
-        "num_spin_sectors": 6,
         "support_variables": [
-            "spin_sector_pairs",
             "energy_table",
             "acquisition_time_per_step",
-            "RGFO_half_spin",
-            "NGO_half_spin",
-            "SW_bias_gain_mode",
-            "ST_bias_gain_mode",
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
             "data_quality",
             "spin_period",
         ],
-        "variable_names": LO_INST_COUNTS_SINGLES_VARIABLE_NAMES,
+        "variable_names": LO_COUNTERS_SINGLES_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_SW_ANGULAR_COUNTS: {
-        "coords": ["epoch", "energy_label", "esa_step", "inst_az", "spin_sector"],
         "dataset_name": "imap_codice_l1a_lo-sw-angular",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"inst_az": 5, "spin_sector": 12, "esa_step": 128},
         "instrument": "lo",
         "num_counters": 4,
-        "num_energy_steps": 128,
-        "num_positions": 5,
-        "num_spin_sectors": 12,
         "support_variables": [
             "energy_table",
             "acquisition_time_per_step",
-            "RGFO_half_spin",
-            "NGO_half_spin",
-            "SW_bias_gain_mode",
-            "ST_bias_gain_mode",
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
             "data_quality",
             "spin_period",
         ],
         "variable_names": LO_SW_ANGULAR_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_NSW_ANGULAR_COUNTS: {
-        "coords": ["epoch", "inst_az", "spin_sector", "esa_step", "energy_label"],
         "dataset_name": "imap_codice_l1a_lo-nsw-angular",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"inst_az": 19, "spin_sector": 12, "esa_step": 128},
         "instrument": "lo",
         "num_counters": 1,
-        "num_energy_steps": 128,
-        "num_positions": 19,
-        "num_spin_sectors": 12,
         "support_variables": [
             "energy_table",
             "acquisition_time_per_step",
-            "RGFO_half_spin",
-            "NGO_half_spin",
-            "SW_bias_gain_mode",
-            "ST_bias_gain_mode",
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
             "data_quality",
             "spin_period",
         ],
         "variable_names": LO_NSW_ANGULAR_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_SW_PRIORITY_COUNTS: {
-        "coords": ["epoch", "inst_az", "spin_sector", "esa_step", "energy_label"],
         "dataset_name": "imap_codice_l1a_lo-sw-priority",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"spin_sector": 12, "esa_step": 128},
         "instrument": "lo",
         "num_counters": 5,
-        "num_energy_steps": 128,
-        "num_positions": 1,
-        "num_spin_sectors": 12,
         "support_variables": [
             "energy_table",
             "acquisition_time_per_step",
-            "RGFO_half_spin",
-            "NGO_half_spin",
-            "SW_bias_gain_mode",
-            "ST_bias_gain_mode",
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
             "data_quality",
             "spin_period",
         ],
         "variable_names": LO_SW_PRIORITY_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_NSW_PRIORITY_COUNTS: {
-        "coords": ["epoch", "inst_az", "spin_sector", "esa_step", "energy_label"],
         "dataset_name": "imap_codice_l1a_lo-nsw-priority",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"spin_sector": 12, "esa_step": 128},
         "instrument": "lo",
         "num_counters": 2,
-        "num_energy_steps": 128,
-        "num_positions": 1,
-        "num_spin_sectors": 12,
         "support_variables": [
             "energy_table",
             "acquisition_time_per_step",
-            "RGFO_half_spin",
-            "NGO_half_spin",
-            "SW_bias_gain_mode",
-            "ST_bias_gain_mode",
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
             "data_quality",
             "spin_period",
         ],
         "variable_names": LO_NSW_PRIORITY_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_SW_SPECIES_COUNTS: {
-        "coords": ["epoch", "inst_az", "spin_sector", "esa_step", "energy_label"],
         "dataset_name": "imap_codice_l1a_lo-sw-species",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"spin_sector": 1, "esa_step": 128},
         "instrument": "lo",
         "num_counters": 16,
-        "num_energy_steps": 128,
-        "num_positions": 1,
-        "num_spin_sectors": 1,
         "support_variables": [
             "energy_table",
             "acquisition_time_per_step",
-            "RGFO_half_spin",
-            "NGO_half_spin",
-            "SW_bias_gain_mode",
-            "ST_bias_gain_mode",
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
             "data_quality",
             "spin_period",
         ],
         "variable_names": LO_SW_SPECIES_VARIABLE_NAMES,
     },
     CODICEAPID.COD_LO_NSW_SPECIES_COUNTS: {
-        "coords": ["epoch", "inst_az", "spin_sector", "esa_step", "energy_label"],
         "dataset_name": "imap_codice_l1a_lo-nsw-species",
-        "dims": ["epoch", "esa_step", "inst_az", "spin_sector"],
+        "dims": {"spin_sector": 1, "esa_step": 128},
         "instrument": "lo",
         "num_counters": 8,
-        "num_energy_steps": 128,
-        "num_positions": 1,
-        "num_spin_sectors": 1,
         "support_variables": [
             "energy_table",
             "acquisition_time_per_step",
-            "RGFO_half_spin",
-            "NGO_half_spin",
-            "SW_bias_gain_mode",
-            "ST_bias_gain_mode",
+            "rgfo_half_spin",
+            "nso_half_spin",
+            "sw_bias_gain_mode",
+            "st_bias_gain_mode",
             "data_quality",
             "spin_period",
         ],