Changeset 5405


Ignore:
Timestamp:
08/19/10 10:01:16 (15 years ago)
Author:
seroussi
Message:

added plots in test 1101

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/test/NightlyRun/test1101.m

    r5098 r5405  
    4040        results{i}=md.results.DiagnosticSolution;
    4141
    42         plotmodel(md,'data',vx,'data',vy,'data',vz,'layer#all',md.numlayers)
     42        %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');
     52
     53        if(L==5000),
     54                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP5000.exp','layer',md.numlayers,...
     55                        'resolution',[10 10],'ylim',[10 18],'xlim',[0 5000],'title','','xlabel','')
     56        elseif(L==10000),
     57                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP10000.exp','layer',md.numlayers,...
     58                        'resolution',[10 10],'ylim',[10 30],'xlim',[0 10000],'title','','xlabel','')
     59        elseif(L==20000),
     60                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP20000.exp','layer',md.numlayers,...
     61                        'resolution',[10 10],'ylim',[0 50],'xlim',[0 20000],'title','','xlabel','')
     62        elseif(L==40000),
     63                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP40000.exp','layer',md.numlayers,...
     64                        'resolution',[10 10],'ylim',[0 80],'xlim',[0 40000],'title','','xlabel','')
     65        elseif(L==80000),
     66                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP80000.exp','layer',md.numlayers,...
     67                        'resolution',[10 10],'ylim',[0 100],'xlim',[0 80000],'title','','xlabel','')
     68        elseif(L==160000),
     69                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP160000.exp','layer',md.numlayers,...
     70                        'resolution',[10 10],'ylim',[0 120],'xlim',[0 160000],'title','','xlabel','')
     71        end
     72        set(gcf,'Color','w')
     73        printmodel(['ismipapattynvxsec' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',2,'hardcopy','off');
    4374end
    4475
Note: See TracChangeset for help on using the changeset viewer.