From 3e08c026c51a585508adc3d4a05c373651ec8787 Mon Sep 17 00:00:00 2001
From: Jason Altekruse <jaltekruse@contractor.usgs.gov>
Date: Tue, 12 Jan 2021 15:59:44 -0700
Subject: [PATCH] initial implementation of PSHAB20. calc mean needs to be
fixed, uncertainty not yet implemented. tests do not pass.
---
.../gov/usgs/earthquake/nshmp/gmm/Gmm.java | 84 +++
.../earthquake/nshmp/gmm/ParkerEtAl_2020.java | 496 +++++++++++++++---
.../gmm/coeffs/PSHAB20_interface.csv | 27 +
.../resources/gmm/coeffs/PSHAB20_slab.csv | 27 +
.../gmm/coeffs/ParkerEtAl_2020_interface.csv | 27 -
.../gmm/coeffs/ParkerEtAl_2020_slab.csv | 27 -
.../usgs/earthquake/nshmp/gmm/PSHAB20.java | 42 ++
src/test/resources/gmm/PSHAB20_inputs.csv | 2 +
src/test/resources/gmm/PSHAB20_results.csv | 23 +
9 files changed, 625 insertions(+), 130 deletions(-)
create mode 100644 src/main/resources/gmm/coeffs/PSHAB20_interface.csv
create mode 100644 src/main/resources/gmm/coeffs/PSHAB20_slab.csv
delete mode 100644 src/main/resources/gmm/coeffs/ParkerEtAl_2020_interface.csv
delete mode 100644 src/main/resources/gmm/coeffs/ParkerEtAl_2020_slab.csv
create mode 100644 src/test/java/gov/usgs/earthquake/nshmp/gmm/PSHAB20.java
create mode 100644 src/test/resources/gmm/PSHAB20_inputs.csv
create mode 100644 src/test/resources/gmm/PSHAB20_results.csv
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java
index 74fc8ae0..21abfa7b 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java
@@ -288,6 +288,90 @@ public enum Gmm {
NgaSubductionUsgs_2018.COEFFS,
NgaSubductionUsgs_2018.CONSTRAINTS),
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_INTERFACE(
+ ParkerEtAl_2020.Interface.class,
+ ParkerEtAl_2020.Interface.NAME,
+ ParkerEtAl_2020.COEFFS_INTERFACE,
+ ParkerEtAl_2020.CONSTRAINTS_INTERFACE),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_SLAB(
+ ParkerEtAl_2020.Slab.class,
+ ParkerEtAl_2020.Slab.NAME,
+ ParkerEtAl_2020.COEFFS_SLAB,
+ ParkerEtAl_2020.CONSTRAINTS_SLAB),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_CASCADIA_INTERFACE(
+ ParkerEtAl_2020.CascadiaInterface.class,
+ ParkerEtAl_2020.CascadiaInterface.NAME,
+ ParkerEtAl_2020.COEFFS_INTERFACE,
+ ParkerEtAl_2020.CONSTRAINTS_INTERFACE),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_CASCADIA_SLAB(
+ ParkerEtAl_2020.CascadiaSlab.class,
+ ParkerEtAl_2020.CascadiaSlab.NAME,
+ ParkerEtAl_2020.COEFFS_SLAB,
+ ParkerEtAl_2020.CONSTRAINTS_SLAB),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_CASCADIA_SEATTLE_INTERFACE(
+ ParkerEtAl_2020.CascadiaSeattleInterface.class,
+ ParkerEtAl_2020.CascadiaSeattleInterface.NAME,
+ ParkerEtAl_2020.COEFFS_INTERFACE,
+ ParkerEtAl_2020.CONSTRAINTS_INTERFACE),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_CASCADIA_SEATTLE_SLAB(
+ ParkerEtAl_2020.CascadiaSeattleSlab.class,
+ ParkerEtAl_2020.CascadiaSeattleSlab.NAME,
+ ParkerEtAl_2020.COEFFS_SLAB,
+ ParkerEtAl_2020.CONSTRAINTS_SLAB),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_CASCADIA_PACNW_INTERFACE(
+ ParkerEtAl_2020.CascadiaPacNWInterface.class,
+ ParkerEtAl_2020.CascadiaPacNWInterface.NAME,
+ ParkerEtAl_2020.COEFFS_INTERFACE,
+ ParkerEtAl_2020.CONSTRAINTS_INTERFACE),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_CASCADIA_PACNW_SLAB(
+ ParkerEtAl_2020.CascadiaPacNWSlab.class,
+ ParkerEtAl_2020.CascadiaPacNWSlab.NAME,
+ ParkerEtAl_2020.COEFFS_SLAB,
+ ParkerEtAl_2020.CONSTRAINTS_SLAB),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_ALASKA_INTERFACE(
+ ParkerEtAl_2020.AlaskaInterface.class,
+ ParkerEtAl_2020.AlaskaInterface.NAME,
+ ParkerEtAl_2020.COEFFS_INTERFACE,
+ ParkerEtAl_2020.CONSTRAINTS_INTERFACE),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_ALASKA_SLAB(
+ ParkerEtAl_2020.AlaskaSlab.class,
+ ParkerEtAl_2020.AlaskaSlab.NAME,
+ ParkerEtAl_2020.COEFFS_SLAB,
+ ParkerEtAl_2020.CONSTRAINTS_SLAB),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_ALEUTIANS_INTERFACE(
+ ParkerEtAl_2020.AleutiansInterface.class,
+ ParkerEtAl_2020.AleutiansInterface.NAME,
+ ParkerEtAl_2020.COEFFS_INTERFACE,
+ ParkerEtAl_2020.CONSTRAINTS_INTERFACE),
+
+ /** @see ParkerEtAl_2020 NGA-Subduction */
+ PSHAB20_ALEUTIANS_SLAB(
+ ParkerEtAl_2020.AleutiansSlab.class,
+ ParkerEtAl_2020.AleutiansSlab.NAME,
+ ParkerEtAl_2020.COEFFS_SLAB,
+ ParkerEtAl_2020.CONSTRAINTS_SLAB),
+
/** @see BcHydro_2012 */
BCHYDRO_12_INTERFACE(
BcHydro_2012.Interface.class,
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java
index b16c1c78..b00ea58c 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/ParkerEtAl_2020.java
@@ -4,10 +4,13 @@ import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.MW;
import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.RRUP;
import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.VS30;
import static gov.usgs.earthquake.nshmp.gmm.GmmInput.Field.ZHYP;
+import static gov.usgs.earthquake.nshmp.gmm.Imt.PGA;
import static java.lang.Math.exp;
import static java.lang.Math.log;
+import static java.lang.Math.log10;
import static java.lang.Math.min;
import static java.lang.Math.pow;
+import static java.lang.Math.sqrt;
import java.util.Map;
@@ -17,7 +20,14 @@ import gov.usgs.earthquake.nshmp.Maths;
import gov.usgs.earthquake.nshmp.gmm.GmmInput.Constraints;
/**
- * Implementation of the PEER NGA-Subduction GMM, Parker et al., 2020 *
+ * Implementation of the PEER NGA-Subduction GMM, Parker et al., 2020 GMM
+ *
+ * Notes:
+ *
+ * This model is only applicable to sites in the forearc, future residuals
+ * analysis is planned to determine if there are attenuation differences between
+ * forearc and backarc zones.
+ *
* <p><b>Reference:</b> Parker, G.A, Stewart, J.P., Boore, D.M., Atkinson, G.M.,
* Hassani, B. (2020) NGA-Subduction global ground-motion models with regional
* adjustment factors. Report no. 2020/03. Berkeley, CA: PEER, 131 p.
@@ -27,6 +37,14 @@ import gov.usgs.earthquake.nshmp.gmm.GmmInput.Constraints;
* @author U.S. Geological Survey
*/
+/*
+ * Regions: global, Alaska, Cascadia. Unused Regions: CAM (Central America),
+ * Japan, SA (South America), Taiwan
+ *
+ * Saturation Regions: global, Aleutian, Alaska, Cascadia. Unused Regions:
+ * Central America S, Central America N, Japan Pacific, Japan Phillipines, South
+ * America N, South America S, Taiwan W, Taiwan E
+ */
public abstract class ParkerEtAl_2020 implements GroundMotionModel {
static final String NAME = "Parker et al. (2020) NGA-Subduction";
@@ -37,58 +55,173 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
.set(ZHYP, Range.closed(0.0, 40.0))
.set(VS30, Range.closed(150.0, 2000.0))
.build();
- static final Constraints CONSTRAINTS_INTRASLAB = Constraints.builder()
+ static final Constraints CONSTRAINTS_SLAB = Constraints.builder()
.set(MW, Range.closed(4.5, 8.5))
.set(RRUP, Range.closed(35.0, 1000.0))
.set(ZHYP, Range.closed(0.0, 200.0))
.set(VS30, Range.closed(150.0, 2000.0))
.build();
- // different for interface and slab?
- // COEFFS_COMMON for shared coefs
- // abstract CoefficientContainer COEFFS for interface/slab-specific coefs set?
- static final CoefficientContainer COEFFS = new CoefficientContainer("TBD.csv");
+ static final class Region {
+ static final String GLOBAL = "Global";
+ static final String ALASKA = "Alaska";
+ static final String ALEUTIANS = "Aleutians";
+ static final String CASCADIA = "Cascadia";
+ }
+
+ static final class Basin {
+ static final String NONE = "None";
+ static final String PACNW_NONE = "PacNW-no-basin";
+ static final String PACNW = "PacNW";
+ static final String SEATTLE = "Seattle";
+ }
+
+ // aleatory uncertainty coefficients are duplicated in interface and interface
+ // CSV files (as are some other coefficients that are also independent of
+ // event type)
+ private static final String INTERFACE_CSV = "PSHAB20_interface.csv";
+ private static final String SLAB_CSV = "PSHAB20_slab.csv";
+ // 4 epistemic coefficients tied to region and interface/slab
+ private static final String EPISTEMIC_INTERFACE_CSV = "PSHAB20_epistemic_interface.csv";
+ private static final String EPISTEMIC_SLAB_CSV = "PSHAB20_epistemic_slab.csv";
+
+ static final CoefficientContainer COEFFS_INTERFACE;
+ static final CoefficientContainer COEFFS_SLAB;
+
+ static {
+ COEFFS_INTERFACE = new CoefficientContainer(INTERFACE_CSV);
+ COEFFS_SLAB = new CoefficientContainer(SLAB_CSV);
+ }
- private static final double B4 = 0.1;
- private static final double V1 = 270.0; // m/s
- private static final double VREF = 760.0; // m/s
+ // private static final double B4 = 0.1;
+ // private static final double V1 = 270.0; // m/s
+ // private static final double VREF = 760.0; // m/s
private static final double DB1 = 20.0; // km
private static final double DB2 = 67.0; // km
private static final class Coefficients {
final Imt imt;
- final double c0, c1, b4, a0, c4, c5, c6, d, m, db, v2, s1, s2, f4, f5, e1, e2, e3;
-
- Coefficients(Imt imt, CoefficientContainer cc) {
+ final double c0; // model constant
+ final double c1, b4, a0; // path scaling
+ final double mc; // saturation corner magnitude
+ final double c4, c5, c6; // magnitude scaling
+ final double d, m, db; // source depth scaling (slab only)
+ final double v1, v2, vRef, s1, s2; // linear site amplification
+ final double f4, f5; // non-linear site amplification
+ final double c_e1, c_e2, c_e3, del_None, del_Seattle; // Cascadia basin
+ final double Ï„, Ï•21, Ï•22, Ï•2v, vm, Ï•2s2s0, a1, Ï•2ss1, Ï•2ss2, a2; // alleatory
+ // uncertainty
+ final double ðœŽEp1, ðœŽEp2, t1, t2;
+
+ Coefficients(Imt imt, CoefficientContainer cc, boolean isSlab, String region) {
this.imt = imt;
Map<String, Double> coeffs = cc.get(imt);
- c0 = 0.0;
- c1 = 0.0;
- b4 = 0.0; // period-independent
- a0 = 0.0;
- c4 = 0.0;
- c5 = 0.0;
- c6 = 0.0;
- d = 0.0;
- m = 0.0;
- db = 0.0;
- v2 = 0.0;
- s1 = 0.0; // s1 = s2 except for Japan and Taiwan
- s2 = 0.0;
- f4 = 0.0;
- f5 = 0.0;
- e1 = 0.0;
- e2 = 0.0;
- e3 = 0.0;
+ // Model constant
+ c0 = coeffs.get(region + "_c0");
+ // Path
+ c1 = coeffs.get("c1");
+ b4 = coeffs.get("b4"); // period-independent
+ a0 = coeffs.get(getA0Col(isSlab, region));
+ // a0 = region.equals(Region.ALASKA) || (region.equals(Region.CASCADIA) &&
+ // isSlab)
+ // ? coeffs.get(region + "_a0") : coeffs.get(Region.GLOBAL + "_a0");
+ // Magnitude saturation corner magnitude
+ mc = getSaturationMagnitude(region, isSlab);
+ // Magnitude-scaling
+ c4 = coeffs.get("c4");
+ c5 = coeffs.get("c5");
+ c6 = coeffs.get("c6");
+ // Source depth scaling (slab only)
+ d = isSlab ? coeffs.get("d") : 0.0;
+ m = isSlab ? coeffs.get("m") : 0.0;
+ db = isSlab ? coeffs.get("db") : 0.0;
+ // Linear site amplification
+ v1 = coeffs.get("V1");
+ v2 = coeffs.get("V2");;
+ vRef = coeffs.get("Vref");
+ s2 = (region.equals(Region.ALEUTIANS)) ? coeffs.get(Region.GLOBAL + "_s2")
+ : coeffs.get(region + "_s2");
+ s1 = s2; // s1 = s2 except for Japan and Taiwan
+ // Non-linear site amplification
+ f4 = coeffs.get("f4");
+ f5 = coeffs.get("f5");
+ // Cascadia basin
+ c_e1 = coeffs.get("C_e1");
+ c_e2 = coeffs.get("C_e2");
+ c_e3 = coeffs.get("C_e3");
+ del_None = coeffs.get("del_None");
+ del_Seattle = coeffs.get("del_Seattle");
+ // aleatory model
+ Ï„ = coeffs.get("Tau");
+ Ï•21 = coeffs.get("phi21");
+ Ï•22 = coeffs.get("phi22");
+ Ï•2v = coeffs.get("phi2V");
+ vm = coeffs.get("VM");
+ Ï•2s2s0 = coeffs.get("phi2S2S0");
+ a1 = coeffs.get("a1");
+ Ï•2ss1 = coeffs.get("phi2SS1");
+ Ï•2ss2 = coeffs.get("phi2SS2");
+ a2 = coeffs.get("a2");
+ // epistemic model
+ ðœŽEp1 = 0.4; // GLOBAL INTERFACE
+ ðœŽEp2 = 0.4; // GLOBAL INTERFACE
+ t1 = 0.2; // GLOBAL INTERFACE
+ t2 = 0.2; // GLOBAL INTERFACE
}
+
+ private static final String getA0Col(boolean isSlab, String region) {
+ String a0Col;
+ if (region.equals(Region.ALASKA)) {
+ a0Col = Region.ALASKA + "_a0";
+ } else if (region.equals(Region.CASCADIA) && isSlab) {
+ a0Col = Region.CASCADIA + "_a0";
+ } else {
+ a0Col = Region.GLOBAL + "_a0";
+ }
+ return a0Col;
+ }
+
+ private static final double getSaturationMagnitude(String region, boolean isSlab) {
+ // TODO? Read from coeff table?
+ // Parket et al. (2020) Table 4.1
+ // Regional saturation magnitudes for interface events computed using
+ // seismogenic fault width [Campbell 2020] and for intraslab events
+ // computed
+ // using slab thickness [Ji and Archuleta 2018].
+ switch (region) {
+ case Region.GLOBAL:
+ return isSlab ? 7.6 : 7.9;
+ case Region.ALASKA:
+ return isSlab ? 7.20 : 8.6;
+ case Region.ALEUTIANS:
+ return isSlab ? 7.98 : 8.0;
+ case Region.CASCADIA:
+ return isSlab ? 7.20 : 7.7;
+ default:
+ throw new IllegalArgumentException("Region [" + region + "] not supported");
+ }
+ }
+
+ // set/get epistemic uncertainty parameters
+ // Region,SigEp1,SigEp2,T1,T2
+ // Global,0.4,0.4,0.2,0.4
+ // Alaska,0.15,0.1,1,4
+ // Aleutian,0.15,0.1,1,4
+ // Cascadia,0.43,0.33,0.2,0.5
+
}
private final Coefficients coeffs;
private final Coefficients coeffsPGA;
ParkerEtAl_2020(final Imt imt) {
- coeffs = new Coefficients(imt, COEFFS);
- coeffsPGA = new Coefficients(imt, COEFFS);
+ coeffs = initCoeffs(imt, isSlab(), getRegion());
+ coeffsPGA = initCoeffs(PGA, isSlab(), getRegion());
+ }
+
+ private static Coefficients initCoeffs(final Imt imt, final boolean isSlab, final String region) {
+ CoefficientContainer coeffs = isSlab ? COEFFS_SLAB : COEFFS_INTERFACE;
+ return new Coefficients(imt, coeffs, isSlab, region);
}
@Override
@@ -97,11 +230,10 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
* Equation 4.1 μ_lny = c0 + Fp + Fm + Fd + Fs
*/
- double pgaRef =
- exp(calcMean(coeffsPGA, isSlab(), false, getMc(), 0.0, in.Mw, in.rRup, in.zHyp, in.vs30,
- in.z2p5));
+ // double Mc = getSaturationMagnitude(getRegion(), isSlab());
+ double pgaRef = exp(calcMean(coeffsPGA, isSlab(), in.Mw, in.rRup, in.zHyp));
- double μ = calcMean(coeffs, isSlab(), useBasinTerm(), getMc(), pgaRef, in.Mw, in.rRup, in.zHyp,
+ double μ = calcMean(coeffs, isSlab(), getBasin(), in.Mw, in.rRup, in.zHyp, pgaRef,
in.vs30, in.z2p5);
double σ = 0.0;
@@ -110,37 +242,49 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
abstract boolean isSlab();
- // Note: mc, corner magnitude, is a regional saturation magnitude
- abstract double getMc();
+ abstract String getRegion();
- abstract boolean useBasinTerm();
+ abstract String getBasin();
+
+ /* calc median at 760 reference conditions */
+ private static final double calcMean(
+ final Coefficients c,
+ final boolean isSlab,
+ final double Mw,
+ final double rRup,
+ final double zHyp) {
+ double h = nearSourceSaturationTerm(Mw, c.mc, isSlab);
+ return c.c0 +
+ getPathTerm(c, rRup, Mw, h) +
+ getMagnitudeTerm(c, Mw) +
+ getDepthTerm(c, isSlab, zHyp);
+ }
/* Equation 4.1 */
private static final double calcMean(
final Coefficients c,
- final boolean slab,
- final boolean useBasinTerm,
- final double mc,
- final double pgaRef,
+ final boolean isSlab,
+ final String basin,
final double Mw,
final double rRup,
final double zHyp,
+ final double pgaRef,
final double vs30,
final double z2p5) {
/* Equation 4.7: Fs = Flin + Fnl + Fb */
double siteTerm = getSiteTermLinear(c, vs30) + getSiteTermNonLinear(c, pgaRef, vs30) +
- getSiteBasinTerm(c, useBasinTerm, z2p5);
+ getSiteBasinTerm(c, basin, vs30, z2p5);
- double h = nearSourceSaturationTerm(Mw, mc, slab);
+ double h = nearSourceSaturationTerm(Mw, c.mc, isSlab);
double μ = c.c0 +
getPathTerm(c, rRup, Mw, h) +
- getMagnitudeTerm(c, mc, Mw) +
- getDepthTerm(c, slab, zHyp) +
+ getMagnitudeTerm(c, Mw) +
+ getDepthTerm(c, isSlab, zHyp) +
siteTerm;
- return 0.0;
+ return μ;
}
/* Equations 4.2 and 4.3 */
@@ -152,9 +296,7 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
private static final double nearSourceSaturationTerm(double Mw, double mc, boolean slab) {
if (slab) {
/* Equation 4.4b */
- return (Mw <= mc)
- ? pow(10, (1.050 / (mc - 4.0)) * (Mw - mc) + 1.544)
- : 35;
+ return (Mw <= mc) ? pow(10, (1.050 / (mc - 4.0)) * (Mw - mc) + 1.544) : 35;
} else {
/* Equation 4.4a */
return pow(10.0, -0.82 + 0.252 * Mw);
@@ -162,11 +304,9 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
}
/* Equation 4.5 */
- private static final double getMagnitudeTerm(final Coefficients c, double mc, double Mw) {
- double dM = Mw - mc;
- return (Mw <= mc)
- ? c.c4 * dM + c.c5 * dM * dM
- : c.c6 * dM;
+ private static final double getMagnitudeTerm(final Coefficients c, double Mw) {
+ double dM = Mw - c.mc;
+ return (Mw <= c.mc) ? c.c4 * dM + c.c5 * dM * dM : c.c6 * dM;
}
/* Equation 4.6 */
@@ -190,12 +330,12 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
/* Equation 4.8 */
private static final double getSiteTermLinear(final Coefficients c, double vs30) {
double Flin;
- if (vs30 <= V1) {
- Flin = c.s1 * log(vs30 / V1) + c.s2 * log(V1 / VREF);
+ if (vs30 <= c.v1) {
+ Flin = c.s1 * log(vs30 / c.v1) + c.s2 * log(c.v1 / c.vRef);
} else if (vs30 > c.v2) {
- Flin = c.s2 * log(c.v2 / VREF);
+ Flin = c.s2 * log(c.v2 / c.vRef);
} else {
- Flin = c.s2 * log(vs30 / VREF);
+ Flin = c.s2 * log(vs30 / c.vRef);
}
return Flin;
}
@@ -203,19 +343,54 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
/* Equation 4.9 */
private static final double getSiteTermNonLinear(final Coefficients c, double pgaRef,
double vs30) {
+ double f3 = 0.05;
double f2 = c.f4 * (exp(c.f5 * (min(vs30, 760) - 200)) - exp(c.f5 * (760 - 200)));
- return 0.0 + f2 * log((pgaRef + 0.05) / 0.05);
+ return 0.0 + f2 * log((pgaRef + f3) / f3);
}
- /* Equations 4.11 - 4.13 */
- private static final double getSiteBasinTerm(final Coefficients c, boolean basin, double z2p5) {
- if (!basin) {
+ /* Equations 4.11/5.5 - 4.13/5.7 */
+ private static final double getSiteBasinTerm(final Coefficients c, String basin, double vs30,
+ double z2p5) {
+ if (basin.equals(Basin.NONE) || Double.isNaN(z2p5)) {
return 0.0;
+ }
+ // Cascadia coefficients
+ double θ0 = 3.94;
+ double θ1 = -0.42;
+ double νμ = 200.0;
+ double νσ = 0.2;
+
+ double e1 = c.c_e1;
+ double e2;
+ double e3;
+
+ if (basin.equals(Basin.PACNW_NONE)) {
+ e2 = c.c_e2 + c.del_None;
+ e3 = c.c_e3 + c.del_None;
+ } else if (basin.equals(Basin.SEATTLE)) {
+ e2 = c.c_e2 + c.del_Seattle;
+ e3 = c.c_e3 + c.del_Seattle;
+ } else if (basin.equals(Basin.PACNW)) {
+ e2 = c.c_e2;
+ e3 = c.c_e3;
} else {
- // TODO: implement for Cascadia - this term is 0 except in Japan and
- // Cascadia
- return 0.0;
+ throw new IllegalArgumentException("Basin [" + basin + "] not supported");
}
+
+ // Eq. 4.13/5.7
+ double μz2p5 = exp(
+ log(10) * θ1 * (1 + Maths.erf((log10(vs30) - log10(νμ)) / νσ / sqrt(2)) + log(10) * θ0));
+ double δz2p5 = log(z2p5) - log(μz2p5);
+
+ double Fb;
+ if (δz2p5 <= (e1 / e3)) {
+ Fb = e1;
+ } else if (δz2p5 >= (e2 / e3)) {
+ Fb = e2;
+ } else {
+ Fb = e3 * δz2p5;
+ }
+ return Fb;
}
/*
@@ -234,15 +409,14 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
}
@Override
- final boolean useBasinTerm() {
- return false;
+ String getRegion() {
+ return Region.GLOBAL;
}
@Override
- final double getMc() {
- return 7.9;
+ String getBasin() {
+ return Basin.NONE;
}
-
}
/*
@@ -261,13 +435,183 @@ public abstract class ParkerEtAl_2020 implements GroundMotionModel {
}
@Override
- final boolean useBasinTerm() {
- return false;
+ String getRegion() {
+ return Region.GLOBAL;
+ }
+
+ @Override
+ String getBasin() {
+ return Basin.NONE;
+ }
+ }
+
+ /*
+ * Cascadia interface model
+ */
+ static class CascadiaInterface extends Interface {
+ static final String NAME = Interface.NAME + " : Cascadia";
+
+ CascadiaInterface(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getRegion() {
+ return Region.CASCADIA;
+ }
+
+ @Override
+ String getBasin() {
+ return Basin.PACNW_NONE;
+ }
+ }
+
+ /*
+ * Cascadia intraslab model
+ */
+ static class CascadiaSlab extends Slab {
+ static final String NAME = Slab.NAME + " : Cascadia";
+
+ CascadiaSlab(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getRegion() {
+ return Region.CASCADIA;
+ }
+
+ @Override
+ String getBasin() {
+ return Basin.PACNW_NONE;
+ }
+ }
+
+ /*
+ * Cascadia interface model, Seattle basin
+ */
+ static class CascadiaSeattleInterface extends CascadiaInterface {
+ static final String NAME = CascadiaInterface.NAME + " : Seattle";
+
+ CascadiaSeattleInterface(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getBasin() {
+ return Basin.SEATTLE;
+ }
+ }
+
+ /*
+ * Cascadia intraslab model, Seattle basin
+ */
+ static class CascadiaSeattleSlab extends CascadiaSlab {
+ static final String NAME = CascadiaSlab.NAME + " : Seattle";
+
+ CascadiaSeattleSlab(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getBasin() {
+ return Basin.SEATTLE;
+ }
+ }
+
+ /*
+ * Cascadia interface model, Other PacNW basin
+ */
+ static class CascadiaPacNWInterface extends CascadiaInterface {
+ static final String NAME = CascadiaInterface.NAME + " : PacNW";
+
+ CascadiaPacNWInterface(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getBasin() {
+ return Basin.PACNW;
+ }
+ }
+
+ /*
+ * Cascadia intraslab model, Other PacNW basin
+ */
+ static class CascadiaPacNWSlab extends CascadiaSlab {
+ static final String NAME = CascadiaSlab.NAME + " : PacNW";
+
+ CascadiaPacNWSlab(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getBasin() {
+ return Basin.PACNW;
+ }
+ }
+
+ /*
+ * Alaska interface model
+ */
+ static class AlaskaInterface extends Interface {
+ static final String NAME = Interface.NAME + " : Alaska";
+
+ AlaskaInterface(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getRegion() {
+ return Region.ALASKA;
+ }
+ }
+
+ /*
+ * Alaska intraslab model
+ */
+ static class AlaskaSlab extends Slab {
+ static final String NAME = Slab.NAME + " : Alaska";
+
+ AlaskaSlab(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getRegion() {
+ return Region.ALASKA;
+ }
+ }
+
+ /*
+ * Aleutians interface model
+ */
+ static class AleutiansInterface extends Interface {
+ static final String NAME = Interface.NAME + " : Aleutians";
+
+ AleutiansInterface(Imt imt) {
+ super(imt);
+ }
+
+ @Override
+ String getRegion() {
+ return Region.ALEUTIANS;
+ }
+ }
+
+ /*
+ * Aleutians intraslab model
+ */
+ static class AleutiansSlab extends Slab {
+ static final String NAME = Slab.NAME + " : Aleutians";
+
+ AleutiansSlab(Imt imt) {
+ super(imt);
}
@Override
- final double getMc() {
- return 7.6;
+ String getRegion() {
+ return Region.ALEUTIANS;
}
}
diff --git a/src/main/resources/gmm/coeffs/PSHAB20_interface.csv b/src/main/resources/gmm/coeffs/PSHAB20_interface.csv
new file mode 100644
index 00000000..d8528531
--- /dev/null
+++ b/src/main/resources/gmm/coeffs/PSHAB20_interface.csv
@@ -0,0 +1,27 @@
+Period (s),Global_c0,Alaska_c0,Aleutian_c0,Cascadia_c0,c1,b4,Global_a0,Alaska_a0,c4,c5,c6,V1,V2,Vref,Global_s2,Alaska_s2,Cascadia_s2,f4,f5,C_e1,C_e2,C_e3,del_None,del_Seattle,Tau,phi21,phi22,phi2V,VM,phi2S2S0,a1,phi2SS1,phi2SS2,a2
+PGV,8.097,9.283796298,8.374796298,7.728,-1.661,0.1,-3.95E-03,-4.04E-03,1.336,-0.039,1.336,270,850,760,-0.601,-1.031,-0.671,-0.31763,-0.0052,0,0.115,0.068,-0.115,0,0.477,0.348,0.288,-0.179,423,0.142,0.047,0.153,0.166,0.011
+PGA,4.082,4.458796298,3.652796298,3.856,-1.662,0.1,-6.57E-03,-5.41E-03,1.246,-0.021,1.128,270,1350,760,-0.498,-0.785,-0.572,-0.44169,-0.0052,0,0,1,0,0,0.48,0.396,0.565,-0.18,423,0.221,0.093,0.149,0.327,0.068
+0.01,3.714,4.094796298,3.288796298,3.488,-1.587,0.1,-6.57E-03,-5.41E-03,1.246,-0.021,1.128,270,1300,760,-0.498,-0.803,-0.571,-0.4859,-0.0052,0,0,1,0,0,0.476,0.397,0.56,-0.18,423,0.223,0.098,0.148,0.294,0.071
+0.02,3.762,4.132796298,3.338796298,3.536,-1.593,0.1,-6.57E-03,-5.41E-03,1.227,-0.021,1.128,270,1225,760,-0.478,-0.785,-0.575,-0.4859,-0.00518,0,0,1,0,0,0.482,0.401,0.563,-0.181,423,0.227,0.105,0.149,0.294,0.073
+0.025,3.859,4.246796298,3.392796298,3.633,-1.607,0.1,-6.57E-03,-5.41E-03,1.221,-0.021,1.128,270,1200,760,-0.464,-0.745,-0.573,-0.4859,-0.00515,0,0,1,0,0,0.49,0.405,0.575,-0.183,423,0.231,0.12,0.15,0.31,0.076
+0.03,4.014,4.386796298,3.535796298,3.788,-1.63,0.1,-6.57E-03,-5.41E-03,1.215,-0.021,1.128,270,1200,760,-0.446,-0.69,-0.565,-0.4908,-0.00511,0,0,1,0,0,0.5,0.413,0.589,-0.188,423,0.239,0.145,0.153,0.313,0.077
+0.04,4.223,4.553796298,3.747796298,3.997,-1.657,0.1,-6.57E-03,-5.41E-03,1.207,-0.021,1.128,270,1200,760,-0.431,-0.636,-0.546,-0.49569,-0.00505,0,0,1,0,0,0.515,0.439,0.616,-0.205,423,0.261,0.177,0.159,0.322,0.08
+0.05,4.456,4.745796298,3.959796298,4.23,-1.687,0.1,-6.57E-03,-5.41E-03,1.201,-0.021,1.128,270,1225,760,-0.42,-0.594,-0.519,-0.49823,-0.00497,0.1,-0.1,-0.063,-0.05,0,0.528,0.473,0.653,-0.23,423,0.285,0.2,0.167,0.33,0.077
+0.075,4.742,4.972796298,4.231796298,4.516,-1.715,0.1,-6.57E-03,-5.41E-03,1.19,-0.021,1.128,270,1350,760,-0.442,-0.586,-0.497,-0.49724,-0.00489,0.3,-0.34,-0.2,-0.075,0.078,0.53,0.529,0.722,-0.262,423,0.339,0.205,0.184,0.299,0.063
+0.1,4.952,5.160796298,4.471796298,4.726,-1.737,0.1,-6.57E-03,-5.41E-03,1.182,-0.021,1.128,270,1450,760,-0.485,-0.629,-0.486,-0.49471,-0.00478,0.333,-0.377,-0.222,-0.081,0.075,0.524,0.517,0.712,-0.239,423,0.347,0.185,0.176,0.31,0.061
+0.15,5.08,5.285796298,4.665796298,4.848,-1.745,0.1,-6.57E-03,-5.41E-03,1.171,-0.021,1.162,270,1500,760,-0.546,-0.729,-0.499,-0.48583,-0.0046,0.29,-0.29,-0.193,-0.091,0.064,0.51,0.457,0.644,-0.185,423,0.313,0.123,0.164,0.307,0.076
+0.2,5.035,5.277796298,4.661796298,4.798,-1.732,0.1,-6.57E-03,-5.41E-03,1.163,-0.021,1.163,270,1425,760,-0.612,-0.867,-0.533,-0.47383,-0.00434,0.177,-0.192,-0.148,-0.092,0.075,0.501,0.432,0.64,-0.138,423,0.277,0.11,0.163,0.301,0.07
+0.25,4.859,5.154796298,4.503796298,4.618,-1.696,0.1,-6.57E-03,-5.41E-03,1.156,-0.021,1.156,270,1350,760,-0.688,-1.011,-0.592,-0.47696,-0.00402,0.1,-0.035,-0.054,0,0,0.492,0.45,0.633,-0.185,423,0.26,0.119,0.169,0.233,0.077
+0.3,4.583,4.910796298,4.276796298,4.34,-1.643,0.1,-6.57E-03,-5.41E-03,1.151,-0.021,1.151,270,1250,760,-0.748,-1.133,-0.681,-0.4845,-0.0037,0,0,1,0,0,0.492,0.436,0.584,-0.158,423,0.254,0.092,0.159,0.22,0.065
+0.4,4.18,4.548796298,3.919796298,3.935,-1.58,0.1,-6.57E-03,-5.41E-03,1.143,-0.022,1.143,270,1150,760,-0.802,-1.238,-0.772,-0.48105,-0.00342,0,0.05,0.2,0,0,0.492,0.433,0.556,-0.19,423,0.23,0.044,0.158,0.222,0.064
+0.5,3.752,4.168796298,3.486796298,3.505,-1.519,0.1,-6.57E-03,-5.41E-03,1.143,-0.023,1.143,270,1025,760,-0.845,-1.321,-0.838,-0.46492,-0.00322,0,0.1,0.2,0,0,0.492,0.428,0.51,-0.186,423,0.225,0.038,0.16,0.243,0.044
+0.75,3.085,3.510796298,2.710796298,2.837,-1.44,0.1,-6.35E-03,-4.78E-03,1.217,-0.026,1.217,270,900,760,-0.911,-1.383,-0.922,-0.43439,-0.00312,0,0.2,0.125,-0.2,0.012,0.492,0.448,0.471,-0.177,422,0.218,0.04,0.175,0.241,0.04
+1,2.644,3.067796298,2.238796298,2.396,-1.419,0.1,-5.80E-03,-4.15E-03,1.27,-0.028,1.24,270,800,760,-0.926,-1.414,-0.932,-0.38484,-0.0031,0,0.245,0.153,-0.245,0.037,0.492,0.43,0.43,-0.166,422,0.227,0.015,0.195,0.195,0.043
+1.5,2.046,2.513796298,1.451796298,1.799,-1.4,0.1,-5.05E-03,-3.42E-03,1.344,-0.031,1.237,270,760,760,-0.888,-1.43,-0.814,-0.32318,-0.0031,0,0.32,0.2,-0.32,0.064,0.492,0.406,0.406,-0.111,422,0.244,-0.047,0.204,0.204,-0.034
+2,1.556,2.061796298,0.906796298,1.31,-1.391,0.1,-4.29E-03,-2.90E-03,1.396,-0.034,1.232,270,760,760,-0.808,-1.421,-0.725,-0.26577,-0.0031,0,0.37,0.239,-0.28,0.14,0.492,0.393,0.393,0,422,0.231,-0.036,0.196,0.196,-0.036
+2.5,1.167,1.709796298,0.392796298,0.922,-1.394,0.1,-3.69E-03,-2.50E-03,1.437,-0.036,1.227,270,760,760,-0.743,-1.391,-0.632,-0.21236,-0.0031,0,0.4,0.264,-0.313,0.19,0.492,0.381,0.381,0,421,0.222,-0.025,0.169,0.169,-0.029
+3,0.92,1.456796298,0.099796298,0.675,-1.416,0.1,-3.21E-03,-2.17E-03,1.47,-0.038,1.223,270,760,760,-0.669,-1.343,-0.57,0,-0.0031,0,0.43,0.287,-0.355,0.165,0.492,0.367,0.367,0,419,0.199,-0.03,0.177,0.177,-0.011
+4,0.595,1.207796298,-0.356203702,0.352,-1.452,0.1,-2.44E-03,-1.65E-03,1.523,-0.044,1.216,270,760,760,-0.585,-1.297,-0.489,0,-0.0031,0,0.44,0.303,-0.417,0.163,0.492,0.33,0.33,0,416,0.191,-0.042,0.158,0.158,0.033
+5,0.465,1.131796298,-0.601203702,0.223,-1.504,0.1,-1.60E-03,-1.25E-03,1.564,-0.048,1.21,270,760,760,-0.506,-1.233,-0.421,0,-0.0031,0,0.45,0.321,-0.45,0.132,0.492,0.298,0.298,0,415,0.181,0.005,0.132,0.132,0.014
+7.5,0.078,0.758796298,-1.137203702,-0.162,-1.569,0.1,-7.66E-04,-5.19E-04,1.638,-0.059,1.2,270,760,760,-0.418,-1.147,-0.357,0,-0.0031,0,0.406,0.312,-0.35,0.15,0.492,0.254,0.254,0,419,0.181,-0.016,0.113,0.113,0.016
+10,0.046,0.708796298,-1.290203702,-0.193,-1.676,0.1,0.00E+00,0.00E+00,1.69,-0.067,1.194,270,760,760,-0.321,-1.06,-0.302,0,-0.0031,0,0.345,0.265,-0.331,0.117,0.492,0.231,0.231,0,427,0.181,0.04,0.11,0.11,0.017
diff --git a/src/main/resources/gmm/coeffs/PSHAB20_slab.csv b/src/main/resources/gmm/coeffs/PSHAB20_slab.csv
new file mode 100644
index 00000000..78a98065
--- /dev/null
+++ b/src/main/resources/gmm/coeffs/PSHAB20_slab.csv
@@ -0,0 +1,27 @@
+T,Global_c0,Alaska_c0,Aleutian_c0,Cascadia_c0,c1,b4,Global_a0,Alaska_a0,Cascadia_a0,c4,c5,c6,d,m,db,V1,V2,Vref,Global_s2,Alaska_s2,Cascadia_s2,f4,f5,C_e1,C_e2,C_e3,del_None,del_Seattle,Tau,phi21,phi22,phi2V,VM,phi2S2S0,a1,phi2SS1,phi2SS2,a2
+PGV,13.194,12.79,13.6,12.874,-2.422,0.1,-1.90E-03,-2.38E-03,-1.09E-03,1.84,-0.05,0.8,0.2693,0.0252,67,270,850,760,-0.601,-1.031,-0.671,-0.31763,-0.0052,0,0.115,0.068,-0.115,0,0.477,0.348,0.288,-0.179,423,0.142,0.047,0.153,0.166,0.011
+PGV,9.907,9.404,9.912,9.6,-2.543,0.1,-2.55E-03,-2.27E-03,-3.54E-03,1.84,-0.05,0.4,0.3004,0.0314,67,270,1350,760,-0.498,-0.785,-0.572,-0.44169,-0.0052,0,0,1,0,0,0.48,0.396,0.565,-0.18,423,0.221,0.093,0.149,0.327,0.068
+0.01,9.962,9.451,9.954,9.802,-2.554,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.84,-0.05,0.4,0.2839,0.0296,67,270,1300,760,-0.498,-0.803,-0.571,-0.4859,-0.0052,0,0,1,0,0,0.476,0.397,0.56,-0.18,423,0.223,0.098,0.148,0.294,0.071
+0.02,10.099,9.587,10.086,9.933,-2.566,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.05,0.415,0.2854,0.0298,67,270,1225,760,-0.478,-0.785,-0.575,-0.4859,-0.00518,0,0,1,0,0,0.482,0.401,0.563,-0.181,423,0.227,0.105,0.149,0.294,0.073
+0.025,10.181,9.667,10.172,10.009,-2.578,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.05,0.43,0.2891,0.0302,67,270,1200,760,-0.464,-0.745,-0.573,-0.4859,-0.00515,0,0,1,0,0,0.49,0.405,0.575,-0.183,423,0.231,0.12,0.15,0.31,0.076
+0.03,10.311,9.808,10.302,10.133,-2.594,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.05,0.445,0.2932,0.0306,67,270,1200,760,-0.446,-0.69,-0.565,-0.4908,-0.00511,0,0,1,0,0,0.5,0.413,0.589,-0.188,423,0.239,0.145,0.153,0.313,0.077
+0.04,10.588,10.086,10.602,10.404,-2.629,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.05,0.46,0.3004,0.0313,67,270,1200,760,-0.431,-0.636,-0.546,-0.49569,-0.00505,0,0,1,0,0,0.515,0.439,0.616,-0.205,423,0.261,0.177,0.159,0.322,0.08
+0.05,10.824,10.379,10.862,10.634,-2.649,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.05,0.475,0.3048,0.0316,67,270,1225,760,-0.42,-0.594,-0.519,-0.49823,-0.00497,0.1,-0.1,-0.063,-0.05,0,0.528,0.473,0.653,-0.23,423,0.285,0.2,0.167,0.33,0.077
+0.075,11.084,10.65,11.184,10.888,-2.65,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.05,0.49,0.2992,0.0321,67,270,1350,760,-0.442,-0.586,-0.497,-0.49724,-0.00489,0.3,-0.34,-0.2,-0.075,0.078,0.53,0.529,0.722,-0.262,423,0.339,0.205,0.184,0.299,0.063
+0.1,11.232,10.816,11.304,11.03,-2.647,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.05,0.505,0.2854,0.032,67,270,1450,760,-0.485,-0.629,-0.486,-0.49471,-0.00478,0.333,-0.377,-0.222,-0.081,0.075,0.524,0.517,0.712,-0.239,423,0.347,0.185,0.176,0.31,0.061
+0.15,11.311,10.883,11.402,11.103,-2.634,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.06,0.52,0.2814,0.0325,67,270,1500,760,-0.546,-0.729,-0.499,-0.48583,-0.0046,0.29,-0.29,-0.193,-0.091,0.064,0.51,0.457,0.644,-0.185,423,0.313,0.123,0.164,0.307,0.076
+0.2,11.055,10.633,11.183,10.841,-2.583,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.068,0.535,0.291,0.0306,67,270,1425,760,-0.612,-0.867,-0.533,-0.47383,-0.00434,0.177,-0.192,-0.148,-0.092,0.075,0.501,0.432,0.64,-0.138,423,0.277,0.11,0.163,0.301,0.07
+0.25,10.803,10.322,10.965,10.583,-2.539,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.075,0.55,0.2758,0.0306,67,270,1350,760,-0.688,-1.011,-0.592,-0.47696,-0.00402,0.1,-0.035,-0.054,0,0,0.492,0.45,0.633,-0.185,423,0.26,0.119,0.169,0.233,0.077
+0.3,10.669,10.116,10.87,10.443,-2.528,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.082,0.565,0.2719,0.0323,67,270,1250,760,-0.748,-1.133,-0.681,-0.4845,-0.0037,0,0,1,0,0,0.492,0.436,0.584,-0.158,423,0.254,0.092,0.159,0.22,0.065
+0.4,10.116,9.561,10.411,9.884,-2.452,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.091,0.58,0.2539,0.0302,67,270,1150,760,-0.802,-1.238,-0.772,-0.48105,-0.00342,0,0.05,0.2,0,0,0.492,0.433,0.556,-0.19,423,0.23,0.044,0.158,0.222,0.064
+0.5,9.579,8.973,9.901,9.341,-2.384,0.1,-2.55E-03,-2.19E-03,-4.01E-03,1.884,-0.1,0.595,0.2482,0.0295,67,270,1025,760,-0.845,-1.321,-0.838,-0.46492,-0.00322,0,0.1,0.2,0,0,0.492,0.428,0.51,-0.186,423,0.225,0.038,0.16,0.243,0.044
+0.75,8.837,8.246,9.335,8.593,-2.338,0.1,-2.11E-03,-1.89E-03,-3.47E-03,1.884,-0.115,0.61,0.2227,0.0266,67,270,900,760,-0.911,-1.383,-0.922,-0.43439,-0.00312,0,0.2,0.125,-0.2,0.012,0.492,0.448,0.471,-0.177,422,0.218,0.04,0.175,0.241,0.04
+1,8.067,7.507,8.68,7.817,-2.267,0.1,-1.87E-03,-1.68E-03,-3.09E-03,1.884,-0.134,0.625,0.1969,0.0231,67,270,800,760,-0.926,-1.414,-0.932,-0.38484,-0.0031,0,0.245,0.153,-0.245,0.037,0.492,0.43,0.43,-0.166,422,0.227,0.015,0.195,0.195,0.043
+1.5,6.829,6.213,7.581,6.573,-2.166,0.1,-1.54E-03,-1.39E-03,-2.54E-03,1.884,-0.154,0.64,0.1452,0.0118,67,270,760,760,-0.888,-1.43,-0.814,-0.32318,-0.0031,0,0.32,0.2,-0.32,0.064,0.492,0.406,0.406,-0.111,422,0.244,-0.047,0.204,0.204,-0.034
+2,5.871,5.206,6.671,5.609,-2.077,0.1,-1.31E-03,-1.18E-03,-2.16E-03,1.884,-0.154,0.655,0.06,0.007,67,270,760,760,-0.808,-1.421,-0.725,-0.26577,-0.0031,0,0.37,0.239,-0.28,0.14,0.492,0.393,0.393,0,422,0.231,-0.036,0.196,0.196,-0.036
+2.5,5.2,4.594,6.047,4.932,-2.015,0.1,-1.13E-03,-1.01E-03,-1.86E-03,1.884,-0.154,0.67,0,0,0,270,760,760,-0.743,-1.391,-0.632,-0.21236,-0.0031,0,0.4,0.264,-0.313,0.19,0.492,0.381,0.381,0,421,0.222,-0.025,0.169,0.169,-0.029
+3,4.83,4.206,5.667,4.556,-2.012,0.1,-9.79E-04,-8.80E-04,-1.61E-03,1.949,-0.154,0.685,0,0,0,270,760,760,-0.669,-1.343,-0.57,-0.17807,-0.0031,0,0.43,0.287,-0.355,0.165,0.492,0.367,0.367,0,419,0.199,-0.03,0.177,0.177,-0.011
+4,4.173,3.517,4.97,3.893,-1.989,0.1,-7.45E-04,-6.70E-04,-1.23E-03,2.031,-0.154,0.7,0,0,0,270,760,760,-0.585,-1.297,-0.489,-0.13729,-0.0031,0,0.44,0.303,-0.417,0.163,0.492,0.33,0.33,0,416,0.191,-0.042,0.158,0.158,0.033
+5,3.833,3.142,4.592,3.547,-1.998,0.1,-5.64E-04,-5.07E-04,-9.29E-04,2.131,-0.154,0.715,0,0,0,270,760,760,-0.506,-1.233,-0.421,-0.07733,-0.0031,0,0.45,0.321,-0.45,0.132,0.492,0.298,0.298,0,415,0.181,0.005,0.132,0.132,0.014
+7.5,3.132,2.391,3.65,2.84,-2.019,0.1,-2.34E-04,-2.10E-04,-3.85E-04,2.185,-0.154,0.73,0,0,0,270,760,760,-0.418,-1.147,-0.357,-0.05443,-0.0031,0,0.406,0.312,-0.35,0.15,0.492,0.254,0.254,0,419,0.181,-0.016,0.113,0.113,0.016
+10,2.72,2.031,2.95,2.422,-2.047,0.1,0.00E+00,0.00E+00,0.00E+00,2.35,-0.154,0.745,0,0,0,270,760,760,-0.321,-1.06,-0.302,-0.03313,-0.0031,0,0.345,0.265,-0.331,0.117,0.492,0.231,0.231,0,427,0.181,0.04,0.11,0.11,0.017
diff --git a/src/main/resources/gmm/coeffs/ParkerEtAl_2020_interface.csv b/src/main/resources/gmm/coeffs/ParkerEtAl_2020_interface.csv
deleted file mode 100644
index e7854a32..00000000
--- a/src/main/resources/gmm/coeffs/ParkerEtAl_2020_interface.csv
+++ /dev/null
@@ -1,27 +0,0 @@
-T,c0_global,c0_alaska,c0_aleutian,c0_cascadia,c0_cam_n,c0_cam_s,c0_japan_pac,c0_japan_phi,c0_sa_n,c0_sa_s,c0_taiwan_e,c0_taiwan_w,c1,b4,a0_global,a0_alaska,a0_cascadia,a0_cam,a0_japan,a0_sa,a0_taiwan,c4,c5,c6,v1,v2,vref,s1_japan,s1_taiwan,s2_global,s2_alaska,s2_cascadia,s2_japan,s2_sa,s2_taiwan,f4,f5,e1_japan,e2_japan,e3_japan,e1_cascadia,e2_cascadia,e3_cascadia,del_none,del_seattle
-PGV,8.097,9.283796298,8.374796298,7.728,7.046899908,7.046899908,8.772125851,7.579125851,8.528671414,8.679671414,7.559846279,7.559846279,-1.661,0.1,-3.95E-03,-4.04E-03,-3.95E-03,-1.53E-03,-2.39E-03,-3.11E-04,-5.14E-03,1.336,-0.039,1.336,270,850,760,-0.738,-0.454,-0.601,-1.031,-0.671,-0.738,-0.681,-0.59,-0.31763,-0.0052,-0.137,0.137,0.091,0,0.115,0.068,-0.115,0
-PGA,4.082,4.458796298,3.652796298,3.856,2.875899908,2.875899908,5.373125851,4.309125851,5.064671414,5.198671414,3.032846279,3.032846279,-1.662,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.246,-0.021,1.128,270,1350,760,-0.586,-0.44,-0.498,-0.785,-0.572,-0.586,-0.333,-0.44,-0.44169,-0.0052,0,0,1,0,0,1,0,0
-0.01,3.714,4.094796298,3.288796298,3.488,2.564899908,2.564899908,5.022125851,3.901125851,4.673671414,4.807671414,2.636846279,2.636846279,-1.587,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.246,-0.021,1.128,270,1300,760,-0.604,-0.44,-0.498,-0.803,-0.571,-0.604,-0.333,-0.44,-0.4859,-0.0052,0,0,1,0,0,1,0,0
-0.02,3.762,4.132796298,3.338796298,3.536,2.636899908,2.636899908,5.066125851,3.935125851,4.694671414,4.827671414,2.698846279,2.698846279,-1.593,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.227,-0.021,1.128,270,1225,760,-0.593,-0.458,-0.478,-0.785,-0.575,-0.593,-0.345,-0.458,-0.4859,-0.00518,0,0,1,0,0,1,0,0
-0.025,3.859,4.246796298,3.392796298,3.633,2.731899908,2.731899908,5.140125851,4.094125851,4.779671414,4.911671414,2.800846279,2.800846279,-1.607,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.221,-0.021,1.128,270,1200,760,-0.569,-0.454,-0.464,-0.745,-0.573,-0.579,-0.362,-0.459,-0.4859,-0.00515,0,0,1,0,0,1,0,0
-0.03,4.014,4.386796298,3.535796298,3.788,2.890899908,2.890899908,5.317125851,4.278125851,4.935671414,5.066671414,2.926846279,2.926846279,-1.63,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.215,-0.021,1.128,270,1200,760,-0.539,-0.455,-0.446,-0.69,-0.565,-0.561,-0.38,-0.464,-0.4908,-0.00511,0,0,1,0,0,1,0,0
-0.04,4.223,4.553796298,3.747796298,3.997,3.075899908,3.075899908,5.564125851,4.531125851,5.182671414,5.312671414,3.069846279,3.069846279,-1.657,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.207,-0.021,1.128,270,1200,760,-0.468,-0.453,-0.431,-0.636,-0.546,-0.508,-0.403,-0.466,-0.49569,-0.00505,0,0,1,0,0,1,0,0
-0.05,4.456,4.745796298,3.959796298,4.23,3.287899908,3.287899908,5.843125851,4.816125851,5.457671414,5.586671414,3.236846279,3.236846279,-1.687,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.201,-0.021,1.128,270,1225,760,-0.403,-0.452,-0.42,-0.594,-0.519,-0.461,-0.427,-0.468,-0.49823,-0.00497,0,0,1,0.1,-0.1,-0.063,-0.05,0
-0.075,4.742,4.972796298,4.231796298,4.516,3.560899908,3.560899908,6.146125851,5.126125851,5.788671414,5.917671414,3.446846279,3.446846279,-1.715,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.19,-0.021,1.128,270,1350,760,-0.325,-0.456,-0.442,-0.586,-0.497,-0.452,-0.458,-0.473,-0.49724,-0.00489,0.05,-0.043,-0.025,0.3,-0.34,-0.2,-0.075,0.078
-0.1,4.952,5.160796298,4.471796298,4.726,3.788899908,3.788899908,6.346125851,5.333125851,5.998671414,6.126671414,3.643846279,3.643846279,-1.737,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.182,-0.021,1.128,270,1450,760,-0.264,-0.468,-0.485,-0.629,-0.486,-0.498,-0.49,-0.482,-0.49471,-0.00478,0.1,-0.085,-0.05,0.333,-0.377,-0.222,-0.081,0.075
-0.15,5.08,5.285796298,4.665796298,4.848,3.945899908,3.945899908,6.425125851,5.420125851,6.103671414,6.230671414,3.798846279,3.798846279,-1.745,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.171,-0.021,1.162,270,1500,760,-0.25,-0.484,-0.546,-0.729,-0.499,-0.568,-0.536,-0.499,-0.48583,-0.0046,0.164,-0.139,-0.082,0.29,-0.29,-0.193,-0.091,0.064
-0.2,5.035,5.277796298,4.661796298,4.798,3.943899908,3.943899908,6.288125851,5.289125851,6.013671414,6.140671414,3.827846279,3.827846279,-1.732,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.163,-0.021,1.163,270,1425,760,-0.288,-0.498,-0.612,-0.867,-0.533,-0.667,-0.584,-0.522,-0.47383,-0.00434,0.164,-0.139,-0.082,0.177,-0.192,-0.148,-0.092,0.075
-0.25,4.859,5.154796298,4.503796298,4.618,3.800899908,3.800899908,5.972125851,4.979125851,5.849671414,5.974671414,3.765846279,3.765846279,-1.696,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.156,-0.021,1.156,270,1350,760,-0.36,-0.511,-0.688,-1.011,-0.592,-0.781,-0.654,-0.555,-0.47696,-0.00402,0.08,-0.08,-0.053,0.1,-0.035,-0.054,0,0
-0.3,4.583,4.910796298,4.276796298,4.34,3.491899908,3.491899908,5.582125851,4.592125851,5.603671414,5.728671414,3.602846279,3.602846279,-1.643,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.151,-0.021,1.151,270,1250,760,-0.455,-0.514,-0.748,-1.133,-0.681,-0.867,-0.725,-0.596,-0.4845,-0.0037,0,0,1,0,0,1,0,0
-0.4,4.18,4.548796298,3.919796298,3.935,3.128899908,3.128899908,5.091125851,4.089125851,5.151671414,5.277671414,3.343846279,3.343846279,-1.58,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.143,-0.022,1.143,270,1150,760,-0.617,-0.51,-0.802,-1.238,-0.772,-0.947,-0.801,-0.643,-0.48105,-0.00342,-0.13,0.113,0.087,0,0.05,0.2,0,0
-0.5,3.752,4.168796298,3.486796298,3.505,2.640899908,2.640899908,4.680125851,3.571125851,4.719671414,4.848671414,3.028846279,3.028846279,-1.519,0.1,-6.57E-03,-5.41E-03,-6.57E-03,-3.87E-03,-8.62E-03,-3.97E-03,-7.87E-03,1.143,-0.023,1.143,270,1025,760,-0.757,-0.506,-0.845,-1.321,-0.838,-1.003,-0.863,-0.689,-0.46492,-0.00322,-0.2,0.176,0.118,0,0.1,0.2,0,0
-0.75,3.085,3.510796298,2.710796298,2.837,1.987899908,1.987899908,3.906125851,2.844125851,3.995671414,4.129671414,2.499846279,2.499846279,-1.44,0.1,-6.35E-03,-4.78E-03,-6.35E-03,-3.42E-03,-7.63E-03,-3.51E-03,-6.80E-03,1.217,-0.026,1.217,270,900,760,-0.966,-0.5,-0.911,-1.383,-0.922,-1.052,-0.942,-0.745,-0.43439,-0.00312,-0.401,0.284,0.167,0,0.2,0.125,-0.2,0.012
-1,2.644,3.067796298,2.238796298,2.396,1.553899908,1.553899908,3.481125851,2.371125851,3.512671414,3.653671414,2.140846279,2.140846279,-1.419,0.1,-5.80E-03,-4.15E-03,-5.80E-03,-2.97E-03,-6.63E-03,-3.05E-03,-6.05E-03,1.27,-0.028,1.24,270,800,760,-0.986,-0.49,-0.926,-1.414,-0.932,-1.028,-0.96,-0.777,-0.38484,-0.0031,-0.488,0.346,0.203,0,0.245,0.153,-0.245,0.037
-1.5,2.046,2.513796298,1.451796298,1.799,0.990899908,0.990899908,2.870125851,1.779125851,2.875671414,3.023671414,1.645846279,1.645846279,-1.4,0.1,-5.05E-03,-3.42E-03,-5.05E-03,-2.45E-03,-5.46E-03,-2.52E-03,-4.98E-03,1.344,-0.031,1.237,270,760,760,-0.966,-0.486,-0.888,-1.43,-0.814,-0.971,-0.942,-0.79,-0.32318,-0.0031,-0.578,0.48,0.24,0,0.32,0.2,-0.32,0.064
-2,1.556,2.061796298,0.906796298,1.31,0.534899908,0.534899908,2.507125851,1.293125851,2.327671414,2.481671414,1.217846279,1.217846279,-1.391,0.1,-4.29E-03,-2.90E-03,-4.29E-03,-2.08E-03,-4.63E-03,-2.14E-03,-4.23E-03,1.396,-0.034,1.232,270,760,760,-0.901,-0.475,-0.808,-1.421,-0.725,-0.901,-0.891,-0.765,-0.26577,-0.0031,-0.645,0.579,0.254,0,0.37,0.239,-0.28,0.14
-2.5,1.167,1.709796298,0.392796298,0.922,0.186899908,0.186899908,2.160125851,0.895125851,1.950671414,2.111671414,0.871846279,0.871846279,-1.394,0.1,-3.69E-03,-2.50E-03,-3.69E-03,-1.79E-03,-3.99E-03,-1.84E-03,-3.64E-03,1.437,-0.036,1.227,270,760,760,-0.822,-0.453,-0.743,-1.391,-0.632,-0.822,-0.842,-0.724,-0.21236,-0.0031,-0.678,0.609,0.267,0,0.4,0.264,-0.313,0.19
-3,0.92,1.456796298,0.099796298,0.675,-0.087100092,-0.087100092,1.969125851,0.607125851,1.766671414,1.932671414,0.596846279,0.596846279,-1.416,0.1,-3.21E-03,-2.17E-03,-3.21E-03,-1.56E-03,-3.47E-03,-1.60E-03,-3.16E-03,1.47,-0.038,1.223,270,760,760,-0.751,-0.428,-0.669,-1.343,-0.57,-0.751,-0.787,-0.675,0,-0.0031,-0.772,0.635,0.265,0,0.43,0.287,-0.355,0.165
-4,0.595,1.207796298,-0.356203702,0.352,-0.353100092,-0.353100092,1.675125851,0.303125851,1.524671414,1.698671414,0.268846279,0.268846279,-1.452,0.1,-2.44E-03,-1.65E-03,-2.44E-03,-1.18E-03,-2.64E-03,-1.22E-03,-2.41E-03,1.523,-0.044,1.216,270,760,760,-0.68,-0.396,-0.585,-1.297,-0.489,-0.68,-0.706,-0.613,0,-0.0031,-0.699,0.709,0.259,0,0.44,0.303,-0.417,0.163
-5,0.465,1.131796298,-0.601203702,0.223,-0.491100092,-0.491100092,1.601125851,0.183125851,1.483671414,1.665671414,0.014846279,0.014846279,-1.504,0.1,-1.60E-03,-1.25E-03,-1.60E-03,-8.95E-04,-2.00E-03,-9.19E-04,-1.82E-03,1.564,-0.048,1.21,270,760,760,-0.592,-0.353,-0.506,-1.233,-0.421,-0.592,-0.621,-0.536,0,-0.0031,-0.642,0.63,0.215,0,0.45,0.321,-0.45,0.132
-7.5,0.078,0.758796298,-1.137203702,-0.162,-0.837100092,-0.837100092,1.270125851,-0.143874149,1.175671414,1.366671414,-0.446153721,-0.446153721,-1.569,0.1,-7.66E-04,-5.19E-04,-7.66E-04,-3.71E-04,-8.28E-04,-3.82E-04,-7.55E-04,1.638,-0.059,1.2,270,760,760,-0.494,-0.311,-0.418,-1.147,-0.357,-0.52,-0.52,-0.444,0,-0.0031,-0.524,0.306,0.175,0,0.406,0.312,-0.35,0.15
-10,0.046,0.708796298,-1.290203702,-0.193,-0.864100092,-0.864100092,1.364125851,-0.195874149,1.271671414,1.462671414,-0.473153721,-0.473153721,-1.676,0.1,0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,1.69,-0.067,1.194,270,760,760,-0.395,-0.261,-0.321,-1.06,-0.302,-0.395,-0.42,-0.352,0,-0.0031,-0.327,0.182,0.121,0,0.345,0.265,-0.331,0.117
\ No newline at end of file
diff --git a/src/main/resources/gmm/coeffs/ParkerEtAl_2020_slab.csv b/src/main/resources/gmm/coeffs/ParkerEtAl_2020_slab.csv
deleted file mode 100644
index d36e0b98..00000000
--- a/src/main/resources/gmm/coeffs/ParkerEtAl_2020_slab.csv
+++ /dev/null
@@ -1,27 +0,0 @@
-T,c0_global,c0_alaska,c0_aleutian,c0_cascadia,c0_cam,c0_japan,c0_sa_n,c0_sa_s,c0_taiwan,c1,b4,a0_global,a0_alaska,a0_cascadia,a0_cam,a0_japan,a0_sa,a0_taiwan,c4,c5,c6,d,m,db (km),v1,v2,vref,s1_japan,s1_taiwan,s2_global,s2_alaska,s2_cascadia,s2_japan,s2_sa,s2_taiwan,f4,f5,e1_japan,e2_japan,e3_japan,e1_cascadia,e2_cascadia,e3_cascadia,del_none,del_seattle
-PGV,13.194,12.79,13.6,12.874,12.81,13.248,12.754,12.927,13.516,-2.422,0.1,-1.90E-03,-2.38E-03,-1.09E-03,-1.92E-03,-2.15E-03,-1.92E-03,-3.66E-03,1.84,-0.05,0.8,0.2693,0.0252,67,270,850,760,-0.738,-0.454,-0.601,-1.031,-0.671,-0.738,-0.681,-0.59,-0.31763,-0.0052,-0.137,0.137,0.091,0,0.115,0.068,-0.115,0
-PGA,9.907,9.404,9.912,9.6,9.58,10.145,9.254,9.991,10.071,-2.543,0.1,-2.55E-03,-2.27E-03,-3.54E-03,-2.38E-03,-3.35E-03,-2.38E-03,-3.62E-03,1.84,-0.05,0.4,0.3004,0.0314,67,270,1350,760,-0.586,-0.44,-0.498,-0.785,-0.572,-0.586,-0.333,-0.44,-0.44169,-0.0052,0,0,1,0,0,1,0,0
-0.01,9.962,9.451,9.954,9.802,9.612,10.162,9.293,9.994,10.174,-2.554,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.84,-0.05,0.4,0.2839,0.0296,67,270,1300,760,-0.604,-0.44,-0.498,-0.803,-0.571,-0.604,-0.333,-0.44,-0.4859,-0.0052,0,0,1,0,0,1,0,0
-0.02,10.099,9.587,10.086,9.933,9.771,10.306,9.403,10.152,10.273,-2.566,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.05,0.415,0.2854,0.0298,67,270,1225,760,-0.593,-0.458,-0.478,-0.785,-0.575,-0.593,-0.345,-0.458,-0.4859,-0.00518,0,0,1,0,0,1,0,0
-0.025,10.181,9.667,10.172,10.009,9.85,10.387,9.481,10.292,10.329,-2.578,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.05,0.43,0.2891,0.0302,67,270,1200,760,-0.569,-0.454,-0.464,-0.745,-0.573,-0.579,-0.362,-0.459,-0.4859,-0.00515,0,0,1,0,0,1,0,0
-0.03,10.311,9.808,10.302,10.133,9.993,10.498,9.592,10.459,10.451,-2.594,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.05,0.445,0.2932,0.0306,67,270,1200,760,-0.539,-0.455,-0.446,-0.69,-0.565,-0.561,-0.38,-0.464,-0.4908,-0.00511,0,0,1,0,0,1,0,0
-0.04,10.588,10.086,10.602,10.404,10.317,10.744,9.834,10.818,10.678,-2.629,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.05,0.46,0.3004,0.0313,67,270,1200,760,-0.468,-0.453,-0.431,-0.636,-0.546,-0.508,-0.403,-0.466,-0.49569,-0.00505,0,0,1,0,0,1,0,0
-0.05,10.824,10.379,10.862,10.634,10.563,10.981,10.027,11.102,10.86,-2.649,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.05,0.475,0.3048,0.0316,67,270,1225,760,-0.403,-0.452,-0.42,-0.594,-0.519,-0.461,-0.427,-0.468,-0.49823,-0.00497,0,0,1,0.1,-0.1,-0.063,-0.05,0
-0.075,11.084,10.65,11.184,10.888,10.785,11.25,10.265,11.424,11.093,-2.65,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.05,0.49,0.2992,0.0321,67,270,1350,760,-0.325,-0.456,-0.442,-0.586,-0.497,-0.452,-0.458,-0.473,-0.49724,-0.00489,0.05,-0.043,-0.025,0.3,-0.34,-0.2,-0.075,0.078
-0.1,11.232,10.816,11.304,11.03,10.841,11.466,10.467,11.49,11.283,-2.647,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.05,0.505,0.2854,0.032,67,270,1450,760,-0.264,-0.468,-0.485,-0.629,-0.486,-0.498,-0.49,-0.482,-0.49471,-0.00478,0.1,-0.085,-0.05,0.333,-0.377,-0.222,-0.081,0.075
-0.15,11.311,10.883,11.402,11.103,10.809,11.619,10.566,11.32,11.503,-2.634,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.06,0.52,0.2814,0.0325,67,270,1500,760,-0.25,-0.484,-0.546,-0.729,-0.499,-0.568,-0.536,-0.499,-0.48583,-0.0046,0.164,-0.139,-0.082,0.29,-0.29,-0.193,-0.091,0.064
-0.2,11.055,10.633,11.183,10.841,10.519,11.351,10.33,10.927,11.32,-2.583,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.068,0.535,0.291,0.0306,67,270,1425,760,-0.288,-0.498,-0.612,-0.867,-0.533,-0.667,-0.584,-0.522,-0.47383,-0.00434,0.164,-0.139,-0.082,0.177,-0.192,-0.148,-0.092,0.075
-0.25,10.803,10.322,10.965,10.583,10.268,11.063,10.124,10.555,11.147,-2.539,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.075,0.55,0.2758,0.0306,67,270,1350,760,-0.36,-0.511,-0.688,-1.011,-0.592,-0.781,-0.654,-0.555,-0.47696,-0.00402,0.08,-0.08,-0.053,0.1,-0.035,-0.054,0,0
-0.3,10.669,10.116,10.87,10.443,10.134,10.878,10.077,10.328,11.079,-2.528,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.082,0.565,0.2719,0.0323,67,270,1250,760,-0.455,-0.514,-0.748,-1.133,-0.681,-0.867,-0.725,-0.596,-0.4845,-0.0037,0,0,1,0,0,1,0,0
-0.4,10.116,9.561,10.411,9.884,9.598,10.296,9.539,9.639,10.547,-2.452,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.091,0.58,0.2539,0.0302,67,270,1150,760,-0.617,-0.51,-0.802,-1.238,-0.772,-0.947,-0.801,-0.643,-0.48105,-0.00342,-0.13,0.113,0.087,0,0.05,0.2,0,0
-0.5,9.579,8.973,9.901,9.341,9.097,9.711,9.03,9.03,10.049,-2.384,0.1,-2.55E-03,-2.19E-03,-4.01E-03,-2.17E-03,-3.11E-03,-2.17E-03,-3.55E-03,1.884,-0.1,0.595,0.2482,0.0295,67,270,1025,760,-0.757,-0.506,-0.845,-1.321,-0.838,-1.003,-0.863,-0.689,-0.46492,-0.00322,-0.2,0.176,0.118,0,0.1,0.2,0,0
-0.75,8.837,8.246,9.335,8.593,8.324,8.934,8.258,8.258,9.327,-2.338,0.1,-2.11E-03,-1.89E-03,-3.47E-03,-1.88E-03,-2.69E-03,-1.88E-03,-3.07E-03,1.884,-0.115,0.61,0.2227,0.0266,67,270,900,760,-0.966,-0.5,-0.911,-1.383,-0.922,-1.052,-0.942,-0.745,-0.43439,-0.00312,-0.401,0.284,0.167,0,0.2,0.125,-0.2,0.012
-1,8.067,7.507,8.68,7.817,7.557,8.164,7.467,7.417,8.504,-2.267,0.1,-1.87E-03,-1.68E-03,-3.09E-03,-1.67E-03,-2.39E-03,-1.67E-03,-2.73E-03,1.884,-0.134,0.625,0.1969,0.0231,67,270,800,760,-0.986,-0.49,-0.926,-1.414,-0.932,-1.028,-0.96,-0.777,-0.38484,-0.0031,-0.488,0.346,0.203,0,0.245,0.153,-0.245,0.037
-1.5,6.829,6.213,7.581,6.573,6.35,6.896,6.22,6.18,7.204,-2.166,0.1,-1.54E-03,-1.39E-03,-2.54E-03,-1.38E-03,-1.97E-03,-1.38E-03,-2.25E-03,1.884,-0.154,0.64,0.1452,0.0118,67,270,760,760,-0.966,-0.486,-0.888,-1.43,-0.814,-0.971,-0.942,-0.79,-0.32318,-0.0031,-0.578,0.48,0.24,0,0.32,0.2,-0.32,0.064
-2,5.871,5.206,6.671,5.609,5.434,5.935,5.261,5.161,6.227,-2.077,0.1,-1.31E-03,-1.18E-03,-2.16E-03,-1.17E-03,-1.67E-03,-1.17E-03,-1.91E-03,1.884,-0.154,0.655,0.06,0.007,67,270,760,760,-0.901,-0.475,-0.808,-1.421,-0.725,-0.901,-0.891,-0.765,-0.26577,-0.0031,-0.645,0.579,0.254,0,0.37,0.239,-0.28,0.14
-2.5,5.2,4.594,6.047,4.932,4.773,5.234,4.567,4.517,5.517,-2.015,0.1,-1.13E-03,-1.01E-03,-1.86E-03,-1.01E-03,-1.44E-03,-1.01E-03,-1.64E-03,1.884,-0.154,0.67,0,0,0,270,760,760,-0.822,-0.453,-0.743,-1.391,-0.632,-0.822,-0.842,-0.724,-0.21236,-0.0031,-0.678,0.609,0.267,0,0.4,0.264,-0.313,0.19
-3,4.83,4.206,5.667,4.556,4.441,4.849,4.176,4.076,5.157,-2.012,0.1,-9.79E-04,-8.80E-04,-1.61E-03,-8.73E-04,-1.25E-03,-8.73E-04,-1.43E-03,1.949,-0.154,0.685,0,0,0,270,760,760,-0.751,-0.428,-0.669,-1.343,-0.57,-0.751,-0.787,-0.675,-0.17807,-0.0031,-0.772,0.635,0.265,0,0.43,0.287,-0.355,0.165
-4,4.173,3.517,4.97,3.893,3.849,4.074,3.495,3.445,4.55,-1.989,0.1,-7.45E-04,-6.70E-04,-1.23E-03,-6.64E-04,-9.52E-04,-6.64E-04,-1.09E-03,2.031,-0.154,0.7,0,0,0,270,760,760,-0.68,-0.396,-0.585,-1.297,-0.489,-0.68,-0.706,-0.613,-0.13729,-0.0031,-0.699,0.709,0.259,0,0.44,0.303,-0.417,0.163
-5,3.833,3.142,4.592,3.547,3.502,3.814,3.038,3.038,4.229,-1.998,0.1,-5.64E-04,-5.07E-04,-9.29E-04,-5.03E-04,-7.20E-04,-5.03E-04,-8.22E-04,2.131,-0.154,0.715,0,0,0,270,760,760,-0.592,-0.353,-0.506,-1.233,-0.421,-0.592,-0.621,-0.536,-0.07733,-0.0031,-0.642,0.63,0.215,0,0.45,0.321,-0.45,0.132
-7.5,3.132,2.391,3.65,2.84,2.821,3.152,2.368,2.368,3.554,-2.019,0.1,-2.34E-04,-2.10E-04,-3.85E-04,-2.09E-04,-2.99E-04,-2.09E-04,-3.41E-04,2.185,-0.154,0.73,0,0,0,270,760,760,-0.494,-0.311,-0.418,-1.147,-0.357,-0.52,-0.52,-0.444,-0.05443,-0.0031,-0.524,0.306,0.175,0,0.406,0.312,-0.35,0.15
-10,2.72,2.031,2.95,2.422,2.408,2.791,1.939,1.939,3.166,-2.047,0.1,0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,0.00E+00,2.35,-0.154,0.745,0,0,0,270,760,760,-0.395,-0.261,-0.321,-1.06,-0.302,-0.395,-0.42,-0.352,-0.03313,-0.0031,-0.327,0.182,0.121,0,0.345,0.265,-0.331,0.117
\ No newline at end of file
diff --git a/src/test/java/gov/usgs/earthquake/nshmp/gmm/PSHAB20.java b/src/test/java/gov/usgs/earthquake/nshmp/gmm/PSHAB20.java
new file mode 100644
index 00000000..20816998
--- /dev/null
+++ b/src/test/java/gov/usgs/earthquake/nshmp/gmm/PSHAB20.java
@@ -0,0 +1,42 @@
+package gov.usgs.earthquake.nshmp.gmm;
+
+import static gov.usgs.earthquake.nshmp.gmm.Gmm.PSHAB20_INTERFACE;
+
+import java.io.IOException;
+import java.util.EnumSet;
+import java.util.Set;
+import java.util.stream.Stream;
+
+import org.junit.jupiter.api.extension.ExtensionContext;
+import org.junit.jupiter.params.ParameterizedTest;
+import org.junit.jupiter.params.provider.Arguments;
+import org.junit.jupiter.params.provider.ArgumentsProvider;
+import org.junit.jupiter.params.provider.ArgumentsSource;
+
+class PSHAB20 implements ArgumentsProvider {
+
+ private static String GMM_INPUTS = "/gmm/PSHAB20_inputs.csv";
+ private static String GMM_RESULTS = "/gmm/PSHAB20_results.csv";
+ // private static String GMM_INPUTS = "/gmm/nga-west-inputs.csv";
+ // private static String GMM_RESULTS = "/gmm/nga-west2-results.csv";
+
+ @ParameterizedTest(name = "({index}) {0} {2} {1}")
+ @ArgumentsSource(PSHAB20.class)
+ void test(int index, Gmm gmm, Imt imt, double exMedian, double exSigma, String inputs) {
+ GmmTest.test(index, gmm, imt, exMedian, exSigma, inputs);
+ }
+
+ @Override
+ public Stream<? extends Arguments> provideArguments(ExtensionContext context) throws Exception {
+ return GmmTest.loadResults(GMM_RESULTS, GMM_INPUTS);
+ }
+
+ /* Result generation sets */
+ private static Set<Gmm> gmms = EnumSet.of(PSHAB20_INTERFACE);
+ private static Set<Imt> imts = PSHAB20_INTERFACE.supportedIMTs();
+
+ public static void main(String[] args) throws IOException {
+ GmmTest.generateResults(gmms, imts, GMM_INPUTS, GMM_RESULTS);
+ }
+
+}
diff --git a/src/test/resources/gmm/PSHAB20_inputs.csv b/src/test/resources/gmm/PSHAB20_inputs.csv
new file mode 100644
index 00000000..5786a8bc
--- /dev/null
+++ b/src/test/resources/gmm/PSHAB20_inputs.csv
@@ -0,0 +1,2 @@
+#Mw,rjB,rRup,rX,dip,width,zTop,zHyp,rake,vs30,vsInf,z1p0,z2p5
+8,0,200,0,90,0,0,55,90,760,0,550,NaN
diff --git a/src/test/resources/gmm/PSHAB20_results.csv b/src/test/resources/gmm/PSHAB20_results.csv
new file mode 100644
index 00000000..d1c97d3b
--- /dev/null
+++ b/src/test/resources/gmm/PSHAB20_results.csv
@@ -0,0 +1,23 @@
+0-PSHAB20_INTERFACE-SA0P01,0.020897,0.885198
+0-PSHAB20_INTERFACE-SA0P02,0.021238,0.890687
+0-PSHAB20_INTERFACE-SA0P03,0.022458,0.907193
+0-PSHAB20_INTERFACE-SA0P05,0.025828,0.954874
+0-PSHAB20_INTERFACE-SA0P075,0.029636,0.984835
+0-PSHAB20_INTERFACE-SA0P1,0.032537,0.975488
+0-PSHAB20_INTERFACE-SA0P15,0.035565,0.936536
+0-PSHAB20_INTERFACE-SA0P2,0.036517,0.918151
+0-PSHAB20_INTERFACE-SA0P25,0.037127,0.923073
+0-PSHAB20_INTERFACE-SA0P3,0.037353,0.915458
+0-PSHAB20_INTERFACE-SA0P4,0.034904,0.913818
+0-PSHAB20_INTERFACE-SA0P5,0.031459,0.911078
+0-PSHAB20_INTERFACE-SA0P75,0.025667,0.921989
+0-PSHAB20_INTERFACE-SA1P0,0.020612,0.912175
+0-PSHAB20_INTERFACE-SA1P5,0.014567,0.898924
+0-PSHAB20_INTERFACE-SA2P0,0.010897,0.891664
+0-PSHAB20_INTERFACE-SA3P0,0.006269,0.876963
+0-PSHAB20_INTERFACE-SA4P0,0.004365,0.855607
+0-PSHAB20_INTERFACE-SA5P0,0.003441,0.836698
+0-PSHAB20_INTERFACE-SA7P5,0.001955,0.809978
+0-PSHAB20_INTERFACE-SA10P0,0.001252,0.795653
+0-PSHAB20_INTERFACE-PGA,0.020287,0.886792
+0-PSHAB20_INTERFACE-PGV,2.053897,0.857630
--
GitLab