from pymoab import core
from pymoab.rng import Range

# create moab instance
mb = core.Core()

# create a vertices
# verts = mb.create_vertices([[.5,2.,.3],[.4,2,3]])
# ^^ does not work (bug in code)
verts = mb.create_vertices([.5,2.,.3,.4,2.,3.])
coords = [[.5,2.,.3],[.4,2.,3.]]

# create a meshset with these vertices
vert_set = mb.create_meshset()
mb.add_entities(vert_set, verts)
# mb.write_file("test-vert-set.h5m", vert_set)
# ^^ this did not work

# write out all meshsets
mb.write_file("test-vert-set.h5m")
mb.write_file("test-vert-set.vtk")

# write out a singe meshet (need to first convert to a range)
r = Range(vert_set)
mb.write_file("test-vert-set.vtk", r)