Skip to content

Commit

Permalink
Changes
Browse files Browse the repository at this point in the history
  • Loading branch information
bjmarfito committed Jul 19, 2023
1 parent 100a025 commit 0fdef9b
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions NA/NA_EQ_inverter.m
Original file line number Diff line number Diff line change
Expand Up @@ -83,7 +83,7 @@
if(old)
id = find(misfit==min(misfit));
oldmod = models(:,id);
[oldres,oldslip,oldresnorm,out,oldsynth]=forward_EQ(models(:,id),xytype,resampstruct);
[~,oldslip,oldresnorm,~,oldsynth]=forward_EQ(models(:,id),xytype,resampstruct);
oldres = data'-oldsynth;
end

Expand All @@ -103,7 +103,7 @@
disp(['Setting up ' num2str(ns) ' initial models'])
for j=1:ns
tic
[res,slip] = forward_EQ(models(:,j),xytype,resampstruct);
[res,~] = forward_EQ(models(:,j),xytype,resampstruct);
misfit(j) = sum(res(1:np).^2)/np;
update_time
EQt=EQt+tmp;
Expand All @@ -130,7 +130,7 @@
close(h)

id=find(misfit==min(misfit));
[res,slip,resnorm,out,synth]=forward_EQ(models(:,id),xytype,resampstruct);
[~,slip,resnorm,out,synth]=forward_EQ(models(:,id),xytype,resampstruct);
res=data'-synth;
bigid = find(slip(1:Npatch)==max(slip(1:Npatch)));

Expand All @@ -151,7 +151,7 @@
mu = 3.3e10;
mom = sum(abs(slip(1:end-nramp)))*out(1).L*out(1).W*mu;
M0 = mom*1e7;
Mw = log10(M0)/1.5-10.73
Mw = log10(M0)/1.5-10.73;


%return
Expand Down

0 comments on commit 0fdef9b

Please sign in to comment.