Index: /issm/trunk/src/m/solutions/ice/ModelProcessorThermal.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/ModelProcessorThermal.m	(revision 993)
+++ /issm/trunk/src/m/solutions/ice/ModelProcessorThermal.m	(revision 994)
@@ -52,6 +52,6 @@
 %the formulation into trias: 
 el6pos=find((md.elements_type(:,1)==macayealenum() | md.elements_type(:,1)==pattynenum()) & (element_partitioning==labindex));
-pos=find(ismember(el6pos,find(md.elements_type(:,1)==macayealenum())));
-elements(pos)=SetStructureField(elements(pos),'element','collapse',ones(length(pos),1));
+pos2=find(ismember(el6pos,find(md.elements_type(:,1)==macayealenum())));
+elements(pos2)=SetStructureField(elements(pos2),'element','collapse',ones(length(pos2),1));
 
 if strcmpi(solutiontype,'thermalsteady'),