diff --git a/modules/util/plot.py b/modules/util/plot.py
index e4052f74b325ba29e1d0aaa9e01b7eaa9d50bea6..effecb3fff494576e7a4c4a0132c385287fb7969 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='magenta', linewidths=2.0,
+            ax.scatter(lons, lats, s=5000.0, marker='o', facecolors='none', edgecolors='white', 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='magenta', linewidths=2.0,
+            ax.scatter(lons, lats, s=5000.0, marker='o', facecolors='none', edgecolors='white', linewidths=2.0,
                        transform=ccrs.PlateCarree())
             lons = ice_lons_vld[ice_intensity == -1]
             lats = ice_lats_vld[ice_intensity == -1]
@@ -310,7 +310,7 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite,
 
     title = satellite+', '+clvrx_str_time+'  >50% Probability Icing'
     plt.title(title, loc='left', fontweight='bold', fontsize=24)
-    cbar = plt.colorbar(p_im)
+    cbar = plt.colorbar(p_im, label='Probability')
     cbar.ax.tick_params(labelsize=30)
 
     ImageDirAndName = os.path.join(homedir, imagename+'_'+clvrx_str_time)