Index: /issm/trunk-jpl/src/m/partition/partitioner.py
===================================================================
--- /issm/trunk-jpl/src/m/partition/partitioner.py	(revision 25131)
+++ /issm/trunk-jpl/src/m/partition/partitioner.py	(revision 25132)
@@ -131,3 +131,3 @@
 
     # Output
-    return part
+    return part, md
Index: /issm/trunk-jpl/test/NightlyRun/test418.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test418.py	(revision 25131)
+++ /issm/trunk-jpl/test/NightlyRun/test418.py	(revision 25132)
@@ -2,18 +2,19 @@
 # this test may crash
 import numpy as np
+
+from AreaAverageOntoPartition import *
 from model import *
+from parameterize import *
+from partitioner import *
+from setflowequation import *
+from setmask import *
 from socket import gethostname
+from solve import *
 from triangle import *
-from setmask import *
-from parameterize import *
-from setflowequation import *
-from solve import *
-from partitioner import *
-from AreaAverageOntoPartition import *
+
 
 #test partitioning, and partition averaging
 # python cannot handle resolutions greater than 30010
 md = triangle(model(), '../Exp/Square.exp', 30000.)
-#print md.mesh.numberofvertices
 md = setmask(md, '../Exp/SquareShelf.exp', '')
 md = parameterize(md, '../Par/SquareSheetShelf.py')
@@ -36,5 +37,5 @@
 # - Run valgrind and fix the above
 #
-partition = partitioner(md, 'package', 'chaco', 'npart', npart) - 1
+partition, md = partitioner(md, 'package', 'chaco', 'npart', npart) - 1
 
 vector = np.arange(1, 1 + md.mesh.numberofvertices, 1).reshape(-1, 1)