From 0a2452d72789c4a7c6439dedde3db8bb892e6752 Mon Sep 17 00:00:00 2001
From: Peter Powers <pmpowers@usgs.gov>
Date: Tue, 4 Jan 2022 16:19:30 -0700
Subject: [PATCH] added check for NaN sed depth in Chapman CPA GMMs
---
.../gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java
index 81a98414..96939fc1 100644
--- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java
+++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/NgaEastUsgs_2017.java
@@ -427,7 +427,9 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
@Override
public LogicTree<GroundMotion> calc(GmmInput in) {
- double cpa = log(ChapmanGuo_2021.cpaPsaRatio(imt, in.z2p5, in.Mw, in.rJB));
+ double cpa = Double.isNaN(in.z2p5)
+ ? 0.0
+ : log(ChapmanGuo_2021.cpaPsaRatio(imt, in.z2p5, in.Mw, in.rJB));
Position p = tables[0].position(in.rRup, in.Mw);
double[] μs = new double[MODEL_COUNT];
for (int i = 0; i < MODEL_COUNT; i++) {
@@ -558,7 +560,9 @@ public abstract class NgaEastUsgs_2017 implements GroundMotionModel {
@Override
public LogicTree<GroundMotion> calc(GmmInput in) {
GmmInput inRock = GmmInput.builder().fromCopy(in).vs30(3000).build();
- double cpa = log(ChapmanGuo_2021.cpaPsaRatio(imt, in.z2p5, in.Mw, in.rJB));
+ double cpa = Double.isNaN(in.z2p5)
+ ? 0.0
+ : log(ChapmanGuo_2021.cpaPsaRatio(imt, in.z2p5, in.Mw, in.rJB));
Position p = tables.values().iterator().next().position(in.rRup, in.Mw);
double[] μs = new double[GMMS.size()];
for (int i = 0; i < GMMS.size(); i++) {
--
GitLab