diff --git a/modules/icing/pireps.py b/modules/icing/pireps.py
index 3c1da4b4c029619f3dce55ca895d8990af32737c..b6f17afa7ea91fbac716907e328f16ac4f03b94e 100644
--- a/modules/icing/pireps.py
+++ b/modules/icing/pireps.py
@@ -109,6 +109,20 @@ def pirep_icing(filename, lon_range=[-180, 180], lat_range=[-55, 55]):
 
     return ice_dict, no_ice_dict
 
+def get_goes_datasource(timestamp):
+    dt_obj, time_tup = get_time_tuple_utc(timestamp)
+
+    yr_dir = str(dt_obj.timetuple().tm_year)
+    date_dir = dt_obj.strftime(dir_fmt)
+    yr_date = yr_dir+date_dir
+
+    files_path = goes16_directory + '/' + yr_dir + '/' + date_dir + '/abi' + '/L1b' + '/RadC/'
+    print(files_path)
+    ds = goes_ds_dct.get(yr_date)
+    if ds is None:
+        ds = GOESL1B(files_path)
+        goes_ds_dct[yr_date] = ds
+    return ds
 
 def analyze(ice_dct, no_ice_dct):
     ice_times = np.array(list(ice_dct.keys()))