Index: /issm/trunk/test/NightlyRun/test1201.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1201.m	(revision 5421)
+++ /issm/trunk/test/NightlyRun/test1201.m	(revision 5422)
@@ -1,3 +1,4 @@
 %This test is a test from the EISMINT for Ice shelves Vincent Rommelaere 1996
+printingflag=true;
 
 results={};
@@ -57,7 +58,8 @@
 plot(s,h1,'r',s,h2,'b',s,h3,'g',s,hth,'k')
 legend('Art. diff.','No Art. diff.','Discontinue','Theorique')
-
-set(gcf,'Color','w')
-printmodel('eismintmasscthickness','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('eismintmasscthickness','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+end
 
 %Fields and tolerances to track changes
Index: /issm/trunk/test/NightlyRun/test1202.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1202.m	(revision 5421)
+++ /issm/trunk/test/NightlyRun/test1202.m	(revision 5422)
@@ -1,3 +1,4 @@
 %Test on the diagnostic model and the prognostic in 2d
+printingflag=true;
 
 %tests 3 and 4: using Glen's flow law
@@ -15,13 +16,17 @@
 vy=PatchToVec(md.results.DiagnosticSolution.Vy);
 
-plotmodel(md,'data',vx,'contourlevels',{0,20,40,60,60,100,120,140,160,180,-20,-40,-60,-80,-100,-120,-140,-160,-180},...
-	'visible','off','contourcolor','k')
-set(gcf,'Color','w')
-printmodel('eismintdiag1vx','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+plotmodel(md,'data',vx,'contourlevels',{0,20,40,60,60,100,120,140,160,180,-20,-40,-60,-80,-100,-120,-140,-160,-180}, ...
+	'contourcolor','k')
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('eismintdiag1vx','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+end
 
 plotmodel(md,'data',vy,'contourlevels',{-100,-200,-300,-400,-500,-600,-700,-800,-900,-1000},...
-	'visible','off','contourcolor','k')
-set(gcf,'Color','w')
-printmodel('eismintdiag1vy','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+	'contourcolor','k')
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('eismintdiag1vy','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+end
 
 %Fields and tolerances to track changes
Index: /issm/trunk/test/NightlyRun/test1203.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1203.m	(revision 5421)
+++ /issm/trunk/test/NightlyRun/test1203.m	(revision 5422)
@@ -1,3 +1,4 @@
 %Test on the diagnostic model and the prognostic in 2d
+printingflag=true;
 
 %test 5 and 6 : 
@@ -19,12 +20,16 @@
 
 %plot results
-plotmodel(md,'data',vx,'contourlevels',{0,20,40,60,80,100,-20,-40,-60,-80,-100},'visible','off',...
+plotmodel(md,'data',vx,'contourlevels',{0,20,40,60,80,100,-20,-40,-60,-80,-100},...
 	'contourcolor','k')
-set(gcf,'Color','w')
-printmodel('eismintdiag2vx','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('eismintdiag2vx','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+end
 plotmodel(md,'data',vy,'contourlevels',{-100,-200,-300,-400,-500,-600,-700,-800,-900,-1000},...
-	'contourcolor','k','visible','off')
-set(gcf,'Color','w')
-printmodel('eismintdiag2vy','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+	'contourcolor','k')
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('eismintdiag2vy','png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
+end
 
 %Fields and tolerances to track changes
Index: /issm/trunk/test/NightlyRun/test1301.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1301.m	(revision 5421)
+++ /issm/trunk/test/NightlyRun/test1301.m	(revision 5422)
@@ -3,5 +3,5 @@
 % at upper surface and an impose flux at its base. The result must be a linear temperature from the upper to the lower
 % surface with an imposed slope (Geothermal flux). if it is not the case, something is thermal modeling has been changed...
-% Just run this file in Matlab, with a properly setup Ice code. 
+printingflag=true;
 
 md=model;
@@ -34,7 +34,9 @@
 plotmodel(md,'data',comp_melting,'title','Modeled melting','data',melting,'title','Analytical melting',...
 	'data',comp_melting-melting,'title','Absolute error','data',relative,'title','Relative error [%]',...
-	'layer#all',1,'caxis#2',[1.02964 1.02966]*10^-4,'FontSize#all',20,'figposition','mathieu','visible','off')
-set(gcf,'Color','w')
-printmodel('thermalmelting','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+	'layer#all',1,'caxis#2',[1.02964 1.02966]*10^-4,'FontSize#all',20,'figposition','mathieu')
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('thermalmelting','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+end
 
 
Index: /issm/trunk/test/NightlyRun/test1302.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1302.m	(revision 5421)
+++ /issm/trunk/test/NightlyRun/test1302.m	(revision 5422)
@@ -2,5 +2,5 @@
 % There is u=v=0 and w=cst everywhere the only thermal boundary conditions are an imposed temperature
 % at upper surface and an impose flux at its base.
-% Just run this file in Matlab, with a properly setup Ice code. 
+printingflag=true;
 
 md=model;
@@ -36,7 +36,9 @@
 	'title','Analytical temperature [K]','view',3,'data',comp_temp-md.observed_temperature,...
 	'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
-	'visible','off','figposition','mathieu','FontSize#all',20)
-set(gcf,'Color','w')
-printmodel('thermaladvection','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+	'figposition','mathieu','FontSize#all',20)
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('thermaladvection','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+end
 
 %Fields and tolerances to track changes
Index: /issm/trunk/test/NightlyRun/test1303.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1303.m	(revision 5421)
+++ /issm/trunk/test/NightlyRun/test1303.m	(revision 5422)
@@ -3,5 +3,5 @@
 % at the lower and upper surface. The result must be a linear temperature from the upper to the lower
 % surface. if it is not the case, something is thermal modeling has been changed...
-% Just run this file in Matlab, with a properly setup Ice code. 
+printingflag=true;
 
 md=model;
@@ -31,7 +31,9 @@
 	'title','Analytical temperature [K]','view',3,'data',comp_temp-md.observed_temperature,...
 	'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
-	'visible','off','figposition','mathieu','FontSize#all',20)
-set(gcf,'Color','w')
-printmodel('thermalconduction','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+	'figposition','mathieu','FontSize#all',20)
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('thermalconduction','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+end
 
 %Fields and tolerances to track changes
Index: /issm/trunk/test/NightlyRun/test1304.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1304.m	(revision 5421)
+++ /issm/trunk/test/NightlyRun/test1304.m	(revision 5422)
@@ -3,4 +3,5 @@
 % at upper surface and an impose flux at its base. The result must be a linear temperature from the upper to the lower
 % surface with an imposed slope (Geothermal flux). if it is not the case, something is thermal modeling has been changed...
+printingflag=true;
 
 md=model;
@@ -31,7 +32,9 @@
 	'title','Analytical temperature','view',3,'data',comp_temp-md.observed_temperature,...
 	'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
-	'visible','off','figposition','mathieu','FontSize#all',20)
-set(gcf,'Color','w')
-printmodel('thermalgeothermalflux','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+	'figposition','mathieu','FontSize#all',20)
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('thermalgeothermalflux','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');
+end
 
 %Fields and tolerances to track changes