diff --git a/modules/aeolus/aeolus_amv.py b/modules/aeolus/aeolus_amv.py index 8b25ffd1d770eb26d529bb05718bbad8f9d05277..acc7bfc5e41901cd5b4d1f4677c4e0e61524e656 100644 --- a/modules/aeolus/aeolus_amv.py +++ b/modules/aeolus/aeolus_amv.py @@ -14,6 +14,7 @@ from metpy.units import units from util.gfs_reader import get_vert_profile_s from amv.intercompare import get_raob_dict_cdf from util.line_plot import do_plot +import h5py amv_file_duration = 60 # minutes @@ -1745,6 +1746,7 @@ def match_calipso_clavrx_to_amvs(calipso_clavrx_path, calipso_clavrx_file, amv_f calipso_clavrx_ds = CLAVRx_CALIPSO(calipso_clavrx_path) nom_time = calipso_clavrx_ds.get_datetime(calipso_clavrx_file).timestamp() + calipso_clavrx_nc4 = Dataset(calipso_clavrx_file) fname, ftime, f_idx = amv_files.get_file_containing_time(nom_time) if f_idx is None: @@ -1775,8 +1777,8 @@ def match_calipso_clavrx_to_amvs(calipso_clavrx_path, calipso_clavrx_file, amv_f amv_ds.close() - xs = calipso_clavrx_ds['x'][:] - ys = calipso_clavrx_ds['y'][:] + xs = calipso_clavrx_nc4['x'][:] + ys = calipso_clavrx_nc4['y'][:] clvr_xy_s = zip(xs, ys) for clvr_xy in clvr_xy_s: