Changeset 5418
- Timestamp:
- 08/19/10 11:35:16 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/NightlyRun/test1101.m
r5405 r5418 1 1 %This test is a test from the ISMP-HOM Intercomparison project 2 2 %Pattyn and Payne 2006 3 printingflag=true; 3 4 4 5 L_list={5000,10000,20000,40000,80000,160000}; 5 6 results={}; 7 minvx=[]; 8 maxvx=[]; 6 9 7 10 for i=1:length(L_list), … … 39 42 vz=PatchToVec(md.results.DiagnosticSolution.Vz); 40 43 results{i}=md.results.DiagnosticSolution; 44 minvx(i)=min(vx(end-md.numberofgrids2d+1:end)); 45 maxvx(i)=max(vx(end-md.numberofgrids2d+1:end)); 41 46 42 47 %Now plot vx, vy, vz and vx on a cross section 43 plotmodel(md,'data',vx,'layer#all',md.numlayers,'xlim',[0 L],'ylim',[0 L]) 44 set(gcf,'Color','w') 45 printmodel(['ismipapattynvx' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off'); 46 plotmodel(md,'data',vy,'layer#all',md.numlayers,'xlim',[0 L],'ylim',[0 L]) 47 set(gcf,'Color','w') 48 printmodel(['ismipapattynvy' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off'); 49 plotmodel(md,'data',vz,'layer#all',md.numlayers,'xlim',[0 L],'ylim',[0 L]) 50 set(gcf,'Color','w') 51 printmodel(['ismipapattynvz' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off'); 48 plotmodel(md,'data',vx,'layer#all',md.numlayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km') 49 if printingflag, 50 set(gcf,'Color','w') 51 printmodel(['ismipapattynvx' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off'); 52 system(['mv ismipapattynvx' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']); 53 end 54 plotmodel(md,'data',vy,'layer#all',md.numlayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km') 55 if printingflag, 56 set(gcf,'Color','w') 57 printmodel(['ismipapattynvy' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off'); 58 system(['mv ismipapattynvy' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']); 59 end 60 plotmodel(md,'data',vz,'layer#all',md.numlayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km') 61 if printingflag, 62 set(gcf,'Color','w') 63 printmodel(['ismipapattynvz' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off'); 64 system(['mv ismipapattynvz' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']); 65 end 52 66 53 67 if(L==5000), … … 70 84 'resolution',[10 10],'ylim',[0 120],'xlim',[0 160000],'title','','xlabel','') 71 85 end 72 set(gcf,'Color','w') 73 printmodel(['ismipapattynvxsec' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off'); 86 if printingflag, 87 set(gcf,'Color','w') 88 printmodel(['ismipapattynvxsec' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off'); 89 system(['mv ismipapattynvxsec' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']); 90 end 74 91 end 75 92 93 %Now plot the min and max values of vx for each size of the square 94 plot([5 10 20 40 80 160],minvx);ylim([0 18]) 95 if printingflag, 96 set(gcf,'Color','w') 97 printmodel('ismipapattynminvx','png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off'); 98 system(['mv ismipapattynminvx.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']); 99 end 100 plot([5 10 20 40 80 160],maxvx);ylim([0 120]) 101 if printingflag, 102 set(gcf,'Color','w') 103 printmodel('ismipapattynmaxvx','png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off'); 104 system(['mv ismipapattynmaxvx.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']); 105 end 76 106 %Fields and tolerances to track changes 77 107 field_names ={ ...
Note:
See TracChangeset
for help on using the changeset viewer.