generalize ability to add systems

[sconde]

At the moment, only Maxwell and Euler 1D systems are working. Moreover these are only working for specific integrators (LF and LW). Need to generalize this for all integrators