- added in export_fig codes from file exchange

  for easy figure export
- added exportfig directory to optional installer
- major update to sph_tri_* codes (output format
  now works better with delaunay codes and higher
  order polygon codes)
- dropped sph_tri_latlon & sph_tri_vertices as
  these are trivial with new output format
- fixed decay constant error bounds in slowdecay*
- minor aesthetic changes to plot_cmb_measurements
  & cmb_pdf_mtx output
- name2rgb now includes over 700 color names!
- wedge introduced (like polar but better)
- plotfkmap & plotfkarf now use wedge
	- better zoom support & simpler to use
- mantledv now outputs wavetype of model too
- mantlemap cleanup
- mantleprofile introduced (uses new wedge function)
- minor doc updates elsewhere

Author: g2e

Diff for: Contents.m

@@ -1,5 +1,5 @@
 % SEIZMO - Passive Seismology Toolbox
-% Version 0.6.189 Haleakala 17-Feb-2011
+% Version 0.6.190 Haleakala 25-Feb-2011
 %
 % Help Info
 %  <a href="matlab:help install_seizmo">INSTALLATION</a>

Diff for: alphabetical_list.m

@@ -285,6 +285,7 @@
 %<a href="matlab:help mancor">mancor</a>                         - Returns mantle travel time corrections for a set of raypaths
 %<a href="matlab:help mantledv">mantledv</a>                       - Returns the seismic velocity deviation for a mantle model
 %<a href="matlab:help mantlemap">mantlemap</a>                      - 3D mantle model map (aka depth slice)
+%<a href="matlab:help mantleprofile">mantleprofile</a>                  - 3D mantle model profile (aka radial slice)
 %<a href="matlab:help map_cmb_profiles">map_cmb_profiles</a>               - Map CMB profile measurements
 %<a href="matlab:help mapeventgrid">mapeventgrid</a>                   - Draws range/azimuth grid (relative to an event) on a map
 %<a href="matlab:help mapfeature">mapfeature</a>                     - Map geographic feature
@@ -457,9 +458,7 @@
 %<a href="matlab:help sph_poly_in">sph_poly_in</a>                    - Returns indexing for which polygon some points are in
 %<a href="matlab:help sph_tri_auto">sph_tri_auto</a>                   - Automatically create spherical triangle polyhedra
 %<a href="matlab:help sph_tri_init">sph_tri_init</a>                   - Returns vertices for regular spherical triangle polyhedra
-%<a href="matlab:help sph_tri_latlon">sph_tri_latlon</a>                 - Convert spherical triangle xyz arrays to lat/lon
 %<a href="matlab:help sph_tri_split">sph_tri_split</a>                  - Splits up faces of spherical triangle polyhedra
-%<a href="matlab:help sph_tri_vertices">sph_tri_vertices</a>               - Unique vertices of a spherical triangle polyhedra
 %<a href="matlab:help sphericalfwd">sphericalfwd</a>                   - Finds a point on a sphere relative to another point
 %<a href="matlab:help sphericalinv">sphericalinv</a>                   - Return distance and azimuth between 2 locations on sphere
 %<a href="matlab:help split">split</a>                          - Lightblue-Blue-Black-Red-Lightred colormap
@@ -542,6 +541,7 @@
 %<a href="matlab:help vertcmp">vertcmp</a>                        - Finds vertical component SEIZMO records
 %<a href="matlab:help vincentyfwd">vincentyfwd</a>                    - Find destination point on an ellipsoid relative to a point
 %<a href="matlab:help vincentyinv">vincentyinv</a>                    - Find distance and azimuth between 2 locations on ellipsoid
+%<a href="matlab:help wedge">wedge</a>                          - Polar coordinate wedge plot
 %<a href="matlab:help whiten">whiten</a>                         - Spectral whitening/normalization of SEIZMO data records
 %<a href="matlab:help wlinem">wlinem</a>                         - Linear fit to traveltime data with weighted least squares
 %<a href="matlab:help write_1dmodel_nd">write_1dmodel_nd</a>               - Writes a 1D model struct in .nd format

Diff for: cmb/cmb_pdf_mtx.m

@@ -59,9 +59,10 @@
 %        Feb.  1, 2011 - initial version
 %        Feb. 10, 2011 - reduced plot code redundancy, use period not freq,
 %                        improved probability estimation, works with 1 pf
+%        Feb. 17, 2011 - aesthetic touches
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Feb. 10, 2011 at 13:35 GMT
+%     Last Updated Feb. 17, 2011 at 13:35 GMT
 
 % todo:
 
@@ -156,10 +157,12 @@
     imagesc(p,x,mtx','parent',ax);
     colormap(ax,fire);
     set(ax,'ydir','normal','color','k','xcolor','w','ycolor','w');
+    set(ax,'xminortick','on','yminortick','on');
+    set(ax,'xgrid','on','ygrid','on','linewidth',1);
     xlabel(ax,'Period (s)');
     cb=colorbar('peer',ax);
     ylabel(cb,'% Probability');
-    grid(ax,'on');
+    set(cb,'linewidth',1);
     switch field
         case 'slow'
             ylabel(ax,'Ray Parameter (s/^o)');

Diff for: cmb/plot_cmb_measurements.m

@@ -33,9 +33,10 @@
 %        Feb.  1, 2011 - initial version
 %        Feb. 10, 2011 - reusing axes works, better labeling, doc update
 %        Feb. 12, 2011 - altered line style code
+%        Feb. 17, 2011 - aesthetic touches
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Feb. 12, 2011 at 13:35 GMT
+%     Last Updated Feb. 17, 2011 at 13:35 GMT
 
 % todo:
 
@@ -78,7 +79,7 @@
 ax=axescheck(varargin{:});
 if(isempty(ax))
     % new plot
-    fh=figure();
+    fh=figure('color','w');
     ax=axes('parent',fh);
     oldpts=0;
 else
@@ -102,7 +103,8 @@
 
     % plot measurements
     [h(i),hyerr(i),hxerr(i)]=ploterr(ax,p0(a),x0(a),{lp(a) up(a)},e0(a),...
-        [l m],'color',c,'markerfacecolor',mc,'markersize',10,varargin{:});
+        [l m],'color',c,'markerfacecolor',mc,'markersize',6,...
+        'linewidth',1,varargin{:});
 end
 
 % change names if corrected
@@ -151,7 +153,12 @@
 % legend (regrab points objects to get those from previous plotting)
 drawnow; % legends can have issues if plot not drawn yet
 lh=legend(findobj(ax,'tag','points'));
-set(lh,'interpreter','none');
+set(lh,'interpreter','none','linewidth',1);
+
+% polish
+box(ax,'on');
+set(ax,'xminortick','on','yminortick','on');
+set(ax,'xgrid','on','ygrid','on','linewidth',1);
 
 % reset hold state
 if(~held); hold(ax,'off'); end

Diff for: cmb/slowdecaypairs.m

@@ -85,9 +85,10 @@
 %        Jan. 31, 2011 - allow no output, odir input, better checks
 %        Feb.  5, 2011 - fix bug when no output specified
 %        Feb. 12, 2011 - include snr-based arrival time error
+%        Feb. 17, 2011 - fixed decay constant error
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Feb. 12, 2011 at 13:35 GMT
+%     Last Updated Feb. 17, 2011 at 13:35 GMT
 
 % todo:
 
@@ -263,16 +264,16 @@
         tmp(b).cslowerr=sqrt(covm(2,2));
         [m,covm]=wlinem(delaz([idx1(b) idx2(b)],1),...
             log(rampl([idx1(b) idx2(b)])),1,...
-            diag(log(rampl([idx1(b) idx2(b)])...
-            +ramplerr([idx1(b) idx2(b)]).^2)....
-            -log(rampl([idx1(b) idx2(b)]))));
+            diag((log(rampl([idx1(b) idx2(b)])...
+            +ramplerr([idx1(b) idx2(b)]))....
+            -log(rampl([idx1(b) idx2(b)]))).^2));
         tmp(b).decay=m(2);
         tmp(b).decayerr=sqrt(covm(2,2));
         [m,covm]=wlinem(delaz([idx1(b) idx2(b)],1),...
             log(crampl([idx1(b) idx2(b)])),1,...
-            diag(log(crampl([idx1(b) idx2(b)])...
-            +ramplerr([idx1(b) idx2(b)]).^2)...
-            -log(crampl([idx1(b) idx2(b)]))));
+            diag((log(crampl([idx1(b) idx2(b)])...
+            +ramplerr([idx1(b) idx2(b)]))...
+            -log(crampl([idx1(b) idx2(b)]))).^2));
         tmp(b).cdecay=m(2);
         tmp(b).cdecayerr=sqrt(covm(2,2));
 

Diff for: cmb/slowdecayprofiles.m

@@ -88,9 +88,10 @@
 %        Jan. 31, 2011 - allow no output, odir input, better checks
 %        Feb.  5, 2011 - fix bug when no output specified
 %        Feb. 12, 2011 - include snr-based arrival time error
+%        Feb. 17, 2011 - fixed decay constant error
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Feb. 12, 2011 at 13:35 GMT
+%     Last Updated Feb. 17, 2011 at 13:35 GMT
 
 % todo:
 
@@ -252,11 +253,11 @@
         tmp(cnt).cslow=m(2);
         tmp(cnt).cslowerr=sqrt(covm(2,2));
         [m,covm]=wlinem(delaz(idx,1),log(rampl(idx)),1,...
-            diag(log(rampl(idx)+ramplerr(idx).^2)-log(rampl(idx))));
+            diag((log(rampl(idx)+ramplerr(idx))-log(rampl(idx))).^2));
         tmp(cnt).decay=m(2);
         tmp(cnt).decayerr=sqrt(covm(2,2));
         [m,covm]=wlinem(delaz(idx,1),log(crampl(idx)),1,...
-            diag(log(crampl(idx)+ramplerr(idx).^2)-log(crampl(idx))));
+            diag((log(crampl(idx)+ramplerr(idx))-log(crampl(idx))).^2));
         tmp(cnt).cdecay=m(2);
         tmp(cnt).cdecayerr=sqrt(covm(2,2));
     end

Diff for: docs/seizmo_to_do.txt

@@ -29,19 +29,19 @@ make noise processing more flexible (file length & overlap)
 
 cmt subtoolbox
     - main thing is to plot CMTS!!!
-plotting profiles for models
-    - need a decent polar plotting program to plot slices of the pie
-      with inner & outer radii specified.  Not even mmpolar does it.
-    - need to add raypath plotting support too
 
 more cmb models (rjk)
     - we need to allow parameter setting (while loop rather than for)
         - cell of numbers  {depths/heights  velo%s}
         - inputs will change per model or series ...
         - can never really do everything here ... just go for basics
 expand sph tri stuff
-	- use delaunay codes
     - function to make g matrix for sph tri stuff
+      - distance from polygon center to pnts on line
+        - simple logical for flat weighting
+        - distance based gaussian weighting
+      - inversion is % from prem?
+        - have to do synthetics for all events
 
 
 %%%%%%%

Diff for: exportfig/.ignore/ghostscript.txt

@@ -0,0 +1 @@
+gs

Diff for: exportfig/eps2pdf.m

@@ -0,0 +1,104 @@
+%EPS2PDF  Convert an eps file to pdf format using ghostscript
+%
+% Examples:
+%   eps2pdf source dest
+%   eps2pdf(source, dest, crop)
+%   eps2pdf(source, dest, crop, append)
+%   eps2pdf(source, dest, crop, append, gray)
+%   eps2pdf(source, dest, crop, append, gray, quality)
+%
+% This function converts an eps file to pdf format. The output can be
+% optionally cropped and also converted to grayscale. If the output pdf
+% file already exists then the eps file can optionally be appended as a new
+% page on the end of the eps file. The level of bitmap compression can also
+% optionally be set.
+%
+% This function requires that you have ghostscript installed on your
+% system. Ghostscript can be downloaded from: http://www.ghostscript.com
+%
+%IN:
+%   source - filename of the source eps file to convert. The filename is
+%            assumed to already have the extension ".eps".
+%   dest - filename of the destination pdf file. The filename is assumed to
+%          already have the extension ".pdf".
+%   crop - boolean indicating whether to crop the borders off the pdf.
+%          Default: true.
+%   append - boolean indicating whether the eps should be appended to the
+%            end of the pdf as a new page (if the pdf exists already).
+%            Default: false.
+%   gray - boolean indicating whether the output pdf should be grayscale or
+%          not. Default: false.
+%   quality - scalar indicating the level of image bitmap quality to
+%             output. A larger value gives a higher quality. quality > 100
+%             gives lossless output. Default: ghostscript prepress default.
+
+% Copyright (C) Oliver Woodford 2009-2010
+
+% Suggestion of appending pdf files provided by Matt C at:
+% http://www.mathworks.com/matlabcentral/fileexchange/23629
+
+% Thank you to Fabio Viola for pointing out compression artifacts, leading
+% to the quality setting.
+% Thank you to Scott for pointing out the subsampling of very small images,
+% which was fixed for lossless compression settings.
+
+function eps2pdf(source, dest, crop, append, gray, quality)
+% Intialise the options string for ghostscript
+options = ['-q -dNOPAUSE -dBATCH -sDEVICE=pdfwrite -dPDFSETTINGS=/prepress -sOutputFile="' dest '"'];
+% Set crop option
+if nargin < 3 || crop
+    options = [options ' -dEPSCrop'];
+end
+% Set the grayscale option
+if nargin > 4 && gray
+    options = [options ' -sColorConversionStrategy=Gray -dProcessColorModel=/DeviceGray'];
+end
+% Set the bitmap quality
+if nargin > 5 && ~isempty(quality)
+    options = [options ' -dAutoFilterColorImages=false -dAutoFilterGrayImages=false'];
+    if quality > 100
+        options = [options ' -dColorImageFilter=/FlateEncode -dGrayImageFilter=/FlateEncode -c ".setpdfwrite << /ColorImageDownsampleThreshold 10 /GrayImageDownsampleThreshold 10 >> setdistillerparams"'];
+    else
+        options = [options ' -dColorImageFilter=/DCTEncode -dGrayImageFilter=/DCTEncode'];
+        v = 1 + (quality < 80);
+        quality = 1 - quality / 100;
+        s = sprintf('<< /QFactor %.2f /Blend 1 /HSample [%d 1 1 %d] /VSample [%d 1 1 %d] >>', quality, v, v, v, v);
+        options = sprintf('%s -c ".setpdfwrite << /ColorImageDict %s /GrayImageDict %s >> setdistillerparams"', options, s, s);
+    end
+end
+% Check if the output file exists
+if nargin > 3 && append && exist(dest, 'file') == 2
+    % File exists - append current figure to the end
+    tmp_nam = tempname;
+    % Copy the file
+    copyfile(dest, tmp_nam);
+    % Add the output file names
+    options = [options ' -f "' tmp_nam '" "' source '"'];
+    try
+        % Convert to pdf using ghostscript
+        [status message] = ghostscript(options);
+    catch
+        % Delete the intermediate file
+        delete(tmp_nam);
+        rethrow(lasterror);
+    end
+    % Delete the intermediate file
+    delete(tmp_nam);
+else
+    % File doesn't exist or should be over-written
+    % Add the output file names
+    options = [options ' -f "' source '"'];
+    % Convert to pdf using ghostscript
+    [status message] = ghostscript(options);
+end
+% Check for error
+if status
+    % Report error
+    if isempty(message)
+        error('Unable to generate pdf. Check destination directory is writable.');
+    else
+        error(message);
+    end
+end
+return
+