diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 55a4b54ca64e2bcb4f2079415fd7b6f3eebf97a4..05486b84f76315f840328b684a03fbc6da2bd00a 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -379,6 +379,30 @@ def run(pirep_dct, platform, outfile=None, outfile_l1b=None, dt_str_start=None, create_file(outfile_l1b, data_dct, l1b_ds_list, l1b_ds_types, lon_c, lat_c, time_s, fl_alt_s, ice_int_s, unq_ids) +def pirep_info(pirep_dct): + time_keys = list(pirep_dct.keys()) + + lat_s = [] + lon_s = [] + flt_lvl_s = [] + ice_intensity_s = [] + for tkey in time_keys: + reports = pirep_dct[tkey] + for tup in reports: + lat, lon, fl, I, uid, rpt_str = tup + lat_s.append(lat) + lon_s.append(lon) + flt_lvl_s.append(fl) + ice_intensity_s.append(I) + + lat_s = np.array(lat_s) + lon_s = np.array(lon_s) + flt_lvl_s = np.array(flt_lvl_s) + ice_intensity_s = np.array(ice_intensity_s) + + return flt_lvl_s, ice_intensity_s, lat_s, lon_s + + def analyze(ice_dct, no_ice_dct): last_file = None