function d_MEDIAN, vec
;
;   replacement for IDL's MEDIAN, called only by median_filtered_mean.pr0
;
;   this replacement is a cut down version designed to do just what I need it to do...
;
;   DPC   20051023  original code
;
nv = size(vec)
nv2 = round(nv[1]/2)
ordered_subscripts = sort(vec)
med = vec[ordered_subscripts[nv2]]
;
return, med
end