Added a new function, gsm_noise_covar(), to global_sky_model.py that will calculate the noise covariance for galactic forground emission. It's not quite done b/c I haven't figured out how to go from healpix coords to Img coords.