Changeset 21705


Ignore:
Timestamp:
05/04/17 14:24:09 (8 years ago)
Author:
Mathieu Morlighem
Message:

CHG: allow for layer plot BC

Location:
issm/trunk-jpl/src/m/plot
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified issm/trunk-jpl/src/m/plot/plot_BC.m

    r19203 r21705  
    77
    88[x y z elements is2d isplanet]=processmesh(md,[],options);
     9spcvx=processdata(md,md.stressbalance.spcvx,options);
     10spcvy=processdata(md,md.stressbalance.spcvy,options);
     11spcvz=processdata(md,md.stressbalance.spcvz,options);
     12nbv = numel(x);
    913
    1014%plot dirichlets
     
    1216if strcmpi(dirichleton,'on'),
    1317        h1=plot3(...
    14                 x(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
    15                 y(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
    16                 z(find(~isnan(md.stressbalance.spcvx(1:md.mesh.numberofvertices,1)))),...
     18                x(find(~isnan(spcvx(1:nbv,1)))),...
     19                y(find(~isnan(spcvx(1:nbv,1)))),...
     20                z(find(~isnan(spcvx(1:nbv,1)))),...
    1721                'ro','MarkerSize',14,'MarkerFaceColor','r');
    1822        h2=plot3(...
    19                 x(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
    20                 y(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
    21                 z(find(~isnan(md.stressbalance.spcvy(1:md.mesh.numberofvertices,1)))),...
     23                x(find(~isnan(spcvy(1:nbv,1)))),...
     24                y(find(~isnan(spcvy(1:nbv,1)))),...
     25                z(find(~isnan(spcvy(1:nbv,1)))),...
    2226                'bo','MarkerSize',10,'MarkerFaceColor','b');
    2327        h3=plot3(...
    24                 x(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
    25                 y(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
    26                 z(find(~isnan(md.stressbalance.spcvz(1:md.mesh.numberofvertices,1)))),...
     28                x(find(~isnan(spcvz(1:nbv,1)))),...
     29                y(find(~isnan(spcvz(1:nbv,1)))),...
     30                z(find(~isnan(spcvz(1:nbv,1)))),...
    2731                'yo','MarkerSize',6 ,'MarkerFaceColor','y');
    2832end
  • TabularUnified issm/trunk-jpl/src/m/plot/plot_icefront.m

    r19203 r21705  
    1919elementzeroice=sum(zeroice(md.mesh.elements),2);
    2020
    21 if dimension(md.mesh)==2,
     21if is2d,
    2222        icefront=(elementice & elementnoice) & ~(elementice==2 & elementzeroice);
    2323
  • TabularUnified issm/trunk-jpl/src/m/plot/plot_manager.m

    r17718 r21705  
    2121                        plot_boundaries(md,options,subplotwidth,i);
    2222                        return;
    23                 case 'BC',
     23                case {'BC','bc'},
    2424                        plot_BC(md,options,subplotwidth,i,data);
    2525                        return;
Note: See TracChangeset for help on using the changeset viewer.