plot_drivingstress

PURPOSE ^

PLOT_DRIVINGSTRESS - plot driving stress

SYNOPSIS ^

function plot_drivingstress(md,options_structure,width,i);

DESCRIPTION ^

PLOT_DRIVINGSTRESS - plot driving stress

   Usage:
      plot_drivingstress(md,options_structure,width,i);

   See also: PLOTMODEL

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function plot_drivingstress(md,options_structure,width,i);
0002 %PLOT_DRIVINGSTRESS - plot driving stress
0003 %
0004 %   Usage:
0005 %      plot_drivingstress(md,options_structure,width,i);
0006 %
0007 %   See also: PLOTMODEL
0008 
0009 %get driving stress
0010 [sx sy s]=drivingstress(md);
0011 
0012 %process data and model
0013 [x y z elements is2d]=processmesh(md,options_structure);
0014 [dstress isongrid]=processdata(md,s,options_structure);
0015 dstress=dstress/1000;
0016 
0017 %plot mesh quivervel
0018 subplot(width,width,i); 
0019 
0020 %edgecolor?
0021 if ~isnan(options_structure.edgecolor),
0022     edgecolor=options_structure.edgecolor;
0023 else
0024     edgecolor='none';
0025 end
0026 
0027 %element data
0028 if ~isongrid
0029     if is2d
0030         A=elements(:,1); B=elements(:,2); C=elements(:,3); 
0031         patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0032     else
0033         A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
0034         patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0035         patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0036         patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0037         patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0038         patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0039     end
0040 %grid data
0041 elseif isongrid
0042     if is2d
0043         A=elements(:,1); B=elements(:,2); C=elements(:,3); 
0044         patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0045     else
0046         if options_structure.layer>=1,
0047             A=elements(:,1); B=elements(:,2); C=elements(:,3); 
0048             patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0049         else
0050             A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
0051             patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0052             patch( 'Faces', [D E F], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0053             patch( 'Faces', [A B E D], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0054             patch( 'Faces', [B E F C ], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0055             patch( 'Faces', [C A D F ], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0056         end
0057     end
0058 end
0059 
0060 %apply options
0061 if isnan(options_structure.title)
0062     options_structure.title='Driving stress [kPa]';
0063 end 
0064 applyoptions(md,dstress,options_structure);

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003