Index: /issm/trunk/src/m/model/modelextract.m
===================================================================
--- /issm/trunk/src/m/model/modelextract.m	(revision 9435)
+++ /issm/trunk/src/m/model/modelextract.m	(revision 9436)
@@ -86,14 +86,14 @@
 
 		%get field
-		field=md1.(model_fields(i));
+		field=md1.(model_fields{i});
 		fieldsize=size(field);
 
 		%size = number of nodes * n
 		if fieldsize(1)==numberofnodes1
-			md2.(model_fields(i))=field(pos_node,:);
+			md2.(model_fields{i})=field(pos_node,:);
 
 		%size = number of elements * n
 		elseif fieldsize(1)==numberofelements1
-			md2.(model_fields(i))=field(pos_elem,:);
+			md2.(model_fields{i})=field(pos_elem,:);
 		end
 	end
@@ -211,12 +211,4 @@
 	if ~isnan(md1.spctemperature),
 		md2.spctemperature(nodestoflag2,1)=1;
-		if ~isnan(md1.observed_temperature)
-			md2.spctemperature(nodestoflag2,2)=md2.observed_temperature(nodestoflag2); 
-		else
-			md2.spctemperature(nodestoflag2,2)=zeros(length(nodestoflag2),2);
-			disp(' ')
-			disp('!! modelextract warning: spc values should be checked !!')
-			disp(' ')
-		end
 	end
 
@@ -253,13 +245,4 @@
 	end
 
-	%reinitialize output parameters
-	md2.viscousheating=NaN;
-	md2.pressure_elem=NaN;
-	md2.stress=NaN;
-	md2.stress_surface=NaN;
-	md2.stress_bed=NaN;
-	md2.deviatoricstress=NaN;
-	md2.strainrate=NaN;
-
 %Keep track of pos_node and pos_elem
 md2.extractednodes=pos_node;