;+ ;PROCEDURE makegap,dg,x,y ;PURPOSE: ; Creates data gaps (by inserting NaN) when the time between data points is ; larger than a value either passed in by the user or calculated to a ; default. ;INPUT: ; dg: If dg is positive, it is the maximum allowed time gap. Any time gaps ; greater than dg will be treated as data gaps. If dg is negative, ; the procedure will calculate a default value for dg of 20 times the ; the smallest time gap in the time series. ; x: The time array. ; y: The data array. ;KEYWORDS: ; v: Optional third dimension array. ; dy: Optional uncertainty in y. ;CREATED BY: Peter Schroeder ;LAST MODIFIED: @(#)makegap.pro 1.2 98/02/18 ;- pro makegap,dg,x,y,dy=dy,v=v dty = size(/type,y) case dty of 4: filly = !values.f_nan 5: filly = !values.d_nan else: return endcase if n_elements(x) le 1 then return diffarray = [x[1:*]-x[0:n_elements(x)-2]] if dg lt 0 then begin posindx = where(diffarray gt 0,poscnt) if poscnt gt 0 then dg = 20d*min(diffarray[posindx]) else return endif gapindx = where(diffarray gt dg,gapcnt) for i=gapcnt-1,0,-1 do begin indx = gapindx[i] x = [x[0:indx],x[indx]+dg,x[indx+1:*]] if size(/n_dimen,y) eq 1 then begin y = [y[0:indx],filly,y[indx+1:*]] if keyword_set(dy) then dy = [dy[0:indx],0.,dy[indx+1:*]] endif else begin y = [y[0:indx,*],replicate(filly,1,dimen2(y)),y[indx+1:*,*]] if keyword_set(dy) then dy = [dy[0:indx,*],$ replicate(0.,1,dimen2(dy)),dy[indx+1:*,*]] endelse if keyword_set(v) then if size(/n_dimen,v) eq 2 then $ v = [v[0:indx,*],v[indx,*],v[indx+1:*,*]] endfor return end