Changeset 9850
- Timestamp:
- 09/20/11 08:10:37 (13 years ago)
- Location:
- issm/trunk/test/NightlyRun
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test1208.m
r9766 r9850 22 22 md.settings.output_frequency=20; 23 23 md.timestepping.final_time=50000; 24 md.prognostic.stabilization= 2;24 md.prognostic.stabilization=1; 25 25 md.thermal.stabilization=2; 26 26 -
issm/trunk/test/NightlyRun/test1302.m
r9766 r9850 34 34 relative=abs((comp_temp-md.initialization.temperature)./md.initialization.temperature)*100; 35 35 relative(find(comp_temp==md.initialization.temperature))=0; 36 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data', 'temperature','view',3,...36 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data',md.initialization.temperature,'view',3,... 37 37 'title','Analytical temperature [K]','view',3,'data',comp_temp-md.initialization.temperature,... 38 38 'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,... -
issm/trunk/test/NightlyRun/test1303.m
r9734 r9850 29 29 relative=abs((comp_temp-md.initialization.temperature)./md.initialization.temperature)*100; 30 30 relative(find(comp_temp==md.initialization.temperature))=0; 31 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data', 'temperature','view',3,...31 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data',md.initialization.temperature,'view',3,... 32 32 'title','Analytical temperature [K]','view',3,'data',comp_temp-md.initialization.temperature,... 33 33 'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,... -
issm/trunk/test/NightlyRun/test1304.m
r9766 r9850 30 30 relative=abs((comp_temp-md.initialization.temperature)./md.initialization.temperature)*100; 31 31 relative(find(comp_temp==md.initialization.temperature))=0; 32 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data', 'temperature','view',3,...32 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data',md.initialization.temperature,'view',3,... 33 33 'title','Analytical temperature','view',3,'data',comp_temp-md.initialization.temperature,... 34 34 'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,... -
issm/trunk/test/NightlyRun/test1401.m
r9734 r9850 12 12 md=squaremesh(md,L,L,nx,ny); 13 13 md.inversion.vel_obs=exp(-(sqrt((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)-0.75).^2*10^6)+((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)/2; 14 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);14 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); 15 15 if printingflag, 16 16 set(gcf,'Color','w') … … 21 21 md=YamsCall(md,md.inversion.vel_obs,0.001,0.3,1.3,10^-4); 22 22 md.inversion.vel_obs=exp(-(sqrt((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)-0.75).^2*10^6)+((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)/2; 23 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);23 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); 24 24 if printingflag, 25 25 set(gcf,'Color','w') … … 30 30 md=YamsCall(md,md.inversion.vel_obs,0.001,0.3,2.5,0.008); 31 31 md.inversion.vel_obs=exp(-(sqrt((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)-0.75).^2*10^6)+((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)/2; 32 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);32 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); 33 33 if printingflag, 34 34 set(gcf,'Color','w') … … 42 42 md=squaremesh(md,L,L,nx,ny); 43 43 md.inversion.vel_obs=exp(-(sqrt((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)-0.75).^2*10^6)+((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)/2; 44 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);44 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); 45 45 if printingflag, 46 46 set(gcf,'Color','w') … … 49 49 end 50 50 51 md. bamg=NaN;51 md.private.bamg=NaN; 52 52 md=bamg(md,'field',md.inversion.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',1.3,'err',10^-4); 53 53 md.inversion.vel_obs=exp(-(sqrt((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)-0.75).^2*10^6)+((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)/2; 54 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);54 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); 55 55 if printingflag, 56 56 set(gcf,'Color','w') … … 59 59 end 60 60 61 md. bamg=NaN;61 md.private.bamg=NaN; 62 62 md=bamg(md,'field',md.inversion.vel_obs,'hmin',0.001,'hmax',0.3,'gradation',2.5,'err',0.008); 63 63 md.inversion.vel_obs=exp(-(sqrt((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)-0.75).^2*10^6)+((md.mesh.x+0.1).^2+(md.mesh.y+0.1).^2)/2; 64 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);64 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); 65 65 if printingflag, 66 66 set(gcf,'Color','w') -
issm/trunk/test/NightlyRun/test1402.m
r9734 r9850 14 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 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 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);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 17 if printingflag, 18 18 set(gcf,'Color','w') … … 26 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 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 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);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 29 if printingflag, 30 30 set(gcf,'Color','w') … … 38 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 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 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);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 41 if printingflag, 42 42 set(gcf,'Color','w') … … 53 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 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 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);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 56 if printingflag, 57 57 set(gcf,'Color','w') … … 60 60 end 61 61 62 md. bamg=NaN;62 md.private.bamg=NaN; 63 63 md=bamg(md,'field',md.inversion.vel_obs,'hmin',0.005,'hmax',0.3,'gradation',2.3,'err',10^-2); 64 64 u=4*md.mesh.x-2; v=4*md.mesh.y-2; … … 66 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 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 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);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 69 if printingflag, 70 70 set(gcf,'Color','w') … … 73 73 end 74 74 75 md. bamg=NaN;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 77 u=4*md.mesh.x-2; v=4*md.mesh.y-2; … … 79 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 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 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);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 82 if printingflag, 83 83 set(gcf,'Color','w') … … 86 86 end 87 87 88 md. bamg=NaN;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 90 u=4*md.mesh.x-2; v=4*md.mesh.y-2; … … 92 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 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 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);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 95 if printingflag, 96 96 set(gcf,'Color','w')
Note:
See TracChangeset
for help on using the changeset viewer.