diff --git a/modules/util/plot.py b/modules/util/plot.py
index a8210844d7682b298468a03b92aa845d3b950e83..e4052f74b325ba29e1d0aaa9e01b7eaa9d50bea6 100644
--- a/modules/util/plot.py
+++ b/modules/util/plot.py
@@ -295,13 +295,13 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite,
             lats = ice_lats_vld[ice_intensity == 1]
             ax.scatter(lons, lats, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree())
             ax.scatter(lons, lats, s=200.0, marker='o', color='black', transform=ccrs.PlateCarree())
-            ax.scatter(lons, lats, s=5000.0, marker='o', facecolors='none', edgecolors='black', linewidths=2.0,
+            ax.scatter(lons, lats, s=5000.0, marker='o', facecolors='none', edgecolors='magenta', linewidths=2.0,
                        transform=ccrs.PlateCarree())
             lons = ice_lons_vld[ice_intensity > 1]
             lats = ice_lats_vld[ice_intensity > 1]
             ax.scatter(lons, lats, s=500.0, marker='^', color='white', transform=ccrs.PlateCarree())
             ax.scatter(lons, lats, s=200.0, marker='^', color='black', transform=ccrs.PlateCarree())
-            ax.scatter(lons, lats, s=5000.0, marker='o', facecolors='none', edgecolors='black', linewidths=2.0,
+            ax.scatter(lons, lats, s=5000.0, marker='o', facecolors='none', edgecolors='magenta', linewidths=2.0,
                        transform=ccrs.PlateCarree())
             lons = ice_lons_vld[ice_intensity == -1]
             lats = ice_lats_vld[ice_intensity == -1]