median_.pro
504 Bytes
function median_, Data, T_AVER
common InternalCommon, FinalTime, num
Npoints = N_elements(FinalTime);
Nreal = n_elements(num);
Timy = FinalTime(num);
FillValue = !Values.F_NAN
outPut = replicate(FillValue, Nreal); outPut = fltarr(Nreal);
for tt = FinalTime[0], FinalTime[Npoints-1]+T_AVER, T_AVER do begin
nummi = where(abs(tt - Timy) LE T_AVER/2.0, Index)
if (Index GT 0) then OutPut[nummi] = median(Data[nummi])
endfor
return, OutPut
end