%test partitioning, and partition averaging
md=triangle(model(),'../Exp/Square.exp',30000);
md=setmask(md,'../Exp/SquareShelf.exp','');
md=parameterize(md,'../Par/SquareSheetShelf.par');
md=setflowequation(md,'macayeal','all');
md.cluster=generic('name',oshostname(),'np',3);

%partitioning
md.qmu.numberofpartitions=100;
md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions);
md.qmu.partition=md.qmu.partition-1;

vector=(1:1:md.mesh.numberofvertices)';
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,...
	};