function select_peak,y,N

;	Returns subscripts corresponding to closest local minimums to
;	a given peak in a curve.

all_peaks=find_peaks(y,/all)

z1=find_peaks(-y,/all)

IF z1(0) gt all_peaks(0) then z1=[0,z1]
IF z1(n_elements(z1)-1) lt all_peaks(n_elements(all_peaks)-1) $
	 then z1=[z1,n_elements(y)-1]


index0=where(z1 lt N)
index0=z1(n_elements(index0)-1)
ind=(where(z1 gt N))(0)

if ind ge 0 then index1=z1(ind) else index1=n_elements(y)-1

return,[index0,index1]

end