Binning time data

        Function FillGaps(delnom,twave,wlist) 
	// Inserts NaN-points for all missing times
	// in a list of waves, using time-wave twave
	// The waves in the list are then x-scaled.
	// Delnom is the nominal delta-t in seconds.
	// Example: FillGaps(3600.0,tsecs,"ww1;ww2;ww3")
	// REMOVES ALL TIME STEPBACKS AND REPEATS
	// use WaveList("*",";","") to get list.
	// fillgaps(3600.,w1,"w2")
                    
	Variable/D delnom                   // Nominal delta-t (seconds)
	Wave/D twave                        // Time wave, must be double-precision real
	String wlist                        // List of waves, semicolon-separated

	Variable i=0, ii=0, iii=1           // Various internal counters
	Variable numinsert=0, pt    
	Variable/D delta, newtlen           // and variables.
	String wv                           // Will hold wavenames
	// Let user know how much
	// the waves will expand
	pt =  numpnts(twave)
	newtlen = round((twave[pt-1]-twave[0])/delnom)+1
	Print "All waves will grow from ",pt,"to",newtlen,"points."
	Print "Time range:",secs2date(twave[0],0),secs2time(twave[0],3),"to",secs2date(twave[pt-1],0),secs2time(twave[pt-1],3)
	do
		if (i+1==numpnts(twave))			// See if done
			break                                		// and leave
		endif
		delta = twave[i+1]-twave[i]              // Size of delta-t in seconds
    
		do
			if(delta<0.5*delnom)                   // Stepback or repeat?
				print/D twave[i], secs2date(twave[i],0), secs2time(twave[i],0),delta, "StepBack!"
				DeletePoints i+1,1,twave
				ii=0
				do                                   // For each wave in list:
					wv=GetStrFromList(wlist,ii,";")
					if (strlen(wv) == 0)
						break                            // Way to leave when list exhausted
					endif
					wave w=$wv                         // Since this is in a function, we need
					// a wave, not a string
					DeletePoints i+1, 1, w             // Put in the point
					w[i+1]=NaN                         // and set it to NaN
					ii+=1
				while(1)
				delta = twave[i+1] - twave[i]
			else
				break
			endif
		while(delta<0.5*delnom)

		if(delta>1.5*delnom)                     // is there a gap?
			if (delta>5*delnom)
				print/D twave[i], secs2date(twave[i],0), secs2time(twave[i],0),delta, "Step Forward!"
			endif
			ii = 0
			numinsert=round(delta/delnom)-1        // Number of missing points in gap
			InsertPoints i+1,numinsert,twave       // Stick them into twave
			iii=1
			do                                     // and calculate the values 
				twave[i+iii]=twave[i]+delnom*iii
				iii+=1
			while(iii<=numinsert)
			do                                     // For each wave in list:
				wv=GetStrFromList(wlist,ii,";")
				if (strlen(wv) == 0)
					break                              // Way to leave when list exhausted
				endif
				wave w=$wv                           		// Since this is in a function, we need
				// a wave, not a string
				InsertPoints i+1, numinsert, w		// Put in the points
				w[i+1,i+numinsert]=NaN 				// and set them to NaN
				ii+=1
			while(1)
		else
			numinsert=0                            		// Cleaning
		endif                                     			// up and
		i=i+numinsert+1                           	// getting out
	while(1)
	ii=0                                       				// Now set up the wave X-scaling
	do                                          				// by going through the wavelist
		wv=GetStrFromList(wlist,ii,";")		// again
		if (strlen(wv) == 0)
			break
		endif
		wave w=$wv
		SetScale/P x,twave[0],delnom,"dat",w
		ii+=1
	while(1)
end
<\igor>

Proc AverageScaled (win,n,ph,wout, nmin)		// Averages n points in wave win,
	// time from x-scaling, Starting at point ph.
	// Output is to wave wout, with times
	// in x-scaling.   KRP 15 October 1993
	// Wavestats version 9 June 1999
	variable n, ph, nmin
	string win, wout
	prompt win,  "Wave to Average:", popup, Wavelist("*",";","")
	prompt wout, "Output Average Wave:"
	prompt n,    "Averaging Interval (# of points):"
	prompt ph,   "Point Number with which to start"
	prompt nmin,   "Minimum # of points accepted"
	
	variable in,be,en,i,np, ii
	silent 1; pauseupdate
	Make /O /N=(numpnts($win)/n) $wout
	$wout = NaN
	DoAvgScaled($win,n,ph,$wout,nmin)
endMacro

Function DoAvgScaled(win,n,ph,wout, nmin)
	wave win, wout
	variable n, ph, nmin
	
	variable  i=0, ii=0, np=numpnts(win), in, be, en
	do 
		in=i/n;be = in*n+ph;en=min(be+n-1,np-1)
		wout[in] = PSW(win,be,en, nmin)
		i += n
		ii += 1
	while(en<(np-1))
	SetScale/P x pnt2x(win,ph),n*deltax(win),"dat", wout
End

// The following function returns the sum of the y values
// of a wave.  The first works from point number pStart up to
//  and including point number pStop. Uses Wavestats.
//	nmin version, 9 June 1999. RP

Function PSW(w,p1,p2,nmin)
	wave w
	variable p1, p2, nmin

	variable V_npnts, V_avg, V_numNaNs
	Wavestats/q/r=[p1,p2] w
	if (V_numNaNs <= (p2 - p1 + 1 - nmin))
		return (V_Avg)
	endif
	return (NaN)
end