Parametric update - Ragnhild is working on finding the right mask for the infilling in model 90.91.
When she fits she gets good fits over the full sky including the Galactic plane. Then she masks, but the power spectra after infilling the mask don't agree with the input. Why is she getting good chi^2 before cutting the sky, but poor power spectra after in-filling?
Question: why not estimate spectra on cut sky? Because of EB leakage affects the B-mode spectra.
Question: does the problem exist only with Foregrounds? No, same problem with/out foregrounds, indicating that something is broken in the code.
Suggestion: perhaps use MASTER for power spectrum estimation, not Anafast?