![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 14310 for issm/trunk/test/NightlyRun/test1402.m
- Timestamp:
- 02/04/13 08:01:04 (12 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 14068-14070,14073-14074,14077-14095,14097-14112,14117,14134-14135,14138-14142,14144,14149-14151,14153,14156-14208,14210-14308
- Property svn:mergeinfo changed
-
issm/trunk/test
- Property svn:mergeinfo changed
/issm/trunk-jpl/test merged: 14068,14091,14102,14106-14107,14134-14135,14138-14139,14150,14156-14158,14169,14179-14181,14190,14198,14212,14269-14270
- Property svn:mergeinfo changed
-
issm/trunk/test/NightlyRun/test1402.m
r13975 r14310 3 3 4 4 %create square mesh 5 L=1 ; %in m5 L=1.; %in m 6 6 nx=30; %numberof nodes in x direction 7 7 ny=30; … … 10 10 %mesh adaptation loop YAMS 11 11 md=squaremesh(md,L,L,nx,ny); 12 u=4 *md.mesh.x-2; v=4*md.mesh.y-2;13 md.inversion.vel_obs=tanh(30 *(u.^2+v.^2-0.25)) ...14 +tanh(30 *((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...15 +tanh(30 *((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25));12 u=4.*md.mesh.x-2.; v=4.*md.mesh.y-2.; 13 md.inversion.vel_obs=tanh(30.*(u.^2+v.^2-0.25)) ... 14 +tanh(30.*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u-0.75).^2+(v+0.75).^2-0.25)) ... 15 +tanh(30.*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u+0.75).^2+(v+0.75).^2-0.25)); 16 16 plotmodel(md,'data',md.inversion.vel_obs,'data',md.inversion.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); 17 if printingflag, 17 if printingflag, 18 18 set(gcf,'Color','w') 19 19 printmodel('mesh2_yams1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off'); … … 21 21 end 22 22 23 md=YamsCall(md,md.inversion.vel_obs,0.005,0.3,2.3,10 ^-2);24 u=4 *md.mesh.x-2; v=4*md.mesh.y-2;25 md.inversion.vel_obs=tanh(30 *(u.^2+v.^2-0.25)) ...26 +tanh(30 *((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...27 +tanh(30 *((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25));23 md=YamsCall(md,md.inversion.vel_obs,0.005,0.3,2.3,10.^-2); 24 u=4.*md.mesh.x-2.; v=4.*md.mesh.y-2.; 25 md.inversion.vel_obs=tanh(30.*(u.^2+v.^2-0.25)) ... 26 +tanh(30.*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u-0.75).^2+(v+0.75).^2-0.25)) ... 27 +tanh(30.*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u+0.75).^2+(v+0.75).^2-0.25)); 28 28 plotmodel(md,'data',md.inversion.vel_obs,'data',md.inversion.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); 29 if printingflag, 29 if printingflag, 30 30 set(gcf,'Color','w') 31 31 printmodel('mesh2_yams2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off'); … … 34 34 35 35 md=YamsCall(md,md.inversion.vel_obs,0.005,0.3,3,0.005); 36 u=4 *md.mesh.x-2; v=4*md.mesh.y-2;37 md.inversion.vel_obs=tanh(30 *(u.^2+v.^2-0.25)) ...38 +tanh(30 *((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...39 +tanh(30 *((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25));36 u=4.*md.mesh.x-2.; v=4.*md.mesh.y-2.; 37 md.inversion.vel_obs=tanh(30.*(u.^2+v.^2-0.25)) ... 38 +tanh(30.*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u-0.75).^2+(v+0.75).^2-0.25)) ... 39 +tanh(30.*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u+0.75).^2+(v+0.75).^2-0.25)); 40 40 plotmodel(md,'data',md.inversion.vel_obs,'data',md.inversion.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); 41 if printingflag, 41 if printingflag, 42 42 set(gcf,'Color','w') 43 43 printmodel('mesh2_yams3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off'); … … 49 49 %mesh adaptation loop BAMG 50 50 md=squaremesh(md,L,L,nx,ny); 51 u=4 *md.mesh.x-2; v=4*md.mesh.y-2;52 md.inversion.vel_obs=tanh(30 *(u.^2+v.^2-0.25)) ...53 +tanh(30 *((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...54 +tanh(30 *((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25));51 u=4.*md.mesh.x-2.; v=4.*md.mesh.y-2.; 52 md.inversion.vel_obs=tanh(30.*(u.^2+v.^2-0.25)) ... 53 +tanh(30.*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u-0.75).^2+(v+0.75).^2-0.25)) ... 54 +tanh(30.*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u+0.75).^2+(v+0.75).^2-0.25)); 55 55 plotmodel(md,'data',md.inversion.vel_obs,'data',md.inversion.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); 56 if printingflag, 56 if printingflag, 57 57 set(gcf,'Color','w') 58 58 printmodel('mesh2_bamg1','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off'); … … 61 61 62 62 md.private.bamg=NaN; 63 md=bamg(md,'field',md.inversion.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',2.3,'err',10 ^-2);64 u=4 *md.mesh.x-2; v=4*md.mesh.y-2;65 md.inversion.vel_obs=tanh(30 *(u.^2+v.^2-0.25)) ...66 +tanh(30 *((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...67 +tanh(30 *((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25));63 md=bamg(md,'field',md.inversion.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',2.3,'err',10.^-2); 64 u=4.*md.mesh.x-2.; v=4.*md.mesh.y-2.; 65 md.inversion.vel_obs=tanh(30.*(u.^2+v.^2-0.25)) ... 66 +tanh(30.*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u-0.75).^2+(v+0.75).^2-0.25)) ... 67 +tanh(30.*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u+0.75).^2+(v+0.75).^2-0.25)); 68 68 plotmodel(md,'data',md.inversion.vel_obs,'data',md.inversion.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); 69 if printingflag, 69 if printingflag, 70 70 set(gcf,'Color','w') 71 71 printmodel('mesh2_bamg2','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off'); … … 75 75 md.private.bamg=NaN; 76 76 md=bamg(md,'field',md.inversion.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',3,'err',0.005); 77 u=4 *md.mesh.x-2; v=4*md.mesh.y-2;78 md.inversion.vel_obs=tanh(30 *(u.^2+v.^2-0.25)) ...79 +tanh(30 *((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...80 +tanh(30 *((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25));77 u=4.*md.mesh.x-2.; v=4.*md.mesh.y-2.; 78 md.inversion.vel_obs=tanh(30.*(u.^2+v.^2-0.25)) ... 79 +tanh(30.*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u-0.75).^2+(v+0.75).^2-0.25)) ... 80 +tanh(30.*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u+0.75).^2+(v+0.75).^2-0.25)); 81 81 plotmodel(md,'data',md.inversion.vel_obs,'data',md.inversion.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); 82 if printingflag, 82 if printingflag, 83 83 set(gcf,'Color','w') 84 84 printmodel('mesh2_bamg3','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off'); … … 88 88 md.private.bamg=NaN; 89 89 md=bamg(md,'field',md.inversion.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',1.5,'err',0.003,'anisomax',1); 90 u=4 *md.mesh.x-2; v=4*md.mesh.y-2;91 md.inversion.vel_obs=tanh(30 *(u.^2+v.^2-0.25)) ...92 +tanh(30 *((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...93 +tanh(30 *((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25));90 u=4.*md.mesh.x-2.; v=4.*md.mesh.y-2.; 91 md.inversion.vel_obs=tanh(30.*(u.^2+v.^2-0.25)) ... 92 +tanh(30.*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u-0.75).^2+(v+0.75).^2-0.25)) ... 93 +tanh(30.*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30.*((u+0.75).^2+(v+0.75).^2-0.25)); 94 94 plotmodel(md,'data',md.inversion.vel_obs,'data',md.inversion.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); 95 if printingflag, 95 if printingflag, 96 96 set(gcf,'Color','w') 97 97 printmodel('mesh2_bamgiso','png','margin','on','marginsize',25,'frame','off','resolution',1,'hardcopy','off');
Note:
See TracChangeset
for help on using the changeset viewer.
