Added functions from qgea to load Q and gsm files so that I wouldn't have to import qgea, which requires healpy, on nersc. Also added a function to compute MQdagNinv so that I would only compute it once for each array instead of for every array and frequency. Also made the plot look prettier.