
function [Xp,Up] = clean_basilisk_surface_velocity(UU,XX)

    nU = length(UU);

    iall = [1:nU];
    
    ii= find(abs(UU)<2.5);
    UU = UU(ii);
    Xp = XX(ii);
    %    dU = diff(UU);
    %ii = find(dU>0.4);
    Up = UU;
    return
    
    
    
end

