Changeset 5418


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

finished ismipa pattyn

File:
1 edited

Legend:

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

    r5405 r5418  
    11%This test is a test from the ISMP-HOM Intercomparison project
    22%Pattyn and Payne 2006
     3printingflag=true;
    34
    45L_list={5000,10000,20000,40000,80000,160000};
    56results={};
     7minvx=[];
     8maxvx=[];
    69
    710for i=1:length(L_list),
     
    3942        vz=PatchToVec(md.results.DiagnosticSolution.Vz);
    4043        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));
    4146
    4247        %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
    5266
    5367        if(L==5000),
     
    7084                        'resolution',[10 10],'ylim',[0 120],'xlim',[0 160000],'title','','xlabel','')
    7185        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
    7491end
    7592
     93%Now plot the min and max values of vx for each size of the square
     94plot([5 10 20 40 80 160],minvx);ylim([0 18])
     95if 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 ']);
     99end
     100plot([5 10 20 40 80 160],maxvx);ylim([0 120])
     101if 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 ']);
     105end
    76106%Fields and tolerances to track changes
    77107field_names     ={ ...
Note: See TracChangeset for help on using the changeset viewer.