11 using namespace
13 string version
= "$Id: rfvel.cpp 1839 2010-03-05 17:04:34Z mmoorkamp $"
<< " This is rfvel: Calculate absolute velocties from receiver functions"
cout << " Reads in radial and vertical components in SAC format"
<< " Writes a file with ending .vel containing apparent velocities and corresponding periods"
<< " The settings for the calculation method, omega, sigma and shift are taken from calcrec.conf"
cout << " This is Version: "
<< endl << endl;
rfmethod = RecCalc::iterdecon;
radfilename = AskFilename("Radial Component: "
verfilename = AskFilename("Vertical Component: "
outfilename = radfilename + ".vel"
cout << "Slowness in s/km: "
RFVel.CalcRFVel(slowness, Radial, Vertical, AppVel);
cerr << e.what() << endl;
There are several ways to calculate receiver functions.
int ReadData(const std::string &filename, tseismicdataformat format=sac)
Read in data from a file, as we cannot determine the type from the ending we have to provide it...
void GetData(std::string filename)
This class implements the method to calculate absolute S-Wave velocities from Receiver function data ...
The basic exception class for all errors that arise in gplib.