Changeset 5369

Timestamp:

08/18/10 11:44:39 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added validation test for Bamg

File:

• issm/trunk/test/NightlyRun/test1401.m (modified) (2 diffs)

issm/trunk/test/NightlyRun/test1401.m

@@ -1 +1 @@
 %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 +11 @@
 %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;