NAME: jbt_extrema (function) PURPOSE: Find extrema in a numerical array and return their indices. CATEGORIES: CALLING SEQUENCE: result = jbt_extrema(array, interp_nan = interp_nan, min_only = min_only, $ max_only = max_only, threshold = threshold) ARGUMENTS: array: (In, required) The array to find extrema in. KEYWORDS: /interp_nan: If set, remove NaNs by linear interpolation before searching for extrema. /min_only: If set, only return minima. /max_only: If set, only return maxima. threshold: (In, optional) Threshold for changing sense. For example, if threshold is 10 and A[i] and A[i+1] are two adjacent points in a local segment that generally has positive slope, then the segment will be treated as a full positive-slope segment if A[i+1]-A[i] > -10. Default = 0. COMMON BLOCKS: EXAMPLES: ; IDL code example npt = 100 a = randomn(seed, npt) x = findgen(npt) ind = jbt_extrema(a) plot, x, a oplot, x[ind], a[ind], psym = 2, color = 6 SEE ALSO: HISTORY: 2012-11-10: Created by Jianbao Tao (JBT), SSL, UC Berkley. 2012-11-12: Initial release in TDAS. JBT, SSL/UCB. VERSION: $LastChangedBy: jianbao_tao $ $LastChangedDate: 2012-11-12 08:36:20 -0800 (Mon, 12 Nov 2012) $ $LastChangedRevision: 11219 $ $URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/trunk/general/missions/rbsp/efw/utils/jbt_extrema.pro $