diff --git a/modules/util/plot.py b/modules/util/plot.py
index a981fbb45287f26e91a18a4b981e5ae3872d36ba..329d4f2f08abc2d7035a57d365d4cd3683277ca1 100644
--- a/modules/util/plot.py
+++ b/modules/util/plot.py
@@ -14,6 +14,7 @@ import h5py
 from util.util import get_grid_values_all, get_cartopy_crs, homedir
 import csv
 from matplotlib.colors import LinearSegmentedColormap
+from matplotlib.markers import MarkerStyle
 
 from util.setup import home_dir
 import mpl_scatter_density
@@ -255,6 +256,7 @@ def plot_image3(image, bt, cmap='Greys'):
 
 def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite, domain, ice_lons_vld=None, ice_lats_vld=None,
                      ice_intensity=None, imagename='icing', extent=[-110, -60, 10, 55]):
+    mrkr_styl = MarkerStyle('o', fillstyle=None)
     geos, xlen, xmin, xmax, ylen, ymin, ymax = get_cartopy_crs(satellite, domain)
     bg_image = None
     if h5f is not None:
@@ -292,12 +294,12 @@ 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=1000.0, fillstyle=None, marker='o', color='black', transform=ccrs.PlateCarree())
+            ax.scatter(lons, lats, s=1000.0, marker=mrkr_styl, color='black', 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=1000.0, fillstyle=None, marker='o', color='black', transform=ccrs.PlateCarree())
+            ax.scatter(lons, lats, s=1000.0, marker=mrkr_styl, color='black', 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='o', color='white', transform=ccrs.PlateCarree())