diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 9e8a7a4048da7727dc611673b39fc7ee1a793eeb..f0007983d64d2397054ce8296d7583bb4329d318 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -4,7 +4,7 @@ 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, \
-    check_oblique, get_timestamp, homedir, write_icing_file, make_for_full_domain_predict, \
+    check_oblique, get_timestamp, homedir, write_icing_file, write_icing_file_nc4, make_for_full_domain_predict, \
     make_for_full_domain_predict2, get_indexes_within_threshold
 from util.plot import make_icing_image
 from util.geos_nav import get_navigation, get_lon_lat_2d_mesh
@@ -2224,7 +2224,8 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
             np.where(keep, preds, -1)
             np.where(keep, probs, -1.0)
 
-        write_icing_file(clvrx_str_time, output_dir, preds_2d_dct, probs_2d_dct, x_rad, y_rad, lons_2d, lats_2d, cc, ll)
+        # write_icing_file(clvrx_str_time, output_dir, preds_2d_dct, probs_2d_dct, x_rad, y_rad, lons_2d, lats_2d, cc, ll)
+        write_icing_file_nc4(clvrx_str_time, output_dir, preds_2d_dct, probs_2d_dct, x_rad, y_rad, lons_2d, lats_2d, cc, ll)
 
         print('Done: ', clvrx_str_time)
         h5f.close()