Create a scatter plot on HEALPix projected axes.
    Inputs: lon (deg), lat (deg)

def projScatter(lon, lat, **kwargs):
    """
    Create a scatter plot on HEALPix projected axes.
    Inputs: lon (deg), lat (deg)
    """
    hp.projscatter(lon, lat, lonlat=True, **kwargs)