00001 #ifndef C1DMTSYNTHDATA_H 00002 #define C1DMTSYNTHDATA_H 00003 00004 #include "MTStation.h" 00005 #include "types.h" 00006 #include <vector> 00007 #include <string> 00008 #include "VecMat.h" 00009 00010 class C1DMTSynthData : public MTStation 00011 { 00012 private: 00013 tcompdata Z; 00014 void Calc(); 00015 trealdata calc_frequencies; 00016 trealdata resistivity; // 00017 trealdata thickness; // Thickness for each layer 00018 trealdata reserrors; // the errors of the model parameters, only used for plotting 00019 trealdata thickerrors; 00020 public: 00021 const trealdata &GetResistivities() 00022 { 00023 return resistivity; 00024 } 00025 const trealdata &GetThicknesses() 00026 { 00027 return thickness; 00028 } 00029 void SetResistivities(const trealdata &res) 00030 { 00031 resistivity = res; 00032 } 00033 void SetResistivityErrors(const trealdata &re) 00034 { 00035 reserrors = re; 00036 } 00037 void SetThicknesses(const trealdata &thick) 00038 { 00039 thickness = thick; 00040 } 00041 void SetThicknessErrors(const trealdata &te) 00042 { 00043 thickerrors = te; 00044 } 00045 gplib::rvec GetModelVector(); 00046 void WriteModel(std::string filename); //write model into file for cagniard algorithm 00047 void ReadModel(std::string filename); 00048 void WritePlot(std::string filename); 00049 using MTStation::GetData; 00050 //! Calculate the synthetic data given the previously set parameters 00051 virtual void GetData(); 00052 virtual C1DMTSynthData *clone() const 00053 { 00054 return new C1DMTSynthData(*this); 00055 } 00056 C1DMTSynthData(const C1DMTSynthData &old); 00057 C1DMTSynthData(); 00058 virtual ~C1DMTSynthData(); 00059 00060 }; 00061 00062 #endif
1.5.5