Harrison - LFAR/PMCL - WIP adding lat/lon to glows l3e

Author: Menlo Innovations - CAVA Project

imap_l3_processing/cdf/config/imap_glows_l3e_survival-probability-ul-hf_variable_attrs.yaml

@@ -198,4 +198,30 @@ elongation_excluded:
   FILLVAL: 65535
   VALIDMAX: 0
   VALIDMIN: 180
-  RECORD_VARYING: RV
+  RECORD_VARYING: RV
+pixel_longitude:
+  NAME: pixel_longitude
+  DATA_TYPE: CDF_FLOAT
+  CATDESC: longitude of the pixel in ecliptic coordinate
+  VAR_TYPE: metadata
+  FIELDNAM: Pixel longitude
+  LABLAXIS: longitude
+  RECORD_VARYING: RV
+  FORMAT: F9.3
+  FILLVAL: -1.0000E+31
+  UNITS: degrees
+  VALIDMIN: 0
+  VALIDMAX: 360
+pixel_latitude:
+  NAME: pixel_latitude
+  DATA_TYPE: CDF_FLOAT
+  CATDESC: latitude of the pixel in ecliptic coordinate
+  VAR_TYPE: metadata
+  FIELDNAM: Pixel latitude
+  LABLAXIS: latitude
+  RECORD_VARYING: RV
+  FORMAT: F9.3
+  FILLVAL: -1.0000E+31
+  UNITS: degrees
+  VALIDMIN: -90
+  VALIDMAX: 90

imap_l3_processing/cdf/config/imap_glows_l3e_survival-probability-ul-sf_variable_attrs.yaml

@@ -198,4 +198,30 @@ elongation_excluded:
   FILLVAL: 65535
   VALIDMAX: 0
   VALIDMIN: 180
-  RECORD_VARYING: RV
+  RECORD_VARYING: RV
+pixel_longitude:
+  NAME: pixel_longitude
+  DATA_TYPE: CDF_FLOAT
+  CATDESC: longitude of the pixel in ecliptic coordinate
+  VAR_TYPE: metadata
+  FIELDNAM: Pixel longitude
+  LABLAXIS: longitude
+  RECORD_VARYING: RV
+  FORMAT: F9.3
+  FILLVAL: -1.0000E+31
+  UNITS: degrees
+  VALIDMIN: 0
+  VALIDMAX: 360
+pixel_latitude:
+  NAME: pixel_latitude
+  DATA_TYPE: CDF_FLOAT
+  CATDESC: latitude of the pixel in ecliptic coordinate
+  VAR_TYPE: metadata
+  FIELDNAM: Pixel latitude
+  LABLAXIS: latitude
+  RECORD_VARYING: RV
+  FORMAT: F9.3
+  FILLVAL: -1.0000E+31
+  UNITS: degrees
+  VALIDMIN: -90
+  VALIDMAX: 90

imap_l3_processing/glows/l3e/glows_l3e_ultra_model.py

@@ -42,6 +42,7 @@ class GlowsL3EUltraData(DataProduct):
     spacecraft_velocity_y: np.ndarray
     spacecraft_velocity_z: np.ndarray
     elongation_excluded: np.ndarray
+    pixel_latitude: np.ndarray
 
     @classmethod
     def convert_dat_to_glows_l3e_ul_product(cls, input_metadata: InputMetadata, file_path: Path,

tests/glows/l3e/test_glows_l3e_ultra_model.py

@@ -70,7 +70,7 @@ def test_convert_dat_to_glows_l3e_ul_product(self):
              25.4914514, 33.1831812, 43.1957952, 56.2295915, 73.1961745, 95.2822137, 124.0324417, 161.4576950,
              210.1755551, 273.5934261, 356.1468544])
 
-        row_1_probability_of_survival = np.array(
+        row_0_probability_of_survival = np.array(
             [0.84827693E+00, 0.86517360E+00, 0.88086645E+00, 0.89551288E+00, 0.90925039E+00, 0.92211605E+00,
              0.93428862E+00, 0.94576000E+00, 0.95660359E+00, 0.96660933E+00, 0.97533752E+00, 0.98226689E+00,
              0.98714107E+00, 0.99025210E+00, 0.99218849E+00, 0.99346478E+00, 0.99439230E+00, 0.99512632E+00,
@@ -123,7 +123,7 @@ def test_convert_dat_to_glows_l3e_ul_product(self):
         np.testing.assert_equal(expected_heal_pix, l3e_ul_product.healpix_index, strict=True)
         np.testing.assert_equal(expected_survival_probability_shape, l3e_ul_product.probability_of_survival.shape,
                                 strict=True)
-        np.testing.assert_equal(l3e_ul_product.probability_of_survival[0].T[0, :], row_1_probability_of_survival,
+        np.testing.assert_equal(l3e_ul_product.probability_of_survival[0].T[0, :], row_0_probability_of_survival,
                                 strict=True)
         np.testing.assert_equal(l3e_ul_product.probability_of_survival[0].T[804, :],
                                 row_804_probability_of_survival, strict=True)
@@ -144,3 +144,10 @@ def test_convert_dat_to_glows_l3e_ul_product(self):
         np.testing.assert_equal(np.array([2.3]), l3e_ul_product.spacecraft_velocity_z, strict=True)
 
         np.testing.assert_equal(np.array([30.0]), l3e_ul_product.elongation_excluded, strict=True)
+
+        np.testing.assert_equal(87.07582, l3e_ul_product.pixel_latitude[0])
+        np.testing.assert_equal(45.00000, l3e_ul_product.pixel_longitude[0])
+
+        np.testing.assert_equal(np.nan, l3e_ul_product.pixel_latitude[804])
+        np.testing.assert_equal(np.nan, l3e_ul_product.pixel_longitude[804])
+

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-01-26T13:05:08.238429",
+    "test_ultra_all_sp_maps": "2026-03-06T13:53:31.090028",
     "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"