Harrison 3009 - EYEA/DLAU - GLOWS: Read pixel lat/lon values from Ultra L3e output

Author: Menlo Innovations - CAVA Project

imap_l3_processing/glows/l3e/glows_l3e_ultra_model.py

@@ -24,6 +24,8 @@
 SPACECRAFT_VELOCITY_Y_VAR_NAME = "spacecraft_velocity_y"
 SPACECRAFT_VELOCITY_Z_VAR_NAME = "spacecraft_velocity_z"
 ELONGATION_EXCLUDED_VAR_NAME = "elongation_excluded"
+PIXEL_LATITUDE_VAR_NAME = "pixel_latitude"
+PIXEL_LONGITUDE_VAR_NAME = "pixel_longitude"
 
 
 @dataclass
@@ -43,6 +45,7 @@ class GlowsL3EUltraData(DataProduct):
     spacecraft_velocity_z: np.ndarray
     elongation_excluded: np.ndarray
     pixel_latitude: np.ndarray
+    pixel_longitude: np.ndarray
 
     @classmethod
     def convert_dat_to_glows_l3e_ul_product(cls, input_metadata: InputMetadata, file_path: Path,
@@ -63,11 +66,18 @@ def convert_dat_to_glows_l3e_ul_product(cls, input_metadata: InputMetadata, file
 
         existing_healpix = data_table[:, 0]
         probability_of_survival = data_table[:, 3:]
+        pixel_latitude = data_table[:, 1]
+        pixel_longitude = data_table[:, 2]
+
 
         probability_of_survival_to_return = np.full((len(energies), len(healpix_indexes)), np.nan, dtype=float)
+        pixel_latitude_to_return = np.full(len(healpix_indexes), np.nan, dtype=float)
+        pixel_longitude_to_return = np.full(len(healpix_indexes), np.nan, dtype=float)
 
-        for healpix, prob_sur in zip(existing_healpix, probability_of_survival):
+        for healpix, prob_sur, pixel_lat, pixel_lon in zip(existing_healpix, probability_of_survival, pixel_latitude, pixel_longitude):
             probability_of_survival_to_return[:, int(healpix)] = prob_sur
+            pixel_latitude_to_return[int(healpix)] = pixel_lat
+            pixel_longitude_to_return[int(healpix)] = pixel_lon
 
         transposed_prob_sur = np.array([probability_of_survival_to_return])
 
@@ -87,6 +97,8 @@ def convert_dat_to_glows_l3e_ul_product(cls, input_metadata: InputMetadata, file
             spacecraft_velocity_y=np.array([args.spacecraft_velocity_y]),
             spacecraft_velocity_z=np.array([args.spacecraft_velocity_z]),
             elongation_excluded=np.array([args.elongation]),
+            pixel_latitude=pixel_latitude_to_return,
+            pixel_longitude=pixel_longitude_to_return,
         )
 
     def to_data_product_variables(self) -> list[DataProductVariable]:
@@ -109,4 +121,6 @@ def to_data_product_variables(self) -> list[DataProductVariable]:
             DataProductVariable(SPACECRAFT_VELOCITY_Y_VAR_NAME, self.spacecraft_velocity_y),
             DataProductVariable(SPACECRAFT_VELOCITY_Z_VAR_NAME, self.spacecraft_velocity_z),
             DataProductVariable(ELONGATION_EXCLUDED_VAR_NAME, self.elongation_excluded),
+            DataProductVariable(PIXEL_LATITUDE_VAR_NAME, self.pixel_latitude),
+            DataProductVariable(PIXEL_LONGITUDE_VAR_NAME, self.pixel_longitude),
         ]

tests/glows/l3e/test_glows_l3e_ultra_model.py

@@ -27,7 +27,9 @@ def test_l3e_ultra_model_to_data_product_variables(self):
             sentinel.spacecraft_velocity_x,
             sentinel.spacecraft_velocity_y,
             sentinel.spacecraft_velocity_z,
-            sentinel.elongation_excluded
+            sentinel.elongation_excluded,
+            sentinel.pixel_latitude,
+            sentinel.pixel_longitude,
         )
 
         expected_energy_labels = ['Energy Label 1', 'Energy Label 2', 'Energy Label 3', 'Energy Label 4',
@@ -57,6 +59,8 @@ def test_l3e_ultra_model_to_data_product_variables(self):
             DataProductVariable("spacecraft_velocity_y", sentinel.spacecraft_velocity_y),
             DataProductVariable("spacecraft_velocity_z", sentinel.spacecraft_velocity_z),
             DataProductVariable("elongation_excluded", sentinel.elongation_excluded),
+            DataProductVariable("pixel_latitude", sentinel.pixel_latitude),
+            DataProductVariable("pixel_longitude", sentinel.pixel_longitude),
         ]
 
         self.assertEqual(expected_data_products, data_products)

tests/periodically_run_tests.json

@@ -2,7 +2,7 @@
     "test_l3bcde_first_time_processing": "2026-03-05T09:59:10.392514",
     "test_spectral_fit_against_validation_data": "2026-01-23T10:39:49.234684",
     "test_calculate_pickup_ions_with_minimize": "2026-03-05T13:18:02.645143",
-    "test_ultra_all_sp_maps": "2026-03-06T13:53:31.090028",
+    "test_ultra_all_sp_maps": "2026-03-12T09:24:59.822141",
     "test_ultra_combined_maps": null,
     "test_ultra_combined_nsp_maps": "2026-03-06T09:59:07.820402",
     "test_ultra_combined_sp_maps": "2026-03-05T15:39:34.795512"