Various updates to noise, xc, response & fsspec codes as well

as many bugfixes (too many to list in detail here so look at
the individual files for those)
- noise codes are likely broken at this point because of new
  additions (noise_c3 in particular)
- xc codes underwent many edits for freq-domain and coherency
- fsspec changes related to backprojection
- response code edits as needed for annotated sacpz files
- bugfixes to changeheader & getheader (of minor concern)
- minor updates elsewhere

Author: g2e

Diff for: README

@@ -35,6 +35,15 @@ seismology research.  See the project website for more details:
    may be issues during the uninstall phase with git-based updates so you
    might need to clean your path occasionally.
 
+3.1 If you can't download through Matlab and are on Linux/Unix/OSX, you may
+    run the bash script setup_seizmo.bash via a terminal to download and
+    unpackage external tools for seizmo.
+
+3.2 If you can't or don't want to save seizmo to your paths, you can use
+    startup_seizmo and shutdown_seizmo to dynamically use seizmo when you
+    want.  This is also useful when the functions in seizmo have the same
+    name as functions in another package that you are using.
+
 4. Getting Started
 ------------------
 4a In Matlab: Type 'help seizmo' in the Matlab command window.  Follow

Diff for: docs/seizmo_to_do.txt

@@ -14,6 +14,8 @@ how about installing the latest octave stable branch instead
 taup Sn/Pn no refracted path bug
 taup strange extra P phases bug
 
+
+
 1. install_seizmo improvements
     - extras
         - iscfmdb
@@ -47,29 +49,42 @@ taup strange extra P phases bug
         - ww3struct
           - how to detect jar files are not in classpath?
 2. fss codes
- - fss & fss_capon?
- - compute spectra not centered on 0,0?
-   - complicates cart2pol & pol2cart
-   - input style: [nsmin nsmax nspts],[esmin esmax espts]
-                  [bazmin bazmax bazpts],[smin smax spts]
-                  smax,[espts nspts]
-                  smax,[spts bazpts]
-                  espts,nspts,spts default to 101
-                  bazpts default to 181
- - amplitudes of unwhitened, tdss & capon
-   - 1/(2*pi) factor?
-   - power spectral density is scaled by 2*delta/npts (td measures)
- - averaging across frequency (mainly for capon method)
- - test again!
+   fss fssxc fsshorz fsshorzxc
+    - nearest frequency option
+    - add in block extensions (fhwidth, ntiles)
+      - fhwidth broken due to slowness magnitude streaking
+        - caused by assuming f_center*s==k for all frequencies?
+    - add in range extensions
+    - test
+   tss tssxc tsshorz tsshorzxc
+    - match fss*
+   arf arfhorz
+    - add in range extensions
+    - test
+   isfss chkfss
+   fsssub fssidx fssgrids
+    - fsssub needs slowness ranging
+   fssavg fssbackproject fsscart2pol fsspol2cart
+    - pol/cart need range extensions, better defaults and use spline by default
+    - fssbackproject from fk/fs_whateveritiscalled
+   fsscorrcoef fssdbinfo
+   plotarf plotfss plotfssupdate
+    - plot option passthrough
+   fssframeslide fssfreqslide
+
+ - an arrays function (yellowknife/lasa)
+ - fssphasecorr should be arrayphasecorr
+ - horz capon formula
+   - invert component matrices and rotate
+   - invert each rotated matrix
 3. geofss codes (use mmap image option)
  - rewrites
    - geo(fssxc,tdssxc)
      - antipodal issue?
  - write
    - geo(fsshorz,fsshorzxc,tdsshorzxc,arfhorz)
      - much of this should come from non-geo versions
- - test!
-4. test & tutorials for fss codes
+4. test & tutorials for (geo)fss codes
 5. 2d inversion on a sphere
     - cubed sphere
         - book keeping for adjacent blocks

Diff for: fk/fkxchorzmap.m

@@ -80,7 +80,7 @@
 %     % One way to perform horizontal fk analysis:
 %     data=rotate(data,'to',0,'kcmpnm1','N','kcmpnm2','E');
 %     xc=correlate(data,'mcxc');
-%     [axc,xc]=split_auto_correlations(xc,false);
+%     [axc,xc]=split_auto_correlations(xc);
 %     [b,e,delta]=getheader(xc,'b','e','delta');
 %     xc=interpolate(xc,delta(1),'spline',max(b),min(e));
 %     xc=rotate_correlations(xc);
@@ -98,9 +98,10 @@
 %        Jan. 29, 2013 - example update for new correlate functions
 %        Sep. 23, 2013 - example update for new correlate functions,
 %                        adjust for rotate_correlations change
+%        May  30, 2014 - adjusted example for split_auto_correlations fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Sep. 23, 2013 at 15:25 GMT
+%     Last Updated May  30, 2014 at 15:25 GMT
 
 % todo:
 

Diff for: fk/fkxchorzvolume.m

@@ -80,7 +80,7 @@
 %     % One way to perform horizontal fk analysis:
 %     data=rotate(data,'to',0,'kcmpnm1','N','kcmpnm2','E');
 %     xc=correlate(data,'mcxc');
-%     [axc,xc]=split_auto_correlations(xc,false);
+%     [axc,xc]=split_auto_correlations(xc);
 %     [b,e,delta]=getheader(xc,'b','e','delta');
 %     xc=interpolate(xc,delta(1),'spline',max(b),min(e));
 %     xc=rotate_correlations(xc);
@@ -107,9 +107,10 @@
 %        Jan. 29, 2013 - example update for new correlate functions
 %        Sep. 23, 2013 - example update for new correlate functions,
 %                        adjust for rotate_correlations change
+%        May  30, 2014 - adjusted example for split_auto_correlations fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Sep. 23, 2013 at 15:25 GMT
+%     Last Updated May  30, 2014 at 15:25 GMT
 
 % todo: