Index: /issm/trunk-jpl/test/NightlyRun/test417.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test417.py	(revision 23257)
+++ /issm/trunk-jpl/test/NightlyRun/test417.py	(revision 23258)
@@ -88,5 +88,5 @@
 
 #Fields and tolerances to track changes
-md.qmu.results = md.results.dakota;
+md.qmu.results = md.results.dakota
 
 #ok, mass flux of 3 profiles should be -3 Gt/yr -3 Gt/yr and the sum, which is -6 Gt/yr
@@ -104,3 +104,3 @@
 field_names      = ['moments']
 field_tolerances = [1e-11]
-field_values = [md.results.dakota.montecarlo]
+field_values     = [md.results.dakota.montecarlo]
Index: /issm/trunk-jpl/test/NightlyRun/test418.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test418.py	(revision 23258)
+++ /issm/trunk-jpl/test/NightlyRun/test418.py	(revision 23258)
@@ -0,0 +1,43 @@
+#Test Name: SquareSheetShelfDiadSSA3dDakotaAreaAverage
+
+# this test may crash
+
+import numpy as np
+from model import *
+from socket import gethostname
+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')
+md = setflowequation(md,'SSA','all')
+md.cluster = generic('name',gethostname(),'np',3)
+
+#partitioning
+md.qmu.numberofpartitions = 100
+
+#corrupted size vs. prev_size
+#Aborted (core dumped)
+md = partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions)
+
+
+md.qmu.partition=md.qmu.partition-1
+
+vector = np.arange(1,1+md.mesh.numberofvertices,1).reshape(-1,1)
+# double check this before committing:
+#print 'before AreaAverageOntoPartition'
+vector_on_partition = AreaAverageOntoPartition(md,vector)
+vector_on_nodes = vector_on_partition[md.qmu.partition+1]
+
+field_names      = ['vector_on_nodes']
+field_tolerances = [1e-11]
+field_values     = [vector_on_nodes]