New Initialisation function for Complex Geometry

With this MR a new function is added to set up the block forest with meshes. This should allow for easier benchmarking of weak and strong scaling.