Just to iterate:
If I have a graph
edge_list=[[0,1,1],[0,2,1],[0,3,1],[1,2,1],[1,3,1],[2,3,1],[0,4,2],[0,5,2],[0,6,2],[1,7,2],[1,8,2],[1,9,2],[2,10,2],[2,11,2],[2,12,2],[3,13,2],[3,14,2],[3,15,2],[4,5,2],[4,6,2],[5,6,2],[7,8,2],[7,9,2],[8,9,2],[10,11,2],[10,12,2],[11,12,2],[13,14,2],[13,15,2],[14,15,2]]
g = Graph(directed=False)
weight = g.new_edge_property('int')
g.add_edge_list(edge_list, eprops=[weight])
and a nested blockmodel
state = minimize_nested_blockmodel_dl(g, state_args=dict(eweight=weight))
then I can either draw the graph directly with a partition from a given level l
state.project_level(l=0).draw()
or store the partition in a property map
group = state.project_partition(l=0, j=0)
and then use this map in conventional drawing
graph_draw(g, vertex_fill_color=group)
Best
Haiko