Index: /issm/trunk/src/m/utils/BC/SetIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetIceSheetBC.m	(revision 1776)
+++ /issm/trunk/src/m/utils/BC/SetIceSheetBC.m	(revision 1777)
@@ -38,6 +38,5 @@
 if (length(md.observed_temperature)==md.numberofgrids),
 	displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
-	md.spctemperature=ones(md.numberofgrids,2); %observed temperauture everywhere
-	md.spctemperature(:,2)=md.observed_temperature; %values of the spcs
+	md.spctemperature=[md.gridonsurface md.observed_temperature]; %impose observed temperature on surface
 	if (length(md.geothermalflux)~=md.numberofgrids),
 		md.geothermalflux=50*10^-3*ones(md.numberofgrids,1); %50 mW/m^2
Index: /issm/trunk/src/m/utils/BC/SetIceShelfBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 1776)
+++ /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 1777)
@@ -51,6 +51,5 @@
 if (length(md.observed_temperature)==md.numberofgrids),
 	displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
-	md.spctemperature=ones(md.numberofgrids,2); %surface temperature
-	md.spctemperature(:,2)=md.observed_temperature; %values of the spcs
+	md.spctemperature=[md.gridonsurface md.observed_temperature]; %impose observed temperature on surface
 	if (length(md.geothermalflux)~=md.numberofgrids),
 		md.geothermalflux=zeros(md.numberofgrids,1);
Index: /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 1776)
+++ /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 1777)
@@ -63,6 +63,5 @@
 if (length(md.observed_temperature)==md.numberofgrids),
 	displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
-	md.spctemperature=ones(md.numberofgrids,2); %surface temperature
-	md.spctemperature(:,2)=md.observed_temperature; %values of the spcs
+	md.spctemperature=[md.gridonsurface md.observed_temperature]; %impose observed temperature on surface
 	if (length(md.geothermalflux)~=md.numberofgrids),
 		md.geothermalflux=zeros(md.numberofgrids,1);