Index: /issm/trunk/src/m/utils/BC/SetIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetIceSheetBC.m	(revision 1274)
+++ /issm/trunk/src/m/utils/BC/SetIceSheetBC.m	(revision 1275)
@@ -35,5 +35,5 @@
 end
 
-disp('      boundary conditions for prognostic model initialization ');
+displaystring(md.debug,'%s',['      boundary conditions for prognostic model initialization']);
 md.gridondirichlet_prog=zeros(md.numberofgrids,1);
 md.dirichletvalues_prog=zeros(md.numberofgrids,1);
@@ -42,5 +42,5 @@
 
 if (length(md.observed_temperature)==md.numberofgrids),
-	disp('      boundary conditions for thermal model');
+	displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
 	md.gridondirichlet_thermal=ones(md.numberofgrids,1); %surface temperature
 	md.dirichletvalues_thermal=md.observed_temperature;
Index: /issm/trunk/src/m/utils/BC/SetIceShelfBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 1274)
+++ /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 1275)
@@ -48,5 +48,5 @@
 end
 
-disp('      boundary conditions for prognostic model initialization ');
+displaystring(md.debug,'%s',['      boundary conditions for prognostic model initialization']);
 md.gridondirichlet_prog=zeros(md.numberofgrids,1);
 md.dirichletvalues_prog=zeros(md.numberofgrids,1);
@@ -55,5 +55,5 @@
 
 if (length(md.observed_temperature)==md.numberofgrids),
-	disp('      boundary conditions for thermal model');
+	displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
 	md.gridondirichlet_thermal=ones(md.numberofgrids,1); %surface temperature
 	md.dirichletvalues_thermal=md.observed_temperature;
Index: /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 1274)
+++ /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 1275)
@@ -60,5 +60,5 @@
 end
 
-disp('      boundary conditions for prognostic model initialization ');
+displaystring(md.debug,'%s',['      boundary conditions for prognostic model initialization']);
 md.gridondirichlet_prog=zeros(md.numberofgrids,1);
 md.dirichletvalues_prog=zeros(md.numberofgrids,1);
@@ -67,5 +67,5 @@
 
 if (length(md.observed_temperature)==md.numberofgrids),
-	disp('      boundary conditions for thermal model');
+	displaystring(md.debug,'%s',['      boundary conditions for thermal model']);
 	md.gridondirichlet_thermal=ones(md.numberofgrids,1); %surface temperature
 	md.dirichletvalues_thermal=md.observed_temperature;