%this is matrix that relates shear stress on every fault patch to slip on
%all fault patches
stressG=make_stressG(pm);

%this matrix relates slip on patches to surface displacements through
%   d=dispG*s where d is data (east,north,up,east,north,up,...) and s is slip
%   vector [s(1:end/2) is strike-slip and s(end/2+1:end) is dip-slip]
%xystats is m x 2 vector of m station coordinates in xy 
dispG=make_dispG(pm,xysites);

