This class stores and writes model for the respktn 1D seismic code that we use for receiver function calculations. More...

#include <ResPkModel.h>

Inheritance diagram for gplib::ResPkModel:

[legend]

Public Types
enum	WaveType { SWave, PWave }

Public Types inherited from gplib::SeismicModel
enum	tArrivalType { DirectS, DirectP }
	Do we want to calculate the arrival of a direct S-Wave or a P-wave. More...

Public Member Functions
void	SetInputWave (WaveType w)

double	GetSlowness () const
	Get the slowness in s/km for the synthetic forward calculation. More...

void	SetSlowness (const double s)
	Set the slowness in s/km for the synthetic forward calculation. More...

const trealdata &	GetStrike () const
	Read-only access to the vector of Strikes in each layer. More...

trealdata &	SetStrike ()
	Read-write access to the vector of Strikes in each layer. More...

const trealdata &	GetDip () const
	Read-only access to the vector of Dips in each layer. More...

trealdata &	SetDip ()
	Read-write access to the vector of Dips in each layer. More...

virtual void	ReadModel (const std::string filename)
	Read the model in its native format from a file. More...

virtual void	WriteModel (const std::string filename)
	Write the model in its native format to a file. More...

virtual void	WriteRunFile (const std::string &filename)
	Write out a script that performs a forward calculation for the current model. More...

	ResPkModel ()

	ResPkModel (const int nlayers)

	ResPkModel (const SeismicModel &source)

	ResPkModel (const ResPkModel &source)

ResPkModel &	operator= (const ResPkModel &source)

virtual	~ResPkModel ()

Public Member Functions inherited from gplib::SeismicModel
unsigned int	GetNpts () const
	Get the number of points for synthetic seismogram calculation. More...

void	SetNpts (const unsigned int s)
	Set the number of points for synthetic seismogram calculation. More...

double	GetSourceDepth () const
	Get the depth to the seismic source that generates the wavefield. More...

void	SetSourceDepth (const double s)
	Set the depth to the seismic source. More...

double	GetDt () const
	Get the time between two samples in s, this is for synthetic forward calculations. More...

void	SetDt (const double s)
	Set the time between two samples in s, this is for synthetic forward calculations. More...

const trealdata &	GetPVelocity () const
	Read-only access to the vector of P-velocities in km/s in each layer. More...

trealdata &	SetPVelocity ()
	Read-write access to the vector of P-velocities in km/s in each layer. More...

const trealdata &	GetSVelocity () const
	Read-only access to the vector of S-velocities in km/s in each layer. More...

trealdata &	SetSVelocity ()
	Read-write access to the vector of S-velocities in km/s in each layer. More...

const trealdata &	GetDensity () const
	Read-only access to the vector of densities in g/cm^3 in each layer. More...

trealdata &	SetDensity ()
	Read-write access to the vector of densities in g/cm^3 in each layer. More...

const trealdata &	GetThickness () const
	Read-only access to the vector of thicknesses in km in each layer. More...

trealdata &	SetThickness ()
	Read-write access to the vector of thicknesses in km in each layer. More...

const trealdata &	GetQp () const
	Read-only access to the vector of P quality factors for each layer. More...

trealdata &	SetQp ()
	Read-write access to the vector of P quality factors for each layer. More...

const trealdata &	GetQs () const
	Read-only access to the vector of S quality factors for each layer. More...

trealdata &	SetQs ()
	Read-write access to the vector of S quality factors for each layer. More...

void	SetSVelErrors (const trealdata &sve)
	Set error bars on S-velocities for plotting. More...

void	SetThickErrors (const trealdata &te)
	Set error bars on Thicknesses for plotting. More...

int	FindLayer (const double depth)
	Returns the index of the layer that correponds to depth in km. More...

double	MatchSlowness (const double slowness, const tArrivalType mode)
	Given a slowness in s/km and a wave type calculate the distance that matches this slowness. More...

double	CalcArrival (const tArrivalType mode, const double recdist)

double	CalcTravelTime (const tArrivalType mode, const double sdepth, const double rdepth, const double p)

void	WritePlot (const std::string &filename)
	Write out an ascii file for plotting with xmgrace etc. More...

void	PlotVelWithErrors (const std::string &filename)
	Write out an ascii file with error bars for plotting with xmgrace etc. More...

void	Init (const int nlayers)
	Init provides a convenient way to allocate memory in all structures for a given number of layers. More...

	SeismicModel (const int nlayers=0)

	SeismicModel (const SeismicModel &source)

SeismicModel &	operator= (const SeismicModel &source)

virtual	~SeismicModel ()

Detailed Description

This class stores and writes model for the respktn 1D seismic code that we use for receiver function calculations.

Our joint inversion uses respktn from C. Ammon to calculate synthetics for receiver function inversion. This class stores the description of this model and provides methods to read such models, write such models and generate scripts that calculate synthetic for such a model.

Definition at line 18 of file ResPkModel.h.