diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020.png
index c398e0c..264f6bc 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I.png
index 9daa9dd..99697eb 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_IQU.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_IQU.png
index 2c167b9..bd4e459 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_IQU.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I_err.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I_err.png
index c7c4363..1f0f0ac 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I_err.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P.png
index 7de026f..f75cdb4 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_PA.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_PA.png
index 4299998..a934f7e 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_PA.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_err.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_err.png
index d4828f5..e291e78 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_err.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_flux.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_flux.png
index ac6f55f..11316dd 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_flux.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRi.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRi.png
index cc8df07..e9536d9 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRi.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRp.png b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRp.png
index b04cff4..42f8f31 100644
Binary files a/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRp.png and b/plots/3C109_x3mc010/3C109_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_errors_POL0_background_location.png b/plots/3C109_x3mc010/3C109_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..fe7e5b9
Binary files /dev/null and b/plots/3C109_x3mc010/3C109_FOC_errors_POL0_background_location.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_errors_background_flux.png b/plots/3C109_x3mc010/3C109_FOC_errors_background_flux.png
new file mode 100644
index 0000000..7888f05
Binary files /dev/null and b/plots/3C109_x3mc010/3C109_FOC_errors_background_flux.png differ
diff --git a/plots/3C109_x3mc010/3C109_FOC_errors_histograms.png b/plots/3C109_x3mc010/3C109_FOC_errors_histograms.png
new file mode 100644
index 0000000..384eae6
Binary files /dev/null and b/plots/3C109_x3mc010/3C109_FOC_errors_histograms.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020.png
new file mode 100644
index 0000000..a156d83
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_I.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_I.png
new file mode 100644
index 0000000..c8e3976
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_I.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_IQU.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_IQU.png
new file mode 100644
index 0000000..59b6768
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_I_err.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_I_err.png
new file mode 100644
index 0000000..ca0b050
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P.png
new file mode 100644
index 0000000..e711622
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_PA.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_PA.png
new file mode 100644
index 0000000..317705d
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P_err.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P_err.png
new file mode 100644
index 0000000..f94299f
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P_flux.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P_flux.png
new file mode 100644
index 0000000..d974771
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_SNRi.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_SNRi.png
new file mode 100644
index 0000000..73f6eb3
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_SNRp.png b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_SNRp.png
new file mode 100644
index 0000000..c1a77aa
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_errors_POL0_background_location.png b/plots/3C273_x0u20/3C273_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..2f178a2
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_errors_POL0_background_location.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_errors_background_flux.png b/plots/3C273_x0u20/3C273_FOC_errors_background_flux.png
new file mode 100644
index 0000000..6d5cea7
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_errors_background_flux.png differ
diff --git a/plots/3C273_x0u20/3C273_FOC_errors_histograms.png b/plots/3C273_x0u20/3C273_FOC_errors_histograms.png
new file mode 100644
index 0000000..7ee0d91
Binary files /dev/null and b/plots/3C273_x0u20/3C273_FOC_errors_histograms.png differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020.png
deleted file mode 100644
index d39ce58..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_I.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_I.png
deleted file mode 100644
index f037897..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_IQU.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 3a69b47..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_I_err.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index c091dfa..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P.png
deleted file mode 100644
index e19625b..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_PA.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index d4d53c6..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P_err.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index baae0fc..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P_flux.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index e353220..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_SNRi.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 940dca5..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_SNRp.png b/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index b1081b1..0000000
Binary files a/plots/3C273_x0u20/3C273_jet_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020.png
index bf77853..800c47a 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I.png
index 30077fb..7eb83cf 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_IQU.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_IQU.png
index ae47077..d7608ac 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_IQU.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I_err.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I_err.png
index 0ffc3fd..1bd609c 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I_err.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P.png
index 8994e97..2bf11e7 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_PA.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_PA.png
index 3dd2f0f..5132e5a 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_PA.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_err.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_err.png
index 4f56194..def2e12 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_err.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_flux.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_flux.png
index 690b812..80489a1 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_flux.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRi.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRi.png
index f744861..e38235e 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRi.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRp.png b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRp.png
index 1a129e8..d3358c3 100644
Binary files a/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRp.png and b/plots/3C405_x136060/3C405_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_errors_POL0_background_location.png b/plots/3C405_x136060/3C405_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..baffb04
Binary files /dev/null and b/plots/3C405_x136060/3C405_FOC_errors_POL0_background_location.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_errors_background_flux.png b/plots/3C405_x136060/3C405_FOC_errors_background_flux.png
new file mode 100644
index 0000000..22d0f79
Binary files /dev/null and b/plots/3C405_x136060/3C405_FOC_errors_background_flux.png differ
diff --git a/plots/3C405_x136060/3C405_FOC_errors_histograms.png b/plots/3C405_x136060/3C405_FOC_errors_histograms.png
new file mode 100644
index 0000000..3110b77
Binary files /dev/null and b/plots/3C405_x136060/3C405_FOC_errors_histograms.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020.png
deleted file mode 100644
index 3bfe474..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_I.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 0f95c6d..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_IQU.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 165ecfb..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_I_err.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 2157407..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P.png
deleted file mode 100644
index a0ebc9e..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_PA.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 60eb4b5..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P_err.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index f50a84d..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P_flux.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 1c6de57..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_SNRi.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 01b9ee7..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_SNRp.png b/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 1e868a2..0000000
Binary files a/plots/CygnusA_x43w0/CygnusA_F275W_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020.png
index 6b54a45..239293f 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I.png
index 93ea52c..93bc1cb 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_IQU.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_IQU.png
index 14a9038..c49e6e8 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_IQU.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I_err.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I_err.png
index 379e4cf..846f87f 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I_err.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P.png
index 761c349..0bd54d1 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_PA.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_PA.png
index de06c55..f332ec1 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_PA.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_err.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_err.png
index aa47fa0..22417c0 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_err.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_flux.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_flux.png
index 702be1f..46d8c31 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_flux.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRi.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRi.png
index d1ffadb..a08f14b 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRi.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRp.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRp.png
index 49d371a..424f6c2 100644
Binary files a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRp.png and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_POL0_background_location.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..e9dd0ac
Binary files /dev/null and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_POL0_background_location.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_background_flux.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_background_flux.png
new file mode 100644
index 0000000..57448f8
Binary files /dev/null and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_background_flux.png differ
diff --git a/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_histograms.png b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_histograms.png
new file mode 100644
index 0000000..aecb5d4
Binary files /dev/null and b/plots/CygnusA_x43w0/CygnusA_F342W_FOC_errors_histograms.png differ
diff --git a/plots/IC5063_x3nl030/18GHz_overplot.png b/plots/IC5063_x3nl030/18GHz_overplot.png
new file mode 100644
index 0000000..1a4a586
Binary files /dev/null and b/plots/IC5063_x3nl030/18GHz_overplot.png differ
diff --git a/plots/IC5063_x3nl030/IC5063.png b/plots/IC5063_x3nl030/IC5063.png
new file mode 100644
index 0000000..72da45b
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..e957ab7
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_I.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..fc69d8b
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_I.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_IQU.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..fbd696e
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_I_err.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..8318bb3
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..0629b08
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_PA.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..b81564b
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P_err.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..d02b714
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P_flux.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..82860a3
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_SNRi.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..8c23009
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_SNRp.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..7665147
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020.png
index b5f3f9f..8a2a873 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I.png
index 3a77ad4..908cdca 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_IQU.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_IQU.png
index 8f0cdbe..d609d7b 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_IQU.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I_err.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I_err.png
index 75e4488..1008a5a 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I_err.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P.png
index 2ea322d..eae84ff 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_PA.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_PA.png
index 6931e2d..0a4d9f7 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_PA.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_err.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_err.png
index 0405934..52039b1 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_err.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_flux.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_flux.png
index 468659c..e8b316b 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_flux.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRi.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRi.png
index 18f3ec2..7c1aee8 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRi.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRp.png b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRp.png
index 94a6529..c1c5e9c 100644
Binary files a/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRp.png and b/plots/IC5063_x3nl030/IC5063_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_errors_POL0_background_location.png b/plots/IC5063_x3nl030/IC5063_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..ac4df42
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_errors_POL0_background_location.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_errors_background_flux.png b/plots/IC5063_x3nl030/IC5063_FOC_errors_background_flux.png
new file mode 100644
index 0000000..d75751a
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_errors_background_flux.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_errors_background_location.png b/plots/IC5063_x3nl030/IC5063_FOC_errors_background_location.png
new file mode 100644
index 0000000..ff5cf06
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_errors_background_location.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_FOC_errors_histograms.png b/plots/IC5063_x3nl030/IC5063_FOC_errors_histograms.png
new file mode 100644
index 0000000..0d62c2f
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_FOC_errors_histograms.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_polang.png b/plots/IC5063_x3nl030/IC5063_polang.png
new file mode 100644
index 0000000..b34883c
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_polang.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_poldeg.png b/plots/IC5063_x3nl030/IC5063_poldeg.png
new file mode 100644
index 0000000..64ff5ea
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_poldeg.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_polflux.png b/plots/IC5063_x3nl030/IC5063_polflux.png
new file mode 100644
index 0000000..227b5dc
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_polflux.png differ
diff --git a/plots/IC5063_x3nl030/IC5063_totalflux.png b/plots/IC5063_x3nl030/IC5063_totalflux.png
new file mode 100644
index 0000000..ff683a4
Binary files /dev/null and b/plots/IC5063_x3nl030/IC5063_totalflux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..d2f434c
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_I.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..13e8b0e
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_I.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_IQU.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..c122eea
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_I_err.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..6626fcd
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..6c82158
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_PA.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..041f793
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P_err.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..8b24c6f
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P_flux.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..45826eb
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_SNRi.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..41886fd
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_SNRp.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..9f5bd24
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020.png
deleted file mode 100644
index 22e43b5..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_I.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_I.png
deleted file mode 100644
index d610112..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_IQU.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 73126b8..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_I_err.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 387722b..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P.png
deleted file mode 100644
index 8aaf4ec..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_PA.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index d4d715c..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P_err.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index bb06953..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 12fd9d6..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 92cc9cb..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 99f1964..0000000
Binary files a/plots/M87/POS1/M87_POS1_1995_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_errors_POL0_background_location.png b/plots/M87/POS1/M87_POS1_1995_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..370c0ed
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_errors_POL0_background_location.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_errors_background_flux.png b/plots/M87/POS1/M87_POS1_1995_FOC_errors_background_flux.png
new file mode 100644
index 0000000..5fd1236
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_errors_background_flux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1995_FOC_errors_histograms.png b/plots/M87/POS1/M87_POS1_1995_FOC_errors_histograms.png
new file mode 100644
index 0000000..06edfc5
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1995_FOC_errors_histograms.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..ae909fe
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_I.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..38c9045
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_I.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_IQU.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..404ba85
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_I_err.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..e4cab65
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..e101e82
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_PA.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..0f7d908
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P_err.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..5b1d348
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P_flux.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..7b98b3d
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_SNRi.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..15ff0ed
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_SNRp.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..fb0a5b5
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020.png
deleted file mode 100644
index 352c66f..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_I.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_I.png
deleted file mode 100644
index c920688..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_IQU.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 17c5d6c..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_I_err.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 6558ab2..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P.png
deleted file mode 100644
index 79bd91e..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_PA.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index fa4cf4f..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P_err.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 9ff8222..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 016dafb..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 1826fc6..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 76a3804..0000000
Binary files a/plots/M87/POS1/M87_POS1_1996_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..c2bacce
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_I.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..24d2379
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_I.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_IQU.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..e85a18c
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_I_err.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..c6d8137
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..8f7d398
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_PA.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..0278704
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P_err.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..3f39e1f
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P_flux.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..cd04bb9
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_SNRi.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..60819a8
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_SNRp.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..859efc4
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020.png
deleted file mode 100644
index 65104c5..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_I.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 1a341e7..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_IQU.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 495989f..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_I_err.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 00034b5..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P.png
deleted file mode 100644
index ad4f95c..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_PA.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index f1e1a17..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P_err.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 06e42e2..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 3dfdf6b..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index ab1869c..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 7128686..0000000
Binary files a/plots/M87/POS1/M87_POS1_1997_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..9d9c9fe
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_I.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..ed22828
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_I.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_IQU.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..fb287fb
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_I_err.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..aa4c199
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..4e85971
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_PA.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..17b85ba
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P_err.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..07f528f
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P_flux.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..f1cc0df
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_SNRi.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..a192333
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_SNRp.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..db2e86d
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020.png
deleted file mode 100644
index dd698cb..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_I.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 86046e9..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_IQU.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 03dd03c..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_I_err.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index e508b44..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P.png
deleted file mode 100644
index 3285bd0..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_PA.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 5f930ee..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P_err.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index d43e29b..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 90b6b01..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 12348bb..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 980f5da..0000000
Binary files a/plots/M87/POS1/M87_POS1_1998_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..2f040ae
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_I.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..6079e01
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_I.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_IQU.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..28923a5
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_I_err.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..e22ac57
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..38f2e45
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_PA.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..669c7d1
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P_err.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..b511b7d
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P_flux.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..1399b5b
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_SNRi.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..82101df
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_SNRp.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..4022e15
Binary files /dev/null and b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020.png
deleted file mode 100644
index 030c023..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_I.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 51db7c2..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_IQU.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 0b31848..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_I_err.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index ecd6419..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P.png
deleted file mode 100644
index e9358c7..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_PA.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index c7cdc15..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P_err.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 658f0d4..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index d659f83..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 08f743f..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 6c219d3..0000000
Binary files a/plots/M87/POS1/M87_POS1_1999_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS1/animated_loop/M87_020_log_core.gif b/plots/M87/POS1/animated_loop/M87_020_log_core.gif
deleted file mode 100644
index 58a27da..0000000
Binary files a/plots/M87/POS1/animated_loop/M87_020_log_core.gif and /dev/null differ
diff --git a/plots/M87/POS1/animated_loop/M87_020_log_core_0.png b/plots/M87/POS1/animated_loop/M87_020_log_core_0.png
deleted file mode 100644
index bc84953..0000000
Binary files a/plots/M87/POS1/animated_loop/M87_020_log_core_0.png and /dev/null differ
diff --git a/plots/M87/POS1/animated_loop/M87_020_log_core_1.png b/plots/M87/POS1/animated_loop/M87_020_log_core_1.png
deleted file mode 100644
index 657984a..0000000
Binary files a/plots/M87/POS1/animated_loop/M87_020_log_core_1.png and /dev/null differ
diff --git a/plots/M87/POS1/animated_loop/M87_020_log_core_2.png b/plots/M87/POS1/animated_loop/M87_020_log_core_2.png
deleted file mode 100644
index c77417f..0000000
Binary files a/plots/M87/POS1/animated_loop/M87_020_log_core_2.png and /dev/null differ
diff --git a/plots/M87/POS1/animated_loop/M87_020_log_core_3.png b/plots/M87/POS1/animated_loop/M87_020_log_core_3.png
deleted file mode 100644
index 70b22b7..0000000
Binary files a/plots/M87/POS1/animated_loop/M87_020_log_core_3.png and /dev/null differ
diff --git a/plots/M87/POS1/animated_loop/M87_020_log_core_4.png b/plots/M87/POS1/animated_loop/M87_020_log_core_4.png
deleted file mode 100644
index dce49b8..0000000
Binary files a/plots/M87/POS1/animated_loop/M87_020_log_core_4.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020.png
deleted file mode 100644
index d1db95b..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_I.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_I.png
deleted file mode 100644
index f5f6acb..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_IQU.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 75f0255..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_I_err.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 47f07c3..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P.png
deleted file mode 100644
index abd66ab..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_PA.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index d142116..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P_err.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 36e9102..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 46b93b8..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 6ba215a..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index b77f409..0000000
Binary files a/plots/M87/POS3/M87_POS3_1995_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020.png
deleted file mode 100644
index 3a96787..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_I.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_I.png
deleted file mode 100644
index f5296a0..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_IQU.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index a9765e1..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_I_err.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index c228a2e..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P.png
deleted file mode 100644
index 5c626cd..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_PA.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 0f2ce64..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P_err.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 5fd5dc9..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index d1d8796..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index f30c8d0..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index f4f26f7..0000000
Binary files a/plots/M87/POS3/M87_POS3_1996_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020.png
deleted file mode 100644
index 63224ce..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_I.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 4943623..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_IQU.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index ae7c403..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_I_err.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 96b0969..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P.png
deleted file mode 100644
index 942ed0f..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_PA.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 0b50378..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P_err.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index fd0146d..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 7502a05..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 8c7032d..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index da33e94..0000000
Binary files a/plots/M87/POS3/M87_POS3_1997_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020.png
deleted file mode 100644
index e8ac230..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_I.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 4028f12..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_IQU.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index a0c76a6..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_I_err.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 3c65e14..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P.png
deleted file mode 100644
index 83f4d07..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_PA.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index c778eb8..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P_err.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index e10f4ff..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index ab72264..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index b4f75c1..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 4241c4e..0000000
Binary files a/plots/M87/POS3/M87_POS3_1998_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020.png
deleted file mode 100644
index 4bb1fd1..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_I.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 8e41446..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_IQU.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 628ed0c..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_I_err.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 8a99453..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P.png
deleted file mode 100644
index 421c0ef..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_PA.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 7f3dbf3..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P_err.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 53077e9..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P_flux.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 77bc355..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_SNRi.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 997ffac..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_SNRp.png b/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 39225da..0000000
Binary files a/plots/M87/POS3/M87_POS3_1999_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/M87/POS3/animated_loop/M87_020_log.gif b/plots/M87/POS3/animated_loop/M87_020_log.gif
deleted file mode 100644
index 31dd440..0000000
Binary files a/plots/M87/POS3/animated_loop/M87_020_log.gif and /dev/null differ
diff --git a/plots/M87/POS3/animated_loop/M87_020_log_0.png b/plots/M87/POS3/animated_loop/M87_020_log_0.png
deleted file mode 100644
index b347d3e..0000000
Binary files a/plots/M87/POS3/animated_loop/M87_020_log_0.png and /dev/null differ
diff --git a/plots/M87/POS3/animated_loop/M87_020_log_1.png b/plots/M87/POS3/animated_loop/M87_020_log_1.png
deleted file mode 100644
index 058ddff..0000000
Binary files a/plots/M87/POS3/animated_loop/M87_020_log_1.png and /dev/null differ
diff --git a/plots/M87/POS3/animated_loop/M87_020_log_2.png b/plots/M87/POS3/animated_loop/M87_020_log_2.png
deleted file mode 100644
index bc0361b..0000000
Binary files a/plots/M87/POS3/animated_loop/M87_020_log_2.png and /dev/null differ
diff --git a/plots/M87/POS3/animated_loop/M87_020_log_3.png b/plots/M87/POS3/animated_loop/M87_020_log_3.png
deleted file mode 100644
index 0b80823..0000000
Binary files a/plots/M87/POS3/animated_loop/M87_020_log_3.png and /dev/null differ
diff --git a/plots/M87/POS3/animated_loop/M87_020_log_4.png b/plots/M87/POS3/animated_loop/M87_020_log_4.png
deleted file mode 100644
index 6e1a5ec..0000000
Binary files a/plots/M87/POS3/animated_loop/M87_020_log_4.png and /dev/null differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..fccc445
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_I.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..fd166ff
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_I.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_IQU.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..c36e6b3
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_I_err.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..cf2ddc1
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..8a55756
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_PA.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..9c8e6f4
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P_err.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..6347e8b
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P_flux.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..2e4888f
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_SNRi.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..262872e
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_SNRp.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..7c657fd
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020.png
index 833fb7f..6affaa3 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I.png
index 3ff806c..1c04847 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_IQU.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_IQU.png
index 4aa23a1..d52c5e7 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_IQU.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I_err.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I_err.png
index b50cd19..776180c 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I_err.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P.png
index e5e3ddc..112390e 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_PA.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_PA.png
index f597654..6ab4948 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_PA.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_err.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_err.png
index a0d7d39..66c76e8 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_err.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_flux.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_flux.png
index 857cc5b..f86d856 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_flux.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRi.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRi.png
index 4418f74..f858e54 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRi.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRp.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRp.png
index 98164f6..3bc66a0 100644
Binary files a/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRp.png and b/plots/MKN3_x3md010/MKN3_F275W_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_POL0_background_location.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..07f0234
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_POL0_background_location.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_background_flux.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_background_flux.png
new file mode 100644
index 0000000..bc3ad86
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_background_flux.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_histograms.png b/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_histograms.png
new file mode 100644
index 0000000..00a10d3
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F275W_FOC_errors_histograms.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..9d2a5d0
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_I.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..184a3b9
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_I.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_IQU.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..8fca642
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_I_err.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..205860f
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..b794d4a
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_PA.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..a3e51f6
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P_err.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..1cf097d
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P_flux.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..9168d8f
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_SNRi.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..24c90c2
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_SNRp.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..14af637
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020.png
index 7a69ea2..c94db92 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I.png
index 8e52501..4c23073 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_IQU.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_IQU.png
index 66bfc9a..dd8789c 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_IQU.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I_err.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I_err.png
index cddf905..fecd649 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I_err.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P.png
index ae31b54..cff525d 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_PA.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_PA.png
index 6230ce3..690c148 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_PA.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_err.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_err.png
index d8bdd47..a585863 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_err.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_flux.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_flux.png
index 75cd872..4e7e051 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_flux.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRi.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRi.png
index 0449c14..c0417f7 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRi.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRp.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRp.png
index 172676a..c78baf1 100644
Binary files a/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRp.png and b/plots/MKN3_x3md010/MKN3_F342W_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_POL120_background_location.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_POL120_background_location.png
new file mode 100644
index 0000000..66a47eb
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_POL120_background_location.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_background_flux.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_background_flux.png
new file mode 100644
index 0000000..e47e23c
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_background_flux.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_background_location.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_background_location.png
new file mode 100644
index 0000000..78bb051
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_background_location.png differ
diff --git a/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_histograms.png b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_histograms.png
new file mode 100644
index 0000000..3889114
Binary files /dev/null and b/plots/MKN3_x3md010/MKN3_F342W_FOC_errors_histograms.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..cfed5e4
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_I.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..1548e75
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_I.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_IQU.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..aca3b2d
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_I_err.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..717b868
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..3e48f07
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_PA.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..39dceea
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P_err.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..fc87c7e
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P_flux.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..6311ee5
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_SNRi.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..b727d99
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_SNRp.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..27ff070
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020.png
index ee9edeb..8dc11fb 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I.png
index 13be95b..fb14511 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_IQU.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_IQU.png
index c4a543a..200c6b4 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_IQU.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I_err.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I_err.png
index e3efb25..d1a4b21 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I_err.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P.png
index 5668bb7..3a5e31f 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_PA.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_PA.png
index 679d725..e690b31 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_PA.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_err.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_err.png
index 213529c..eefc9bf 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_err.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_flux.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_flux.png
index 342e302..7f85408 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_flux.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRi.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRi.png
index 920f067..32bb5e4 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRi.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRp.png b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRp.png
index 7430c78..c5eefef 100644
Binary files a/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRp.png and b/plots/MKN3_x3nl010/MKN3_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_errors_POL0_background_location.png b/plots/MKN3_x3nl010/MKN3_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..c74dd5f
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_errors_POL0_background_location.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_errors_background_flux.png b/plots/MKN3_x3nl010/MKN3_FOC_errors_background_flux.png
new file mode 100644
index 0000000..0db6b4c
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_errors_background_flux.png differ
diff --git a/plots/MKN3_x3nl010/MKN3_FOC_errors_histograms.png b/plots/MKN3_x3nl010/MKN3_FOC_errors_histograms.png
new file mode 100644
index 0000000..0f5b7de
Binary files /dev/null and b/plots/MKN3_x3nl010/MKN3_FOC_errors_histograms.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..f241f44
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_I.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..41342f1
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_I.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_IQU.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..d0d107d
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_I_err.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..372fc9e
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..6c627f2
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_PA.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..b71dc6a
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P_err.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..53e927a
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P_flux.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..65fcc9f
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_SNRi.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..6862100
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_SNRp.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..455da18
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020.png
deleted file mode 100644
index 3cf858d..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_I.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 419d292..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_IQU.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 68d5e57..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_I_err.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 6eeff46..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P.png
deleted file mode 100644
index c447ec2..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_PA.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 15ca10a..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P_err.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 884a392..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P_flux.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index d6bf62c..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_SNRi.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 97dc2ae..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_SNRp.png b/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 9f62156..0000000
Binary files a/plots/MKN463_x2rp030/MKN463_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_errors_POL0_background_location.png b/plots/MKN463_x2rp030/MKN463_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..7c76d09
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_errors_POL0_background_location.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_errors_background_flux.png b/plots/MKN463_x2rp030/MKN463_FOC_errors_background_flux.png
new file mode 100644
index 0000000..074fb36
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_errors_background_flux.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_errors_background_location.png b/plots/MKN463_x2rp030/MKN463_FOC_errors_background_location.png
new file mode 100644
index 0000000..a35ae65
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_errors_background_location.png differ
diff --git a/plots/MKN463_x2rp030/MKN463_FOC_errors_histograms.png b/plots/MKN463_x2rp030/MKN463_FOC_errors_histograms.png
new file mode 100644
index 0000000..002f0ef
Binary files /dev/null and b/plots/MKN463_x2rp030/MKN463_FOC_errors_histograms.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..30dbdd9
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_I.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..47aa8d3
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_I.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_IQU.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..8666088
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_I_err.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..4367f95
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..3085e7a
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_PA.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..20f2f1b
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P_err.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..6f77247
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P_flux.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..b26edaa
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_SNRi.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..757c232
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_SNRp.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..82f6241
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020.png
deleted file mode 100644
index 5529a75..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_I.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 6de9019..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_IQU.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 019ae27..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_I_err.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index 16659e4..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P.png
deleted file mode 100644
index d5a2850..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_PA.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 4facb0d..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P_err.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 1698f88..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P_flux.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index f20ea56..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_SNRi.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 83af244..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_SNRp.png b/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 5094277..0000000
Binary files a/plots/MKN78_x3nl020/MKN78_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_errors_POL0_background_location.png b/plots/MKN78_x3nl020/MKN78_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..948d8a9
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_errors_POL0_background_location.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_errors_background_flux.png b/plots/MKN78_x3nl020/MKN78_FOC_errors_background_flux.png
new file mode 100644
index 0000000..9b1b055
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_errors_background_flux.png differ
diff --git a/plots/MKN78_x3nl020/MKN78_FOC_errors_histograms.png b/plots/MKN78_x3nl020/MKN78_FOC_errors_histograms.png
new file mode 100644
index 0000000..2552e88
Binary files /dev/null and b/plots/MKN78_x3nl020/MKN78_FOC_errors_histograms.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010.png
new file mode 100644
index 0000000..06fece9
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_I.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_I.png
new file mode 100644
index 0000000..ad6148d
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_I.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_IQU.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_IQU.png
new file mode 100644
index 0000000..9b9cb2a
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_I_err.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_I_err.png
new file mode 100644
index 0000000..1b884b9
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P.png
new file mode 100644
index 0000000..4b5838e
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_PA.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_PA.png
new file mode 100644
index 0000000..c7676e1
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P_err.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P_err.png
new file mode 100644
index 0000000..f135e5c
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P_flux.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P_flux.png
new file mode 100644
index 0000000..4b97be7
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_SNRi.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_SNRi.png
new file mode 100644
index 0000000..9fb9716
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_SNRp.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_SNRp.png
new file mode 100644
index 0000000..32d9658
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020.png
deleted file mode 100644
index ce73fac..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_I.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_I.png
deleted file mode 100644
index 4f8a010..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_I.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_IQU.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_IQU.png
deleted file mode 100644
index 9bcfa4f..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_IQU.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_I_err.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_I_err.png
deleted file mode 100644
index cdfa4da..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_I_err.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P.png
deleted file mode 100644
index e299da3..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_PA.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_PA.png
deleted file mode 100644
index 465ba26..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_PA.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P_err.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P_err.png
deleted file mode 100644
index 6b718f8..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P_err.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P_flux.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P_flux.png
deleted file mode 100644
index 20f1e7d..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_P_flux.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_SNRi.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_SNRi.png
deleted file mode 100644
index 184d229..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_SNRi.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_SNRp.png b/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_SNRp.png
deleted file mode 100644
index 3aa182c..0000000
Binary files a/plots/NGC1068_x14w010/NGC1068_FOC_combine_FWHM020_SNRp.png and /dev/null differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_errors_POL0_background_location.png b/plots/NGC1068_x14w010/NGC1068_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..84f083e
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_errors_POL0_background_location.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_errors_background_flux.png b/plots/NGC1068_x14w010/NGC1068_FOC_errors_background_flux.png
new file mode 100644
index 0000000..f9a38fc
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_errors_background_flux.png differ
diff --git a/plots/NGC1068_x14w010/NGC1068_FOC_errors_histograms.png b/plots/NGC1068_x14w010/NGC1068_FOC_errors_histograms.png
new file mode 100644
index 0000000..ac3dc7d
Binary files /dev/null and b/plots/NGC1068_x14w010/NGC1068_FOC_errors_histograms.png differ
diff --git a/plots/NGC1068_x274020/NGC1068.png b/plots/NGC1068_x274020/NGC1068.png
index 791565b..4c5a774 100644
Binary files a/plots/NGC1068_x274020/NGC1068.png and b/plots/NGC1068_x274020/NGC1068.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_POL0_crop_region.png b/plots/NGC1068_x274020/NGC1068_FOC_POL0_crop_region.png
index e9558b3..4382a54 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_POL0_crop_region.png and b/plots/NGC1068_x274020/NGC1068_FOC_POL0_crop_region.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_center_image.png b/plots/NGC1068_x274020/NGC1068_FOC_center_image.png
index dd4279d..2e92d50 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_center_image.png and b/plots/NGC1068_x274020/NGC1068_FOC_center_image.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010.png
index e79c567..08e7b20 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I.png
index 69fe46b..06bc0e4 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_IQU.png
index dd7c477..51fdca7 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_IQU.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_IQU.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I_err.png
index 42d9495..059925b 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_I_err.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P.png
index a185435..30b4550 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_PA.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_PA.png
index b6070f9..7c16e93 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_PA.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_PA.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_err.png
index a2aa310..e76024c 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_err.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_flux.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_flux.png
index bf468f7..dcb0abf 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_flux.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_P_flux.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRi.png
index 3d35e85..918dd64 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRi.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRp.png
index 6638190..e1324d2 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM010_SNRp.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020.png
index aa4b93c..066e594 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I.png
index 6290d67..8911008 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_IQU.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_IQU.png
index 3000e93..7621233 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_IQU.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_IQU.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I_err.png
index 1a4cccc..4b4b302 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_I_err.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P.png
index 2e9ce19..c221167 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_PA.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_PA.png
index 3c29ccc..e4f68f3 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_PA.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_PA.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_err.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_err.png
index 65d22b8..c051f14 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_err.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_err.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_flux.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_flux.png
index 37d822f..590ae8b 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_flux.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_P_flux.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRi.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRi.png
index b44e597..c59d9ec 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRi.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRi.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRp.png b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRp.png
index ab95ad1..546a110 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRp.png and b/plots/NGC1068_x274020/NGC1068_FOC_combine_FWHM020_SNRp.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_crop_region.png b/plots/NGC1068_x274020/NGC1068_FOC_crop_region.png
index b0514ef..cec3dd3 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_crop_region.png and b/plots/NGC1068_x274020/NGC1068_FOC_crop_region.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_errors_POL0_background_location.png b/plots/NGC1068_x274020/NGC1068_FOC_errors_POL0_background_location.png
index 21cb987..6b47622 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_errors_POL0_background_location.png and b/plots/NGC1068_x274020/NGC1068_FOC_errors_POL0_background_location.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_errors_background_flux.png b/plots/NGC1068_x274020/NGC1068_FOC_errors_background_flux.png
index 32e79c2..2002723 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_errors_background_flux.png and b/plots/NGC1068_x274020/NGC1068_FOC_errors_background_flux.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_errors_background_location.png b/plots/NGC1068_x274020/NGC1068_FOC_errors_background_location.png
index 174a7fd..235c214 100644
Binary files a/plots/NGC1068_x274020/NGC1068_FOC_errors_background_location.png and b/plots/NGC1068_x274020/NGC1068_FOC_errors_background_location.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_FOC_errors_histograms.png b/plots/NGC1068_x274020/NGC1068_FOC_errors_histograms.png
new file mode 100644
index 0000000..003b43c
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_FOC_errors_histograms.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px.png
new file mode 100644
index 0000000..fb501ab
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_I.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_I.png
new file mode 100644
index 0000000..a8a79d0
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_I.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_IQU.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_IQU.png
new file mode 100644
index 0000000..4f663a0
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_IQU.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_I_err.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_I_err.png
new file mode 100644
index 0000000..2a925fc
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_I_err.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P.png
new file mode 100644
index 0000000..0d834aa
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_PA.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_PA.png
new file mode 100644
index 0000000..6c002d5
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_PA.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P_err.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P_err.png
new file mode 100644
index 0000000..b269ba5
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P_err.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P_flux.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P_flux.png
new file mode 100644
index 0000000..4e2ae83
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_P_flux.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_SNRi.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_SNRi.png
new file mode 100644
index 0000000..a23816b
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_SNRi.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_SNRp.png b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_SNRp.png
new file mode 100644
index 0000000..8c6117d
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_bin10px_SNRp.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_errors_POL0_background_location.png b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_POL0_background_location.png
new file mode 100644
index 0000000..6b47622
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_POL0_background_location.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_errors_background_flux.png b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_background_flux.png
new file mode 100644
index 0000000..2002723
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_background_flux.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_errors_background_location.png b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_background_location.png
new file mode 100644
index 0000000..e1e7a3f
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_background_location.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_FOC_errors_histograms.png b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_histograms.png
new file mode 100644
index 0000000..003b43c
Binary files /dev/null and b/plots/NGC1068_x274020/NGC1068_K_FOC_errors_histograms.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_K_comparison.png b/plots/NGC1068_x274020/NGC1068_K_comparison.png
index 092e763..c4bea69 100644
Binary files a/plots/NGC1068_x274020/NGC1068_K_comparison.png and b/plots/NGC1068_x274020/NGC1068_K_comparison.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_polang.png b/plots/NGC1068_x274020/NGC1068_polang.png
index 1b96784..e540509 100644
Binary files a/plots/NGC1068_x274020/NGC1068_polang.png and b/plots/NGC1068_x274020/NGC1068_polang.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_poldeg.png b/plots/NGC1068_x274020/NGC1068_poldeg.png
index 7f78e36..f73b67b 100644
Binary files a/plots/NGC1068_x274020/NGC1068_poldeg.png and b/plots/NGC1068_x274020/NGC1068_poldeg.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_polflux.png b/plots/NGC1068_x274020/NGC1068_polflux.png
index 1f33da6..5de6f8b 100644
Binary files a/plots/NGC1068_x274020/NGC1068_polflux.png and b/plots/NGC1068_x274020/NGC1068_polflux.png differ
diff --git a/plots/NGC1068_x274020/NGC1068_totalflux.png b/plots/NGC1068_x274020/NGC1068_totalflux.png
index 64cbd95..80272ad 100644
Binary files a/plots/NGC1068_x274020/NGC1068_totalflux.png and b/plots/NGC1068_x274020/NGC1068_totalflux.png differ
diff --git a/src/FOC_reduction.py b/src/FOC_reduction.py
index 6939b91..5871f37 100755
--- a/src/FOC_reduction.py
+++ b/src/FOC_reduction.py
@@ -40,11 +40,11 @@ globals()['plots_folder'] = "../plots/NGC1068_x274020/"
 #globals()['plots_folder'] = "../plots/3C405_x136060/"
 
 #globals()['data_folder'] = "../data/CygnusA_x43w0/"
-##globals()['infiles'] = ['x43w0101r_c0f.fits', 'x43w0102r_c0f.fits', 'x43w0103r_c0f.fits',
-##   'x43w0104r_c0f.fits', 'x43w0105r_c0f.fits', 'x43w0106r_c0f.fits',
-##   'x43w0107r_c0f.fits', 'x43w0108r_c0f.fits', 'x43w0109r_c0f.fits']
-#globals()['infiles'] = ['x43w0201r_c0f.fits', 'x43w0202r_c0f.fits', 'x43w0203r_c0f.fits',
-#   'x43w0204r_c0f.fits', 'x43w0205r_c0f.fits', 'x43w0206r_c0f.fits']
+#globals()['infiles'] = ['x43w0101r_c0f.fits', 'x43w0102r_c0f.fits', 'x43w0103r_c0f.fits',
+#   'x43w0104r_c0f.fits', 'x43w0105r_c0f.fits', 'x43w0106r_c0f.fits',
+#   'x43w0107r_c0f.fits', 'x43w0108r_c0f.fits', 'x43w0109r_c0f.fits'] #F342W
+##globals()['infiles'] = ['x43w0201r_c0f.fits', 'x43w0202r_c0f.fits', 'x43w0203r_c0f.fits',
+##   'x43w0204r_c0f.fits', 'x43w0205r_c0f.fits', 'x43w0206r_c0f.fits'] #F275W
 #globals()['plots_folder'] = "../plots/CygnusA_x43w0/"
 
 #globals()['data_folder'] = "../data/3C109_x3mc010/"
@@ -72,8 +72,8 @@ globals()['plots_folder'] = "../plots/NGC1068_x274020/"
 #globals()['plots_folder'] = "../plots/MKN3_x3nl010/"
 
 #globals()['data_folder'] = "../data/MKN3_x3md010/"
-##globals()['infiles'] = ['x3md0101r_c0f.fits', 'x3md0102r_c0f.fits', 'x3md0103r_c0f.fits']
-#globals()['infiles'] = ['x3md0104r_c0f.fits', 'x3md0105r_c0f.fits', 'x3md0106r_c0f.fits']
+#globals()['infiles'] = ['x3md0101r_c0f.fits', 'x3md0102r_c0f.fits', 'x3md0103r_c0f.fits'] #F275W
+##globals()['infiles'] = ['x3md0104r_c0f.fits', 'x3md0105r_c0f.fits', 'x3md0106r_c0f.fits'] #F342W
 #globals()['plots_folder'] = "../plots/MKN3_x3md010/"
 
 #globals()['data_folder'] = "../data/MKN78_x3nl020/"
@@ -91,19 +91,19 @@ globals()['plots_folder'] = "../plots/NGC1068_x274020/"
 
 #BEWARE: 5 observations separated by 1 year each (1995, 1996, 1997, 1998, 1999)
 #globals()['data_folder'] = "../data/M87/POS1/"
-#globals()['infiles'] = ['x2py010ct_c0f.fits','x2py010dt_c0f.fits','x2py010et_c0f.fits','x2py010ft_c0f.fits'] #1995
-#globals()['infiles'] = ['x3be010ct_c0f.fits','x3be010dt_c0f.fits','x3be010et_c0f.fits','x3be010ft_c0f.fits'] #1996
-#globals()['infiles'] = ['x43r010km_c0f.fits','x43r010mm_c0f.fits','x43r010om_c0f.fits','x43r010rm_c0f.fits'] #1997
-#globals()['infiles'] = ['x43r110kr_c0f.fits','x43r110mr_c0f.fits','x43r110or_c0f.fits','x43r110rr_c0f.fits'] #1998
+##globals()['infiles'] = ['x2py010ct_c0f.fits','x2py010dt_c0f.fits','x2py010et_c0f.fits','x2py010ft_c0f.fits'] #1995
+##globals()['infiles'] = ['x3be010ct_c0f.fits','x3be010dt_c0f.fits','x3be010et_c0f.fits','x3be010ft_c0f.fits'] #1996
+##globals()['infiles'] = ['x43r010km_c0f.fits','x43r010mm_c0f.fits','x43r010om_c0f.fits','x43r010rm_c0f.fits'] #1997
+##globals()['infiles'] = ['x43r110kr_c0f.fits','x43r110mr_c0f.fits','x43r110or_c0f.fits','x43r110rr_c0f.fits'] #1998
 #globals()['infiles'] = ['x43r210kr_c0f.fits','x43r210mr_c0f.fits','x43r210or_c0f.fits','x43r210rr_c0f.fits'] #1999
 #globals()['plots_folder'] = "../plots/M87/POS1/"
 
 #BEWARE: 5 observations separated by 1 year each (1995, 1996, 1997, 1998, 1999)
 #globals()['data_folder'] = "../data/M87/POS3/"
-#globals()['infiles'] = ['x2py030at_c0f.fits','x2py030bt_c0f.fits','x2py030ct_c0f.fits','x2py0309t_c0f.fits'] #1995
-#globals()['infiles'] = ['x3be030at_c0f.fits','x3be030bt_c0f.fits','x3be030ct_c0f.fits','x3be0309t_c0f.fits'] #1996
-#globals()['infiles'] = ['x43r030em_c0f.fits','x43r030gm_c0f.fits','x43r030im_c0f.fits','x43r030lm_c0f.fits'] #1997
-#globals()['infiles'] = ['x43r130er_c0f.fits','x43r130fr_c0f.fits','x43r130ir_c0f.fits','x43r130lr_c0f.fits'] #1998
+##globals()['infiles'] = ['x2py030at_c0f.fits','x2py030bt_c0f.fits','x2py030ct_c0f.fits','x2py0309t_c0f.fits'] #1995
+##globals()['infiles'] = ['x3be030at_c0f.fits','x3be030bt_c0f.fits','x3be030ct_c0f.fits','x3be0309t_c0f.fits'] #1996
+##globals()['infiles'] = ['x43r030em_c0f.fits','x43r030gm_c0f.fits','x43r030im_c0f.fits','x43r030lm_c0f.fits'] #1997
+##globals()['infiles'] = ['x43r130er_c0f.fits','x43r130fr_c0f.fits','x43r130ir_c0f.fits','x43r130lr_c0f.fits'] #1998
 #globals()['infiles'] = ['x43r230er_c0f.fits','x43r230fr_c0f.fits','x43r230ir_c0f.fits','x43r230lr_c0f.fits'] #1999
 #globals()['plots_folder'] = "../plots/M87/POS3/"
 
@@ -123,20 +123,20 @@ def main():
     # Initial crop
     display_crop = False
     # Error estimation
-    error_sub_shape = (10,10)
+    error_sub_shape = (80,80)
     display_error = False
     # Data binning
     rebin = True
     if rebin:
-        pxsize = 10
-        px_scale = 'pixel'         #pixel, arcsec or full
+        pxsize = 0.05
+        px_scale = 'arcsec'         #pixel, arcsec or full
         rebin_operation = 'sum'     #sum or average
     # Alignement
     align_center = 'image'          #If None will align image to image center
     display_data = False
     # Smoothing
     smoothing_function = 'combine'  #gaussian_after, weighted_gaussian_after, gaussian, weighted_gaussian or combine
-    smoothing_FWHM = None           #If None, no smoothing is done
+    smoothing_FWHM = 0.10           #If None, no smoothing is done
     smoothing_scale = 'arcsec'      #pixel or arcsec
     # Rotation
     rotate_stokes = True            #rotation to North convention can give erroneous results
@@ -146,7 +146,7 @@ def main():
     final_display = False
     # Polarization map output
     figname = 'NGC1068_FOC'         #target/intrument name
-    figtype = '_bin10px'    #additionnal informations
+    figtype = '_combine_FWHM010'    #additionnal informations
     SNRp_cut = 5.    #P measurments with SNR>3
     SNRi_cut = 50.   #I measurments with SNR>30, which implies an uncertainty in P of 4.7%.
     step_vec = 1    #plot all vectors in the array. if step_vec = 2, then every other vector will be plotted
@@ -163,15 +163,13 @@ def main():
     # Deconvolve data using Richardson-Lucy iterative algorithm with a gaussian PSF of given FWHM.
     if deconvolve:
         data_array = proj_red.deconvolve_array(data_array, headers, psf=psf, FWHM=psf_FWHM, scale=psf_scale, shape=psf_shape, iterations=iterations, algo=algo)
-    # Rotate data to have North up
-    if rotate_data:
-        data_mask = np.ones(data_array.shape[1:]).astype(bool)
-        alpha = headers[0]['orientat']
-        mrot = np.array([[np.cos(-alpha), -np.sin(-alpha)], [np.sin(-alpha), np.cos(-alpha)]])
-        data_array, error_array, data_mask, headers = proj_red.rotate_data(data_array, error_array, data_mask, headers, -alpha)
+    # Estimate error from data background, estimated from sub-image of desired sub_shape.
+    if px_scale.lower() not in ['full','integrate']:
+        #data_array, error_array, headers, background = proj_red.get_error(data_array, headers, error_array, sub_shape=error_sub_shape, display=display_error, savename=figname+"_errors", plots_folder=plots_folder, return_background=True)
+        data_array, error_array, headers, background = proj_red.get_error2(data_array, headers, error_array, display=display_error, savename=figname+"_errors", plots_folder=plots_folder, return_background=True)
 
     # Align and rescale images with oversampling.
-    data_array, error_array, headers, data_mask = proj_red.align_data(data_array, headers, error_array=error_array, upsample_factor=10, ref_center=align_center, return_shifts=False)
+    data_array, error_array, headers, data_mask = proj_red.align_data(data_array, headers, error_array=error_array, background=background, upsample_factor=10, ref_center=align_center, return_shifts=False)
 
     # Rebin data to desired pixel size.
     if rebin:
@@ -179,9 +177,11 @@ def main():
             data_array, error_array, headers = proj_red.get_error(data_array, headers, error_array, data_mask, sub_shape=error_sub_shape, display=display_error, savename=figname+"_errors", plots_folder=plots_folder)
         data_array, error_array, headers, Dxy, data_mask = proj_red.rebin_array(data_array, error_array, headers, pxsize=pxsize, scale=px_scale, operation=rebin_operation, data_mask=data_mask)
 
-    # Estimate error from data background, estimated from sub-image of desired sub_shape.
-    if px_scale.lower() not in ['full','integrate']:
-        data_array, error_array, headers = proj_red.get_error(data_array, headers, error_array, data_mask, sub_shape=error_sub_shape, display=display_error, savename=figname+"_errors", plots_folder=plots_folder)
+    # Rotate data to have North up
+    if rotate_data:
+        data_mask = np.ones(data_array.shape[1:]).astype(bool)
+        alpha = headers[0]['orientat']
+        data_array, error_array, data_mask, headers = proj_red.rotate_data(data_array, error_array, data_mask, headers, -alpha)
 
     #Plot array for checking output
     if display_data and px_scale.lower() not in ['full','integrate']:
diff --git a/src/comparison_Kishimoto.py b/src/comparison_Kishimoto.py
index 87bd47c..d408784 100755
--- a/src/comparison_Kishimoto.py
+++ b/src/comparison_Kishimoto.py
@@ -44,8 +44,8 @@ for d in data_K:
 for d in [data_S, data_K]:
     for i in d:
         d[i][np.isnan(d[i])] = 0.
-    d['P'] = np.where(d['I']>0.,np.sqrt(d['Q']**2+d['U']**2)/d['I'],0.)
-    d['sP'] = np.where(d['I']>0.,np.sqrt((d['Q']**2*d['sQ']**2+d['U']**2*d['sU']**2)/(d['Q']**2+d['U']**2)+((d['Q']/d['I'])**2+(d['U']/d['I'])**2)*d['sI']**2)/d['I'],0.)
+    d['P'] = np.where(np.logical_and(np.isfinite(d['I']),d['I']>0.),np.sqrt(d['Q']**2+d['U']**2)/d['I'],0.)
+    d['sP'] = np.where(np.logical_and(np.isfinite(d['I']),d['I']>0.),np.sqrt((d['Q']**2*d['sQ']**2+d['U']**2*d['sU']**2)/(d['Q']**2+d['U']**2)+((d['Q']/d['I'])**2+(d['U']/d['I'])**2)*d['sI']**2)/d['I'],0.)
     d['PA'] = princ_angle((90./np.pi)*np.arctan2(d['U'],d['Q'])+180.)
     d['SNRp'] = np.zeros(d['P'].shape)
     d['SNRp'][d['sP']>0.] = d['P'][d['sP']>0.]/d['sP'][d['sP']>0.]
diff --git a/src/lib/plots.py b/src/lib/plots.py
index e19b8fe..08575fc 100755
--- a/src/lib/plots.py
+++ b/src/lib/plots.py
@@ -337,7 +337,7 @@ def polarization_map(Stokes, data_mask=None, rectangle=None, SNRp_cut=3., SNRi_c
         # Display polarization degree map
         display='pa'
         vmin, vmax = 0., 360.
-        im = ax.imshow(pang.data, vmin=vmin, vmax=vmax, aspect='equal', cmap='inferno', alpha=1.)
+        im = ax.imshow(princ_angle(pang.data), vmin=vmin, vmax=vmax, aspect='equal', cmap='inferno', alpha=1.)
         cbar = plt.colorbar(im, cax=cbar_ax, label=r"$\theta_P$ [°]")
     elif display.lower() in ['s_p','pol_err','pol_deg_err']:
         # Display polarization degree error map
@@ -1743,7 +1743,7 @@ class pol_map(object):
             vmin, vmax = 0., np.max(self.data[self.data > 0.])
             label = r"$P$ [%]"
         elif self.display_selection.lower() in ['pol_ang']:
-            self.data = self.PA
+            self.data = princ_angle(self.PA)
             vmin, vmax = 0, 360.
             label = r"$\theta_{P}$ [°]"
         elif self.display_selection.lower() in ['snri']:
@@ -1833,7 +1833,7 @@ class pol_map(object):
             P_cut = np.sqrt(Q_cut**2+U_cut**2)/I_cut
             P_cut_err = np.sqrt((Q_cut**2*Q_cut_err**2 + U_cut**2*U_cut_err**2 + 2.*Q_cut*U_cut*QU_cut_err)/(Q_cut**2 + U_cut**2) + ((Q_cut/I_cut)**2 + (U_cut/I_cut)**2)*I_cut_err**2 - 2.*(Q_cut/I_cut)*IQ_cut_err - 2.*(U_cut/I_cut)*IU_cut_err)/I_cut
 
-            PA_cut = princ_angle(np.degrees((1./2.)*np.arctan2(U_cut,Q_cut))+180.)
+            PA_cut = 360.-princ_angle(np.degrees((1./2.)*np.arctan2(U_cut,Q_cut)))
             PA_cut_err = princ_angle(np.degrees((1./(2.*(Q_cut**2+U_cut**2)))*np.sqrt(U_cut**2*Q_cut_err**2 + Q_cut**2*U_cut_err**2 - 2.*Q_cut*U_cut*QU_cut_err)))
 
         else:
@@ -1858,7 +1858,7 @@ class pol_map(object):
             P_reg = np.sqrt(Q_reg**2+U_reg**2)/I_reg
             P_reg_err = np.sqrt((Q_reg**2*Q_reg_err**2 + U_reg**2*U_reg_err**2 + 2.*Q_reg*U_reg*QU_reg_err)/(Q_reg**2 + U_reg**2) + ((Q_reg/I_reg)**2 + (U_reg/I_reg)**2)*I_reg_err**2 - 2.*(Q_reg/I_reg)*IQ_reg_err - 2.*(U_reg/I_reg)*IU_reg_err)/I_reg
 
-            PA_reg = princ_angle((90./np.pi)*np.arctan2(U_reg,Q_reg)+180.)
+            PA_reg = 360.-princ_angle((90./np.pi)*np.arctan2(U_reg,Q_reg))
             PA_reg_err = (90./(np.pi*(Q_reg**2+U_reg**2)))*np.sqrt(U_reg**2*Q_reg_err**2 + Q_reg**2*U_reg_err**2 - 2.*Q_reg*U_reg*QU_reg_err)
 
             new_cut = np.logical_and(self.region, self.cut)
@@ -1875,7 +1875,7 @@ class pol_map(object):
             P_cut = np.sqrt(Q_cut**2+U_cut**2)/I_cut
             P_cut_err = np.sqrt((Q_cut**2*Q_cut_err**2 + U_cut**2*U_cut_err**2 + 2.*Q_cut*U_cut*QU_cut_err)/(Q_cut**2 + U_cut**2) + ((Q_cut/I_cut)**2 + (U_cut/I_cut)**2)*I_cut_err**2 - 2.*(Q_cut/I_cut)*IQ_cut_err - 2.*(U_cut/I_cut)*IU_cut_err)/I_cut
 
-            PA_cut = princ_angle((90./np.pi)*np.arctan2(U_cut,Q_cut)+180.)
+            PA_cut = 360.-princ_angle((90./np.pi)*np.arctan2(U_cut,Q_cut))
             PA_cut_err = (90./(np.pi*(Q_cut**2+U_cut**2)))*np.sqrt(U_cut**2*Q_cut_err**2 + Q_cut**2*U_cut_err**2 - 2.*Q_cut*U_cut*QU_cut_err)
 
         if hasattr(self, 'cont'):
diff --git a/src/lib/reduction.py b/src/lib/reduction.py
index e1ba275..d58f57c 100755
--- a/src/lib/reduction.py
+++ b/src/lib/reduction.py
@@ -407,6 +407,200 @@ def deconvolve_array(data_array, headers, psf='gaussian', FWHM=1., scale='px',
     return deconv_array
 
 
+def get_error2(data_array, headers, error_array=None, data_mask=None,
+            display=False, savename=None, plots_folder="",
+            return_background=False):
+    """
+    ----------
+    Inputs:
+    data_array : numpy.ndarray
+        Array containing the data to study (2D float arrays).
+    headers : header list
+        Headers associated with the images in data_array.
+    error_array : numpy.ndarray, optional
+        Array of images (2D floats, aligned and of the same shape) containing
+        the error in each pixel of the observation images in data_array.
+        If None, will be initialized to zeros.
+        Defaults to None.
+    data_mask : numpy.ndarray, optional
+        2D boolean array delimiting the data to work on.
+        If None, will be initialized with a full true mask.
+        Defaults to None.
+    display : boolean, optional
+        If True, data_array will be displayed with a rectangle around the
+        sub-image selected for background computation.
+        Defaults to False.
+    savename : str, optional
+        Name of the figure the map should be saved to. If None, the map won't
+        be saved (only displayed). Only used if display is True.
+        Defaults to None.
+    plots_folder : str, optional
+        Relative (or absolute) filepath to the folder in wich the map will
+        be saved. Not used if savename is None.
+        Defaults to current folder.
+    return_background : boolean, optional
+        If True, the pixel background value for each image in data_array is
+        returned.
+        Defaults to False.
+    ----------
+    Returns:
+    data_array : numpy.ndarray
+        Array containing the data to study minus the background.
+    headers : header list
+        Updated headers associated with the images in data_array.
+    error_array : numpy.ndarray
+        Array containing the background values associated to the images in
+        data_array.
+    background : numpy.ndarray
+        Array containing the pixel background value for each image in
+        data_array.
+        Only returned if return_background is True.
+    """
+    # Crop out any null edges
+    if error_array is None:
+        error_array = np.zeros(data_array.shape)
+    n_data_array, n_error_array = deepcopy(data_array), deepcopy(error_array)
+    if not data_mask is None:
+        data, mask = n_data_array, deepcopy(data_mask)
+    else:
+        data, _, _ = crop_array(n_data_array, headers, step=5, null_val=0., inside=False)
+        data_mask = np.ones(data_array[0].shape, dtype=bool)
+        mask = np.ones(data[0].shape, dtype=bool)
+
+    error_bkg = np.ones(n_data_array.shape)
+    std_bkg = np.zeros((data.shape[0]))
+    background = np.zeros((data.shape[0]))
+    
+    if display:
+        fig_h, ax_h = plt.subplots(figsize=(10,6), constrained_layout=True)
+        date_time = np.array([headers[i]['date-obs']+';'+headers[i]['time-obs']
+            for i in range(len(headers))])
+        date_time = np.array([datetime.strptime(d,'%Y-%m-%d;%H:%M:%S')
+            for d in date_time])
+    for i, image in enumerate(data):
+        #Compute the Count-rate histogram for the image
+        n_mask = np.logical_and(mask,image>0.)
+        
+        #n_bins = np.fix(np.sqrt(image[n_mask].size)).astype(int) # Square-root
+        #n_bins = np.ceil(np.log2(image[n_mask].size)).astype(int)+1 # Sturges
+        #n_bins = 2*np.fix(np.power(image[n_mask].size,1/3)).astype(int) # Rice
+        #n_bins = np.fix((image[n_mask].max()-image[n_mask].min())/(3.5*image[n_mask].std()/np.power(image[n_mask].size,1/3))).astype(int) # Scott
+        n_bins = np.fix((image[n_mask].max()-image[n_mask].min())/(2*np.subtract(*np.percentile(image[n_mask], [75, 25]))/np.power(image[n_mask].size,1/3))).astype(int) # Freedman-Diaconis
+        
+        hist, bin_edges = np.histogram(np.log(image[n_mask]),bins=n_bins)
+        bin_centers = np.exp((bin_edges[:-1]+bin_edges[1:])/2)
+        #Take the background as the count-rate with the maximum number of pixels
+        hist_max = bin_centers[np.argmax(hist)]
+        bkg = np.sqrt(np.sum(image[np.abs(image-hist_max)/hist_max<0.5]**2)/image[np.abs(image-hist_max)/hist_max<0.5].size)
+        #bkg = np.sum(image[n_mask]*1/np.abs(image[n_mask]-hist_max))/np.sum(1/np.abs(image[n_mask]-hist_max))
+        #bkg = np.percentile(image[image<hist_max],25.)
+        #bkg = 0.95*hist_max
+        if display:
+            ax_h.plot(bin_centers,hist,'+',color="C{0:d}".format(i),alpha=0.8,label=headers[i]['filtnam1']+' ('+str(date_time[i])+") with n_bins = {0:d}".format(n_bins))
+            ax_h.plot([bkg,bkg],[hist.min(), hist.max()],'x--',color="C{0:d}".format(i),alpha=0.8)
+        error_bkg[i] *= bkg
+       
+        # Quadratically add uncertainties in the "correction factors" (see Kishimoto 1999)
+        #wavelength dependence of the polariser filters
+        #estimated to less than 1%
+        err_wav = data_array[i]*0.01
+        #difference in PSFs through each polarizers
+        #estimated to less than 3%
+        err_psf = data_array[i]*0.03
+        #flatfielding uncertainties
+        #estimated to less than 3%
+        err_flat = data_array[i]*0.03
+
+        n_error_array[i] = np.sqrt(n_error_array[i]**2 + error_bkg[i]**2 + err_wav**2 + err_psf**2 + err_flat**2)
+        
+        #Substract background
+        n_data_array[i][data_mask] = n_data_array[i][data_mask] - bkg
+        n_data_array[i][np.logical_and(data_mask,n_data_array[i] <= 0.01*bkg)] = 0.01*bkg#n_data_array[i][np.logical_and(data_mask,n_data_array[i] > 0.)].min()
+ 
+        std_bkg[i] = image[image<2*bkg].std()
+        background[i] = bkg
+
+        if (data_array[i] < 0.).any():
+            print(data_array[i])
+        #if i==0:
+            #np.savetxt("output/s_bg.txt",error_bkg[i])
+            #np.savetxt("output/s_wav.txt",err_wav)
+            #np.savetxt("output/s_psf.txt",err_psf)
+            #np.savetxt("output/s_flat.txt",err_flat)
+    if display:   
+        ax_h.set_xscale('log')
+        ax_h.set_xlim([background.mean()*1e-2,background.mean()*1e2])
+        ax_h.set_xlabel(r"Count rate [$s^{-1}$]")
+        #ax_h.set_yscale('log')
+        ax_h.set_ylabel(r"Number of pixels in bin")
+        ax_h.set_title("Histogram for each observation")
+        plt.legend()
+
+        plt.rcParams.update({'font.size': 15})
+        convert_flux = np.array([head['photflam'] for head in headers])
+        filt = np.array([headers[i]['filtnam1'] for i in range(len(headers))])
+        dict_filt = {"POL0":'r', "POL60":'g', "POL120":'b'}
+        c_filt = np.array([dict_filt[f] for f in filt])
+        
+        fig,ax = plt.subplots(figsize=(10,6), constrained_layout=True)
+        for f in np.unique(filt):
+            mask = [fil==f for fil in filt]
+            ax.scatter(date_time[mask], background[mask]*convert_flux[mask],
+                    color=dict_filt[f],label="{0:s}".format(f))
+        ax.errorbar(date_time, background*convert_flux,
+                yerr=std_bkg*convert_flux, fmt='+k',
+                markersize=0, ecolor=c_filt)
+        # Date handling
+        locator = mdates.AutoDateLocator()
+        formatter = mdates.ConciseDateFormatter(locator)
+        ax.xaxis.set_major_locator(locator)
+        ax.xaxis.set_major_formatter(formatter)
+        ax.set_xlabel("Observation date and time")
+        ax.set_ylabel(r"Flux [$ergs \cdot cm^{-2} \cdot s^{-1} \cdot \AA^{-1}$]")
+        #ax.set_title("Background flux and error computed for each image")
+        plt.legend()
+
+        plt.rcParams.update({'font.size': 15})
+        fig2, ax2 = plt.subplots(figsize=(10,10))
+        data0 = data[0]*convert_flux[0]
+        bkg_data0 = data0 <= background[0]*convert_flux[0]
+        instr = headers[0]['instrume']
+        rootname = headers[0]['rootname']
+        exptime = headers[0]['exptime']
+        filt = headers[0]['filtnam1']
+        #plots
+        #im = ax2.imshow(data0, vmin=data0.min(), vmax=data0.max(), origin='lower', cmap='gray')
+        im = ax2.imshow(data0, norm=LogNorm(data0[data0>0.].mean()/10.,data0.max()), origin='lower', cmap='gray')
+        bkg_im = ax2.imshow(bkg_data0, origin='lower', cmap='Reds', alpha=0.7)
+        ax2.annotate(instr+":"+rootname, color='white', fontsize=10,
+                xy=(0.02, 0.95), xycoords='axes fraction')
+        ax2.annotate(filt, color='white', fontsize=14, xy=(0.02, 0.02),
+                xycoords='axes fraction')
+        ax2.annotate(str(exptime)+" s", color='white', fontsize=10, xy=(0.80, 0.02),
+                xycoords='axes fraction')
+        ax2.set(#title="Location of background computation.",
+                xlabel='pixel offset',
+                ylabel='pixel offset')
+
+        fig2.subplots_adjust(hspace=0, wspace=0, right=0.85)
+        cbar_ax = fig2.add_axes([0.9, 0.12, 0.02, 0.75])
+        fig2.colorbar(im, cax=cbar_ax, label=r"Flux [$ergs \cdot cm^{-2} \cdot s^{-1} \cdot \AA^{-1}$]")
+
+        if not(savename is None):
+            fig.savefig(plots_folder+savename+"_background_flux.png",
+                    bbox_inches='tight')
+            fig2.savefig(plots_folder+savename+'_'+filt+'_background_location.png',
+                    bbox_inches='tight')
+            fig_h.savefig(plots_folder+savename+'_histograms.png', bbox_inches='tight')
+
+        plt.show()
+
+    if return_background:
+        return n_data_array, n_error_array, headers, background #np.array([n_error_array[i][0,0] for i in range(n_error_array.shape[0])])
+    else:
+        return n_data_array, n_error_array, headers
+
+
 def get_error(data_array, headers, error_array=None, data_mask=None,
             sub_shape=None, display=False, savename=None, plots_folder="",
             return_background=False):
@@ -465,10 +659,12 @@ def get_error(data_array, headers, error_array=None, data_mask=None,
     # Crop out any null edges
     if error_array is None:
         error_array = np.zeros(data_array.shape)
+    n_data_array, n_error_array = deepcopy(data_array), deepcopy(error_array)
     if not data_mask is None:
-        data, error, mask = data_array, error_array, data_mask
+        data, error, mask = n_data_array, n_error_array, deepcopy(data_mask)
     else:
-        data, error, _ = crop_array(data_array, headers, error_array, step=5, null_val=0., inside=False)
+        data, _, _ = crop_array(n_data_array, headers, n_error_array, step=5, null_val=0., inside=False)
+        data_mask = np.ones(n_data_array[0].shape, dtype=bool)
         mask = np.ones(data[0].shape, dtype=bool)
     if sub_shape is None:
         sub_shape = (np.array(data_array.shape[1:])/10).astype(int)
@@ -481,8 +677,9 @@ def get_error(data_array, headers, error_array=None, data_mask=None,
     diff = (sub_shape-1).astype(int)
     temp = np.zeros((shape[0],shape[1]-diff[0],shape[2]-diff[1]))
     error_bkg = np.ones(data_array.shape)
-    rectangle = []
+    std_bkg = np.zeros((data.shape[0]))
     background = np.zeros((shape[0]))
+    rectangle = []
 
     for i,image in enumerate(data):
         # Find the sub-image of smallest integrated flux (suppose no source)
@@ -518,9 +715,12 @@ def get_error(data_array, headers, error_array=None, data_mask=None,
         error_array[i] = np.sqrt(error_array[i]**2 + error_bkg[i]**2 + err_wav**2 + err_psf**2 + err_flat**2)
         
         #Substract background
-        data_array[i] = np.abs(data_array[i] - sub_image.mean())
+        n_data_array[i][data_mask] = n_data_array[i][data_mask] - bkg
+        n_data_array[i][np.logical_and(data_mask,n_data_array[i] <= 0.)] = n_data_array[i][np.logical_and(data_mask,n_data_array[i] > 0.)].min()
  
-        background[i] = sub_image.sum()
+        std_bkg[i] = image[image<2*bkg].std()
+        background[i] = bkg
+
         if (data_array[i] < 0.).any():
             print(data_array[i])
         #if i==0:
@@ -531,7 +731,7 @@ def get_error(data_array, headers, error_array=None, data_mask=None,
 
     if display:
         plt.rcParams.update({'font.size': 15})
-        convert_flux = headers[0]['photflam']
+        convert_flux = np.array([head['photflam'] for head in headers])
         date_time = np.array([headers[i]['date-obs']+';'+headers[i]['time-obs']
             for i in range(len(headers))])
         date_time = np.array([datetime.strptime(d,'%Y-%m-%d;%H:%M:%S')
@@ -543,10 +743,10 @@ def get_error(data_array, headers, error_array=None, data_mask=None,
         fig,ax = plt.subplots(figsize=(10,6), constrained_layout=True)
         for f in np.unique(filt):
             mask = [fil==f for fil in filt]
-            ax.scatter(date_time[mask], background[mask]*convert_flux,
+            ax.scatter(date_time[mask], background[mask]*convert_flux[mask],
                     color=dict_filt[f],label="{0:s}".format(f))
         ax.errorbar(date_time, background*convert_flux,
-                yerr=error_array[:,0,0]*convert_flux, fmt='+k',
+                yerr=std_bkg*convert_flux, fmt='+k',
                 markersize=0, ecolor=c_filt)
         # Date handling
         locator = mdates.AutoDateLocator()
@@ -560,7 +760,7 @@ def get_error(data_array, headers, error_array=None, data_mask=None,
 
         plt.rcParams.update({'font.size': 15})
         fig2, ax2 = plt.subplots(figsize=(10,10))
-        data0 = data[0]*convert_flux
+        data0 = data[0]*convert_flux[0]
         instr = headers[0]['instrume']
         rootname = headers[0]['rootname']
         exptime = headers[0]['exptime']
@@ -589,16 +789,12 @@ def get_error(data_array, headers, error_array=None, data_mask=None,
                     bbox_inches='tight')
             fig2.savefig(plots_folder+savename+'_'+filt+'_background_location.png',
                     bbox_inches='tight')
-            vmin = np.min(np.log10(data[data > 0.]))
-            vmax = np.max(np.log10(data[data > 0.]))
             plot_obs(data, headers, vmin=data.min(), vmax=data.max(),
                     rectangle=rectangle,
                     savename=savename+"_background_location",
                     plots_folder=plots_folder)
         else:
-            vmin = np.min(np.log10(data[data > 0.]))
-            vmax = np.max(np.log10(data[data > 0.]))
-            plot_obs(np.log10(data), headers, vmin=vmin, vmax=vmax,
+            plot_obs(data, headers, vmin=vmin, vmax=vmax,
                     rectangle=rectangle)
 
         plt.show()
@@ -696,13 +892,17 @@ def rebin_array(data_array, error_array, headers, pxsize, scale,
             mask = sum_image > 0.
             new_error = np.zeros(rms_image.shape)
             if operation.lower() in ["mean", "average", "avg"]:
-                new_error[mask] = np.sqrt(bin_ndarray(error**2*image,
-                    new_shape=new_shape, operation='average')[mask]/sum_image[mask])
+                new_error = np.sqrt(bin_ndarray(error**2,
+                    new_shape=new_shape, operation='average'))
+                #new_error[mask] = np.sqrt(bin_ndarray(error**2*image,
+                #    new_shape=new_shape, operation='average')[mask]/sum_image[mask])
                 #new_error[mask] = np.sqrt(bin_ndarray(error**2,
                 #    new_shape=new_shape, operation='average')[mask])
             else:
-                new_error[mask] = np.sqrt(bin_ndarray(error**2*image,
-                    new_shape=new_shape, operation='sum')[mask]/sum_image[mask])
+                new_error = np.sqrt(bin_ndarray(error**2,
+                    new_shape=new_shape, operation='sum'))
+                #new_error[mask] = np.sqrt(bin_ndarray(error**2*image,
+                #    new_shape=new_shape, operation='sum')[mask]/sum_image[mask])
                 #new_error[mask] = np.sqrt(bin_ndarray(error**2,
                 #    new_shape=new_shape, operation='sum')[mask])
             rebinned_error.append(np.sqrt(rms_image**2 + new_error**2))
@@ -724,8 +924,8 @@ def rebin_array(data_array, error_array, headers, pxsize, scale,
         return rebinned_data, rebinned_error, rebinned_headers, Dxy, data_mask
 
 
-def align_data(data_array, headers, error_array=None, upsample_factor=1.,
-        ref_data=None, ref_center=None, return_shifts=False):
+def align_data(data_array, headers, error_array=None, background=None,
+        upsample_factor=1., ref_data=None, ref_center=None, return_shifts=False):
     """
     Align images in data_array using cross correlation, and rescale them to
     wider images able to contain any rotation of the reference image.
@@ -789,10 +989,8 @@ def align_data(data_array, headers, error_array=None, upsample_factor=1.,
     if not same:
         raise ValueError("All images in data_array must have same shape as\
             ref_data")
-    if error_array is None:
-        _, error_array, headers, background = get_error(data_array, headers, return_background=True)
-    else:
-        _, _, headers, background = get_error(data_array, headers, error_array=error_array, return_background=True)
+    if (error_array is None) or (background is None):
+        _, error_array, headers, background = get_error2(data_array, headers, return_background=True)
 
     # Crop out any null edges
     #(ref_data must be cropped as well)
@@ -801,8 +999,8 @@ def align_data(data_array, headers, error_array=None, upsample_factor=1.,
     full_headers.append(headers[0])
     err_array = np.concatenate((error_array,[np.zeros(ref_data.shape)]),axis=0)
 
-    full_array, err_array, full_headers = crop_array(full_array, full_headers, err_array, step=5,
-            inside=False)
+    full_array, err_array, full_headers = crop_array(full_array, full_headers,
+            err_array, step=5, inside=False, null_val=0.)
 
     data_array, ref_data, headers = full_array[:-1], full_array[-1], full_headers[:-1]
     error_array = err_array[:-1]
@@ -833,7 +1031,7 @@ def align_data(data_array, headers, error_array=None, upsample_factor=1.,
     shifts, errors = [], []
     for i,image in enumerate(data_array):
         # Initialize rescaled images to background values
-        rescaled_error[i] *= background[i]
+        rescaled_error[i] *= 0.01*background[i]
         # Get shifts and error by cross-correlation to ref_data
         shift, error, phase_diff = phase_cross_correlation(ref_data/ref_data.max(), image/image.max(),
                 upsample_factor=upsample_factor)
@@ -876,7 +1074,7 @@ def align_data(data_array, headers, error_array=None, upsample_factor=1.,
         headers[i].update(headers_wcs[i].to_header())
 
     data_mask = rescaled_mask.all(axis=0)
-    data_array, error_array, data_mask, headers = crop_array(rescaled_image, headers, rescaled_error, data_mask)
+    data_array, error_array, data_mask, headers = crop_array(rescaled_image, headers, rescaled_error, data_mask, null_val=0.01*background)
 
     if return_shifts:
         return data_array, error_array, headers, data_mask, shifts, errors
@@ -1316,7 +1514,7 @@ def compute_Stokes(data_array, error_array, data_mask, headers,
         P_diluted = np.sqrt(Q_diluted**2+U_diluted**2)/I_diluted
         P_diluted_err = (1./I_diluted)*np.sqrt((Q_diluted**2*Q_diluted_err**2 + U_diluted**2*U_diluted_err**2 + 2.*Q_diluted*U_diluted*QU_diluted_err)/(Q_diluted**2 + U_diluted**2) + ((Q_diluted/I_diluted)**2 + (U_diluted/I_diluted)**2)*I_diluted_err**2 - 2.*(Q_diluted/I_diluted)*IQ_diluted_err - 2.*(U_diluted/I_diluted)*IU_diluted_err)
 
-        PA_diluted = princ_angle((90./np.pi)*np.arctan2(U_diluted,Q_diluted)+180.)
+        PA_diluted = 360.-princ_angle((90./np.pi)*np.arctan2(U_diluted,Q_diluted))
         PA_diluted_err = (90./(np.pi*(Q_diluted**2 + U_diluted**2)))*np.sqrt(U_diluted**2*Q_diluted_err**2 + Q_diluted**2*U_diluted_err**2 - 2.*Q_diluted*U_diluted*QU_diluted_err)
 
         for header in headers:
@@ -1376,7 +1574,7 @@ def compute_pol(I_stokes, Q_stokes, U_stokes, Stokes_cov, headers):
     P = np.zeros(I_stokes.shape)
     P[mask] = I_pol[mask]/I_stokes[mask]
     PA = np.zeros(I_stokes.shape)
-    PA[mask] = princ_angle((90./np.pi)*np.arctan2(U_stokes[mask],Q_stokes[mask])+180.)
+    PA[mask] = (90./np.pi)*np.arctan2(U_stokes[mask],Q_stokes[mask])
 
     if (P>1).any():
         print("WARNING : found {0:d} pixels for which P > 1".format(P[P>1.].size))
@@ -1575,7 +1773,7 @@ def rotate_Stokes(I_stokes, Q_stokes, U_stokes, Stokes_cov, data_mask, headers,
     P_diluted = np.sqrt(Q_diluted**2+U_diluted**2)/I_diluted
     P_diluted_err = (1./I_diluted)*np.sqrt((Q_diluted**2*Q_diluted_err**2 + U_diluted**2*U_diluted_err**2 + 2.*Q_diluted*U_diluted*QU_diluted_err)/(Q_diluted**2 + U_diluted**2) + ((Q_diluted/I_diluted)**2 + (U_diluted/I_diluted)**2)*I_diluted_err**2 - 2.*(Q_diluted/I_diluted)*IQ_diluted_err - 2.*(U_diluted/I_diluted)*IU_diluted_err)
 
-    PA_diluted = princ_angle((90./np.pi)*np.arctan2(U_diluted,Q_diluted)+180.)
+    PA_diluted = 360.-princ_angle((90./np.pi)*np.arctan2(U_diluted,Q_diluted))
     PA_diluted_err = (90./(np.pi*(Q_diluted**2 + U_diluted**2)))*np.sqrt(U_diluted**2*Q_diluted_err**2 + Q_diluted**2*U_diluted_err**2 - 2.*Q_diluted*U_diluted*QU_diluted_err)
 
     for header in new_headers:
diff --git a/src/overplot.py b/src/overplot.py
index 4bb2ae5..5940e1d 100755
--- a/src/overplot.py
+++ b/src/overplot.py
@@ -20,30 +20,30 @@ levelsMorganti = np.array([1.,2.,3.,8.,16.,32.,64.,128.])
 #levels18GHz = np.array([0.6, 1.5, 3, 6, 12, 24, 48, 96])/100.*Stokes_18GHz[0].data.max()
 levels18GHz = levelsMorganti*0.28*1e-3
 A = overplot_radio(Stokes_UV, Stokes_18GHz)
-A.plot(levels=levels18GHz, SNRp_cut=3.0, SNRi_cut=60.0, savename='../plots/IC5063_x3nl030/18GHz_overplot_forced.png')
+A.plot(levels=levels18GHz, SNRp_cut=3.0, SNRi_cut=60.0, savename='../plots/IC5063_x3nl030/18GHz_overplot.png')
 
 #levels24GHz = np.array([1.,1.5, 3, 6, 12, 24, 48, 96])/100.*Stokes_24GHz[0].data.max()
 levels24GHz = levelsMorganti*0.46*1e-3
 B = overplot_radio(Stokes_UV, Stokes_24GHz)
-B.plot(levels=levels24GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/24GHz_overplot_forced.png')
+B.plot(levels=levels24GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/24GHz_overplot.png')
 
 levels103GHz = np.linspace(1,99,11)/100.*np.max(deepcopy(Stokes_103GHz[0].data[Stokes_103GHz[0].data > 0.]))
 C = overplot_radio(Stokes_UV, Stokes_103GHz)
-C.plot(levels=levels103GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/103GHz_overplot_forced.png')
+C.plot(levels=levels103GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/103GHz_overplot.png')
 
 levels229GHz = np.linspace(1,99,11)/100.*np.max(deepcopy(Stokes_229GHz[0].data[Stokes_229GHz[0].data > 0.]))
 D = overplot_radio(Stokes_UV, Stokes_229GHz)
-D.plot(levels=levels229GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/229GHz_overplot_forced.png')
+D.plot(levels=levels229GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/229GHz_overplot.png')
 
 levels357GHz = np.linspace(1,99,11)/100.*np.max(deepcopy(Stokes_357GHz[0].data[Stokes_357GHz[0].data > 0.]))
 E = overplot_radio(Stokes_UV, Stokes_357GHz)
-E.plot(levels=levels357GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/357GHz_overplot_forced.png')
+E.plot(levels=levels357GHz, SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/357GHz_overplot.png')
 
 #F = overplot_pol(Stokes_UV, Stokes_S2)
-#F.plot(SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/S2_overplot_forced.png', norm=LogNorm(vmin=5e-20,vmax=5e-18))
+#F.plot(SNRp_cut=3.0, SNRi_cut=80.0, savename='../plots/IC5063_x3nl030/S2_overplot.png', norm=LogNorm(vmin=5e-20,vmax=5e-18))
 
 G = overplot_pol(Stokes_UV, Stokes_IR, cmap='inferno')
-G.plot(SNRp_cut=3.0, SNRi_cut=60.0, savename='../plots/IC5063_x3nl030/IR_overplot_forced.png', norm=LogNorm(vmin=1e-17,vmax=5e-15), cmap='inferno_r')
+G.plot(SNRp_cut=3.0, SNRi_cut=60.0, savename='../plots/IC5063_x3nl030/IR_overplot.png', norm=LogNorm(vmin=1e-17,vmax=5e-15), cmap='inferno_r')
 
 #data_folder1 = "../data/M87/POS1/"
 #plots_folder1 = "../plots/M87/POS1/"