EQ_ana_gmfitsw2.m

function [ jj,EQsmfit ] = EQ_ana_gmfitsw2 (Files,i,lon1,lat1,d)

%                              EQ_svanalysis

bfitEQFile = Files(i).name;
EQsmfit = textread(bfitEQFile,'%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;
EQsmfit (a,:) = [];
    
vare = EQsmfit(:,11); [ ~,i ] = max(vare); jj = EQsmfit(i,2);
EQsmfit = textread(bfitEQFile,'%f','commentstyle','shell');
EQsmfit = reshape(EQsmfit,11,[])';

end