diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01.png
index 76df9fc..b5b2a7c 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_I.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_I.png
new file mode 100644
index 0000000..9d13f4a
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_I.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_IQU.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_IQU.png
index ffadb81..f14b120 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_IQU.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_IQU.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_I_err.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_I_err.png
index 433ac7a..0abc9e6 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_I_err.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_I_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P.png
index 9c273eb..b81e7b4 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_err.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_err.png
index d02e5a4..3b05273 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_err.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_flux.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_flux.png
index d175c33..f708510 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_flux.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_P_flux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRi.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRi.png
index 2858269..b10aacc 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRi.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRi.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRp.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRp.png
index f0775be..1c543cb 100644
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRp.png and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM01_SNRp.png differ
diff --git a/src/lib/plots.py b/src/lib/plots.py
index c3ffe4b..db0dc20 100755
--- a/src/lib/plots.py
+++ b/src/lib/plots.py
@@ -43,6 +43,7 @@ from matplotlib.patches import Rectangle, Circle
 from matplotlib.path import Path
 from matplotlib.widgets import RectangleSelector, LassoSelector, Button, Slider, TextBox
 from matplotlib.colors import LogNorm
+from matplotlib.ticker import LogFormatter 
 import matplotlib.font_manager as fm
 from mpl_toolkits.axes_grid1.anchored_artists import AnchoredSizeBar, AnchoredDirectionArrows
 from astropy.wcs import WCS
@@ -359,8 +360,10 @@ def polarization_map(Stokes, data_mask=None, rectangle=None, SNRp_cut=3., SNRi_c
         #ax.clabel(cont,inline=True,fontsize=6)
     else:
         # Defaults to intensity map
-        vmin, vmax = 0., np.max(stkI.data[stkI.data > 0.]*convert_flux*2.)
-        im = ax.imshow(stkI.data*convert_flux, vmin=vmin, vmax=vmax, aspect='equal', cmap='inferno', alpha=1.)
+        vmin, vmax = np.min(stkI.data[SNRi > SNRi_cut]*convert_flux)/10., np.max(stkI.data[SNRi > SNRi_cut]*convert_flux)
+        #im = ax.imshow(stkI.data*convert_flux, vmin=vmin, vmax=vmax, aspect='equal', cmap='inferno', alpha=1.)
+        #cbar = plt.colorbar(im, cax=cbar_ax, label=r"$F_{\lambda}$ [$ergs \cdot cm^{-2} \cdot s^{-1} \cdot \AA$]")
+        im = ax.imshow(stkI.data*convert_flux, norm=LogNorm(vmin,vmax), aspect='equal', cmap='inferno', alpha=1.)
         cbar = plt.colorbar(im, cax=cbar_ax, label=r"$F_{\lambda}$ [$ergs \cdot cm^{-2} \cdot s^{-1} \cdot \AA$]")
 
     #Get integrated values from header
@@ -395,7 +398,7 @@ def polarization_map(Stokes, data_mask=None, rectangle=None, SNRp_cut=3., SNRi_c
         if display.lower() == 'default':
             ax.add_artist(px_sc)
             ax.add_artist(north_dir)
-        ax.annotate(r"$F_{{\lambda}}^{{int}}$({0:.0f} $\AA$) = {1} $ergs \cdot cm^{{-2}} \cdot s^{{-1}} \cdot \AA^{{-1}}$".format(pivot_wav,sci_not(I_diluted*convert_flux,I_diluted_err*convert_flux,2)), color='white', xy=(0.01, 0.92), xycoords='axes fraction')
+        ax.annotate(r"$F_{{\lambda}}^{{int}}$({0:.0f} $\AA$) = {1} $ergs \cdot cm^{{-2}} \cdot s^{{-1}} \cdot \AA^{{-1}}$".format(pivot_wav,sci_not(I_diluted*convert_flux,I_diluted_err*convert_flux,2)), color='white', xy=(0.01, 0.97), xycoords='axes fraction')
 
     # Display instrument FOV
     if not(rectangle is None):
@@ -405,7 +408,7 @@ def polarization_map(Stokes, data_mask=None, rectangle=None, SNRp_cut=3., SNRi_c
             edgecolor=color, fill=False))
 
 
-    ax.coords.grid(True, color='white', ls='dotted', alpha=0.5)
+    #ax.coords.grid(True, color='white', ls='dotted', alpha=0.5)
     ax.coords[0].set_axislabel('Right Ascension (J2000)')
     ax.coords[0].set_axislabel_position('t')
     ax.coords[0].set_ticklabel_position('t')