-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathEQ_ana_svbfit.m
executable file
·67 lines (55 loc) · 2.59 KB
/
EQ_ana_svbfit.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
function EQ_ana_svbfit(EQID,data,sites,vel,d)
% EQ_ana_sv
EQID = strrep (EQID,'co','');
regID = data.type(5); mu = data.slip(2);
regname = EQ_create_printregname (regID);
File = [ EQID '_' regname '_' 'Rig' num2str(mu) '_svbfit' ];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PEAK 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
cd GTdef_SV_Misfit; EQmfitFiles = dir('*.txt'); cd ..
[ i,jj,EQsvmfit ] = EQ_ana_bmfitsv1 (EQmfitFiles);
cd GTdef_L_Misfit; EQmfitFiles = dir('*.txt');
[ EQlmfit ] = EQ_ana_mfitl (EQmfitFiles,i); cd ..
cd GTdef_SV_Input; EQinFiles = dir('*.txt');
[ lon1,lat1,EQin ] = EQ_ana_ins (EQinFiles,i,jj); cd ..
if exist('Analysis_ve1','dir') ~= 7; mkdir ('Analysis_ve1'); end
cd Analysis_ve1;
iFile = [ File '1' ];
[ Files ] = EQ_print_bfitFile (iFile,data,EQin,sites,vel);
GTdef(Files.in,0); GTdef_project(Files.out);
[ ~ ] = EQ_print_svbfit (Files.fit,EQsvmfit,data,1);
[ ~ ] = EQ_print_svbfit (Files.fat,EQlmfit,data,1);
[ ~ ] = EQ_print_maxve1 (Files.mve,EQsvmfit(jj,:),data,1);
cd ..
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PEAK 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
cd GTdef_SV_Misfit; EQmfitFiles = dir('*.txt'); cd ..
[ i,jj,EQsvmfit ] = EQ_ana_bmfitsv2 (EQmfitFiles,lon1,lat1,d);
cd GTdef_L_Misfit; EQmfitFiles = dir('*.txt');
[ EQlmfit ] = EQ_ana_mfitl (EQmfitFiles,i); cd ..
cd GTdef_SV_Input; EQinFiles = dir('*.txt');
[ lon2,lat2,EQin ] = EQ_ana_ins (EQinFiles,i,jj); cd ..
if exist('Analysis_ve1','dir') ~= 7; mkdir ('Analysis_ve1'); end
cd Analysis_ve1;
iFile = [ File '2' ];
[ Files ] = EQ_print_bfitFile (iFile,data,EQin,sites,vel);
GTdef(Files.in,0); GTdef_project(Files.out);
[ ~ ] = EQ_print_svbfit (Files.fit,EQsvmfit,data,2);
[ ~ ] = EQ_print_svbfit (Files.fat,EQlmfit,data,2);
[ ~ ] = EQ_print_maxve1 (Files.mve,EQsvmfit(jj,:),data,2);
cd ..
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PEAK 3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
cd GTdef_SV_Misfit; EQmfitFiles = dir('*.txt'); cd ..
[ i,jj,EQsvmfit ] = EQ_ana_bmfitsv3 (EQmfitFiles,lon1,lat1,lon2,lat2,d);
cd GTdef_L_Misfit; EQmfitFiles = dir('*.txt');
[ EQlmfit ] = EQ_ana_mfitl (EQmfitFiles,i); cd ..
cd GTdef_SV_Input; EQinFiles = dir('*.txt');
[ ~,~,EQin ] = EQ_ana_ins (EQinFiles,i,jj); cd ..
if exist('Analysis_ve1','dir') ~= 7; mkdir ('Analysis_ve1'); end
cd Analysis_ve1;
iFile = [ File '3' ];
[ Files ] = EQ_print_bfitFile (iFile,data,EQin,sites,vel);
GTdef(Files.in,0); GTdef_project(Files.out);
[ ~ ] = EQ_print_svbfit (Files.fit,EQsvmfit,data,3);
[ ~ ] = EQ_print_svbfit (Files.fat,EQlmfit,data,3);
[ ~ ] = EQ_print_maxve1 (Files.mve,EQsvmfit(jj,:),data,3);
cd ..
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%