Changeset 15564 for issm/trunk-jpl/src/m/plot
- Timestamp:
- 07/23/13 14:23:07 (12 years ago)
- Location:
- issm/trunk-jpl/src/m/plot
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/plot/plot_elementstype.m
r13730 r15564 18 18 19 19 if is2d 20 % Hutterelements20 %SIA elements 21 21 posH=find(data==1); 22 22 A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); … … 26 26 A=elements(posM,1); B=elements(posM,2); C=elements(posM,3); 27 27 p2=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',2,'FaceColor','flat','EdgeColor',edgecolor); 28 % Pattynelement28 %HO element 29 29 posP=find(data==3); 30 30 A=elements(posP,1); B=elements(posP,2); C=elements(posP,3); 31 31 p3=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',3,'FaceColor','flat','EdgeColor',edgecolor); 32 %MacAyeal Pattynelement32 %MacAyealHO element 33 33 posMP=find(data==5); 34 34 A=elements(posMP,1); B=elements(posMP,2); C=elements(posMP,3); 35 35 p5=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',5,'FaceColor','flat','EdgeColor',edgecolor); 36 % Stokeselements36 %FS elements 37 37 posS=find(data==4); 38 38 A=elements(posS,1); B=elements(posS,2); C=elements(posS,3); 39 39 p6=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',4,'FaceColor','flat','EdgeColor',edgecolor); 40 %MacAyeal Stokeselements40 %MacAyealFS elements 41 41 posMS=find(data==6); 42 42 A=elements(posMS,1); B=elements(posMS,2); C=elements(posMS,3); 43 43 p7=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',6,'FaceColor','flat','EdgeColor',edgecolor); 44 % PattynStokeselements44 %HOFS elements 45 45 posPS=find(data==7); 46 46 A=elements(posPS,1); B=elements(posPS,2); C=elements(posPS,3); … … 52 52 53 53 legend([p1 p2 p3 p5 p6 p7 p8 p9],... 54 ' Hutter''s elements','MacAyeal''s elements','Pattyn''s elements',...55 'MacAyeal Pattyn''s elements','Stokes''s elements','MacAyealStokes''s elements','PattynStokes''s elements','None element');54 'SIA''s elements','MacAyeal''s elements','HO''s elements',... 55 'MacAyealHO''s elements','FS''s elements','MacAyealFS''s elements','HOFS''s elements','None element'); 56 56 57 57 else 58 % Hutterelements58 %SIA elements 59 59 posH=find(data==1); 60 60 A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); D=elements(posH,4); E=elements(posH,5); F=elements(posH,6); … … 72 72 patch( 'Faces', [B E F C],'Vertices', [x y z],'CData', 2,'FaceColor','flat','EdgeColor',edgecolor); 73 73 patch( 'Faces', [C A D F],'Vertices', [x y z],'CData', 2,'FaceColor','flat','EdgeColor',edgecolor); 74 % Pattynelements74 %HO elements 75 75 posP=find(data==3); 76 76 A=elements(posP,1); B=elements(posP,2); C=elements(posP,3); D=elements(posP,4); E=elements(posP,5); F=elements(posP,6); … … 80 80 patch( 'Faces', [B E F C],'Vertices', [x y z],'CData', 3,'FaceColor','flat','EdgeColor',edgecolor); 81 81 patch( 'Faces', [C A D F],'Vertices', [x y z],'CData', 3,'FaceColor','flat','EdgeColor',edgecolor); 82 % Stokeselements82 %FS elements 83 83 posS=find(data==4); 84 84 A=elements(posS,1); B=elements(posS,2); C=elements(posS,3); D=elements(posS,4); E=elements(posS,5); F=elements(posS,6); … … 88 88 patch( 'Faces', [B E F C],'Vertices', [x y z],'CData', 4,'FaceColor','flat','EdgeColor',edgecolor); 89 89 patch( 'Faces', [C A D F],'Vertices', [x y z],'CData', 4,'FaceColor','flat','EdgeColor',edgecolor); 90 %MacAyeal Pattynelements90 %MacAyealHO elements 91 91 posP=find(data==5); 92 92 A=elements(posP,1); B=elements(posP,2); C=elements(posP,3); D=elements(posP,4); E=elements(posP,5); F=elements(posP,6); … … 96 96 patch( 'Faces', [B E F C],'Vertices', [x y z],'CData', 5,'FaceColor','flat','EdgeColor',edgecolor); 97 97 patch( 'Faces', [C A D F],'Vertices', [x y z],'CData', 5,'FaceColor','flat','EdgeColor',edgecolor); 98 % PattynStokeselements98 %HOFS elements 99 99 PosPS=find(data==7); 100 100 A=elements(PosPS,1); B=elements(PosPS,2); C=elements(PosPS,3); D=elements(PosPS,4); E=elements(PosPS,5); F=elements(PosPS,6); … … 104 104 patch( 'Faces', [B E F C],'Vertices', [x y z],'CData', 7,'FaceColor','flat','EdgeColor',edgecolor); 105 105 patch( 'Faces', [C A D F],'Vertices', [x y z],'CData', 7,'FaceColor','flat','EdgeColor',edgecolor); 106 %MacAyeal Stokeselements106 %MacAyealFS elements 107 107 PosMS=find(data==6); 108 108 A=elements(PosMS,1); B=elements(PosMS,2); C=elements(PosMS,3); D=elements(PosMS,4); E=elements(PosMS,5); F=elements(PosMS,6); … … 122 122 123 123 legend([p1 p2 p3 p4 p5 p6 p7 p8],... 124 ' Hutter''s elements','MacAyeal''s elements','Pattyn''s elements','Stokes''s elements',...125 'MacAyeal Pattyn''s elements','PattynStokes''s elements','MacAyealStokes''s elements','None elements');124 'SIA''s elements','MacAyeal''s elements','HO''s elements','FS''s elements',... 125 'MacAyealHO''s elements','HOFS''s elements','MacAyealFS''s elements','None elements'); 126 126 end 127 127 -
issm/trunk-jpl/src/m/plot/plot_penalties.m
r13730 r15564 40 40 P2=plot3(x(md.penalties(i,:)),y(md.penalties(i,:)),z(md.penalties(i,:)),'bo-','LineWidth',2,'MarkerSize',8,'MarkerFaceColor','b'); 41 41 end 42 legend([P1 P2],'MacAyeal''s penalized nodes',' Pattyn''s penalized nodes');42 legend([P1 P2],'MacAyeal''s penalized nodes','HO''s penalized nodes'); 43 43 end 44 44
Note:
See TracChangeset
for help on using the changeset viewer.