Index: /issm/trunk/test/NightlyRun/test1402.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1402.m	(revision 5409)
+++ /issm/trunk/test/NightlyRun/test1402.m	(revision 5410)
@@ -1,76 +1,97 @@
 %test the anisotropic mesh adaptation
+printingflag=true;
 
 %create square mesh
 L=1; %in m
-nx=70; %numberof nodes in x direction
-ny=70;
+nx=30; %numberof nodes in x direction
+ny=30;
 md=model;
 
 %mesh adaptation loop YAMS
 md=squaremesh(md,L,L,nx,ny);
-plotmodel(md,'data','mesh');pause(0.5);
-for i=1:1
-	u=4*md.x-2; v=4*md.y-2;
-	md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
-		+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
-		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
-	md=YamsCall(md,md.vel_obs,0.005,0.3,1.3,10^-4);
-	plotmodel(md,'data','mesh');pause(0.5);
+u=4*md.x-2; v=4*md.y-2;
+md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
+	+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+	+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
+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('mesh2_yams1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh2_yams1.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:1
-	u=4*md.x-2; v=4*md.y-2;
-	md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
-		+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
-		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
-	md=YamsCall(md,md.vel_obs,0.005,0.3,2,10^-4);
-	plotmodel(md,'data','mesh');pause(0.5);
+
+md=YamsCall(md,md.vel_obs,0.005,0.3,2.3,10^-2);
+u=4*md.x-2; v=4*md.y-2;
+md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
+	+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+	+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
+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('mesh2_yams2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh2_yams2.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:3,
-	u=4*md.x-2; v=4*md.y-2;
-	md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
-		+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
-		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
-	md=YamsCall(md,md.vel_obs,0.005,0.3,2,0.001);
-	plotmodel(md,'data','mesh');pause(0.5);
+
+md=YamsCall(md,md.vel_obs,0.005,0.3,3,0.005);
+u=4*md.x-2; v=4*md.y-2;
+md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
+	+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+	+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
+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('mesh2_yams3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh2_yams3.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-md.vel_obs=tanh(30*((4*md.x-2).^2+(4*md.y-2).^2-0.25))+...
-	tanh(30*((4*md.x-2.75).^2+(4*md.y-1.25).^2-2.75))+...
-	tanh(30*((4*md.x-2.75).^2+(4*md.y-2).^2-1.25))+...
-	tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-2.75))+...
-	tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-1.25));
 
 %mesh adaptation loop BAMG
 md=squaremesh(md,L,L,nx,ny);
-plotmodel(md,'data','mesh');pause(0.5);
-for i=1:1
-	u=4*md.x-2; v=4*md.y-2;
-	md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
-		+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
-		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
-	md.bamg=NaN;
-	md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',1.3,'err',10^-4);
-	plotmodel(md,'data','mesh');pause(0.5);
+u=4*md.x-2; v=4*md.y-2;
+md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
+	+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+	+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
+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('mesh2_bamg1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh2_bamg1.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:1
-	u=4*md.x-2; v=4*md.y-2;
-	md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
-		+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
-		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
-	md.bamg=NaN;
-	md=bamg(md,'field',md.vel_obs,'hmin',0.005,'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.005,'hmax',0.3,'gradation',2.3,'err',10^-2);
+u=4*md.x-2; v=4*md.y-2;
+md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
+	+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+	+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
+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('mesh2_bamg2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh2_bamg2.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-for i=1:3,
-	u=4*md.x-2; v=4*md.y-2;
-	md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
-		+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
-		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
-	md.bamg=NaN;
-	md=bamg(md,'field',md.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',2,'err',10^-3);
-	plotmodel(md,'data','mesh');pause(0.5);
+
+md.bamg=NaN;
+md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',3,'err',0.005);
+u=4*md.x-2; v=4*md.y-2;
+md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
+	+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+	+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
+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('mesh2_bamg3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh2_bamg3.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
 end
-md.vel_obs=tanh(30*((4*md.x-2).^2+(4*md.y-2).^2-0.25))+...
-	tanh(30*((4*md.x-2.75).^2+(4*md.y-1.25).^2-2.75))+...
-	tanh(30*((4*md.x-2.75).^2+(4*md.y-2).^2-1.25))+...
-	tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-2.75))+...
-	tanh(30*((4*md.x-1.25).^2+(4*md.y-2).^2-1.25));
+
+md.bamg=NaN;
+md=bamg(md,'field',md.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',1.5,'err',0.003,'anisomax',1);
+u=4*md.x-2; v=4*md.y-2;
+md.vel_obs=tanh(30*(u.^2+v.^2-0.25)) ...
+	+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+	+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;
+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('mesh2_bamgiso','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
+	system(['mv mesh2_bamgiso.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/Mesh/ ']);
+end