Index: /issm/trunk-jpl/src/m/contrib/hack/tres.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/hack/tres.m	(revision 13805)
+++ /issm/trunk-jpl/src/m/contrib/hack/tres.m	(revision 13806)
@@ -22,13 +22,7 @@
 		md.initialization.pressure=md.results.DiagnosticSolution.Pressure;
 	end
-	if md.rifts.numrifts,
+	if ~isempty(md.rifts.riftstruct),
 		if isfield(md.results.DiagnosticSolution,'riftproperties'),
 			md.rifts.riftproperties=md.results.DiagnosticSolution.riftproperties;
-		end
-	end
-	if md.inversion.iscontrol==1,
-		for control_parameters=md.inversion.control_parameters
-			%Will need to be updated... good luck ;)
-			md.(EnumToModelField(control_parameters))=md.results.DiagnosticSolution.(EnumToString(control_parameters));
 		end
 	end
Index: /issm/trunk-jpl/src/m/mesh/rifts/meshprocessoutsiderifts.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/rifts/meshprocessoutsiderifts.m	(revision 13805)
+++ /issm/trunk-jpl/src/m/mesh/rifts/meshprocessoutsiderifts.m	(revision 13806)
@@ -7,5 +7,5 @@
 
 %go through rifts, and figure out which ones touch the domain outline
-for i=1:md.rifts.numrifts,
+for i=1:length(md.rifts.riftstruct),
 
 	%first, flag nodes that belong to the domain outline
@@ -84,5 +84,4 @@
 md.mesh.z=zeros(md.mesh.numberofvertices,1);
 md.mesh.vertexonboundary=zeros(length(md.mesh.x),1); md.mesh.vertexonboundary(md.mesh.segments(:,1:2))=1;
-md.rifts.numrifts=length(md.rifts.riftstruct);
 md.flowequation.element_equation=3*ones(md.mesh.numberofelements,1);
 md.mesh.vertexonbed=ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/src/m/mesh/rifts/rifttipsonmesh.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/rifts/rifttipsonmesh.m	(revision 13805)
+++ /issm/trunk-jpl/src/m/mesh/rifts/rifttipsonmesh.m	(revision 13806)
@@ -4,5 +4,5 @@
 
 %read rifts from outline file
-rifts=expread(riftoutline,1);
+rifts=expread(riftoutline);
 
 tips=[];