def get_water_atom_indices(topology):
    """Get list of atom indices in "WAT" residues"""
    indices = []
    water_residues = [r for r in topology.residues() if r.name == "WAT"]
    for water_residue in water_residues:
        for a in water_residue.atoms():
            indices.append(a.index)
    return indices