diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index c93e29db63fdd5eacd124de753b674c84a9f7680..715bb872f3d69763cd0381be24c8e68027f8476b 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -929,7 +929,7 @@ def analyze2(filename, filename_l1b):
     # cld_frac = np.sum(cld_mask, axis=1) / 256 # need to do this earlier
 
     #  Simple model: Everything below the melting point is icing=True else icing=False  -------------------------
-    preds = cld_top_tmp < 273.15
+    preds = cld_top_tmp <= 273.15
 
     true_ice = ice & preds
     false_ice = no_ice & preds