You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
113 lines
3.7 KiB
C
113 lines
3.7 KiB
C
1 year ago
|
/***************************************************************************
|
||
|
|
* Copyright (C) gempa GmbH *
|
||
|
|
* All rights reserved. *
|
||
|
|
* Contact: gempa GmbH (seiscomp-dev@gempa.de) *
|
||
|
|
* *
|
||
|
|
* GNU Affero General Public License Usage *
|
||
|
|
* This file may be used under the terms of the GNU Affero *
|
||
|
|
* Public License version 3.0 as published by the Free Software Foundation *
|
||
|
|
* and appearing in the file LICENSE included in the packaging of this *
|
||
|
|
* file. Please review the following information to ensure the GNU Affero *
|
||
|
|
* Public License version 3.0 requirements will be met: *
|
||
|
|
* https://www.gnu.org/licenses/agpl-3.0.html. *
|
||
|
|
* *
|
||
|
|
* Other Usage *
|
||
|
|
* Alternatively, this file may be used in accordance with the terms and *
|
||
|
|
* conditions contained in a signed written agreement between you and *
|
||
|
|
* gempa GmbH. *
|
||
|
|
***************************************************************************/
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
#ifndef SEISCOMP_SEISMOLOGY_MAGNITUDES_H
|
||
|
|
#define SEISCOMP_SEISMOLOGY_MAGNITUDES_H
|
||
|
|
|
||
|
|
#include <seiscomp/core.h>
|
||
|
|
|
||
|
|
namespace Seiscomp {
|
||
|
|
namespace Magnitudes {
|
||
|
|
|
||
|
|
/**
|
||
|
|
* Compute Mwp after Tsuboi(1999) with linear correction from GFZ
|
||
|
|
* @param amplitude maximum reading of absolute integrated displacement in m*s
|
||
|
|
* @param delta epicentral distance in degrees
|
||
|
|
* @param mag return value
|
||
|
|
* @param offset taken from fit against Harvard Mw (Mwp=Mw*slope+offset)
|
||
|
|
* @param slope taken from fit against Harvard Mw
|
||
|
|
* @param alpha P-wave velocity along ray path in m/s
|
||
|
|
* @param rho mass density along ray path in kg/m^3
|
||
|
|
* @param fp radiation pattern influence factor
|
||
|
|
* @return false Mwp undefined for given epicentral distance
|
||
|
|
* @return true Mwp calculated successfully
|
||
|
|
*/
|
||
|
|
SC_SYSTEM_CORE_API
|
||
|
|
bool compute_Mwp(
|
||
|
|
double amplitude,
|
||
|
|
double delta,
|
||
|
|
double &mag, // resulting magnitude
|
||
|
|
double offset = 0, double slope = 1,
|
||
|
|
// double offset=2.25,
|
||
|
|
// double slope=1./1.6,
|
||
|
|
double alpha=7900.,
|
||
|
|
double rho=3400.,
|
||
|
|
double fp=0.52);
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
/**
|
||
|
|
* Compute Mjma after Katsumata (1988) PRELIMINARY!!! only valid for depth < 60 km !!!
|
||
|
|
* @param amplitude maximum reading of absolute integrated displacement in m*s
|
||
|
|
* @param delta epicentral distance in degrees
|
||
|
|
* @return false Mjma undefined for given epicentral distance
|
||
|
|
* @return true Mwp calculated successfully
|
||
|
|
*/
|
||
|
|
SC_SYSTEM_CORE_API
|
||
|
|
bool compute_Mjma(
|
||
|
|
double amplitude,
|
||
|
|
double delta, // in degrees
|
||
|
|
double depth, // in kilometers
|
||
|
|
double& Mjma);
|
||
|
|
|
||
|
|
/**
|
||
|
|
* Mjma variant for BMG/Indonesia
|
||
|
|
*/
|
||
|
|
SC_SYSTEM_CORE_API
|
||
|
|
bool compute_Mbmg(
|
||
|
|
double amplitude,
|
||
|
|
double delta, // in degrees
|
||
|
|
double depth, // in kilometers
|
||
|
|
double& Mbmg);
|
||
|
|
|
||
|
|
|
||
|
|
SC_SYSTEM_CORE_API
|
||
|
|
bool
|
||
|
|
compute_mb(
|
||
|
|
double amplitude, // in micrometers
|
||
|
|
double period, // in seconds
|
||
|
|
double delta, // in degrees
|
||
|
|
double depth, // in kilometers
|
||
|
|
double *mag); // resulting magnitude
|
||
|
|
|
||
|
|
|
||
|
|
SC_SYSTEM_CORE_API
|
||
|
|
bool
|
||
|
|
compute_mb_fromVelocity(
|
||
|
|
double amplitude, // in micrometers/second
|
||
|
|
double delta, // in degrees
|
||
|
|
double depth, // in kilometers
|
||
|
|
double *mag); // resulting magnitude
|
||
|
|
|
||
|
|
|
||
|
|
SC_SYSTEM_CORE_API
|
||
|
|
bool
|
||
|
|
compute_ML(
|
||
|
|
double amplitude, // in micrometers
|
||
|
|
double delta, // in degrees
|
||
|
|
double depth, // in kilometers
|
||
|
|
double *mag); // resulting magnitude
|
||
|
|
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
#endif
|