Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 805)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 806)
@@ -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;
@@ -258,5 +258,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;
@@ -329,5 +329,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;
@@ -371,6 +371,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;
@@ -382,6 +382,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;
@@ -393,6 +393,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;
@@ -444,6 +444,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;
@@ -455,6 +455,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;
@@ -466,6 +466,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;