def read_cif(cif):
    (atoms, alllines) = fix_cif(cif)
    (ts, tm, tinds) = get_tsites(cif)
    (os, om, oinds) = get_osites(cif)
    return (atoms, ts, tm, tinds, os, om, oinds)
def read_cif(cif):
    (atoms, alllines) = fix_cif(cif)
    (ts, tm, tinds) = get_tsites(cif)
    (os, om, oinds) = get_osites(cif)
    return (atoms, ts, tm, tinds, os, om, oinds)