diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index bf47973a94ffd425ab04125aabfe84e2172f6695..b9bedf80632fba5f036db1345bc2fd8f8e1a5ebc 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -23,9 +23,15 @@ time_keys = None
 pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv'
 #pirep_file = '/home/rink/data/pireps/pireps_20180101_20200331.csv'
 
-ds_list = ['temp_8_5um_nom', 'temp_10_4um_nom', 'temp_11_0um_nom', 'temp_13_3um_nom', 'cld_height_acha',
+l1b_ds_list = ['temp_10_4um_nom', 'temp_11_0um_nom', 'temp_12_0um_nom', 'temp_13_3um_nom', 'temp_3_75um_nom',
+               'temp_6_2um_nom', 'temp_6_7um_nom', 'temp_7_3_nom', 'temp_8_5um_nom', 'temp_9_7um_nom',
+               'refl_0_47um_nom', 'refl_0_55_nom', 'refl_0_65um_nom', 'refl_0_86um_nom', 'refl_1_38um_nom',
+               'refl_1_60um_nom']
+l1b_grd_dct = {name: [] for name in l1b_ds_list}
+
+ds_list = ['cld_height_acha', 'cld_geo_thick', 'cld_press_acha', 'sensor_zenith_angle', 'cloud_type', 'supercooled_prob_acha',
            'supercooled_cloud_fraction', 'cld_temp_acha', 'cld_opd_acha', 'cloud_phase', 'solar_zenith_angle',
-           'cloud_mask']
+           'cloud_mask', 'cld_reff_acha']
 ds_grd_dct = {name: [] for name in ds_list}
 
 a_clvr_file = '/Users/tomrink/data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc'
@@ -163,6 +169,12 @@ def run(pirep_dct, outfile=None):
                 gvals = get_grid_values(h5f, ds_name, ll[0], cc[0], 20)
                 if gvals is not None:
                     ds_grd_dct[ds_name].append(gvals)
+
+            for didx, ds_name in enumerate(l1b_ds_list):
+                gvals = get_grid_values(h5f, ds_name, ll[0], cc[0], 20)
+                if gvals is not None:
+                    l1b_grd_dct[ds_name].append(gvals)
+
             lon_c.append(lon_s[0])
             lat_c.append(lat_s[0])
             time_s.append(time)
@@ -176,6 +188,12 @@ def run(pirep_dct, outfile=None):
 
     create_file(outfile, data_dct, ds_list, lon_c, lat_c, time_s)
 
+    data_dct = {}
+    for ds_name in ds_list:
+        data_dct[ds_name] = np.array(l1b_grd_dct[ds_name])
+
+    create_file(outfile, data_dct, l1b_ds_list, lon_c, lat_c, time_s)
+
 
 def analyze(ice_dct, no_ice_dct):