calcrec2.cpp

Go to the documentation of this file.

#include "RecCalc.h"
#include "CCalcRecConf.h"
#include "ResPkModel.h"
#include "SeismicDataComp.h"

#include <iostream>
#include <string>
#include <algorithm>
#include <readline/readline.h>

using namespace std;

int main()
{
        
        int npts;
        double dt;
        double slowness;
        
        CCalcRecConf Config;
        
        Config.GetData("calcrec.conf");
                
        //cout << "Model Filename: ";
        //cin >> modelfilename;
        char *inputline = readline("Model Filename without extension (.mod): ");
        string modelfilename(inputline);
        free(inputline);
        cout << "Dt: ";
        cin >> dt;
        cout << "Npts: ";
        cin >> npts;
        cout << "Slowness: ";
        cin >> slowness;
        SeismicDataComp Result;
        RecCalc Receiver(Config.shift,Config.omega,Config.sigma,Config.cc);
        ResPkModel Model;
        Model.GetData(modelfilename+".mod");
        Model.dt = dt;
        Model.npts = npts;
        Model.slowness = slowness;
        Receiver.CalcRecSynth(modelfilename,Model,Result,false);
        Result.WriteAsSac(modelfilename+".rec");
        Model.WritePlot(modelfilename);
}

---