C1DAnisoMTSynthData.h

Go to the documentation of this file.

#ifndef C1DANISOMTSYNTHDATA_H_
#define C1DANISOMTSYNTHDATA_H_

#include "MTStation.h"
#include "VecMat.h"
//! Calculate response of a 1D anisotropic model, code is based on Pek and Santos fortran code
class C1DAnisoMTSynthData : public MTStation
{
private:
        trealdata calc_frequencies;
        trealdata thicknesses;  // Thickness for each layer 
        trealdata strikes;
        trealdata slants;
        trealdata dips;
        trealdata rho1;
        trealdata rho2;
        trealdata rho3;
        trealdata effcond1;
        trealdata effcond2;
        trealdata effstrike;
        //! Transform angles and resistivities into effective horizontal values, has to be called before CalcZ
        void TransformRho();
        //! Calculate the response
        void CalcZ();
public:
        //! Set the anisotropy strike for each layer in degree
        void SetStrikes(const trealdata &a){strikes = a;}
        //! Set the anisotropy slant for each layer in degree
        void SetSlants(const trealdata &a){slants = a;}
        //! Set the anisotropy dip for each layer in degree
        void SetDips(const trealdata &a){dips = a;}
        //! Set the first principal resistivity for each layer in Ohm.m
        void SetRho1(const trealdata &a){rho1 = a;}
        //! Set the second principal resistivity for each layer in Ohm.m
        void SetRho2(const trealdata &a){rho2 = a;}
        //! Set the first principal resistivity for each layer in Ohm.m
        void SetRho3(const trealdata &a){rho3 = a;}
        //! Set the thicknes in km
        void SetThicknesses(const trealdata &thick){thicknesses = thick;}
        const trealdata &GetThicknesses() {return thicknesses;}
        using MTStation::GetData;
        //! Calculate the synthetic data given the previously set parameters
        virtual void GetData(); 
        //! write model into file 
        gplib::rvec GetModelVector(); 
        void WriteModel(std::string filename); 
        void ReadModel(std::string filename);
        void WritePlot(std::string filename);
        virtual C1DAnisoMTSynthData *clone() const{return new C1DAnisoMTSynthData(*this);}
        C1DAnisoMTSynthData(const C1DAnisoMTSynthData &old);
        C1DAnisoMTSynthData();
        virtual ~C1DAnisoMTSynthData();
};

#endif /*C1DANISOMTSYNTHDATA_H_*/

---