changed compare_beam_spherical_harmonics.py. The function now truncates lmax and compared it with the original map. Much faster with the use of built in functions