program mesh_to_vtk

use ESMF
implicit none

type(ESMF_Mesh) :: mesh
integer :: rc, localrc

call ESMF_Initialize(defaultCalKind=ESMF_CALKIND_GREGORIAN, rc=localrc)
mesh = ESMF_MeshCreate("hgrid.nc", ESMF_FILEFORMAT_ESMFMESH)
call ESMF_MeshWrite(mesh, "hgrid", rc=rc)
call ESMF_Finalize()

end program