Index: /issm/trunk/test/Par/SquareThermal.par
===================================================================
--- /issm/trunk/test/Par/SquareThermal.par	(revision 9045)
+++ /issm/trunk/test/Par/SquareThermal.par	(revision 9046)
@@ -42,6 +42,5 @@
 disp('      boundary conditions for thermal model');
 md.temperature=md.observed_temperature;
-md.spctemperature(:,1)=ones(md.numberofnodes,1); %surface temperature
-md.spctemperature(:,2)=md.observed_temperature;
+md.spctemperature(:)=md.observed_temperature;
 md.geothermalflux=zeros(md.numberofnodes,1); 
 pos=find(md.elementonicesheet);md.geothermalflux(md.elements(pos,:))=1*10^-3; %1 mW/m^2