diff --git a/package/FOC_reduction.py b/package/FOC_reduction.py
index f890f05..1d4443d 100755
--- a/package/FOC_reduction.py
+++ b/package/FOC_reduction.py
@@ -274,7 +274,7 @@ def main(target=None, proposal_id=None, infiles=None, output_dir="./data", crop=
         flux_data=np.array([flux_data, flux_error, flux_mask]),
         flux_head=flux_head,
     )
-    outfiles.append("/".join([data_folder, Stokes_hdul["I_STOKES"].header["FILENAME"] + ".fits"]))
+    outfiles.append("/".join([data_folder, Stokes_hdul[0].header["FILENAME"] + ".fits"]))
 
     # Step 5:
     # crop to desired region of interest (roi)
diff --git a/package/lib/plots.py b/package/lib/plots.py
index 0b5ab5d..b98beac 100755
--- a/package/lib/plots.py
+++ b/package/lib/plots.py
@@ -2951,7 +2951,7 @@ class pol_map(object):
 
     @property
     def wcs(self):
-        return WCS(self.Stokes["I_STOKES"].header).celestial.deepcopy()
+        return WCS(self.Stokes[0].header).celestial.deepcopy()
 
     @property
     def Flux(self):