From 564b1dfbc005569618971c34d66d826d915f75f9 Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Mon, 9 Aug 2021 13:16:10 -0500 Subject: [PATCH] new method --- modules/icing/pirep_goes.py | 23 +++++++++++++++++++++++ 1 file changed, 23 insertions(+) diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 5e38a529..2635dd2d 100644 --- a/modules/icing/pirep_goes.py +++ b/modules/icing/pirep_goes.py @@ -1490,6 +1490,29 @@ def spatial_filter(icing_dict): return grd_bins +# dt_str_0: start datetime string in format YYYY-MM-DD_HH:MM +# dt_str_1: end datetime string in format YYYY-MM-DD_HH:MM +# return a flatten list of icing reports +def time_filter(dt_str_0, dt_str_1, icing_dct): + dto_0 = datetime.datetime.strptime(dt_str_0, '%Y-%m-%d_%H:%M').replace(tzinfo=timezone.utc) + ts_0 = dto_0.timestamp() + + dto_1 = datetime.datetime.strptime(dt_str_1, '%Y-%m-%d_%H:%M').replace(tzinfo=timezone.utc) + ts_1 = dto_1.timestamp() + + keep_reports = [] + keep_times = [] + + for ts in list(icing_dct.keys()): + if not (ts_0 <= ts < ts_1): + rpts = icing_dct[ts] + for idx, tup in enumerate(rpts): + keep_reports.append(tup) + keep_times.append(ts) + + return keep_times, keep_reports + + def analyze_moon_phase(icing_dict): ts = api.load.timescale() eph = api.load('de421.bsp') -- GitLab