Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m	(revision 330)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m	(revision 331)
@@ -15,8 +15,5 @@
 	soln(count).u_f=[];
 
-	if m.parameters.debug,
-		disp(sprintf('\n%s',['   starting direct shooting method']));
-	end
-	
+	displaystring(m.parameters.debug,'\n%s',['   starting direct shooting method']);
 	while(~converged),
 		
@@ -49,7 +46,5 @@
 
 		%Solve	
-		if m.parameters.debug,
-			disp(sprintf('%s%g','      condition number of stiffness matrix: ',condest(K_ff)));
-		end
+		displaystring(m.parameters.debug,'\n%s%g','      condition number of stiffness matrix: ',condest(K_ff));
 		[soln(count).u_f]=Solver(K_ff,p_f,[],m.parameters);
 	
@@ -75,8 +70,8 @@
 		%Relative criterion
 		if (ndu/nu<=m.parameters.eps_rel),
-			if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: norm(du)/norm(u)=',ndu/nu,' < ',m.parameters.eps_rel)); end
+			displaystring(m.parameters.debug,'%s%g%s%g','      convergence criterion: norm(du)/norm(u)=',ndu/nu,' < ',m.parameters.eps_rel);
 			converged=1;
 		else
-			if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: norm(du)/norm(u)=',ndu/nu,' > ',m.parameters.eps_rel)); end
+			displaystring(m.parameters.debug,'%s%g%s%g','      convergence criterion: norm(du)/norm(u)=',ndu/nu,' > ',m.parameters.eps_rel);
 			converged=0;
 		end
@@ -85,8 +80,8 @@
 		if ~isnan(m.parameters.eps_abs),
 			if (nduinf<=m.parameters.eps_abs),
-				if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: max(du)=',nduinf,' < ',m.parameters.eps_abs)); end
+				displaystring(m.parameters.debug,'%s%g%s%g','      convergence criterion: max(du)=',nduinf,' < ',m.parameters.eps_abs);
 			else
-				if m.parameters.debug, disp(sprintf('%s %g %s %g','      convergence criterion: max(du)=',nduinf,' > ',m.parameters.eps_abs)); end
-					converged=0;
+				displaystring(m.parameters.debug,'%s%g%s%g','      convergence criterion: max(du)=',nduinf,' > ',m.parameters.eps_abs);
+				converged=0;
 			end
 		end
@@ -117,4 +112,3 @@
 		varargout(2)={K_fs};
 	end
-	
 end