Create skeleton for imposing boundaries for mhd.
Introduce update_us function for updating penalization field.
Move mean-flow stuff to fsi.
Integral output data files now names *.t