Index: /issm/trunk/src/m/classes/public/ReadData.m
===================================================================
--- /issm/trunk/src/m/classes/public/ReadData.m	(revision 4201)
+++ /issm/trunk/src/m/classes/public/ReadData.m	(revision 4202)
@@ -27,5 +27,7 @@
 	elseif type==3,
 		N=fread(fid,1,'int');
-		field=fread(fid,M*N,'double');
+		field=fread(fid,[M N],'double');
+	else
+		error(['cannot read data of type ' num2str(type) ]);
 	end
 
Index: /issm/trunk/src/m/classes/public/isresultconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/isresultconsistent.m	(revision 4201)
+++ /issm/trunk/src/m/classes/public/isresultconsistent.m	(revision 4202)
@@ -26,5 +26,5 @@
 
 %DIAGNOSTIC
-if analysis_type==DiagnosticAnalysisEnum() & md.control_analysis==0,
+if analysis_type==DiagnosticSolutionEnum() & md.control_analysis==0,
 
 	if (md.dim==3)
@@ -59,5 +59,5 @@
 
 %CONTROL
-if analysis_type==DiagnosticAnalysisEnum() & md.control_analysis==1,
+if analysis_type==DiagnosticSolutionEnum() & md.control_analysis==1,
 
 	if ~md.cm_gradient,
@@ -118,5 +118,5 @@
 
 %THERMAL
-if analysis_type==ThermalAnalysisEnum(),
+if analysis_type==ThermalSolutionEnum(),
 
 	for iter=1:length(md.results.ThermalAnalysis)
@@ -154,5 +154,5 @@
 end
 
-if analysis_type==TransientAnalysisEnum(),
+if (analysis_type==Transient2DSolutionEnum() | analysis_type==Transient3DSolutionEnum()),
 
 	for iter=1:length(md.results.TransientAnalysis)
Index: /issm/trunk/src/m/classes/public/loadresultsfromdisk.m
===================================================================
--- /issm/trunk/src/m/classes/public/loadresultsfromdisk.m	(revision 4201)
+++ /issm/trunk/src/m/classes/public/loadresultsfromdisk.m	(revision 4202)
@@ -20,8 +20,8 @@
 	%load results onto model
 	structure=parseresultsfromdisk(filename);
-	md.results.(structure.analysis_type)=structure;
+	md.results.(structure.SolutionType)=structure;
 
 	%recover analysis_type from results
-	md.analysis_type=structure(1).analysis_type;
+	md.analysis_type=structure(1).SolutionType;
 	if isscalar(md.sub_analysis_type),
 		md.sub_analysis_type=EnumAsString(md.sub_analysis_type);
@@ -30,5 +30,5 @@
 	%Check result is consistent, only if it exists
 	disp(sprintf('%s\n','checking result consistency'));
-	if ~isresultconsistent(md,structure(1).analysis_type),
+	if ~isresultconsistent(md,structure(1).SolutionType),
 		%it would be very cruel to use an error, it would kill the computed results (although they are not consistent...)
 		disp('!! results not consistent correct the model !!')