Index: sm/trunk/src/m/solutions/processresults.m
===================================================================
--- /issm/trunk/src/m/solutions/processresults.m	(revision 9362)
+++ 	(revision )
@@ -1,157 +1,0 @@
-function newresults=processresults(models,results)
-%PROCESSRESULTS - process results from core solution
-%
-%   The solution (vel,pressure,...) are on the g-set,
-%   use index to get the value on each node
-%
-%   Usage:
-%      results=processresults(results,models,analysis_type)
-
-%initialize output
-newresults=struct();
-
-%get analysis_type:
-analysis_type=models.analysis_type;
-
-%recover models first
-if (analysis_type==DiagnosticAnalysisEnum | analysis_type==TransientAnalysisEnum | analysis_type==SteadystateAnalysisEnum()),
-	m_dh=models.dh;
-	m_ds=models.ds;
-	m_dhu=models.dhu;
-
-	%some flags needed
-	dim=m_dhu.parameters.Dim;
-	ishutter=m_dhu.parameters.Ishutter;
-	ismacayealpattyn=m_dh.parameters.Ismacayealpattyn;
-	isstokes=m_ds.parameters.Isstokes;
-end
-if (analysis_type==ControlAnalysisEnum()),
-	m_dh=models.dh;
-
-	%some flags needed
-	dim=m_dh.parameters.Dim;
-	ishutter=m_dh.parameters.Ishutter;
-	ismacayealpattyn=m_dh.parameters.Ismacayealpattyn;
-	isstokes=m_dh.parameters.Isstokes;
-end
-if (analysis_type==ThermalAnalysisEnum() | analysis_type==TransientAnalysisEnum() | analysis_type==SteadystateAnalysisEnum()),
-	m_m=models.m;
-end
-
-
-%go through results
-for i=1:length(results),
-
-	%get field names of results
-	resultsname=fieldnames(results(i));
-
-	%go through all field
-	for j=1:length(resultsname),
-
-		if strcmpi(resultsname{j},'u_g'),
-
-			u_g=results(i).u_g;
-			yts=m_dh.parameters.Yts;
-
-			if dim==2,
-				newresults(i).step=results(i).step;
-				newresults(i).time=results(i).time;
-				newresults(i).vx=u_g(1:2:end)*yts;
-				newresults(i).vy=u_g(2:2:end)*yts;
-				newresults(i).vel=sqrt(newresults(i).vx.^2+newresults(i).vy.^2);
-
-			else
-				newresults(i).step=results(i).step;
-				newresults(i).time=results(i).time;
-				if isstokes,
-					newresults(i).vx=u_g(1:4:end)*yts;
-					newresults(i).vy=u_g(2:4:end)*yts;
-					newresults(i).vz=u_g(3:4:end)*yts;
-				else
-					newresults(i).vx=u_g(1:3:end)*yts;
-					newresults(i).vy=u_g(2:3:end)*yts;
-					newresults(i).vz=u_g(3:3:end)*yts;
-				end
-				newresults(i).vel=sqrt(newresults(i).vx.^2+newresults(i).vy.^2+newresults(i).vz.^2);
-			end
-
-		elseif strcmpi(resultsname{j},'p_g'),
-
-			p_g=results(i).p_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).pressure=p_g;
-
-		elseif strcmpi(resultsname{j},'t_g'),
-
-			t_g=results(i).t_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).temperature=t_g;
-
-		elseif strcmpi(resultsname{j},'m_g'),
-
-			m_g=results(i).m_g;
-			yts=m_m.parameters.Yts;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).melting=m_g*yts;
-
-		elseif strcmpi(resultsname{j},'h_g'),
-
-			h_g=results(i).h_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).thickness=h_g;
-
-		elseif strcmpi(resultsname{j},'s_g'),
-
-			s_g=results(i).s_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).surface=s_g;
-
-		elseif strcmpi(resultsname{j},'b_g'),
-
-			b_g=results(i).b_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).bed=b_g;
-
-		elseif strcmpi(resultsname{j},'param_g'),
-
-			param_g=results(i).param_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).parameter=param_g;
-
-		elseif strcmpi(resultsname{j},'v_g'),
-
-			v_g=results(i).v_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).vel=v_g;
-
-		elseif strcmpi(resultsname{j},'sx_g'),
-
-			sy_g=results(i).sy_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).slopey=sy_g;
-
-		elseif strcmpi(resultsname{j},'sy_g'),
-
-			sx_g=results(i).sx_g;
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).slopex=sx_g;
-
-		else
-
-			newresults(i).step=results(i).step;
-			newresults(i).time=results(i).time;
-			newresults(i).(resultsname{j})=results(i).(resultsname{j});
-
-		end
-	end
-end