static char help[] = "blablabla";


/* The PETSc packages we need: */

#include <petscsys.h>
#include <petscviewer.h>
#include <petscdmplex.h>
#include <petscdmlabel.h>

int main(int argc, char **args){
  PetscViewer mshviewer;
  PetscErrorCode ierr;
  DM mesh;  
  /* Firing up Petsc */
 ierr= PetscInitialize(&argc, &args,(char*)0,help);
 /* Creating the mesh from an existing gmsh file */
 ierr =  PetscViewerVTKOpen(PETSC_COMM_WORLD, "mymesh.vtk", FILE_MODE_WRITE, &mshviewer);

ierr =  DMPlexCreateFromFile(PETSC_COMM_WORLD,"testmesh_2D_box_quad.msh", PETSC_FALSE, &mesh);
 

 ierr =   DMView(mesh, mshviewer);

 /* Freein Workspace */
 ierr =  PetscViewerDestroy(&mshviewer);
 ierr =  PetscFinalize();

  return 0;
  
}