CMTStation Class Reference

The class CMTStation is used to store the transfer functions and related information for a MT-site. More...

#include <CMTStation.h>

Inheritance diagram for CMTStation:

[legend]Collaboration diagram for CMTStation:

[legend]List of all members.

Public Types
enum	tmtdataformat {   unknown, mtt, j, edi,   pek }
Public Member Functions
void	AssignAll (const int nfreq)
void	Rotate (const double rotangle)
void	Rotate (void)
	Rotate to zero rotation angle.
trealdata	GetFrequencies () const
	return the available frequencies in a single vector
void	SetFrequencies (const trealdata &freqs)
	Set the frequencies of the tensor elements, invalidates the previously stored impedance data, mainly for creating synthetic data.
double	GetLatitude () const
	access funtion for Latitude
void	SetLatitude (double lat)
double	GetLongitude () const
void	SetLongitude (double lon)
double	GetElevation () const
std::string	GetName ()
double	GetAzimuth () const
const CMTTensor &	at (const unsigned int i) const
	direct acces to a tensor at a given index
const std::vector< CMTTensor > &	GetMTData () const
	Get the full vector of Tensor elements read only.
std::vector< CMTTensor > &	SetMTData ()
	Get the full vector of Tensor elements for reading and writing.
virtual CMTStation *	clone () const
	for parallel runs we need to make a copy of the object and its derived classes
	CMTStation ()
	CMTStation (const CMTStation &old)
	CMTStation (const int size)
	CMTStation (const std::string filename)
virtual	~CMTStation ()
virtual void	GetData (const std::string filename)
	read in data from file, determines format by ending
virtual void	GetData ()
virtual void	WriteData (const std::string filename)
void	WriteAsMtt (const std::string filename)
	Write data in goettingen .mtt format.
void	WriteAsEdi (const std::string filename)
	Write data as edi (no functionality yet).
void	WriteAsJ (const std::string filename)
	Write data to j-file.
void	WriteBack ()
	Write data back in original format, with filename given by station name.
CMTStation &	operator= (const CMTStation &source)
Protected Member Functions
void	Update ()
	Update all derived quantities.
void	Assign (const int nfreq)
	Assign() assigns zero to all derived quantities, this makes the calculation.
Friends
class	C1DMTSynthData

---

Detailed Description

The class CMTStation is used to store the transfer functions and related information for a MT-site.

Definition at line 13 of file CMTStation.h.

---

Member Enumeration Documentation

enum CMTStation::tmtdataformat

Enumerator:

unknown
mtt
j
edi
pek

Definition at line 15 of file CMTStation.h.

{unknown, mtt, j, edi, pek};

---

Constructor & Destructor Documentation

CMTStation::CMTStation

(

)

Definition at line 37 of file CMTStation.cpp.

Referenced by clone().

{
        InitialSetup();
}

CMTStation::CMTStation

(

const CMTStation &

old

)

Definition at line 614 of file CMTStation.cpp.

References Update().

                                           :
        MTData(old.MTData),
        TFData(old.TFData),
        latitude(old.latitude),
        longitude(old.longitude),
        elevation(old.elevation),
        azimuth(old.azimuth),
        name(old.name),
        dataformat(old.dataformat)
{
        Update();
}

Here is the call graph for this function:

CMTStation::CMTStation

(

const int

size

)

Definition at line 51 of file CMTStation.cpp.

References Assign().

                                    :
MTData(size)
{
        InitialSetup();
        Assign(size);
}

Here is the call graph for this function:

CMTStation::CMTStation

(

const std::string

filename

)

Definition at line 45 of file CMTStation.cpp.

References GetData().

{
        InitialSetup();
        GetData(filename);
}

Here is the call graph for this function:

CMTStation::~CMTStation

(

)

[virtual]

Definition at line 42 of file CMTStation.cpp.

{}

---

Member Function Documentation

void CMTStation::Update

(

)

[protected]

Update all derived quantities.

The Update() method first updates all components individually and then recalculates all derived quantities that depend on more than one component

Definition at line 71 of file CMTStation.cpp.

Referenced by CMTStation(), C1DMTSynthData::GetData(), C1DAnisoMTSynthData::GetData(), operator=(), Rotate(), and WriteAsMtt().

{
        
}

void CMTStation::Assign

(

const int

nfreq

)

[protected]

Assign() assigns zero to all derived quantities, this makes the calculation.

Definition at line 61 of file CMTStation.cpp.

References MTData.

Referenced by AssignAll(), CMTStation(), C1DMTSynthData::GetData(), and C1DAnisoMTSynthData::GetData().

{
        MTData.assign(nfreq,CMTTensor());
        TFData.assign(nfreq,CMagneticTF());
}

void CMTStation::AssignAll

(

const int

nfreq

)

Definition at line 57 of file CMTStation.cpp.

References Assign().

Referenced by main().

{
        Assign(nfreq);
}

Here is the call graph for this function:

void CMTStation::Rotate

(

const double

rotangle

)

Given a rotation angle in radian this method rotates the impedance tensor and updates all derived quantities

Definition at line 98 of file CMTStation.cpp.

References MTData, and Update().

{
    for (unsigned int i = 0; i < MTData.size() ; ++i)
    {
           DoRotate(i,rotangle);
    }
    Update();
}

Here is the call graph for this function:

void CMTStation::Rotate

(

void

)

Rotate to zero rotation angle.

Rotate(void) rotates the impedance tensor by the negative of the angle given in the rotangles field

Definition at line 109 of file CMTStation.cpp.

References at(), MTData, CMTTensor::rotangle, and Update().

Referenced by GetData().

{
    for (unsigned int i = 0; i < MTData.size() ; ++i)
    {
                DoRotate(i,-MTData.at(i).rotangle);
    }
    Update();
}

Here is the call graph for this function:

trealdata CMTStation::GetFrequencies

(

)

const

return the available frequencies in a single vector

Return a vector containing the available frequencies; This should be rewritten to be more efficient

Definition at line 119 of file CMTStation.cpp.

References CMTTensor::GetFrequency(), and MTData.

Referenced by C1DMTSynthData::GetData(), C1DAnisoMTSynthData::GetData(), main(), PTensor1DMTObjective::SafeParallel(), and C1DMTObjective::SafeParallel().

{
        trealdata temp(MTData.size());
        transform(MTData.begin(),MTData.end(),temp.begin(),mem_fun_ref(&CMTTensor::GetFrequency));
        return temp;
}

Here is the call graph for this function:

void CMTStation::SetFrequencies

(

const trealdata &

freqs

)

Set the frequencies of the tensor elements, invalidates the previously stored impedance data, mainly for creating synthetic data.

Definition at line 126 of file CMTStation.cpp.

References at(), CMTTensor::frequency, and MTData.

Referenced by main(), PTensor1DMTObjective::SafeParallel(), and C1DMTObjective::SafeParallel().

{
        //assert(MTData.size() == freqs.size());
        const unsigned int nfreqs = freqs.size();
        if (MTData.size() != nfreqs) //this invalidates the data stored there before
        {
                MTData.resize(nfreqs);
                TFData.resize(nfreqs);
        }
        for (unsigned int i = 0; i < nfreqs; ++i)
                MTData.at(i).frequency = freqs.at(i);
        assert(MTData.size() == TFData.size());
        
        
}

Here is the call graph for this function:

double CMTStation::GetLatitude

(

)

const [inline]

access funtion for Latitude

Definition at line 68 of file CMTStation.h.

Referenced by PrintComponent().

{return latitude;}

void CMTStation::SetLatitude

(

double

lat

)

[inline]

Definition at line 69 of file CMTStation.h.

Referenced by CMTStationList::GetData().

{latitude = lat;}

double CMTStation::GetLongitude

(

)

const [inline]

Definition at line 70 of file CMTStation.h.

Referenced by PrintComponent().

{return longitude;}

void CMTStation::SetLongitude

(

double

lon

)

[inline]

Definition at line 71 of file CMTStation.h.

Referenced by CMTStationList::GetData().

{longitude = lon;}

double CMTStation::GetElevation

(

)

const [inline]

Definition at line 72 of file CMTStation.h.

{return elevation;}

std::string CMTStation::GetName

(

)

[inline]

Definition at line 73 of file CMTStation.h.

Referenced by main(), HasSameName::operator()(), and PrintComponent().

{return name;}

double CMTStation::GetAzimuth

(

)

const [inline]

Definition at line 74 of file CMTStation.h.

{return azimuth;};

const CMTTensor& CMTStation::at

(

const unsigned int

i

)

const [inline]

direct acces to a tensor at a given index

Definition at line 76 of file CMTStation.h.

Referenced by C1DMTSynthData::GetData(), main(), PTensorMTStation::operator=(), Rotate(), PTensor1DMTObjective::SafeParallel(), SetFrequencies(), and WriteAsJ().

{return MTData.at(i);}

const std::vector<CMTTensor>& CMTStation::GetMTData

(

)

const [inline]

Get the full vector of Tensor elements read only.

Definition at line 78 of file CMTStation.h.

Referenced by main(), PTensorMTStation::operator=(), PrintComponent(), C1DMTObjective::SafeParallel(), and WriteParameterToFile().

{return MTData;}

std::vector<CMTTensor>& CMTStation::SetMTData

(

)

[inline]

Get the full vector of Tensor elements for reading and writing.

Definition at line 80 of file CMTStation.h.

Referenced by main().

{return MTData;}

virtual CMTStation* CMTStation::clone

(

)

const [inline, virtual]

for parallel runs we need to make a copy of the object and its derived classes

Reimplemented in C1DAnisoMTSynthData, and C1DMTSynthData.

Definition at line 82 of file CMTStation.h.

References CMTStation().

{return new CMTStation(*this);}

Here is the call graph for this function:

void CMTStation::GetData

(

const std::string

filename

)

[virtual]

read in data from file, determines format by ending

Implements CGeneralData.

Definition at line 419 of file CMTStation.cpp.

References edi, j, MTData, mtt, pek, Rotate(), and unknown.

Referenced by CMTStationList::GetData(), init(), and main().

{
        if (!boost::filesystem::exists(filename))
                throw CFatalException("File does not exist : "+filename);
        string ending;
        unsigned int dotpos = filename.find('.',0);
        if (dotpos != string::npos)
                ending = filename.substr(dotpos);
        else
        {
                cerr << "File has no extension ! " << filename << endl;
                dataformat = unknown;
                throw CFatalException("File has no extension ! " + filename);
        }
        if ( ending == ".mtt")
                                        {
                                                ReadMtt(filename);
                                                dataformat = mtt;
                                        }
                                        else
                                                if      (ending == ".dat" || ending == ".j")
                                                {
                                                        ReadJ(filename);
                                                        dataformat = j;
                                                }
                                                else
                                                        if      (filename.substr(filename.size()-4) == ".edi")
                                                        {
                                                                ReadEdi(filename);
                                                                dataformat = edi;
                                                        }
                                                        else
                                                                if (ending == ".pek")
                                                                {
                                                                        ReadPek1D(filename);
                                                                        dataformat = pek;
                                                                }
                                                                else
                                                                {
                                                                        dataformat = unknown;
                                                                        throw CFatalException("File not recognised: "+filename); 
                                                                }
        assert(MTData.size() == TFData.size());
        Rotate(); //rotate data to 0 degree 
                
}

Here is the call graph for this function:

virtual void CMTStation::GetData

(

)

[inline, virtual]

Declaration for GetData() without any parameter. This function is intended for use with objects that read in data that is generated within the program.

Implements CGeneralData.

Reimplemented in C1DAnisoMTSynthData, and C1DMTSynthData.

Definition at line 90 of file CMTStation.h.

Referenced by CMTStation().

{};

void CMTStation::WriteData

(

const std::string

filename

)

[virtual]

The abstract declaration for a method that writes the data to a file. Therefore no version without a filename is present

Implements CGeneralData.

Definition at line 466 of file CMTStation.cpp.

References EDIParser::name, and WriteBack().

Referenced by main().

{
        name = filename;
        WriteBack();
}

Here is the call graph for this function:

void CMTStation::WriteAsMtt

(

const std::string

filename

)

Write data in goettingen .mtt format.

Definition at line 491 of file CMTStation.cpp.

References Update().

Referenced by main(), WriteBack(), PTensor1DMTObjective::WriteData(), Iso1DMTObjective::WriteData(), and Aniso1DMTObjective::WriteData().

{
        Update();
    WriteMtt((filename+".mtt").c_str());
}

Here is the call graph for this function:

void CMTStation::WriteAsEdi

(

const std::string

filename

)

Write data as edi (no functionality yet).

Definition at line 497 of file CMTStation.cpp.

Referenced by WriteBack().

{
        throw CFatalException("WriteAsEdi not implemented yet !");
}

void CMTStation::WriteAsJ

(

const std::string

filename

)

Write data to j-file.

Definition at line 520 of file CMTStation.cpp.

References at(), EDIParser::azimuth, EDIParser::elevation, CMTTensor::frequency, CMTTensor::GetdZxx(), CMTTensor::GetdZxy(), CMTTensor::GetdZyx(), CMTTensor::GetdZyy(), CMTTensor::GetZxx(), CMTTensor::GetZxy(), CMTTensor::GetZyx(), CMTTensor::GetZyy(), EDIParser::latitude, EDIParser::longitude, MTData, and EDIParser::name.

Referenced by main(), and WriteBack().

{
        ofstream outfile;
        outfile.open((filename+".j").c_str());
        int actfreqs = 0;
        const double convfactor = 4. * acos(-1.)*1e-4;
        
        outfile << "# J-File Produced by CJData.cpp" << endl;
        outfile << ">LATITUDE  = " << latitude  << endl; 
        outfile << ">LONGITUDE = " << longitude << endl; 
        outfile << ">ELEVATION = " << elevation << endl;
        outfile << ">AZIMUTH   = " << azimuth << endl;
        outfile << name << endl;
        for (unsigned int i = 0; i < MTData.size(); ++i)
                if (MTData.at(i).frequency != 0)
                        actfreqs++;
        outfile << "ZXX S.I." << endl;
        outfile << actfreqs << endl;
        //WriteJBlock(DataXX,outfile,convfactor);
        WriteJBlock(&CMTTensor::GetZxx,&CMTTensor::GetdZxx,outfile,convfactor);
        outfile << "ZXY S.I." << endl;
        outfile << actfreqs << endl;
        WriteJBlock(&CMTTensor::GetZxy,&CMTTensor::GetdZxy,outfile,convfactor);
        
        outfile << "ZYX S.I." << endl;
        outfile << actfreqs << endl;
        WriteJBlock(&CMTTensor::GetZyx,&CMTTensor::GetdZyx,outfile,convfactor);
        
        outfile << "ZYY S.I." << endl;
        outfile << actfreqs << endl;
        WriteJBlock(&CMTTensor::GetZyy,&CMTTensor::GetdZyy,outfile,convfactor);
        
/*      if (norm(DataZY.Z.at(0)) !=0)
        {
                outfile << "TZX " << endl;
                outfile << actfreqs << endl;
                WriteJBlock(DataZX,outfile,1.);
        }
        
        if (norm(DataZY.Z.at(0)) !=0)
        {
                outfile << "TZY " << endl;
                outfile << actfreqs << endl;
                WriteJBlock(DataZY,outfile,1.);
        }*/
        outfile.close();
}

Here is the call graph for this function:

void CMTStation::WriteBack

(

)

Write data back in original format, with filename given by station name.

Definition at line 472 of file CMTStation.cpp.

References edi, j, mtt, EDIParser::name, WriteAsEdi(), WriteAsJ(), and WriteAsMtt().

Referenced by WriteData().

{
        switch (dataformat)
        {
                case mtt:
                        WriteAsMtt(name.c_str());
                        break;
                case j:
                        WriteAsJ(name.c_str());
                        break;
                case edi:
                        WriteAsEdi(name.c_str());
                        break;
                default:
                        throw CFatalException("Unknown File Format for Writing ! This should not happen");
        }
        
}

Here is the call graph for this function:

CMTStation & CMTStation::operator=

(

const CMTStation &

source

)

Definition at line 626 of file CMTStation.cpp.

References azimuth, EDIParser::azimuth, dataformat, elevation, EDIParser::elevation, latitude, EDIParser::latitude, longitude, EDIParser::longitude, MTData, name, EDIParser::name, TFData, and Update().

{
        if (this == &source) return *this;
        
        copy(source.MTData.begin(),source.MTData.end(),back_inserter(this->MTData));
        copy(source.TFData.begin(),source.TFData.end(),back_inserter(this->TFData));

        
        this->latitude = source.latitude;
        this->longitude = source.longitude;
        this->elevation = source.elevation;
        this->azimuth = source.azimuth;
        this->name = source.name;
        this->dataformat = source.dataformat;
        this->Update();
        return *this;
}

Here is the call graph for this function:

---

Friends And Related Function Documentation

friend class C1DMTSynthData [friend]

Definition at line 93 of file CMTStation.h.

---

The documentation for this class was generated from the following files:

• CMTStation.h
• CMTStation.cpp

---