EQ_ana_bmfitsv3.m

function [ i,jj,EQsmfit ] = EQ_ana_bmfits3 (Files,lon1,lat1,lon2,lat2,d)

%                              EQ_svanalysis

mat_vare = zeros(length(Files),11);
parfor ii = 1 : length(Files)
    
    iiEQFile = Files(ii).name;
    EQsmfit = textread(['./GTdef_SV_Misfit/' iiEQFile],'%f','commentstyle','shell');
    EQsmfit = reshape(EQsmfit,11,[])';
    
    lon = EQsmfit(:,3); lat = EQsmfit(:,4);
    a = lon<lon1+d & lon>lon1-d & lat<lat1+d & lat>lat1-d;
    b = lon<lon2+d & lon>lon2-d & lat<lat2+d & lat>lat2-d;
    EQsmfit (a|b,:) = [];
    
    vare = EQsmfit(:,9); [ ~,i ] = max(vare);
    mat_vare(ii,:) = EQsmfit(i,:);
    
end
vare = mat_vare(:,9); [ ~,i ] = max(vare); jj = mat_vare(i,2);
bfitEQFile = Files(i).name;
EQsmfit = textread(['./GTdef_SV_Misfit/' bfitEQFile],'%f','commentstyle','shell');
EQsmfit = reshape(EQsmfit,11,[])';

end