MS_plot
MS_PLOT - Plot phasevels/anisotropy on pole figures.
Given an elasticity matrix and density, produce pole figures showing the P- and S-wave anisotrpy.
MS_plot(C, rh, ...)
Usage: MS_plot(C, rh) Produce three pole figures showing P-wave velocity, S-wave anisotropy and fast S-wave polarisation direction
MS_plot(..., 'fontsize', f) Set minimum fontsize in plots to f.
MS_plot(..., 'wtitle', S) Set the (window) title to be string S. If a minus sign is prepended to this title, the figure number is also suppressed (the minus sign is not shown).
MS_plot(..., 'cmap', CM) Redefine the colormap. CM can either be a (nx3) matrix containing a colormap, or a string describing a function to generate such a matrix (such as the built-in MATLAB colormap functions). E.g.: MS_plot(C,rh,'cmap','cool(64)') - uses the MATLAB function cool to generate a cyan-to-purple colourmap. The default is 'jet', reversed so blue is fast/high, as is conventional for seismic velocity colorscales.
MS_plot(..., 'reverse') Unreverse the sense of the colour map. Default (no argument) is for blue to be fast/high as is conventional for seismic velocity colorscales. Setting this option results in red being fast/high.
MS_plot(..., 'pcontours', pcvect) MS_plot(..., 'scontours', scvect) Set the contour levels for the P-wave velocity and S-wave anisotropy plots. The values of pcvect and scvect can be scalars indicating the number of contour lines to use or vectors, where each value represents a different contour line. Thus, the minimum, maximum values and spacing of contours can easily be set using the syntax [minval:spacing:maxval] for pcvect and scvect.
MS_plot(..., 'limitsonpol') Include markers for the maximum and minimum values of S-wave anisotropy on the plot of fast S-wave polarisation direction. Not shown by default.
MS_plot(..., 'polsize' s1, s2, w1, w2) Set the size of the markers used to indicate the fast S-wave polarisation direction. There are two markers used a larger white "background" marker with a length set by s1 (defaults to 0.18) and width set by w1 (defaults to 3.0) and an inner black marker with width set by w2 (defaults to 2.0) and length set by s2 (defaults to 0.18). A useful "neet" alternitive to the default is s1 = 0.18, s2 = 0.16, w1 = 2.0 and w2 = 1.0.
MS_plot(..., 'plotmap', pmap) Change the order, or content, of the subplots. This allows a single pole figure (e.g. the P-wave velocity) to be drawn, or the pole figures to be drawn in a different order. The argument pmap must be a cell array of strings. Each string represents a different pole figure and its location in the cell array represents the location of the pole figure. Thus the default {'VP', 'AVS', 'AVSPOL'} draws three pole figures in a row. If the argument were {'VP'; 'AVS'} two pole figures would be drawn, one above the other. Valid strings are 'VP", 'AVS' and 'AVSPOL'. These are not case sensitive.
MS_plot(..., 'pdata', azimuth, inclination, vp) MS_plot(..., 'sdata', azimuth, inclination, polarisation, avs) Add data points to the P-wave velocity ('VP') or S-wave polarisation plot ('AVSPOL'), respectively. This can be used, for example, to compare an elastic model with shear-wave splitting measurements. In each case data points relate to a ray propagating in a direction described by the azimuth and inclination. The point is coloured according to either the P-wave velocity, vp, or S-wave anisotropy, avs, in km/s or. In the case of S-wave anisotropy, the fast polarisation direction is also shown. Each argument (azimuth, inclination, vp, polarisation and avs) is an array and each must be the same length.
MS_plot(..., 'quiet') Don't write isotropic velocities to the terminal.
See also: MS_SPHERE, MS_PHASEVELS