Fun_RebuildGreenModel.m

% =========================================================================
% Rebuild Green's functions after the fault grid has been updated.
%
% Important:
%   This function only rebuilds Green's functions
%   It must not reprocess observations because 'ob' has already been
%   baseline-corrected, rearranged, and integrated into displacement
% =========================================================================
function [g, locasub, locadep, muSubfault]=Fun_RebuildGreenModel( ...
    PATHg, SDR, grid, gridsize, source, loca, epicenter, ...
    srate, IndexEQ, index, depth, dataVelFolder, nTime)
% Read raw Green's functions for the updated fault geometry
[gRaw, locasub, locadep]=ids_getg( ...
    PATHg, SDR, grid, gridsize, source, loca, ...
    epicenter, srate, IndexEQ, index, depth);

% Number of stations is fixed after station screening
nsta=size(loca,1);

% Convert raw velocity Green's functions to formatted displacement Green's
% functions
g=Fun_SetGreenModel3(gRaw, nsta, grid, srate, nTime);

% Apply the new depth-dependent shear-modulus correction
[g, muSubfault]=Fun_SetGreenModel4(g, locadep, dataVelFolder, grid);

end