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 all the vals temporaries here?