Calculate delta (physical time in bond_harmonic_swimmer.cpp)
Add tags in force calculation, add a template for swimmer wave.
Works for one bond
Add n1 and n2 in the definition of the harmonic/swimmer bond