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(..., 'avscontours', ascvect)
MS_plot(..., 'scontours', scvect)
     Set the contour levels for the P-wave velocity, S-wave
     anisotropy and (optional) S-wave velocity plots. The values of
     pcvect, ascvect 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 ascvect or 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