- updated many functions to make discrete time vectors properly

  - probably a bug or two made by making all these changes
- added misc/maketime which shows how to do this
- added misc/freqcount for occurance counting purposes

Author: g2e

Diff for: docs/seizmo_to_do.txt

@@ -5,11 +5,6 @@ way to make movies while you do other things on the screen:
 ## hf = figure(‘visible’, ‘off’)
 ## addframe(aviobj, hf)
 
-times from b+npts+delta:
-   float issues: b+(0:npts-1)*delta (THIS IS WHAT I'VE BEEN USING!)
-no float issues: b+(0:delta:delta*(npts-1))  (MERGE OF UNEVEN USES THIS!)
-- grep 0:npts */*.m
-
 octave_364_fixes in install_seizmo & startup
 mp_tmerc goes from m_map_fixes to m_map_fixes/private in webinstall_mmap
 figure out best way to use scatter in octave
@@ -19,33 +14,31 @@ how about installing the latest octave stable branch instead
 taup Sn/Pn no refracted path bug
 taup strange extra P phases bug
 
-3. isc_fm
-4. sacpz
-5. maketimes issue
+download_extras.bash
 
 1. install_seizmo improvements
     - extras
-        x webinstall_extras
-        x uninstall_extras
         - iscfmdb
           x parse_isc_fm
           - clean_isc_fm
             - no depth -> out
             - no mag -> out
             - fill in remaining fields
           - update_isc_fm
-        x readsacpz_rdseed
-        x writesacpz_rdseed
-        x issacpz_rdseed
-        x fix_old_sacpz
-        - update sacpzdb
+          - compat w/ findcmt findcmts plotmt mapcmts etc ?
+        - sacpz
+          x readsacpz_rdseed
+          x writesacpz_rdseed
+          x issacpz_rdseed
+          x fix_old_sacpz
+          - update sacpzdb
             - makesacpz
             - getsacpz (sacpz ws support)
             - db2sacpz
             - (remove|apply)sacpz
             - iris_sacpzdb_build
             - iris_sacpzdb_fixes
-        - station_check (uses new sacpzdb)
+          - station_check (uses new sacpzdb)
         - add new mapping features?
             - natural earth: http://www.naturalearthdata.com/downloads/
             - gplates stuff
@@ -62,8 +55,6 @@ taup strange extra P phases bug
                 - add an ani flag of some sort (should indicate type too)
                 - dependent functions need to handle flag logic and math
             - add more models shared via iris emc
-    - install_seizmo
-    - uninstall_seizmo
     - any function that requires a jar does
       a check first and tries to add it to
       the dynamic path if not found

Diff for: ftran/omegamultiply.m

@@ -44,9 +44,10 @@
 %        Feb. 11, 2011 - mass nargchk fix
 %        Feb.  4, 2012 - doc update, multiplyomega to omegamultiply
 %        Feb. 14, 2013 - use strcmpi for consistency
+%        Mar.  1, 2014 - maketime fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Feb. 14, 2013 at 15:05 GMT
+%     Last Updated Mar.  1, 2014 at 15:05 GMT
 
 % todo:
 
@@ -99,7 +100,8 @@
 
         % differentiate
         cols=size(data(i).dep,2)/2;
-        omega=2*pi*delta(i)*[0:npts2(i) (npts2(i)-1):-1:1].';
+        omega=2*pi*delta(i);
+        omega=[0:omega:omega*npts2(i) (npts2(i)-1)*omega:-omega:omega].';
         if(strcmpi(iftype(i),'irlim'))
             % rlim %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
             % 0Hz real/imag == 0 else

Diff for: ftran/omegashift.m

@@ -33,9 +33,10 @@
 %        Mar. 15, 2012 - fix example
 %        May  31, 2012 - minor doc update
 %        Feb. 14, 2013 - use strcmpi for consistency
+%        Mar.  1, 2014 - maketime fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Feb. 14, 2013 at 15:05 GMT
+%     Last Updated Mar.  1, 2014 at 15:05 GMT
 
 % todo:
 
@@ -96,7 +97,8 @@
 
         % time shift
         cols=size(data(i).dep,2)/2;
-        wt=shift(i)*2*pi*delta(i)*(0:npts(i)-1).';
+        wt=shift(i)*2*pi*delta(i);
+        wt=(0:wt:wt*(npts(i)-1)).';
         if(strcmpi(iftype(i),'irlim'))
             data(i).dep(:,[1:2:end 2:2:end])=...
                 [data(i).dep(:,1:2:end).*cos(wt(:,ones(1,cols)))...

Diff for: ftran/splitpad.m

@@ -41,9 +41,10 @@
 %                        tolerance to the point at 0 requirement
 %        Dec. 21, 2011 - doc update
 %        May  29, 2012 - pow2pad=0 by default
+%        Mar.  1, 2014 - maketime fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated May  29, 2012 at 10:05 GMT
+%     Last Updated Mar.  1, 2014 at 10:05 GMT
 
 % todo:
 
@@ -96,7 +97,8 @@
     depmen=nan(nrecs,1); depmin=depmen; depmax=depmen;
     for i=1:nrecs
         % check has point at 0 (within some tolerance)
-        idx=find(abs(b(i)+(0:npts(i)-1)*delta(i))<delta(i)/10,1);
+        idx=find(...
+            abs(b(i)+(0:delta(i):delta(i)*(npts(i)-1)))<delta(i)/10,1);
         if(isempty(idx))
             error('seizmo:splitpad:noZeroPoint',...
                 'Record %d does not have a point at zero!',i);

Diff for: hdr/makeuneven.m

@@ -32,9 +32,10 @@
 %        Dec. 30, 2010 - fixed 1-liner description
 %        Feb. 11, 2011 - mass nargchk fix
 %        Jan. 30, 2012 - doc update, better checkheader/getheader usage
+%        Mar.  1, 2014 - maketime fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Jan. 30, 2012 at 15:05 GMT
+%     Last Updated Mar.  1, 2014 at 15:05 GMT
 
 % todo:
 
@@ -71,7 +72,7 @@
         % loop over even, add .ind
         idx=find(even);
         for i=1:sum(even)
-            data(idx(i)).ind=b(i)+(0:npts(i)-1)'*delta(i);
+            data(idx(i)).ind=b(i)+(0:delta(i):delta(i)*(npts(i)-1)).';
         end
     end
 

Diff for: forge/freqcount.m renamed to misc/freqcount.m

@@ -308,12 +308,12 @@
 %                        debugging output, unevenoverlap option, uneven &
 %                        multi-rate merging
 %        Feb. 22, 2014 - overlap blend method, unevenoverlap ignore method
+%        Mar.  1, 2014 - maketime fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Feb. 22, 2014 at 15:05 GMT
+%     Last Updated Mar.  1, 2014 at 15:05 GMT
 
 % todo:
-% - ignore unevenoverlap method
 
 % check nargin
 if(mod(nargin-1,2))
@@ -1320,9 +1320,9 @@
 end
 
 % interpolate gap
-lasttimes=ab(last,2)+(0:(npts(last)-1))*delta;
-gaptimes=sae(2)+(1:nsamples)*delta;
-newtimes=sab(2)+(0:(npts(first)-1))*delta;
+lasttimes=ab(last,2)+(0:delta:delta*(npts(last)-1));
+gaptimes=sae(2)+(delta:delta:delta*nsamples);
+newtimes=sab(2)+(0:delta:delta*(npts(first)-1));
 gapdata=interp1([newtimes lasttimes],...
     [data(first).dep; data(last).dep],...
     gaptimes,option.INTERPOLATE,'extrap');
@@ -1365,9 +1365,9 @@
 end
 
 % interpolate gap
-firsttimes=ab(first,2)+(0:(npts(first)-1))*delta;
-gaptimes=ae(first,2)+(1:nsamples)*delta;
-newtimes=sab(2)+(0:(npts(last)-1))*delta;
+firsttimes=ab(first,2)+(0:delta:delta*(npts(first)-1));
+gaptimes=ae(first,2)+(delta:delta:delta*nsamples);
+newtimes=sab(2)+(0:delta:delta*(npts(last)-1));
 gapdata=interp1([firsttimes newtimes],...
     [data(first).dep; data(last).dep],...
     gaptimes,option.INTERPOLATE,'extrap');
@@ -1410,10 +1410,10 @@
 end
 
 % interpolate first record and gap
-firsttimes=ab(first,2)+(0:(npts(first)-1))*delta;
-lasttimes=ab(last,2)+(0:(npts(last)-1))*delta;
-gaptimes=sae(2)+(1:nsamples)*delta;
-newtimes=sab(2)+(0:(npts(first)-1))*delta;
+firsttimes=ab(first,2)+(0:delta:delta*(npts(first)-1));
+lasttimes=ab(last,2)+(0:delta:delta*(npts(last)-1));
+gaptimes=sae(2)+(delta:delta:delta*nsamples);
+newtimes=sab(2)+(0:delta:delta*(npts(first)-1));
 newdata=interp1([firsttimes lasttimes],...
     [data(first).dep; data(last).dep],...
     [newtimes gaptimes],option.INTERPOLATE,'extrap');
@@ -1456,10 +1456,10 @@
 end
 
 % interpolate last record and gap
-firsttimes=ab(first,2)+(0:(npts(first)-1))*delta;
-lasttimes=ab(last,2)+(0:(npts(last)-1))*delta;
-gaptimes=ae(first,2)+(1:nsamples)*delta;
-newtimes=sab(2)+(0:(npts(last)-1))*delta;
+firsttimes=ab(first,2)+(0:delta:delta*(npts(first)-1));
+lasttimes=ab(last,2)+(0:delta:delta*(npts(last)-1));
+gaptimes=ae(first,2)+(delta:delta:delta*nsamples);
+newtimes=sab(2)+(0:delta:delta*(npts(last)-1));
 newdata=interp1([firsttimes lasttimes],...
     [data(first).dep; data(last).dep],...
     [gaptimes newtimes],option.INTERPOLATE,'extrap');
@@ -1570,8 +1570,8 @@
 sae=ae(first,:)+[0 shift];
 
 % interpolate first record
-oldtimes=ab(first,2)+(0:(npts(first)-1))*delta;
-newtimes=sab(2)+(0:(npts(first)-1))*delta;
+oldtimes=ab(first,2)+(0:delta:delta*(npts(first)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(first)-1));
 data(first).dep=...
     interp1(oldtimes,data(first).dep,newtimes,option.INTERPOLATE,'extrap');
 data(first).dep=data(first).dep.';
@@ -1615,8 +1615,8 @@
 sae=ae(last,:)-[0 shift];
 
 % interpolate last record
-oldtimes=ab(last,2)+(0:(npts(last)-1))*delta;
-newtimes=sab(2)+(0:(npts(last)-1))*delta;
+oldtimes=ab(last,2)+(0:delta:delta*(npts(last)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(last)-1));
 data(last).dep=...
     interp1(oldtimes,data(last).dep,newtimes,option.INTERPOLATE,'extrap');
 data(last).dep=data(last).dep.';
@@ -1735,8 +1735,8 @@
 sae=ae(first,:)+[0 shift];
 
 % interpolate first record
-oldtimes=ab(first,2)+(0:(npts(first)-1))*delta;
-newtimes=sab(2)+(0:(npts(first)-1))*delta;
+oldtimes=ab(first,2)+(0:delta:delta*(npts(first)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(first)-1));
 data(first).dep=...
     interp1(oldtimes,data(first).dep,newtimes,option.INTERPOLATE,'extrap');
 data(first).dep=data(first).dep.';
@@ -1780,8 +1780,8 @@
 sae=ae(last,:)-[0 shift];
 
 % interpolate last record
-oldtimes=ab(last,2)+(0:(npts(last)-1))*delta;
-newtimes=sab(2)+(0:(npts(last)-1))*delta;
+oldtimes=ab(last,2)+(0:delta:delta*(npts(last)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(last)-1));
 data(last).dep=...
     interp1(oldtimes,data(last).dep,newtimes,option.INTERPOLATE,'extrap');
 data(last).dep=data(last).dep.';
@@ -1908,8 +1908,8 @@
 sae=ae(first,:)+[0 shift];
 
 % interpolate first record
-oldtimes=ab(first,2)+(0:(npts(first)-1))*delta;
-newtimes=sab(2)+(0:(npts(first)-1))*delta;
+oldtimes=ab(first,2)+(0:delta:delta*(npts(first)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(first)-1));
 data(first).dep=...
     interp1(oldtimes,data(first).dep,newtimes,option.INTERPOLATE,'extrap');
 data(first).dep=data(first).dep.';
@@ -1957,8 +1957,8 @@
 sae=ae(last,:)-[0 shift];
 
 % interpolate last record
-oldtimes=ab(last,2)+(0:(npts(last)-1))*delta;
-newtimes=sab(2)+(0:(npts(last)-1))*delta;
+oldtimes=ab(last,2)+(0:delta:delta*(npts(last)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(last)-1));
 data(last).dep=...
     interp1(oldtimes,data(last).dep,newtimes,option.INTERPOLATE,'extrap');
 data(last).dep=data(last).dep.';
@@ -2092,8 +2092,8 @@
 sae=ae(first,:)+[0 shift];
 
 % interpolate first record
-oldtimes=ab(first,2)+(0:(npts(first)-1))*delta;
-newtimes=sab(2)+(0:(npts(first)-1))*delta;
+oldtimes=ab(first,2)+(0:delta:delta*(npts(first)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(first)-1));
 data(first).dep=...
     interp1(oldtimes,data(first).dep,newtimes,option.INTERPOLATE,'extrap');
 data(first).dep=data(first).dep.';
@@ -2142,8 +2142,8 @@
 sae=ae(last,:)-[0 shift];
 
 % interpolate last record
-oldtimes=ab(last,2)+(0:(npts(last)-1))*delta;
-newtimes=sab(2)+(0:(npts(last)-1))*delta;
+oldtimes=ab(last,2)+(0:delta:delta*(npts(last)-1));
+newtimes=sab(2)+(0:delta:delta*(npts(last)-1));
 data(last).dep=...
     interp1(oldtimes,data(last).dep,newtimes,option.INTERPOLATE,'extrap');
 data(last).dep=data(last).dep.';

Diff for: forge/maketime.m renamed to misc/maketime.m

@@ -185,9 +185,10 @@
 %        Jan. 30, 2012 - drop SEIZMO global
 %        Feb.  7, 2012 - fix no horizontals bug
 %        Mar. 13, 2012 - use getheader improvements
+%        Mar.  1, 2014 - maketime fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Mar. 13, 2012 at 23:18 GMT
+%     Last Updated Mar.  1, 2014 at 23:18 GMT
 
 % todo:
 % - a 'norotate' would be nice
@@ -607,10 +608,12 @@
             % get independent component
             t1=cell(n1,1); t2=cell(n2,1);
             for j=1:n1
-                t1{j}=cb1(j)+(0:npts(cidx1(j))-1).'*delta(pidx(1));
+                t1{j}=cb1(j)+...
+                    (0:delta(pidx(1)):delta(pidx(1))*(npts(cidx1(j))-1)).';
             end
             for j=1:n2
-                t2{j}=cb2(j)+(0:npts(cidx2(j))-1).'*delta(pidx(1));
+                t2{j}=cb2(j)+...
+                    (0:delta(pidx(1)):delta(pidx(1))*(npts(cidx2(j))-1)).';
             end
 
             % loop over all possible pairings
@@ -846,25 +849,25 @@
 
 % get new times
 if(move1)
-    t=t2(bidx2+(0:npts-1));
+    t=t2(bidx2+(0:(npts-1)));
 else % move 2
-    t=t1(bidx1+(0:npts-1));
+    t=t1(bidx1+(0:(npts-1)));
 end
 
 % get x1,x2
 if(shift<=shiftmax)
     % ok we are just shifting one record to time align
     if(option.DEBUG); disp('Adjust Method: SHIFT'); end
-    x1=x1(bidx1+(0:npts-1));
-    x2=x2(bidx2+(0:npts-1));
+    x1=x1(bidx1+(0:(npts-1)));
+    x2=x2(bidx2+(0:(npts-1)));
 else % interpolate
     % interpolating one record to time align
     if(option.DEBUG); disp('Adjust Method: INTERPOLATE'); end
     if(move1)
         x1=interp1(t1,x1,t,option.INTERPOLATE,'extrap');
-        x2=x2(bidx2+(0:npts-1));
+        x2=x2(bidx2+(0:(npts-1)));
     else % interpolate 2
-        x1=x1(bidx1+(0:npts-1));
+        x1=x1(bidx1+(0:(npts-1)));
         x2=interp1(t2,x2,t,option.INTERPOLATE,'extrap');
     end
 end

Diff for: multi/meld.m

@@ -40,9 +40,10 @@
 %        Apr. 15, 2013 - initial version, no 50% data in subwindow bugfix
 %        Apr. 17, 2013 - use powerspectraldensity to get psd
 %        Jan. 26, 2014 - abs path exist fix
+%        Mar.  1, 2014 - maketime fix
 %
 %     Written by Garrett Euler (ggeuler at wustl dot edu)
-%     Last Updated Jan. 26, 2014 at 13:30 GMT
+%     Last Updated Mar.  1, 2014 at 13:30 GMT
 
 % todo:
 
@@ -204,7 +205,7 @@
         % grab frequency once
         if(~exist('f','var'))
             [npts,delta]=getheader(swdata(1),'npts','delta');
-            f=(0:npts-1)'*delta;
+            f=(0:delta:delta*(npts-1)).';
         end
 
         % append spectra to storage cells (one per record)