def read(filename, HDU=0, return_header=False):
    """Quick function for reading one HDU and header

    HDU is name or index of HDU, default 0
    Returns the data array and a header dictionary"""
    with open(filename) as f:
        if return_header:
            return f[HDU].read_all(), f[HDU].header
        else:
            return f[HDU].read_all()
