First commit. Does it work? Need to better use arguments
PBC lack. To be done now
Works with a kind of time-periodic boundary conditions
Optimizations in compilation flags