snapshot...

modules/contrail/util.py

@@ -119,11 +119,11 @@ def extract(mask_image, image_ts, clavrx_path):
     for key in bins_dict.keys():
         print('working on pressure level: ', bins[key])
         for c_idx in bins_dict[key]:
+            press = contrail_press[c_idx]
+            lat = contrail_lats[c_idx]
             lon = contrail_lons[c_idx]
             if lon < 0:  # Match GFS convention
                 lon += 360.0
-            lat = contrail_lats[c_idx]
-            press = contrail_press[c_idx]
 
             horz_shear_value = horz_shear_3d.interp(Pressure=press, Longitude=lon, Latitude=lat, method='nearest')
             static_value = static_3d.interp(Pressure=press, Longitude=lon, Latitude=lat, method='nearest')