Use some P4EST_F90_TOPIDX in geometry wrap
Divided ReadMeshFromHDF5 into two parts

The second one only constructs the p4est
Added some comments