Cubedsphere

This code implements a cubed sphere mesh, with gnonomic projection. I also left some code there for the dmplex refinement approach, which currently causes dmplex to segfault, since dmplex doesn't expect embedded quads in 3D.
Why not do a panel_numbering.shape = (6, Nx**2) before starting and thereby get rid of all the reshapes?