diff --git a/modules/contrail/util.py b/modules/contrail/util.py index ded285a35e2e632d0f76091af00b1df82bd70aba..beecd27bfd5c7ff9bcc017583f7ab7583eba4f8a 100644 --- a/modules/contrail/util.py +++ b/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')