def compare_gsm_to_old():
g = GlobalSkyModel2016(freq_unit='GHz', resolution='hi', unit='TCMB')
d = g.generate(0.408)
g.view()
g_old = GlobalSkyModel(freq_unit='GHz', basemap='haslam')
d_old = g_old.generate(0.408)
g_old.view()
def compare_gsm_to_old():
g = GlobalSkyModel2016(freq_unit='GHz', resolution='hi', unit='TCMB')
d = g.generate(0.408)
g.view()
g_old = GlobalSkyModel(freq_unit='GHz', basemap='haslam')
d_old = g_old.generate(0.408)
g_old.view()
def test_set_methods():
gsm = GlobalSkyModel()
gsm.generate(40)
#gsm.view()
gsm.set_basemap('haslam')
gsm.view()
gsm.set_basemap('5deg')
gsm.view()
gsm.set_basemap('wmap')
gsm.view()
gsm.set_freq_unit('GHz')
gsm.view()