from netgen.occ import *
import netgen.gui

disk = Circle((0,0), 1).Face()
geo = OCCGeometry(disk, dim=2)

ngmesh = geo.GenerateMesh(maxh=1.)
ngmesh.Save("disk_original.vol")

from ngsPETSc import MeshMapping

petscMap = MeshMapping(ngmesh)
plex = petscMap.petscPlex
plex.setRefinementUniform(True)
rplex = plex.refine()

input("Refine ?")

ngMap = MeshMapping(rplex, geo=geo)
rngmesh = ngMap.ngMesh
rngmesh.SetGeometry(ngmesh.GetGeometry())

rngmesh.Curve(2)
rngmesh.Save("disk.vol")