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