Changeset 14286


Ignore:
Timestamp:
01/25/13 08:01:41 (12 years ago)
Author:
Eric.Larour
Message:

CHG: plot BC without dirichlet on

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/plot/plot_BC.m

    r13730 r14286  
    77
    88%plot dirichlets
    9 h1=plot3(...
    10         md.mesh.x(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
    11         md.mesh.y(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
    12         md.mesh.z(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
    13         'ro','MarkerSize',14,'MarkerFaceColor','r');
    14 h2=plot3(...
    15         md.mesh.x(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
    16         md.mesh.y(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
    17         md.mesh.z(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
    18         'bo','MarkerSize',10,'MarkerFaceColor','b');
    19 h3=plot3(...
    20         md.mesh.x(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
    21         md.mesh.y(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
    22         md.mesh.z(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
    23         'yo','MarkerSize',6 ,'MarkerFaceColor','y');
     9dirichleton=getfieldvalue(options,'dirichlet','on');
     10if strcmpi(dirichleton,'on'),
     11        h1=plot3(...
     12                md.mesh.x(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
     13                md.mesh.y(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
     14                md.mesh.z(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),...
     15                'ro','MarkerSize',14,'MarkerFaceColor','r');
     16        h2=plot3(...
     17                md.mesh.x(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
     18                md.mesh.y(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
     19                md.mesh.z(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),...
     20                'bo','MarkerSize',10,'MarkerFaceColor','b');
     21        h3=plot3(...
     22                md.mesh.x(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
     23                md.mesh.y(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
     24                md.mesh.z(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),...
     25                'yo','MarkerSize',6 ,'MarkerFaceColor','y');
     26end
    2427
    2528%update legend
    2629[legend_h,object_h,plot_h,text_strings]=legend();
    2730legend('off');
    28 text_strings{end+1}='vx Dirichlet';
    29 text_strings{end+1}='vy Dirichlet';
    30 if h3, text_strings{end+1}='vz Dirichlet'; end
    31 plot_h(end+1)=h1;
    32 plot_h(end+1)=h2;
    33 if h3, plot_h(end+1)=h3; end
     31if strcmpi(dirichleton,'on'),
     32        text_strings{end+1}='vx Dirichlet';
     33        text_strings{end+1}='vy Dirichlet';
     34        if h3, text_strings{end+1}='vz Dirichlet'; end
     35        plot_h(end+1)=h1;
     36        plot_h(end+1)=h2;
     37        if h3, plot_h(end+1)=h3; end
     38end
    3439legend(plot_h,text_strings,'location','NorthEast')
    3540
Note: See TracChangeset for help on using the changeset viewer.