imap-processing 1.0.2__py3-none-any.whl → 1.0.3__py3-none-any.whl

imap_processing/_version.py

@@ -1,3 +1,3 @@
 # These version placeholders will be replaced later during substitution.
-__version__ = "1.0.2"
-__version_tuple__ = (1, 0, 2)
+__version__ = "1.0.3"
+__version_tuple__ = (1, 0, 3)

imap_processing/cdf/config/imap_codice_global_cdf_attrs.yaml

@@ -241,6 +241,24 @@ imap_codice_l2_lo-sw-species:
   Logical_source: imap_codice_l2_lo-sw-species
   Logical_source_description: IMAP Mission CoDICE Lo Level-2 Sunward Species Intensity Data.
 
+imap_codice_l2_lo-nsw-species:
+  <<: *instrument_base
+  Data_type: L2_lo-nsw-species>Level-2 Lo Non-Sunward Species Intensities Data
+  Logical_source: imap_codice_l2_lo-nsw-species
+  Logical_source_description: IMAP Mission CoDICE Lo Level-2 Non-Sunward Species Intensity Data.
+
+imap_codice_l2_lo-sw-angular:
+  <<: *instrument_base
+  Data_type: L2_lo-sw-species>Level-2 Lo Sunward Angular Intensities Data
+  Logical_source: imap_codice_l2_lo-sw-angular
+  Logical_source_description: IMAP Mission CoDICE Lo Level-2 Sunward Angular Intensity Data.
+
+imap_codice_l2_lo-nsw-angular:
+  <<: *instrument_base
+  Data_type: L2_lo-nsw-species>Level-2 Lo Non-Sunward Angular Intensities Data
+  Logical_source: imap_codice_l2_lo-nsw-angular
+  Logical_source_description: IMAP Mission CoDICE Lo Level-2 Non-Sunward Angular Intensity Data.
+
 imap_codice_l2_hi-omni:
   <<: *instrument_base
   Data_type: L2_hi-omni>Level-2 Hi Omnidirectional Species Intensities Data

imap_processing/cdf/config/imap_codice_l1a_variable_attrs.yaml

@@ -297,13 +297,13 @@ direct_events_attrs: &direct_events
     VALIDMAX: 10000
 
 voltage_table:
-    CATDESC: ElectroStatic Analyzer Energy Values
+    CATDESC: ElectroStatic Analyzer Voltage Values
     DEPEND_1: esa_step
-    FIELDNAM: Energy Table
+    FIELDNAM: Voltage Table
     FORMAT: F12.6
-    LABLAXIS: eV
+    LABLAXIS: V
     SCALETYP: log
-    UNITS: eV
+    UNITS: V
     VALIDMIN: 1.0
     VALIDMAX: 14100.0
     VAR_TYPE: support_data

imap_processing/cdf/config/imap_hi_variable_attrs.yaml

@@ -432,8 +432,18 @@ hi_de_nominal_bin:
 
 # Define override values for epoch as defined in imap_constant_attrs.yaml
 hi_pset_epoch:
-  CATDESC: Midpoint time of pointing, number of nanoseconds since J2000 with leap seconds included
-  BIN_LOCATION: 0.5
+  CATDESC: Start time of pointing, number of nanoseconds since J2000 with leap seconds included
+  DELTA_PLUS_VAR: epoch_delta
+  BIN_LOCATION: 0
+
+epoch_delta:
+  <<: *default_int64
+  CATDESC: Number of nanoseconds in spacecraft pointing covered by this pointing set product
+  FIELDNAM: epoch delta
+  UNITS: ns
+  VAR_TYPE: support_data
+  DISPLAY_TYPE: no_plot
+  TIME_SCALE: Terrestrial Time
 
 hi_pset_esa_energy_step:
   <<: *default_esa_energy_step

imap_processing/cdf/config/imap_idex_global_cdf_attrs.yaml

@@ -58,16 +58,9 @@ imap_idex_l2b_sci:
   Logical_source: imap_idex_l2b_sci-1mo
   Logical_source_description: IMAP Mission IDEX Instrument Level-2B Monthly Data
 
-imap_idex_l2c_sci-healpix:
-  <<: *instrument_base
-  Data_level: 2C
-  Data_type: L2C_HEALPIX-MAP-1MO>Level-2C HEALPIX Map Monthly Data
-  Logical_source: imap_idex_l2c_healpix-map-1mo
-  Logical_source_description: IMAP Mission IDEX Instrument Level-2C Monthly Data Healpix
-
 imap_idex_l2c_sci-rectangular:
   <<: *instrument_base
   Data_level: 2C
   Data_type: L2C_RECTANGULAR-MAP-1MO>Level-2C Rectangular Map Monthly Data
   Logical_source: imap_idex_l2c_rectangular-map-1mo
-  Logical_source_description: IMAP Mission IDEX Instrument Level-2C Monthly Data Rectangular
+  Logical_source_description: IMAP Mission IDEX Instrument Level-2C Rectangular Map Monthly Data

imap_processing/cdf/config/imap_idex_l1b_variable_attrs.yaml

@@ -3,6 +3,8 @@ double_fillval: &double_fillval -1.0E31
 
 data_min: &data_min 0
 data_max: &data_max 4095
+spice_data_min: &spice_data_min -2.0e7
+spice_data_max: &spice_data_max 2.0e7
 
 # <=== Base Attributes ===>
 string_base_attrs: &string_base
@@ -47,7 +49,8 @@ spice_base: &spice_base
   <<: *l1b_data_base
   VAR_TYPE: data
   DISPLAY_TYPE: time_series
-  UNITS: Degrees
+  VALIDMIN: *spice_data_min
+  VALIDMAX: *spice_data_max
 
 # <=== Instrument Setting Attributes ===>
 trigger_mode:
@@ -342,7 +345,6 @@ current_neg2p5v_bus:
   UNITS: A
 
 # <=== Spice Data Attributes ===>
-# TODO: Get actual validmin and vaildmax for ephemeris attrs
 ephemeris_position_x:
   <<: *spice_base
   FIELDNAM: Position X
@@ -362,7 +364,7 @@ ephemeris_position_y:
 ephemeris_position_z:
   <<: *spice_base
   FIELDNAM: Position Z
-  CATDESC: Cartesian coordinate Z positions for the IMAP spacecraft in the ECLIPJ2000 frame.
+  CATDESC: Cartesian coord Z positions for the IMAP spacecraft in the ECLIPJ2000 frame.
   LABLAXIS: Position Z
   UNITS: km
   DICT_KEY: SPASE>Support>SupportQuantity:Positional
@@ -395,33 +397,42 @@ ephemeris_velocity_z:
 longitude:
   <<: *spice_base
   VALIDMIN: 0
+  VALIDMAX: 360
   FIELDNAM: Longitude
   CATDESC: Longitude of the dust event as observed from Earth in the ECLIPJ2000 frame.
   LABLAXIS: Longitude
+  UNITS: Degrees
   DICT_KEY: SPASE>Support>SupportQuantity:Longitude
 
 latitude:
   <<: *spice_base
+  VALIDMIN: -90
+  VALIDMAX: 90
   FIELDNAM: Latitude
   CATDESC: Latitude of the dust event as observed from Earth in the ECLIPJ2000 frame.
   LABLAXIS: Latitude
+  UNITS: Degrees
   DICT_KEY: SPASE>Support>SupportQuantity:Latitude
 
 spin_phase:
   <<: *spice_base
   VALIDMIN: 0
+  VALIDMAX: 360
   FIELDNAM: Spin Phase
   CATDESC: IMAP Spin Phase
   LABLAXIS: Spin Phase
-  FORMAT: F12.6
-  FILLVAL: *double_fillval
+  FORMAT: I3
+  FILLVAL: *int_fillval
+  UNITS: Degrees
   DICT_KEY: SPASE>Support>SupportQuantity:SpinPhase
 
 solar_longitude:
   <<: *spice_base
+  VALIDMIN: -180
   VALIDMAX: 180
   FIELDNAM: Solar Longitude
   CATDESC: Solar Longitude of the IMAP spacecraft
   LABLAXIS: Solar Longitude
+  UNITS: Degrees
   DICT_KEY: SPASE>Support>SupportQuantity:Longitude
 

imap_processing/cdf/config/imap_idex_l2a_variable_attrs.yaml

@@ -180,7 +180,7 @@ tof_peak_fit_parameters:
   LABL_PTR_2: peak_fit_parameter_labels
   VALIDMIN: 0
   VALIDMAX: *int_maxval
-  LABLAXIS: TOF Peak Fit Parameters
+  LABLAXIS: TOF Fit Parameters
   FORMAT: F20.6
   FILLVAL: *double_fillval
   DISPLAY_TYPE: no_plot
@@ -205,14 +205,16 @@ tof_peak_area_under_fit:
 tof_peak_chi_squared:
   <<: *chi_square_base
   DEPEND_1: mass_index
-  LABL_PTR_1: mass_labels
+  DISPLAY_TYPE: spectrogram
+  LABLAXIS: Chi Square (TOF)
   CATDESC: Chi squared value for the TOF peak emg fits.
-  FIELDNAM: Chi Square (TOF peak fits))
+  FIELDNAM: Chi Square (TOF peak fits)
 
 tof_peak_reduced_chi_squared:
   <<: *chi_square_base
   DEPEND_1: mass_index
-  LABL_PTR_1: mass_labels
+  DISPLAY_TYPE: spectrogram
+  LABLAXIS: Red Chi Square (TOF)
   CATDESC: Reduced chi squared value for the TOF peak emg fits.
   FIELDNAM: Reduced Chi Square (TOF peak fits)
 
@@ -240,37 +242,37 @@ target_low_fit_parameters:
 
 target_low_impact_charge:
   <<: *impact_charge_base
-  LABLAXIS: Target Low Impact Charge
-  CATDESC: Total charge from the dust impact on Target Low channel derived from the fitted curve.
+  LABLAXIS: TL Impact Charge
+  CATDESC: Total charge from the dust impact on Target Low (TL) channel derived from the fitted curve.
   FIELDNAM: Target Low Impact Charge
 
 target_low_dust_mass_estimate:
   <<: *mass_estimate_base
-  LABLAXIS: Target Low Dust Mass Estimate
+  LABLAXIS: TL Dust Mass Est
   CATDESC: Estimated dust mass from Target Low signal based on impact charge.
   FIELDNAM: Estimated Dust Mass (Target Low)
 
 target_low_velocity_estimate:
   <<: *velocity_estimate_base
-  LABLAXIS: Target Low Velocity Estimate
+  LABLAXIS: TL Velocity Est
   CATDESC: Estimated particle velocity from Target Low signal based on impact charge.
   FIELDNAM: Estimated Velocity (Target Low)
 
 target_low_chi_squared:
   <<: *chi_square_base
-  LABLAXIS: Target Low Chi Square
+  LABLAXIS: TL Chi Square
   CATDESC: Chi squared value for the Target Low signal fit.
   FIELDNAM: Chi Square (Target Low)
 
 target_low_reduced_chi_squared:
   <<: *chi_square_base
-  LABLAXIS: Target Low Reduced Chi Square
+  LABLAXIS: TL Red Chi Square
   CATDESC: Reduced chi squared value for the Target Low signal fit.
   FIELDNAM: Reduced Chi Square (Target Low)
 
 target_low_fit_results:
   <<: *fit_results_base
-  LABLAXIS: Target Low Fit Results
+  LABLAXIS: TL Fit Results
   CATDESC: Values of Target Low signal fit evaluated at each time point.
   FIELDNAM: Target Low Fit results
 
@@ -282,37 +284,37 @@ target_high_fit_parameters:
 
 target_high_impact_charge:
   <<: *impact_charge_base
-  LABLAXIS: Target High Impact Charge
-  CATDESC: Charge from the dust impact on Target High derived from the fitted curve.
+  LABLAXIS: TH Impact Charge
+  CATDESC: Charge from the dust impact on Target High (TH) derived from the fitted curve.
   FIELDNAM: Target High Impact Charge
 
 target_high_dust_mass_estimate:
   <<: *mass_estimate_base
-  LABLAXIS: Target High Dust Mass Estimate
+  LABLAXIS: TH Dust Mass Est
   CATDESC: Estimated dust mass from Target High signal based on impact charge.
   FIELDNAM: Estimated Dust Mass (Target High)
 
 target_high_velocity_estimate:
   <<: *velocity_estimate_base
-  LABLAXIS: Target High Velocity Estimate
+  LABLAXIS: TH Velocity Est
   CATDESC: Estimated particle velocity from Target High signal based on impact charge.
   FIELDNAM: Estimated Velocity (Target High)
 
 target_high_chi_squared:
   <<: *chi_square_base
-  LABLAXIS: Target High Chi Square
+  LABLAXIS: TH Chi Square
   CATDESC: Chi squared value for the Target High signal fit.
   FIELDNAM: Chi Square (Target High)
 
 target_high_reduced_chi_squared:
   <<: *chi_square_base
-  LABLAXIS: Target High Reduced Chi Square
+  LABLAXIS: TH Red Chi Square
   CATDESC: Reduced chi squared value for the Target Low signal fit.
   FIELDNAM: Reduced Chi Square (Target High)
 
 target_high_fit_results:
   <<: *fit_results_base
-  LABLAXIS: Target High Fit Results
+  LABLAXIS: TH Fit Results
   CATDESC: Values of Target High signal fit evaluated at each time point.
   FIELDNAM: Target High Fit Results
 
@@ -324,36 +326,36 @@ ion_grid_fit_parameters:
 
 ion_grid_impact_charge:
   <<: *impact_charge_base
-  LABLAXIS: Ion Grid Impact Charge
-  CATDESC: Charge from the dust impact on Ion Grid channel derived from the fitted curve.
+  LABLAXIS: IG Impact Charge
+  CATDESC: Charge from the dust impact on Ion Grid (IG) channel derived from the fitted curve.
   FIELDNAM: Ion Grid Impact Charge
 
 ion_grid_dust_mass_estimate:
   <<: *mass_estimate_base
-  LABLAXIS: Ion Grid Dust Mass Estimate
+  LABLAXIS: IG Dust Mass Est
   CATDESC: Estimated dust mass from Ion Grid signal based on impact charge.
   FIELDNAM: Estimated Dust Mass (Ion Grid)
 
 ion_grid_velocity_estimate:
   <<: *velocity_estimate_base
-  LABLAXIS: Ion Grid Velocity Estimate
+  LABLAXIS: IG Velocity Est
   CATDESC: Estimated particle velocity from Ion Grid signal based on impact charge.
   FIELDNAM: Estimated Velocity (Ion Grid)
 
 ion_grid_chi_squared:
   <<: *chi_square_base
-  LABLAXIS: Ion Grid Chi Square
+  LABLAXIS: IG Chi Square
   CATDESC: Chi squared value for the Ion Grid signal fit.
   FIELDNAM: Chi Square (Ion Grid)
 
 ion_grid_reduced_chi_squared:
   <<: *chi_square_base
-  LABLAXIS: Ion Grid Reduced Chi Square
+  LABLAXIS: IG Red Chi Square
   CATDESC: Reduced chi squared value for the Ion Grid signal fit.
   FIELDNAM: Reduced Chi Square (Ion Grid)
 
 ion_grid_fit_results:
   <<: *fit_results_base
-  LABLAXIS: Ion Grid Fit Results
+  LABLAXIS: IG Fit Results
   CATDESC: Values of Ion Grid signal fit evaluated at each time point.
   FIELDNAM: Ion Grid Fit Results

imap_processing/cdf/config/imap_idex_l2b_variable_attrs.yaml

@@ -4,8 +4,8 @@ double_fillval: &double_fillval -1.0E31
 
 # Label attributes
 mass_labels:
-  CATDESC: Labels for Mass (kg)
-  FIELDNAM: Mass (kg)
+  CATDESC: Labels for Mass
+  FIELDNAM: Mass
   VAR_TYPE: metadata
   FORMAT: A8
   DEPEND_1: mass
@@ -30,23 +30,23 @@ spin_phase_labels:
 # Index attributes
 mass:
   CATDESC: Log-spaced mass
-  FIELDNAM: Mass (kg)
+  FIELDNAM: Mass
   UNITS: kg
   FORMAT: E10.4
   VAR_TYPE: support_data
   SCALETYP: log
-  LABLAXIS: Mass (kg)
+  LABLAXIS: Mass
   VALIDMIN: 0.0
   VALIDMAX: 1.00e-14
   FILLVAL: *double_fillval
   LABL_PTR_1: mass_labels
-  DICT_KEY: SPASE>Support>SupportQuantity:Other
+  DICT_KEY: SPASE>Particle>ParticleType:Dust,ParticleQuantity:Mass
 
 spin_phase:
     CATDESC: The spacecraft spin phase at the time of detection
     VAR_NOTES: This is given in 4 bins [315-45, 45-135, 135-225, 225-315] degrees.
-    FIELDNAM: Spin Phase (deg)
-    LABLAXIS: Spin Phase (deg)
+    FIELDNAM: Spin Phase
+    LABLAXIS: Spin Phase
     VAR_TYPE: support_data
     SCALETYP: linear
     FILLVAL: *int_fillval
@@ -55,12 +55,12 @@ spin_phase:
     VALIDMAX: 360
     LABL_PTR_1: spin_phase_labels
     UNITS: deg
-    DICT_KEY: SPASE>SupportQuantity:SpinPhase,Qualifier:Array
+    DICT_KEY: SPASE>Support>SupportQuantity:SpinPhase,Qualifier:Array
 
 impact_charge:
     CATDESC: Log-spaced impact charge
-    FIELDNAM: Impact Charge (fC)
-    LABLAXIS: Impact Charge (fC)
+    FIELDNAM: Impact Charge
+    LABLAXIS: Impact Charge
     VAR_TYPE: support_data
     SCALETYP: log
     FILLVAL: *int_fillval
@@ -106,18 +106,18 @@ base_mass: &base_mass
 rate_by_charge:
     <<: *base_charge
     CATDESC: Count rate per day by impact charge and spin phase.
-    FIELDNAM: Rate (day^-1) by Charge
+    FIELDNAM: Rate by Charge
     UNITS: day^-1
     FILLVAL: *double_fillval
-    DICT_KEY: SPASE>SupportQuantity:CountRate,Qualifier:Array
+    DICT_KEY: SPASE>Particle>ParticleType:Dust,ParticleQuantity:CountRate,Qualifier:Array
 
 rate_by_mass:
     <<: *base_mass
     CATDESC: Count rate per day by mass and spin phase.
-    FIELDNAM: Rate (day^-1) by Mass
+    FIELDNAM: Rate by Mass
     UNITS: day^-1
     FILLVAL: *double_fillval
-    DICT_KEY: SPASE>SupportQuantity:CountRate,Qualifier:Array
+    DICT_KEY: SPASE>Particle>ParticleType:Dust,ParticleQuantity:CountRate,Qualifier:Array
 
 counts_by_charge:
     <<: *base_charge
@@ -147,7 +147,7 @@ impact_day_of_year:
     VALIDMAX: 366
     VAR_TYPE: data
     UNITS: days
-    DICT_KEY: SPASE>SupportQuantity:Temporal,Qualifier:Array
+    DICT_KEY: SPASE>Support>SupportQuantity:Temporal
 
 rate_calculation_quality_flags:
   CATDESC: Quality flag for rate calculation (1 = good, 0 = insufficient IDEX uptime data)
@@ -161,4 +161,4 @@ rate_calculation_quality_flags:
   VALIDMAX: 1
   VAR_TYPE: data
   UNITS: " "
-  DICT_KEY: SPASE>SupportQuantity:QualityFlag,Qualifier:Array
+  DICT_KEY: SPASE>Support>SupportQuantity:QualityFlag

imap_processing/cdf/config/imap_idex_l2c_variable_attrs.yaml

@@ -90,7 +90,7 @@ rectangular_lat_pixel:
 rate_by_charge_map:
   <<: *base_charge_map
   CATDESC: Count rate per day by impact charge, longitude, and latitude.
-  FIELDNAM: Rate (day^-1) by Charge Map
+  FIELDNAM: Rate by Charge Map
   UNITS: day^-1
   FILLVAL: *double_fillval
   DICT_KEY: SPASE>SupportQuantity:CountRate,Qualifier:Array
@@ -114,7 +114,7 @@ counts_by_mass_map:
 rate_by_mass_map:
   <<: *base_mass_map
   CATDESC: Count rate per day by mass, longitude, and latitude.
-  FIELDNAM: Rate (day^-1) by Mass Map
+  FIELDNAM: Rate by Mass Map
   UNITS: day^-1
   FILLVAL: *double_fillval
   DICT_KEY: SPASE>SupportQuantity:CountRate,Qualifier:Array

imap_processing/cdf/config/imap_ultra_l1c_variable_attrs.yaml

@@ -172,6 +172,18 @@ shcoarse:
   # TODO: come back to format
   UNITS: seconds
 
+epoch_delta:
+  <<: *default
+  CATDESC: Number of nanoseconds in spacecraft pointing covered by this pointing set product
+  FIELDNAM: epoch delta
+  UNITS: ns
+  VAR_TYPE: support_data
+  DISPLAY_TYPE: no_plot
+  TIME_SCALE: Terrestrial Time
+  SCALE_TYPE: linear
+  FORMAT: I20
+  dtype: int64
+
 energy_bin_delta:
   <<: *default_float32
   CATDESC: Difference between the energy bin edges.

imap_processing/codice/codice_l1b.py

@@ -47,6 +47,23 @@ def convert_to_rates(dataset: xr.Dataset, descriptor: str) -> np.ndarray:
         constants, f"{descriptor.upper().replace('-', '_')}_VARIABLE_NAMES"
     )
 
+    if descriptor.startswith("lo-"):
+        # Calculate energy_table using voltage_table and k_factor
+        energy_attrs = dataset["voltage_table"].attrs | {
+            "UNITS": "keV/e",
+            "LABLAXIS": "E/q",
+            "CATDESC": "Energy per charge",
+            "FIELDNAM": "Energy per charge",
+        }
+        # 1e3 is to convert eV to keV
+        dataset["energy_table"] = xr.DataArray(
+            dataset["voltage_table"].values * dataset["k_factor"].values * 1e-3,
+            dims=[
+                "esa_step",
+            ],
+            attrs=energy_attrs,
+        )
+
     if descriptor in [
         "lo-counters-aggregated",
         "lo-counters-singles",
@@ -69,6 +86,8 @@ def convert_to_rates(dataset: xr.Dataset, descriptor: str) -> np.ndarray:
             "sw_bias_gain_mode",
             "st_bias_gain_mode",
             "spin_period",
+            "voltage_table",
+            "acquisition_time_per_step",
         ]
         dataset = dataset.drop_vars(drop_variables)
     elif descriptor in [
@@ -86,7 +105,6 @@ def convert_to_rates(dataset: xr.Dataset, descriptor: str) -> np.ndarray:
 
         # Denominator to convert counts to rates
         denominator = dataset.acquisition_time_per_step * n_sector
-
         # Do not carry these variable attributes from L1a to L1b for above products
         drop_variables = [
             "k_factor",

imap_processing/codice/codice_l2.py

@@ -31,8 +31,9 @@ from imap_processing.codice.constants import (
     LO_POSITION_TO_ELEVATION_ANGLE,
     LO_SW_ANGULAR_VARIABLE_NAMES,
     LO_SW_PICKUP_ION_SPECIES_VARIABLE_NAMES,
-    LO_SW_SPECIES_VARIABLE_NAMES,
+    LO_SW_SOLAR_WIND_SPECIES_VARIABLE_NAMES,
     NSW_POSITIONS,
+    PIXEL_ORIENTATIONS,
     PUI_POSITIONS,
     SOLAR_WIND_POSITIONS,
     SW_POSITIONS,
@@ -354,11 +355,14 @@ def process_lo_angular_intensity(
         positions,
         average_across_positions=False,
     )
+
     # transform positions to elevation angles
     if positions == SW_POSITIONS:
         pos_to_el = LO_POSITION_TO_ELEVATION_ANGLE["sw"]
+        position_index_to_adjust = 0
     elif positions == NSW_POSITIONS:
         pos_to_el = LO_POSITION_TO_ELEVATION_ANGLE["nsw"]
+        position_index_to_adjust = 9
     else:
         raise ValueError("Unknown positions for elevation angle mapping.")
 
@@ -369,6 +373,8 @@ def process_lo_angular_intensity(
             [pos_to_el[pos] for pos in dataset["inst_az"].data],
         )
     )
+    # add uncertainties to species list
+    species_list = species_list + [f"unc_{var}" for var in species_list]
     # Take the mean across elevation angles and restore the original dimension order
     dataset_converted = (
         dataset[species_list]
@@ -383,8 +389,41 @@ def process_lo_angular_intensity(
     dataset = dataset.assign_coords(
         spin_angle=("spin_sector", dataset["spin_sector"].data * 15.0 + 7.5)
     )
-
     dataset = dataset.drop_vars(species_list).merge(dataset_converted)
+    # Positions 0 and 10 only observe half of the 24 spins for each esa step.
+    # To account for this, we replicate the counts observed in position 0 and 10 for
+    # each esa step to either spin angles 0-11 or 12-23, depending on the pixel
+    # orientation (A/B). See section 11.2.2 of the CoDICE algorithm document
+    a_inds = np.array(
+        [pos for pos, orientation in PIXEL_ORIENTATIONS.items() if orientation == "A"]
+    )
+    b_inds = np.array(
+        [pos for pos, orientation in PIXEL_ORIENTATIONS.items() if orientation == "B"]
+    )
+
+    position_index = position_index_to_adjust
+    for species in species_list:
+        # Determine the correct spin indices based on the position
+        spin_sectors = dataset["spin_sector"].data
+        spin_inds_1 = np.where(spin_sectors >= 12)[0]
+        spin_inds_2 = np.where(spin_sectors < 12)[0]
+        # if position_index is 9, swap the spin indices
+        if position_index == 9:
+            spin_inds_1, spin_inds_2 = spin_inds_2, spin_inds_1
+
+        # Assign the values to the correct positions and spin sectors
+        dataset[species].values[
+            :, a_inds[:, np.newaxis], spin_inds_1, position_index
+        ] = dataset[species].values[
+            :, a_inds[:, np.newaxis], spin_inds_2, position_index
+        ]
+
+        dataset[species].values[
+            :, b_inds[:, np.newaxis], spin_inds_2, position_index
+        ] = dataset[species].values[
+            :, b_inds[:, np.newaxis], spin_inds_1, position_index
+        ]
+
     return dataset
 
 
@@ -777,7 +816,6 @@ def process_codice_l2(
     # This should get science files since ancillary or spice doesn't have data_type
     # as data level.
     file_path = dependencies.get_file_paths(descriptor=descriptor)[0]
-
     # Now form product name from descriptor
     descriptor = ScienceFilePath(file_path).descriptor
     dataset_name = f"imap_codice_l2_{descriptor}"
@@ -790,6 +828,9 @@ def process_codice_l2(
         "imap_codice_l2_lo-nsw-angular",
         "imap_codice_l2_lo-sw-angular",
     ]:
+        cdf_attrs = ImapCdfAttributes()
+        cdf_attrs.add_instrument_global_attrs("codice")
+
         l2_dataset = load_cdf(file_path).copy()
 
         geometric_factor_lookup = get_geometric_factor_lut(dependencies)
@@ -797,12 +838,13 @@ def process_codice_l2(
         geometric_factors = compute_geometric_factors(
             l2_dataset, geometric_factor_lookup
         )
+
         if dataset_name == "imap_codice_l2_lo-sw-species":
             # Filter the efficiency lookup table for solar wind efficiencies
             efficiencies = efficiency_lookup[efficiency_lookup["product"] == "sw"]
             # Calculate the pickup ion sunward solar wind intensities using equation
             # described in section 11.2.3 of algorithm document.
-            process_lo_species_intensity(
+            l2_dataset = process_lo_species_intensity(
                 l2_dataset,
                 LO_SW_PICKUP_ION_SPECIES_VARIABLE_NAMES,
                 geometric_factors,
@@ -811,25 +853,31 @@ def process_codice_l2(
             )
             # Calculate the sunward solar wind species intensities using equation
             # described in section 11.2.3 of algorithm document.
-            process_lo_species_intensity(
+            l2_dataset = process_lo_species_intensity(
                 l2_dataset,
-                LO_SW_SPECIES_VARIABLE_NAMES,
+                LO_SW_SOLAR_WIND_SPECIES_VARIABLE_NAMES,
                 geometric_factors,
                 efficiencies,
                 SOLAR_WIND_POSITIONS,
             )
+            l2_dataset.attrs.update(
+                cdf_attrs.get_global_attributes("imap_codice_l2_lo-sw-species")
+            )
         elif dataset_name == "imap_codice_l2_lo-nsw-species":
             # Filter the efficiency lookup table for non-solar wind efficiencies
             efficiencies = efficiency_lookup[efficiency_lookup["product"] == "nsw"]
             # Calculate the non-sunward species intensities using equation
             # described in section 11.2.3 of algorithm document.
-            process_lo_species_intensity(
+            l2_dataset = process_lo_species_intensity(
                 l2_dataset,
                 LO_NSW_SPECIES_VARIABLE_NAMES,
                 geometric_factors,
                 efficiencies,
                 NSW_POSITIONS,
             )
+            l2_dataset.attrs.update(
+                cdf_attrs.get_global_attributes("imap_codice_l2_lo-nsw-species")
+            )
         elif dataset_name == "imap_codice_l2_lo-sw-angular":
             efficiencies = efficiency_lookup[efficiency_lookup["product"] == "sw"]
             # Calculate the sunward solar wind angular intensities using equation
@@ -841,6 +889,9 @@ def process_codice_l2(
                 efficiencies,
                 SW_POSITIONS,
             )
+            l2_dataset.attrs.update(
+                cdf_attrs.get_global_attributes("imap_codice_l2_lo-sw-angular")
+            )
         if dataset_name == "imap_codice_l2_lo-nsw-angular":
             # Calculate the non sunward angular intensities
             efficiencies = efficiency_lookup[efficiency_lookup["product"] == "nsw"]
@@ -851,7 +902,9 @@ def process_codice_l2(
                 efficiencies,
                 NSW_POSITIONS,
             )
-
+            l2_dataset.attrs.update(
+                cdf_attrs.get_global_attributes("imap_codice_l2_lo-nsw-angular")
+            )
     if dataset_name in [
         "imap_codice_l2_hi-counters-singles",
         "imap_codice_l2_hi-counters-aggregated",