Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 801)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 802)
@@ -67,6 +67,6 @@
 	'tolx','np','eps_rel','exclusive','n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
 for i=1:length(fields),
-	if ~isempty(eval(['md.' fields{i}])),
-		if any(isnan(eval(['md.' fields{i}]))),
+	if ~isempty(md.(fields{i})),
+		if any(isnan(md.(fields{i}))),
 			disp(['model ' md.name ' has an NaN value in field ' fields{i} '!']);
 			bool=0; return;
@@ -80,6 +80,6 @@
 	'sparsity','lowmem','n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
 for i=1:length(fields),
-	if ~isempty(eval(['md.' fields{i}])),
-		if any((eval(['md.' fields{i}]))<0),
+	if ~isempty(md.(fields{i})),
+		if any(md.(fields{i})<0),
 			disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
 			bool=0; return;
@@ -97,6 +97,6 @@
 	'sparsity','deltaH','DeltaH','timeacc','timedec'};
 for i=1:length(fields),
-	if ~isempty(eval(['md.' fields{i}])),
-		if any((eval(['md.' fields{i}]))==0),
+	if ~isempty(md.(fields{i})),
+		if any(md.(fields{i})==0),
 			disp(['model ' md.name ' has a =0 value in field ' fields{i} '!']);
 			bool=0; return;
@@ -108,5 +108,5 @@
 fields={'elements','p','q','elementoniceshelf','n','elementonbed'};
 for i=1:size(fields,2),
-	if (size(eval(['md.' fields{i}]),1)~=md.numberofelements),
+	if (size(md.(fields{i}),1)~=md.numberofelements),
 		disp(['model ' md.name ' field ' fields{i} ' should be of size ' num2str(md.numberofelements) '!']);
 		bool=0; return;
@@ -117,5 +117,5 @@
 fields={'x','y','z','B','drag','gridondirichlet_diag','melting','accumulation','surface','thickness','bed','gridonbed','gridonsurface'};
 for i=1:length(fields),
-	if length(eval(['md.' fields{i}]))~=md.numberofgrids,
+	if length(md.(fields{i}))~=md.numberofgrids,
 		disp(['model ' md.name ' field ' fields{i} ' should be of size ' num2str(md.numberofgrids) '!']);
 		bool=0; return;
@@ -249,5 +249,5 @@
 	fields={'dt','ndt'};
 	for i=1:length(fields),
-		if any((eval(['md.' fields{i}]))<0),
+		if any(md.(fields{i})<0),
 			disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
 			bool=0; return;
@@ -320,5 +320,5 @@
 	fields={'vx_obs','vy_obs'};
 	for i=1:length(fields),
-		if any(length(eval(['md.' fields{i}]))~=md.numberofgrids),
+		if any(length(md.(fields{i}))~=md.numberofgrids),
 			disp(['model ' md.name ' field ' fields{i} ' should be of size ' num2str(md.numberofgrids) '!']);
 			bool=0; return;
@@ -362,6 +362,6 @@
 	fields={'sparsity'};
 	for i=1:length(fields),
-		if ~isempty(eval(['md.' fields{i}])),
-			if any(isnan(eval(['md.' fields{i}]))),
+		if ~isempty(md.(fields{i})),
+			if any(isnan(md.(fields{i}))),
 				disp(['model ' md.name ' has an NaN value in field ' fields{i} '!']);
 				bool=0; return;
@@ -373,6 +373,6 @@
 	fields={'sparsity'};
 	for i=1:length(fields),
-		if ~isempty(eval(['md.' fields{i}])),
-			if any((eval(['md.' fields{i}]))<0),
+		if ~isempty(md.(fields{i})),
+			if any(md.(fields{i})<0),
 				disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
 				bool=0; return;
@@ -384,6 +384,6 @@
 	fields={'sparsity'};
 	for i=1:length(fields),
-		if ~isempty(eval(['md.' fields{i}])),
-			if any((eval(['md.' fields{i}]))==0),
+		if ~isempty(md.(fields{i})),
+			if any(md.(fields{i})==0),
 				disp(['model ' md.name ' has a =0 value in field ' fields{i} '!']);
 				bool=0; return;
@@ -435,6 +435,6 @@
 		fields={'time','np'};
 		for i=1:length(fields),
-			if ~isempty(eval(['md.' fields{i}])),
-				if any(isnan(eval(['md.' fields{i}]))),
+			if ~isempty(md.(fields{i})),
+				if any(isnan(md.(fields{i}))),
 					disp(['model ' md.name ' has an NaN value in field ' fields{i} '!']);
 					bool=0; return;
@@ -446,6 +446,6 @@
 		fields={'time','np'};
 		for i=1:length(fields),
-			if ~isempty(eval(['md.' fields{i}])),
-				if any((eval(['md.' fields{i}]))<0),
+			if ~isempty(md.(fields{i})),
+				if any(md.(fields{i})<0),
 					disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
 					bool=0; return;
@@ -457,6 +457,6 @@
 		fields={'time','np'};
 		for i=1:length(fields),
-			if ~isempty(eval(['md.' fields{i}])),
-				if any((eval(['md.' fields{i}]))==0),
+			if ~isempty(md.(fields{i})),
+				if any(md.(fields{i})==0),
 					disp(['model ' md.name ' has a =0 value in field ' fields{i} '!']);
 					bool=0; return;
Index: /issm/trunk/src/m/classes/public/parseresultsfromdisk.m
===================================================================
--- /issm/trunk/src/m/classes/public/parseresultsfromdisk.m	(revision 801)
+++ /issm/trunk/src/m/classes/public/parseresultsfromdisk.m	(revision 802)
@@ -17,7 +17,7 @@
 while ~isempty(result),
 
-	eval(['results(' num2str(result.step) ').step=result.step;']);
-	eval(['results(' num2str(result.step) ').time=result.time;']);
-	eval(['results(' num2str(result.step) ').' result.fieldname '=result.field;']);
+	results(result.step).step=result.step;
+	results(result.step).time=result.time;
+	results(result.step).(result.fieldname)=result.field;
 	result=ReadData(fid);