diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 31b5f166a9d84f788839512488ad59f1286e3b91..1f5c71dc20ae0ec351b2ab9392e7e6f188d97cf4 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -3,7 +3,7 @@ import numpy as np
 import pickle
 import matplotlib.pyplot as plt
 import os
-from util.util import get_time_tuple_utc, GenericException, add_time_range_to_filename, is_night, is_day
+from util.util import get_time_tuple_utc, GenericException, add_time_range_to_filename, is_night, is_day, get_grid_values_all
 from aeolus.datasource import CLAVRx, GOESL1B
 from util.geos_nav import GEOSNavigation
 import h5py
@@ -303,6 +303,44 @@ def run(pirep_dct, outfile=None, outfile_l1b=None, dt_str_start=None, dt_str_end
         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 make_for_full_domain_predict(clvrx_file, name_list=l1b_ds_list):
+    w_x = 16
+    w_y = 16
+    xlen = 2500
+    ylen = 1500
+
+    h5f = h5py.File(clvrx_file, 'r')
+
+    grd_dct = {name: [] for name in name_list}
+
+    cnt_a = 0
+    for didx, ds_name in enumerate(name_list):
+        gvals = get_grid_values_all(h5f, ds_name)
+        if gvals is not None:
+            grd_dct[ds_name].append(gvals)
+            cnt_a += 1
+
+    if cnt_a > 0 and cnt_a != len(ds_list):
+        raise GenericException('weirdness')
+
+    grd_dct_n = {name: [] for name in name_list}
+
+    n_x = int(xlen/w_x)
+    n_y = int(ylen/w_y)
+
+    i_0 = 0
+    j_0 = 0
+
+    for didx, ds_name in enumerate(ds_list):
+        for j in range(4, n_y-4, 1):
+            j_ul = j_0 + j * w_y
+            for i in range(4, n_x-4, 1):
+                i_ul = i_0 + i * w_x
+                grd_dct_n[ds_name].append(grd_dct[ds_name][j_ul:j_ul+w_y, i_ul:i_ul+w_x])
+
+    return grd_dct_n
+
+
 def analyze(ice_dct, no_ice_dct):
 
     last_file = None