Index: /issm/trunk-jpl/src/m/contrib/hack/tres.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/hack/tres.m	(revision 13020)
+++ /issm/trunk-jpl/src/m/contrib/hack/tres.m	(revision 13021)
@@ -105,5 +105,5 @@
 		case MaterialsRheologyBEnum(), string='rheology_B'; return
 		case MaterialsRheologyBbarEnum(), string='rheology_B'; return
-		case BalancethicknessThickeningRateEnum: string='dhdt'; return
+		case BalancethicknessThickeningRateEnum(), string='dhdt'; return
 		case VxEnum(), string='vx'; return
 		case InversionVxObsEnum(), string='vx_obs'; return
Index: /issm/trunk-jpl/src/m/exp/expwrite.m
===================================================================
--- /issm/trunk-jpl/src/m/exp/expwrite.m	(revision 13020)
+++ /issm/trunk-jpl/src/m/exp/expwrite.m	(revision 13021)
@@ -17,5 +17,5 @@
 fid=fopen(filename,'w');
 for n=1:length(a),
-	if(length(contours(n).x)~=length(contours(n).y)),
+	if(length(a(n).x)~=length(a(n).y)),
 		error('contours x and y coordinates must be of identical size');
 	end
Index: /issm/trunk-jpl/src/m/mesh/roundmesh.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/roundmesh.m	(revision 13020)
+++ /issm/trunk-jpl/src/m/mesh/roundmesh.m	(revision 13021)
@@ -20,5 +20,5 @@
 y_list=radius*y_list.*sin(theta);
 A=struct('x',x_list,'y',y_list,'density',1);
-expwrite('RoundDomainOutline.exp',A);
+expwrite(A,'RoundDomainOutline.exp');
 
 %Call Bamg
Index: /issm/trunk-jpl/test/NightlyRun/test218.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test218.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test218.m	(revision 13021)
@@ -84,5 +84,5 @@
 
 %Fields and tolerances to track changes
-md=tres(md,'dakota');
+md.qmu.results=md.results.dakota;
 md.results.dakota.importancefactors=importancefactors(md,'scaled_MaterialsRheologyB','MaxVel')';
 field_names     ={'importancefactors'};
Index: /issm/trunk-jpl/test/NightlyRun/test234.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test234.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test234.m	(revision 13021)
@@ -56,5 +56,5 @@
 %solve
 md=solve(md,TransientSolutionEnum,'overwrite','y');
-md=tres(md,'dakota');
+md.qmu.results=md.results.dakota;
 
 %Fields and tolerances to track changes
Index: /issm/trunk-jpl/test/NightlyRun/test235.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test235.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test235.m	(revision 13021)
@@ -56,5 +56,5 @@
 %solve
 md=solve(md,TransientSolutionEnum,'overwrite','y');
-md=tres(md,'dakota');
+md.qmu.results=md.results.dakota;
 
 %Fields and tolerances to track changes
Index: /issm/trunk-jpl/test/NightlyRun/test412.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test412.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test412.m	(revision 13021)
@@ -36,5 +36,5 @@
 
 %Fields and tolerances to track changes
-md=tres(md,'dakota');
+md.qmu.results=md.results.dakota;
 md.results.dakota.importancefactors=importancefactors(md,'scaled_FrictionCoefficient','MaxVel')';
 field_names     ={'importancefactors'};
Index: /issm/trunk-jpl/test/NightlyRun/test413.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test413.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test413.m	(revision 13021)
@@ -37,5 +37,5 @@
 
 %Fields and tolerances to track changes
-md=tres(md,'dakota');
+md.qmu.results=md.results.dakota;
 md.results.dakota.importancefactors=importancefactors(md,'scaled_FrictionCoefficient','MaxVel')';
 field_names     ={'importancefactors'};
Index: /issm/trunk-jpl/test/NightlyRun/test414.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test414.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test414.m	(revision 13021)
@@ -49,5 +49,5 @@
 %solve
 md=solve(md,DiagnosticSolutionEnum,'overwrite','y');
-md=tres(md,'dakota');
+md.qmu.results=md.results.dakota;
 
 %Fields and tolerances to track changes
Index: /issm/trunk-jpl/test/NightlyRun/test417.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test417.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test417.m	(revision 13021)
@@ -53,5 +53,5 @@
 
 %Fields and tolerances to track changes
-md=tres(md,'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
Index: /issm/trunk-jpl/test/NightlyRun/test420.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test420.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test420.m	(revision 13021)
@@ -33,5 +33,5 @@
 %solve
 md=solve(md,DiagnosticSolutionEnum,'overwrite','y');
-md=tres(md,'dakota');
+md.qmu.results=md.results.dakota;
 
 %test on thickness
Index: /issm/trunk-jpl/test/NightlyRun/test511.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test511.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test511.m	(revision 13021)
@@ -8,5 +8,5 @@
 md=extrude(md,3,1);
 md=setflowequation(md,'stokes','all');
-md=modelextract(md,md.mask.elementonfloatingice);
+md=extract(md,md.mask.elementonfloatingice);
 
 %control parameters
Index: /issm/trunk-jpl/test/NightlyRun/test613.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test613.m	(revision 13020)
+++ /issm/trunk-jpl/test/NightlyRun/test613.m	(revision 13021)
@@ -6,5 +6,5 @@
 
 %Ice sheet only
-md=modelextract(md,md.mask.elementongroundedice);
+md=extract(md,md.mask.elementongroundedice);
 pos=find(md.mesh.vertexonboundary);
 md.balancethickness.spcthickness(pos)=md.geometry.thickness(pos);