A profile through the Earth model that stores interface radius, and interval P or S wave velocity information about the crustal stack below either the source or the receiver.
More...
|
| CrustalProfile () |
| Default constructor. More...
|
|
| CrustalProfile (const CrustalProfile &CrustalProfile) |
| Copy constructor. More...
|
|
virtual | ~CrustalProfile () |
|
CrustalProfile & | operator= (const CrustalProfile &other) |
| Equal operator. More...
|
|
bool | operator== (const CrustalProfile &other) |
| Equality operator. More...
|
|
bool | operator!= (const CrustalProfile &other) |
| Inequality operator. More...
|
|
void | setup (Grid &grid, const int &phase, const double &lat, const double &lon, const double &depth) |
| Parameterized constructor. More...
|
|
const Location & | getLocation () const |
| Retrieve a const reference to the Location associated with this Profile. More...
|
|
double | getEarthRadius () |
|
int | getNIntervals () |
| Retrieve the number of intervals associated with this Profile. More...
|
|
void | getData (vector< double > &depths, vector< double > &velocities) |
| Retrieve all the interval depth and velocity information contained in this CrustalProfile object. More...
|
|
double | getInterfaceRadius (const int &k) |
| Retrieve the radius of the top of the k'th interval, in km. More...
|
|
void | getDepths (double *depths) |
| Retrieve the depths of the top of each intervals, in km. More...
|
|
double * | getVelocities () |
| Retrieve the P or S wave velocity of each interval, in km/sec. More...
|
|
double | getDepth (const int &k) |
| Retrieve the depth of the top of the k'th interval, in km. More...
|
|
double | getVelocity (const int &k) |
| Retrieve the P or S wave velocity of the k'th interval, in km/sec. More...
|
|
double | getPCrit (GreatCircle *greatCircle) |
| Retrieve the critical spherical ray parameter for this profile, in radians/seconds. More...
|
|
void | xtCrust (GreatCircle *greatCircle, const double &rayParameter, double &xTotal, double &tTotal) |
| Compute the horizontal offset (radians) and the travel time (sec) of the ray through the crust. More...
|
|
void | xtCrust (GreatCircle *greatCircle, const double &rayParameter, int layid[], double x[], double r[], double v[], double t[], int &npoints) |
| Compute the horizontal offset (radians) and the travel time (sec) of the ray through the crust. More...
|
|
const int * | getLayid () |
|
int | getNLayid () |
|
int | getTopLayid () |
|
string | toString (GreatCircle *greatCircle, double rayParameter) |
| Returns a formatted string containing detailed information about this Profile. More...
|
|
string | toString () |
|
size_t | memSize () |
|
void | setTauPSite (taup::TauPSite *tps) |
| Set the input TauPSite into this CrustalProfile;. More...
|
|
taup::TauPSite * | getTauPSite () |
| Return this CrustalProfiles TauPSite object;. More...
|
|
bool | isInCrust () |
| true location radius > moho radius. More...
|
|
| InterpolatedProfile () |
| Default constructor. More...
|
|
| InterpolatedProfile (Grid &grid, Location &location) |
| Parameterized constructor. More...
|
|
| InterpolatedProfile (const InterpolatedProfile &other) |
| Copy constructor. More...
|
|
virtual | ~InterpolatedProfile () |
| Destructor. More...
|
|
InterpolatedProfile & | operator= (const InterpolatedProfile &other) |
| Equal operator. More...
|
|
bool | operator== (const InterpolatedProfile &other) |
| Equality operator. More...
|
|
bool | operator!= (const InterpolatedProfile &other) |
| Inequality operator. More...
|
|
int | getNCoefficients () |
|
vector< GridProfile * > & | getNodes () |
| Retrieve a list of pointers to the GridProfile objects upon which this InterpolatedProfile is dependent. More...
|
|
GridProfile * | getNode (const int &i) |
|
vector< int > & | getNodeIds () |
|
int | getNodeId (const int &i) |
|
vector< double > & | getCoefficients () |
| Retrieve the interpolation coefficients that define the dependency of this InterpolatedProfile on its neighbors. More...
|
|
double | getCoefficient (const int &i) |
|
void | getNodeIds (int *nodeIds, int &size) |
| Retrieve the ID numbers of the GridProfiles that contributed to the interpolated values at this InterpolatedProfile object. More...
|
|
void | getCoefficients (double *coeff, int &size) |
| Retrieve the interpolation coefficients that were applied to the neighboring GridProfiles that contributed to the interpolated values at this InterpolatedProfile object. More...
|
|
void | getWeights (int *nodeIds, double *coeff, int &size) |
|
void | interpRadius (const int &k, double &radius) |
| Calculate a single radius value based on the neighboring GridProfile objects and the interpolation coefficients. More...
|
|
void | interpDepth (const int &k, double &depth) |
| Calculate a single depth value based on the neighboring GridProfile objects and the interpolation coefficients. More...
|
|
void | interpVelocity (const int &type, const int &k, double &velocity) |
| Calculate a single velocity value based on the neighboring GridProfile objects and the interpolation coefficients. More...
|
|
void | interpGradient (const int &type, double &gradient) |
| Calculate a single gradient value based on the neighboring GridProfile objects and the interpolation coefficients. More...
|
|
bool | isActiveProfile () |
| Returns true if all of the neighboring GridProfile objects. More...
|
|
The CrustalProfile class represents a Profile based on radius and P or S wave velocity values interpolated from values of nearby GridProfile objects. A CrustalProfile object should be acquired by calling Grid::getCrustalProfile().
CrustalProfile owns a Location object which defines a point in the Earth. For CrustalProfile objects that represent a seismic source, the Location object represents the latitude, longitude and depth of the event. For CrustalProfile objects that represent a seismic receiver, the Location object represents the latitude, longitude and depth (negative elevation) of the receiver.
Each instance of a CrustalProfile is populated with only the information needed to compute crustal travel time for a particular phase. Supported phases are Pn, Sn, Pg and Lg. Since mantle velocity gradient information is not required to compute crustal travel times, CrustalProfile objects do not store mantle gradient information. CrustalProfile objects to support Pn or Pg store only P wave velocities while those instantiated to support Sn or Lg only store S wave velocity information.
The xtCrust() function computes the travel time (in seconds) and the horizontal offset (in radians) of the ray as it traverses the crust from the depth of the Location object down to the head wave interface (Moho or top of middle crust). xtCrust takes a ray parameter as an arguement. The function will always return an offset and travel time because if the ray parameter exceeds the critical ray parameter for any of the velocity interfaces that the ray has to interact with, the function will decrease the ray parameter enough to successfully travel from the Location depth to the deepest interface it knows about, using a constant ray parameter over the distance traveled.
CrustalProfile objects can be members of multiple GreatCircle objects. Grid maintains a map that associates a phase/Location combination with a particular instance of a CruatalProfile. The map can be cleared by calling Grid::clearCrustalProfiles(), which is also accessible from SlbmInterface::clear().
void slbm::CrustalProfile::xtCrust |
( |
GreatCircle * |
greatCircle, |
|
|
const double & |
rayParameter, |
|
|
int |
layid[], |
|
|
double |
x[], |
|
|
double |
r[], |
|
|
double |
v[], |
|
|
double |
t[], |
|
|
int & |
npoints |
|
) |
| |
Computes the travel time (in seconds) and the horizontal offset (in radians) of the ray as it traverses the crust from the depth of the Location object down to the head wave interface (Moho for Pn and Sn, top of middle crust for Pg and Lg).
Unlike the other version of xtCrust(), this version returns information about the horizontal offset and travel time at each interface of the crust that the ray interacts with.
The horizontal offset is the angular separation of the Location object owned by this CrustalProfile and the pierce point of the ray on the head wave interface.
The function will always return offset and travel time information because if the ray parameter exceeds the critical ray parameter for any of the velocity interfaces that the ray has to interact with, the function will decrease the ray parameter enough to successfully travel from the Location depth to the deepest interface it knows about, using a constant ray parameter over the distance traveled.
- Parameters
-
greatCircle | the GreatCircle object, of which this CrustalProfile is a member, that is requesting the information. |
rayParameter | the spherical ray parameter, in radians/sec. |
layid | the layerid number of each interface that the ray interacted with. Only includes layers below the depth of the Location object and excludes the layerids of zero thickness layers. |
x | the horizontal offset of the point where the ray intersected each interface, in radians. The first element is zero, and the last element is the same as the value of xTotal returned by the other version of xtCrust(). |
r | the radii of the interfaces that the ray interacted with, in km. |
v | the velocities of the layers that the ray interacted with, in km/sec. |
t | the cumulative travel time of the ray across all the preceeding layers. The first element is zero, and the last element is the same as the value of tTotal returned by the other version of xtCrust(). |
npoints | the number of elements in arrays layid, x, r, v and t. |