diff --git a/package/lib/fits.py b/package/lib/fits.py
index 549b0b5..57ed7e5 100755
--- a/package/lib/fits.py
+++ b/package/lib/fits.py
@@ -89,7 +89,7 @@ def get_obs_data(infiles, data_folder="", compute_flux=False):
 
     # force WCS for POL60 to have same pixel size as POL0 and POL120
     is_pol60 = np.array([head["filtnam1"].lower() == "pol60" for head in headers], dtype=bool)
-    cdelt = np.round(np.array([WCS(head).wcs.cdelt[:2] for head in headers]), 10)
+    cdelt = np.round(np.abs(np.array([WCS(head).wcs.cdelt[:2] for head in headers])), 10)
     if np.unique(cdelt[np.logical_not(is_pol60)], axis=0).size != 2:
         print(np.unique(cdelt[np.logical_not(is_pol60)], axis=0))
         raise ValueError("Not all images have same pixel size")
diff --git a/package/lib/plots.py b/package/lib/plots.py
index 10fa302..0df38ad 100755
--- a/package/lib/plots.py
+++ b/package/lib/plots.py
@@ -189,7 +189,7 @@ def plot_Stokes(Stokes, savename=None, plots_folder=""):
     for dataset in [stkI, stkQ, stkU]:
         dataset[np.logical_not(data_mask)] = np.nan
 
-    wcs = WCS(Stokes[0]).deepcopy()
+    wcs = WCS(Stokes[0]).celestial.deepcopy()
 
     # Plot figure
     plt.rcParams.update({"font.size": 14})
@@ -370,7 +370,9 @@ def polarization_map(
     ax.coords[0].set_axislabel("Right Ascension (J2000)")
     ax.coords[0].set_axislabel_position("t")
     ax.coords[0].set_ticklabel_position("t")
-    ax.set_ylabel("Declination (J2000)", labelpad=-1)
+    ax.coords[1].set_axislabel("Declination (J2000)")
+    ax.coords[1].set_axislabel_position("l")
+    ax.coords[1].set_ticklabel_position("l")
 
     vmin, vmax = 0.0, stkI.max() * convert_flux
     for key, value in [["cmap", [["cmap", "inferno"]]], ["width", [["width", 0.4]]], ["linewidth", [["linewidth", 0.6]]]]: