diff --git a/src/main/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018.java b/src/main/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018.java index 5c62030ec93160745449378c049a8a17d84a2790..8220873d504c602ee74b22857708d389418774de 100644 --- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018.java +++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018.java @@ -34,14 +34,8 @@ import gov.usgs.earthquake.nshmp.tree.LogicTree; * The prediction of GMIMs for site conditions other than Vs30 = 3000 m/s and k0 * = 0.006s requires the use of appropriate site-amplification factors. * - * Note: it appears as if an "average empirical calibration factor of 0.32" (p. - * 2287, in log10 space) should be included in the value of coefficient c1 in - * Table 4 for the PZCT-2ES model but is not. This factor is added here in the - * Coefficients for the empirical-scaling model in order to match the response - * spectra in Figures 3 and 4. - * - * - Table 4, should C10 at 3 sec be negative? There are other coefs that look - * like they may have the wrong sign + * Note: Figure 5 in the paper is incorrect. Pezeshk sent a corrected figure + * that is used to validate the PZCT18-2ES model. * * <p><b>Note:</b> Direct instantiation of {@code GroundMotionModel}s is * prohibited. Use {@link Gmm#instance(Imt)} to retrieve an instance for a @@ -101,16 +95,12 @@ public abstract class PezeshkEtAl_2018 implements GroundMotionModel { Map<String, Double> coeffs_sigma = cc_sigma.get(imt); /* - * An "average empirical calibration factor" of 0.32 (log10-space) - * discussed on p. 2287 apparently did not make it into coefficient c1 in - * Table 5 (PZCT18_2ES), but is included in Table 4 (PZCT18_1SS). - * **clarification from Pezeshk pending, we are taking Table 5 c1 value as - * correct and Figure 4 is incorrect. + * Pezeshk confirmed that the "average empirical calibration factor" of + * 0.32 (log10-space), discussed on p. 2287, is included in coefficient c1 + * in Table 5 (PZCT18_2ES) and sent a corrected Figure 4 for validation. */ - // double c1delta = (model == Model.PZCT18_2ES) ? -0.32 : 0.0; // median model coefficients - // c1 = coeffs.get("c1") + c1delta; c1 = coeffs.get("c1"); c2 = coeffs.get("c2"); c3 = coeffs.get("c3");