EQ_print_bfitFile.m

function [ fID ] = EQ_print_bfitFile (name,data,EQin,sites,vel)

%                           EQ_print_bfitFile.m
%           EQ Function that defines the names of the bfit files
%                     Nathanael Wong Zhixin, Feng Lujia
%
% This function defines the names of the various EQ files that are used for
% post processing and analysis.  The function then creates an input file
% based on the imported data for the earthquake event.
%
% INPUT:
% -- name  : earthquake event name/ID
% -- data  : information on the earthquake event
% -- EQin  : earthquake fault input information
% -- sites : GPS station input information
% -- vel   : GPS displacement information
%
% OUTPUT:
% -- fID   : structure containing information on file names
%
% FORMAT OF CALL: EQ_print_bfitFile (name,data,EQin,sites,vel)
% 
% OVERVIEW:
% 1) This function extracts information on the GTdef-fault and GPS
%    information type
%
% 2) The function then creates the filenames for the respective earthquake
%
% 3) The function then creates an input file based on the information for
%    the earthquake imported from the parent function.
%
% VERSIONS:
% 1) -- Final version validated and commented on 20190811 by Nathanael Wong

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

flt = data.type(1); pnt = data.type(2);

fID.in  = {}; fID.in  = [ name '.in' ];
fID.out = {}; fID.out = [ name '_fwd.out' ];
fID.pch = {}; fID.pch = [ name '_fwd_patches.out' ];
fID.pnt = {}; fID.pnt = [ name '_fwd_point.out' ];
fID.xsc = {}; fID.xsc = [ name '_fwd_xsection.out' ];
fID.fit = {}; fID.fit = [ name '_fit1.txt' ];
fID.fat = {}; fID.fat = [ name '_fit2.txt' ];
fID.mve = {}; fID.mve = [ name '_mve.txt' ];

fEQin = fopen (fID.in,'w');

if     flt == 1, EQ_print_flt1 (fEQin,data,EQin);
elseif flt == 2, EQ_print_flt2 (fEQin,data,EQin);
elseif flt == 3, EQ_print_flt3 (fEQin,data,EQin);
elseif flt == 4, EQ_print_flt4 (fEQin,data,EQin);
end

if     pnt == 1, EQ_print_pnt1 (fEQin,sites,vel);
elseif pnt == 2, EQ_print_pnt2 (fEQin,sites,vel);
elseif pnt == 3, EQ_print_pnt3 (fEQin,sites,vel);
end

fclose (fEQin);

end