SeismicModel.h

Go to the documentation of this file.

#ifndef CSEISMICMODEL_H
#define CSEISMICMODEL_H

#include "types.h"

/** \addtogroup seistools Seismic data analysis and modeling */
/* @{ */
class SeismicModel
  {
public:
  enum tArrivalType
    { DirectS, DirectP};
private:
    trealdata SVelErrors;
    trealdata ThickErrors;
    bool
        FuzzComp(const double elem1, const double elem2, const double tolerance);
    trealdata PVelocity;
    trealdata SVelocity;
    trealdata Density;
    trealdata Thickness;
    trealdata Qp;
    trealdata Qs;
    trealdata RecDist;
    trealdata Strike;
    trealdata Dip;
    double SourceDepth;
    double RedVelocity;
    double dt;
    unsigned int npts;
    double slowness;
public:
    double GetSourceDepth() const
      {
        return SourceDepth;
      }
    void SetSourceDepth(const double s){SourceDepth =s;}
    double GetRedVelocity() const
      {
        return RedVelocity;
      }
    void SetRedVelocity(const double s){RedVelocity =s;}
    double GetDt() const
      {
        return dt;
      }
    void SetDt(const double s){dt =s;}
    unsigned int GetNpts() const
      {
        return npts;
      }
    void SetNpts(const unsigned int s){npts =s;}
    double GetSlowness() const
      {
        return slowness;
      }
    void SetSlowness(const double s){slowness =s;}
    const trealdata &GetPVelocity() const {return PVelocity;}
    trealdata &SetPVelocity(){return PVelocity;}
    const trealdata &GetSVelocity() const {return SVelocity;}
    trealdata &SetSVelocity() {return SVelocity;}
    const trealdata &GetDensity() const {return Density;}
    trealdata &SetDensity(){return Density;}
    const trealdata &GetThickness() const {return Thickness;}
    trealdata &SetThickness(){return Thickness;}
    const trealdata &GetQp() const {return Qp;}
    trealdata &SetQp(){return Qp;}
    const trealdata &GetQs() const {return Qs;}
    trealdata &SetQs(){return Qs;}
    const trealdata &GetRecDist() const {return RecDist;}
    trealdata &SetRecDist(){return RecDist;}
    const trealdata &GetStrike() const {return Strike;}
    trealdata &SetStrike(){return Strike;}
    const trealdata &GetDip() const {return Dip;}
    trealdata &SetDip(){return Dip;}
    void SetSVelErrors(const trealdata &sve)
      {
        SVelErrors = sve;
      }
    void SetThickErrors(const trealdata &te)
      {
        ThickErrors = te;
      }
    int FindLayer(const double depth);
    double MatchSlowness(const double slowness, const tArrivalType mode);
    double CalcArrival(const tArrivalType mode, const int recindex);
    double CalcTravelTime(const tArrivalType mode, const double sdepth,
        const double rdepth, const double p);
    virtual void GetData()=0;
    virtual void GetData(const std::string filename)=0;
    virtual void WriteData(const std::string filename)=0;
    virtual void WriteRunFile(const std::string &filename)=0;
    void WritePlot(const std::string &filename);
    void PlotVelWithErrors(const std::string &filename);
    void Init(const int nlayers);
    SeismicModel();
    SeismicModel(const int nlayers);
    SeismicModel(const SeismicModel& source);
    SeismicModel& operator=(const SeismicModel& source);
    virtual ~SeismicModel();
  };
/* @} */
#endif // CSEISMICMODEL_H

---