Fun_PlotWave.m

function Fun_PlotWave(Sortob, da1, srate,scale)
% =========================================================================
% Plot strong-motion waveforms sorted by epicentral distance
% The function assumes that each column is one waveform component.
% =========================================================================
if nargin < 3 || isempty(srate)
    srate = 4;
end
% Check data dimensions
[numberOfSamples, numberOfChannels] = size(Sortob);
da1 = da1(:);
if numel(da1) ~= numberOfChannels
    error(['The length of da1 must equal the number of columns ', ...
        'in sortedObservedData.']);
end
% Time vector
time = (0:numberOfSamples - 1)' ./ srate;
% Use the distance vector as vertical positions.
dist = da1;
% Scale waveforms so that neighboring traces do not overlap excessively.


figure('Color', 'w');
hold on;
box off;
% Plot each waveform at its corresponding epicentral distance
cols=slanCM('viridis',ceil(length(da1)/3));
k=1;
for iChannel = 1:numberOfChannels
    if mod(iChannel,3)==0
    waveform = Sortob(:, iChannel);
    % Remove the mean only for display.
    waveform = waveform - mean(waveform, 'omitnan');
    plottedWaveform = dist(iChannel) + scale .* waveform;
    plot(time, plottedWaveform, ...
        'k-','LineWidth', 0.3,'Color',cols(k,:));
    k=k+1;
    end
end
xlabel('Time (s)');
ylabel('Epicentral distance (km)');
set(gca, ...
    'FontName', 'Times New Roman', ...
    'FontSize', 18, ...
    'TickDir', 'out','linewidth',0.8,'YDir', 'reverse');
axis('tight');
set(gcf,'position',[200,50,1000,800]);
hold off;
end