Index: /issm/trunk/test/NightlyRun/test1401.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1401.m	(revision 5368)
+++ /issm/trunk/test/NightlyRun/test1401.m	(revision 5369)
@@ -1,4 +1,5 @@
 %test the anisotropic mesh adaptation
 %function to capture = exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
+printingflag=true;
 
 %create square mesh
@@ -10,42 +11,57 @@
 %mesh adaptation loop YAMS
 md=squaremesh(md,L,L,nx,ny);
-plotmodel(md,'data','mesh');pause(0.5);
-for i=1:1
-	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
-	md=YamsCall(md,md.vel_obs,0.001,0.3,1.3,10^-4);
-	plotmodel(md,'data','mesh');pause(0.5);
+md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
+plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('mesh1_yams1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh1_yams1.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:1
-	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
-	md=YamsCall(md,md.vel_obs,0.001,0.3,2,10^-4);
-	plotmodel(md,'data','mesh');pause(0.5);
+
+md=YamsCall(md,md.vel_obs,0.001,0.3,1.3,10^-4);
+md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
+plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('mesh1_yams2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh1_yams2.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:3,
-	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
-	md=YamsCall(md,md.vel_obs,0.001,0.3,2,0.005);
-	plotmodel(md,'data','mesh');pause(0.5);
+
+md=YamsCall(md,md.vel_obs,0.001,0.3,2.5,0.008);
+md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
+plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('mesh1_yams3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh1_yams3.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
 
 %mesh adaptation loop BAMG
 md=squaremesh(md,L,L,nx,ny);
-plotmodel(md,'data','mesh');pause(0.5);
-for i=1:1
-	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
-	md.bamg=NaN;
-	md=bamg(md,'field',md.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',1.3,'err',10^-4);
-	plotmodel(md,'data','mesh');pause(0.5);
+md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
+plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('mesh1_bamg1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh1_bamg1.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:1
-	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
-	md.bamg=NaN;
-	md=bamg(md,'field',md.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',2,'err',10^-4);
-	plotmodel(md,'data','mesh');pause(0.5);
+
+md.bamg=NaN;
+md=bamg(md,'field',md.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',1.3,'err',10^-4);
+md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
+plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('mesh1_bamg2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh1_bamg2.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:3,
-	md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
-	md.bamg=NaN;
-	md=bamg(md,'field',md.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',2,'err',0.005);
-	plotmodel(md,'data','mesh');pause(0.5);
+
+md.bamg=NaN;
+md=bamg(md,'field',md.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',2.5,'err',0.008);
+md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;
+plotmodel(md,'data','vel_obs','data','vel_obs','nlines',1,'ncols',2,'title','','figposition',[500 500 1000 500],'axis#all','equal','xlim#all',[0 1],'ylim#all',[0 1],'edgecolor#1','w'); pause(0.5);
+if printingflag, 
+	set(gcf,'Color','w')
+	printmodel('mesh1_bamg3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh1_bamg3.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-md.vel_obs=exp(-(sqrt((md.x+0.1).^2+(md.y+0.1).^2)-0.75).^2*10^6)+((md.x+0.1).^2+(md.y+0.1).^2)/2;