Now compiles and runs (without h5 output yet)