;+NAME/ONE LINE DESCRIPTION OF ROUTINE: ; IQR computes the interquartile range of an input array. ; function iqr,dat id = sort(dat) nd = n_elements(dat) iqr = dat(id(3*nd/4)) - dat(id(nd/4)) return,iqr end