The flux container currently supports indexing:
returns a new container, with only specified flavor.
But if you try to access several flavors in a list:
flux[['NU_E','NU_E_BAR']]
it fails.
This is annoying behaviour, as it prevents the user from summing over all neutrinos or antineutrinos etc.