def create_mesh(outfile, topo_dim, geom_dim, x_coord, y_coord, z_coord):
rho_coords, rho_shape = get_coord_array(x_coord, y_coord, z_coord)
print "shape of mesh: %s" % rho_shape
print "number of points in mesh %s" % len(rho_coords)
mesh_example = MeshExample(topo_dim,geom_dim) # topo_dim = 1, geom_dim = 3
num_vertices = len(rho_coords)
num_cells = num_vertices
mesh_example.initializing_empty_grid(num_vertices, num_cells)
mesh_example.create_vertices(rho_coords)
mesh_example.create_cells(num_cells)
mesh_example.close()
fx=File(outfile)
fx << mesh_example.mesh
return mesh_example.mesh
