function[Data]= FilterVx(location,Data,filter)
%Subfunction of readData, Filters away all datapoints that have a meridional wind |maginitude| bigger then "filter"
%%
    a=size(Data.Year,2);
    keep=find((Data.Vx<=filter)&(Data.Vx>=-filter));
    Data.Year    = Data.Year(keep);
    Data.Month   = Data.Month(keep);
    Data.Day     = Data.Day(keep);
    Data.Hour    = Data.Hour(keep);
    Data.NMeteor = Data.NMeteor(keep);
    Data.Frang   = Data.Frang(keep);
    Data.Rsep    = Data.Rsep(keep);
    Data.Vx      = Data.Vx(keep);
    Data.Vy      = Data.Vy(keep);
    Data.Lat     = Data.Lat(keep);
    Data.Lon     = Data.Lon(keep);
    Data.Vm      = Data.Vm(keep);
    Data.Vm_lat  = Data.Vm_lat(keep);
    Data.Vm_lon  = Data.Vm_lon(keep);
    Data.SD_Vx   = Data.SD_Vx(keep);
    Data.SD_Vy   = Data.SD_Vy(keep);
    Data.Filter = filter;
    %Data.Name = [Data.Name,'_FilterVx',int2str(filter)];
    fprintf(['# of filtered data: ' int2str(a),'\n'])
    %cd C:\Users\kleinkne\Documents\PHD\DATA
    %cd(location)
    %save(Data.Name,'Data')

    %cd ..\
end