calcrec.cpp

Go to the documentation of this file.

#include "RecCalc.h"
#include "CCalcRecConf.h"
#include "SeismicDataComp.h"
#include "CFatalException.h"
#include <iostream>
#include <string>
#include <algorithm>
#include "Util.h"

/*! \file 
 * This program calculates a receiver function from single radial and vertical components.
 * When started the program only asks for the name of the files for these two components. 
 * The output receiver function has the same name as the radial component with .rec appended.
 * All other aspects of the calculation are specified in the file calcrec.conf in the directory
 * in which the program is executed.
 * 
 * The parameters for that file are as follows:
 * 
 * <TABLE>
 * <TR> <TD> <B>Name </B></TD> <TD> <B>Description </B></TD></TR>
 * <TR> <TD> recmethod </TD>  <TD> Can be specdiv for spectral division method or iterdecon for iterative deconvolution. </TD></TR>
 * <TR> <TD> omega </TD>  <TD>The center of the gaussian used to filter the receiver function during estimation. </TD> </TR>
 * <TR> <TD> sigma </TD>  <TD>The width of the gaussian used to filter the receiver function during estimation. 
 *    Larger values mean a more broad band reciever function that contains more detailed structure but is also susceptible to noise.</TD> </TR>
 * <TR> <TD> shift </TD>  <TD>The shift in seconds of the initial correlation peak. This is purely used for display purposes. </TD> </TR>
 * <TR> <TD> cc </TD>  <TD>(only specdiv) The waterlevel relative to the maximum autopower of the vertical component. </TD> </TR>
 * </TABLE> 
 * More details about these parameters can be found for example in 
 * Langston, 1979, JGR 84(B9) 4749-4762
 * Liggoria, 1999, BSSA 89(5), 1395-1400
 */
using namespace std;

string version = "$Id: calcrec.cpp 1663 2008-02-13 17:24:46Z mmoorkamp $";

int main(int argc, char *argv[]) {
        cout << " This is datarec: Calculate receiver functions from input data"
                        << endl; // write some info
        cout << " Reads in radial and vertical components in SAC format" << endl;
        cout
                        << " Outputs a file with the name of the radial component and .rec appended"
                        <<endl;
        cout << " Some behaviour can be configured with the file calcrec.conf"
                        << endl;
        cout << " This is Version: " << version << endl << endl;
        string modelfilename;

        CCalcRecConf Config;
        SeismicDataComp Radial, Vertical;
        ResPkModel Model;
        string outfilename;
        try
        {

                Config.GetData("calcrec.conf");
                RecCalc::trfmethod rfmethod = RecCalc::specdiv;
                if (Config.recmethod == "iterdecon")
                {
                        rfmethod = RecCalc::iterdecon;
                }
                string radfilename, verfilename;
                if (argc == 3) {
                        radfilename = argv[1];
                        verfilename = argv[2];
                } else {
                        radfilename = AskFilename("Radial Component: ");
                        verfilename = AskFilename("Vertical Component: ");
                }
                outfilename = radfilename+".rec";
                Radial.GetData(radfilename);
                Vertical.GetData(verfilename);

                RecCalc Receiver(Config.shift,Config.omega,Config.sigma,Config.cc,false,rfmethod);
                Receiver.SetNormalize(true);
                SeismicDataComp RecFunc;
                Receiver.CalcRecData(Radial,Vertical,RecFunc);

                RecFunc.WriteAsSac(outfilename);
        }
        catch (CFatalException &e)
        {
                cerr << e.what() << endl; // if something fails print error
                return -1; // and stop execution
        }
}

---