Changeset 18162


Ignore:
Timestamp:
06/18/14 13:46:23 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: tetras in plot elements

File:
1 edited

Legend:

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

    r17482 r18162  
    1616
    1717                pos=find(~isnan(data)); %needed for element on water
    18                 if is2d,
    19                         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', data(pos),'FaceColor','flat','EdgeColor',edgecolor);
    20                 else
     18                if size(elements,2)==6, %prisms
    2119                        A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); D=elements(pos,4); E=elements(pos,5); F=elements(pos,6);
    2220                        patch( 'Faces', [A B C],  'Vertices', [x y z],'CData', data(pos),'FaceColor','flat','EdgeColor',edgecolor);
     
    2523                        patch( 'Faces', [B E F C],'Vertices', [x y z],'CData', data(pos),'FaceColor','flat','EdgeColor',edgecolor);
    2624                        patch( 'Faces', [C A D F],'Vertices', [x y z],'CData', data(pos),'FaceColor','flat','EdgeColor',edgecolor);
     25                elseif size(elements,2)==4, %tetras
     26                        A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4);
     27                        patch( 'Faces',[A B C],'Vertices', [x y z],'CData',data(pos),'FaceColor','flat','EdgeColor',edgecolor);
     28                        patch( 'Faces',[A B D],'Vertices', [x y z],'CData',data(pos),'FaceColor','flat','EdgeColor',edgecolor);
     29                        patch( 'Faces',[B C D],'Vertices', [x y z],'CData',data(pos),'FaceColor','flat','EdgeColor',edgecolor);
     30                        patch( 'Faces',[C A D],'Vertices', [x y z],'CData',data(pos),'FaceColor','flat','EdgeColor',edgecolor);
     31                else
     32                        A=elements(pos,1); B=elements(pos,2); C=elements(pos,3);
     33                        patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', data(pos),'FaceColor','flat','EdgeColor',edgecolor);
     34                end
     35
     36                if is2d,
    2737                end
    2838
Note: See TracChangeset for help on using the changeset viewer.