snapshot...

modules/contrail/util.py

@@ -48,13 +48,12 @@ def extract(mask_image, image_ts, clavrx_path):
     clvrx_lons = get_grid_values_all(clvrx_h5f, 'longitude').flatten()
     clvrx_lats = get_grid_values_all(clvrx_h5f, 'latitude').flatten()
 
-    # Assuming GOES FD for now -------------------
-    elems, lines = np.meshgrid(np.arange(5424), np.arange(5424))
-    lines, elems = lines.flatten(), elems.flatten()
-
     contrail_idxs = (mask_image == 1).flatten()
     print('number of contrail pixels: ', np.sum(contrail_idxs))
 
+    # Assuming GOES FD for now -------------------
+    # elems, lines = np.meshgrid(np.arange(5424), np.arange(5424))
+    # lines, elems = lines.flatten(), elems.flatten()
     # See note above regarding support for GOES-18
     # contrail_lines, contrail_elems = lines[contrail_idxs], elems[contrail_idxs]
     # contrail_lons, contrail_lats = geos_nav.lc_to_earth(contrail_elems, contrail_lines)
@@ -68,7 +67,7 @@ def extract(mask_image, image_ts, clavrx_path):
     contrail_lats = contrail_lats[keep]
 
     wind = get_point_s(xr_dataset, ['u-wind','v-wind'], contrail_lons, contrail_lats, contrail_press)
-    print(wind)
+