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 3f6aeeefbe03f2a5d857af1e0c95da9165be159d..9c5fdee03e1662f70c037f5a9fb9c482b54094d8 100644 --- a/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java +++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/Gmm.java @@ -1101,6 +1101,20 @@ public enum Gmm { NgaEastUsgs_2017.COEFFS_SIGMA_PANEL, NgaEastUsgs_2017.CONSTRAINTS), + /** @see PezeshkEtAl_2018 */ + PZCT18_1SS( + PezeshkEtAl_2018.PZCT18_1SS.class, + PezeshkEtAl_2018.PZCT18_1SS.NAME, + PezeshkEtAl_2018.COEFFS_1SS, + PezeshkEtAl_2018.CONSTRAINTS), + + /** @see PezeshkEtAl_2018 */ + PZCT18_2ES( + PezeshkEtAl_2018.PZCT18_2ES.class, + PezeshkEtAl_2018.PZCT18_2ES.NAME, + PezeshkEtAl_2018.COEFFS_2ES, + PezeshkEtAl_2018.CONSTRAINTS), + /* Combined: must be declared after any dependent models above. */ /** 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 new file mode 100644 index 0000000000000000000000000000000000000000..8220873d504c602ee74b22857708d389418774de --- /dev/null +++ b/src/main/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018.java @@ -0,0 +1,213 @@ +/** + * + */ +package gov.usgs.earthquake.nshmp.gmm; + +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.GmmUtils.BASE_10_TO_E; +import static java.lang.Math.exp; +import static java.lang.Math.log10; +import static java.lang.Math.max; +import static java.lang.Math.min; +import static java.lang.Math.sqrt; + +import java.util.Map; + +import com.google.common.annotations.Beta; +import com.google.common.collect.Range; + +import gov.usgs.earthquake.nshmp.gmm.GmmInput.Constraints; +import gov.usgs.earthquake.nshmp.gmm.NgaEastUsgs_2017.SiteAmp; +import gov.usgs.earthquake.nshmp.tree.LogicTree; + +/** + * Implementation of the Pezeshk et al. (2018) ground motion model for central + * and eastern North America (CENA). Pezeshk et al. (2018) provide coefficients + * for two hybrid empirical method (HEM) ground-motion prediction equations, one + * based on a stochastic-scaling approach (PZCT-1SS) and one using the + * empirical-scaling approach (PZCT-2ES). The latter model is preferred by the + * authors as the empirical-scaling approach reduces oversaturation of GMIM + * predictions at large magnitudes, short distances, and short periods. + * + * 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: 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 + * desired {@link Imt}. + * + * <p><b>Reference:</b> Pezeshk, S., Zandieh, A., Campbell, K.W., and Tavakoli, + * B., 2018, Ground-motion prediction equations for Central and Eastern North + * America using the hybrid empirical method and NGA-West2 Empirical + * ground-motion models: Bulletin of the Seismological Society of America, v. + * 108, n. 4, p. 2278-2304, August 2018. + * + * + * <p><b>doi:</b> <a href="http://dx.doi.org/10.1785/0120170179"> + * 10.1785/0120170179</a> + * + * <p><b>Component:</b> **** RotD50 ? + * + * @author U.S. Geological Survey + * @see Gmm#PZCT18_1SS + * @see Gmm#PZCT18_2ES + */ +@Beta +public abstract class PezeshkEtAl_2018 implements GroundMotionModel { + + static final String NAME = "Pezeshk et al. (2018)"; + + static final Constraints CONSTRAINTS = Constraints.builder() + .set(MW, Range.closed(4.0, 8.0)) + .set(RRUP, Range.closed(0.0, 1000.0)) + .set(VS30, Range.closed(0.0, 3000.0)) // ? + .build(); + + static final class Model { + static final String PZCT18_1SS = "PZCT18-1SS"; + static final String PZCT18_2ES = "PZCT18-2ES"; + } + + /* + * Note: 0.04 and 0.08 second coefficents are commented out because NGA-East + * SiteAmp doesn't support SA0P04 (also commented out there) or SA0P08. + */ + static final CoefficientContainer COEFFS_1SS = new CoefficientContainer("PZCT18_1SS.csv"); + static final CoefficientContainer COEFFS_2ES = new CoefficientContainer("PZCT18_2ES.csv"); + static final CoefficientContainer COEFFS_SIGMAS = new CoefficientContainer("PZCT18_sigmas.csv"); + + private static final class Coefficients { + final Imt imt; + final double c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11, σReg; + final double c12, c13, c14, c15, c16, c17, c18, c19, c20, c21, c22, c23, c24, c25; + final String model; + + Coefficients(Imt imt, CoefficientContainer cc, CoefficientContainer cc_sigma, String model) { + this.imt = imt; + this.model = model; + + Map<String, Double> coeffs = cc.get(imt); + Map<String, Double> coeffs_sigma = cc_sigma.get(imt); + + /* + * 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. + */ + + // median model coefficients + c1 = coeffs.get("c1"); + c2 = coeffs.get("c2"); + c3 = coeffs.get("c3"); + c4 = coeffs.get("c4"); + c5 = coeffs.get("c5"); + c6 = coeffs.get("c6"); + c7 = coeffs.get("c7"); + c8 = coeffs.get("c8"); + c9 = coeffs.get("c9"); + c10 = coeffs.get("c10"); + c11 = coeffs.get("c11"); + σReg = coeffs.get("σReg"); + + // between-event variability coefficients + c12 = coeffs_sigma.get("c12"); + c13 = coeffs_sigma.get("c13"); + c14 = coeffs_sigma.get("c14"); + c15 = coeffs_sigma.get("c15"); + c16 = coeffs_sigma.get("c16"); + c17 = coeffs_sigma.get("c17"); + c18 = coeffs_sigma.get("c18"); + // within-event variability coefficients + c19 = coeffs_sigma.get("c19"); + c20 = coeffs_sigma.get("c20"); + c21 = coeffs_sigma.get("c21"); + c22 = coeffs_sigma.get("c22"); + c23 = coeffs_sigma.get("c23"); + c24 = coeffs_sigma.get("c24"); + c25 = coeffs_sigma.get("c25"); + } + + } + + private final Coefficients coeffs; + private final Coefficients coeffsPga; + private final SiteAmp siteAmp; + + PezeshkEtAl_2018(Imt imt, String model) { + CoefficientContainer cc = (model == Model.PZCT18_1SS) ? COEFFS_1SS : COEFFS_2ES; + coeffs = new Coefficients(imt, cc, COEFFS_SIGMAS, model); + coeffsPga = new Coefficients(Imt.PGA, cc, COEFFS_SIGMAS, model); + this.siteAmp = new NgaEastUsgs_2017.SiteAmp(imt); + } + + @Override + public LogicTree<GroundMotion> calc(GmmInput in) { + double μPga = exp(calcMean(coeffsPga, in.Mw, in.rRup)); + double μ = calcMean(coeffs, in.Mw, in.rRup); + double σ = calcStdDev(coeffs, in.Mw); + + SiteAmp.Value fSite = siteAmp.calc(μPga, in.vs30); + μ = fSite.apply(μ); + + return GroundMotions.createTree(μ, σ); + } + + private double calcMean(Coefficients c, double Mw, double Rrup) { + /* Eq. 5 */ + double R = sqrt(Rrup * Rrup + c.c11 * c.c11); + + /* Eq. 4 - log10 units */ + double μ = c.c1 + + (c.c2 * Mw) + + (c.c3 * Mw * Mw) + + (c.c4 + c.c5 * Mw) * min(log10(R), log10(60.0)) + + (c.c6 + c.c7 * Mw) * max(min(log10(R / 60.0), log10(2.0)), 0.0) + + (c.c8 + c.c9 * Mw) * max(log10(R / 120.0), 0.0) + + (c.c10 * R); + return μ * BASE_10_TO_E; + } + + private double calcStdDev(Coefficients c, double Mw) { + /* Aleatory standard deviation */ + /* Eq. 6 (between-event term) and Eq. 7 (within-event term) - natural log */ + double τ, φ; + if (Mw <= 4.5) { + τ = c.c12; + φ = c.c19 + c.c20 * Mw; + } else if (Mw <= 5.0) { + τ = c.c13 + c.c14 * Mw; + φ = c.c21 + c.c22 * Mw; + } else if (Mw <= 6.5) { + τ = c.c15 + c.c16 * Mw; + φ = c.c23 + c.c24 * Mw; + } else { + τ = c.c17 + c.c18 * Mw; + φ = c.c25; + } + return sqrt(τ * τ + φ * φ); + } + + static final class PZCT18_1SS extends PezeshkEtAl_2018 { + static final String MODEL = Model.PZCT18_1SS; + static final String NAME = PezeshkEtAl_2018.NAME + " : " + MODEL; + + PZCT18_1SS(Imt imt) { + super(imt, MODEL); + } + } + + static final class PZCT18_2ES extends PezeshkEtAl_2018 { + static final String MODEL = Model.PZCT18_2ES; + static final String NAME = PezeshkEtAl_2018.NAME + " : " + MODEL; + + PZCT18_2ES(Imt imt) { + super(imt, MODEL); + } + } + +} diff --git a/src/main/resources/gmm/coeffs/PZCT18_1SS.csv b/src/main/resources/gmm/coeffs/PZCT18_1SS.csv new file mode 100644 index 0000000000000000000000000000000000000000..150322ddd8ee009469978486eeac2a7a78f16929 --- /dev/null +++ b/src/main/resources/gmm/coeffs/PZCT18_1SS.csv @@ -0,0 +1,24 @@ + T, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11, σReg + PGA, -2.043e-1, 7.805e-1, -8.063e-2, -3.956, 3.901e-1, -3.259e-1, 8.243e-2, -1.928e+0, 2.511e-1, -2.347e-3, -6.368, 6.757e-2 + 0.01, 4.265e-1, 7.075e-1, -7.832e-2, -4.184, 4.091e-1, -6.455e-1, 1.049e-1, -2.056e+0, 2.585e-1, -2.286e-3, 6.653, 6.683e-2 + 0.02, 6.573e-1, 6.803e-1, -7.444e-2, -4.088, 3.909e-1, -6.327e-1, 5.930e-2, -2.645e+0, 2.934e-1, -2.045e-3, 6.723, 7.019e-2 + 0.03, 3.420e-1, 7.205e-1, -7.437e-2, -3.896, 3.729e-1, 6.628e-3, -1.971e-2, -2.848e+0, 3.026e-1, -2.068e-3, 6.523, 7.611e-2 +#0.04, -1.266e-2, 7.652e-1, -7.558e-2, -3.753, 3.644e-1, 5.479e-1, -6.757e-2, -2.577e+0, 2.711e-1, -2.376e-3, -6.332, 8.139e-2 + 0.05, -2.608e-1, 7.963e-1, -7.698e-2, -3.683, 3.639e-1, 7.904e-1, -7.123e-2, -2.094e+0, 2.232e-1, -2.747e-3, -6.190, 8.635e-2 +#0.08, -7.995e-1, 9.087e-1, -8.460e-2, -3.587, 3.591e-1, 6.022e-1, 3.520e-3, -1.075e+0, 1.398e-1, -3.319e-3, 6.105, 9.263e-2 + 0.10, -1.391e+0, 1.060e+0, -9.486e-2, -3.475, 3.439e-1, 2.369e-1, 6.864e-2, -4.517e-1, 8.967e-2, -3.481e-3, 6.130, 9.041e-2 + 0.15, -2.425e+0, 1.322e+0, -1.117e-1, -3.322, 3.177e-1, -3.225e-1, 1.450e-1, 1.825e-1, 8.081e-3, -3.216e-3, 6.229, 8.438e-2 + 0.20, -3.243e+0, 1.537e+0, -1.262e-1, -3.278, 3.062e-1, -5.643e-1, 1.689e-1, 1.641e-1, -4.100e-3, -2.638e-3, 6.366, 8.335e-2 + 0.25, -3.858e+0, 1.685e+0, -1.361e-1, -3.265, 3.024e-1, -6.514e-1, 1.742e-1, 5.760e-2, -3.273e-3, -2.137e-3, 6.331, 8.384e-2 + 0.30, -4.240e+0, 1.772e+0, -1.419e-1, -3.291, 3.053e-1, -6.865e-1, 1.736e-1, -6.002e-2, 2.221e-3, -1.744e-3, 6.318, 8.163e-2 + 0.40, -5.013e+0, 1.932e+0, -1.501e-1, -3.235, 2.919e-1, -6.095e-1, 1.481e-1, -3.715e-1, 2.499e-2, -1.072e-3, 6.493, 7.561e-2 + 0.50, -5.441e+0, 2.005e+0, -1.532e-1, -3.232, 2.888e-1, -5.287e-1, 1.262e-1, -6.771e-1, 5.550e-2, -6.642e-4, 6.558, 6.897e-2 + 0.75, -6.085e+0, 2.070e+0, -1.533e-1, -3.201, 2.863e-1, -4.266e-1, 1.014e-1, -9.856e-1, 8.410e-2, -3.040e-4, 6.303, 6.329e-2 + 1.00, -6.446e+0, 2.077e+0, -1.498e-1, -3.165, 2.814e-1, -3.605e-1, 8.293e-2, -1.113e+0, 9.241e-2, -1.857e-4, 6.187, 6.259e-2 + 1.50, -6.727e+0, 1.992e+0, -1.381e-1, -3.132, 2.848e-1, -2.712e-1, 6.172e-2, -1.272e+0, 1.104e-1, -4.798e-5, 5.818, 6.435e-2 + 2.00, -6.801e+0, 1.907e+0, -1.289e-1, -3.162, 2.935e-1, -2.309e-1, 4.892e-2, -1.438e+0, 1.316e-1, -3.612e-6, 5.655, 6.504e-2 + 3.00, -6.740e+0, 1.729e+0, -1.104e-1, -3.199, 3.019e-1, -1.443e-1, 2.933e-2, -1.619e+0, 1.535e-1, 1.309e-5, -5.447, 6.739e-2 + 4.00, -6.638e+0, 1.589e+0, -9.647e-2, -3.220, 3.074e-1, -1.276e-1, 2.714e-2, -1.575e+0, 1.411e-1, -4.202e-5, -5.592, 6.594e-2 + 5.00, -6.726e+0, 1.536e+0, -8.978e-2, -3.231, 3.089e-1, -1.063e-1, 2.324e-2, -1.572e+0, 1.354e-1, -9.189e-5, 5.726, 6.501e-2 + 7.50, -7.194e+0, 1.514e+0, -8.441e-2, -3.236, 3.170e-1, -1.766e-2, 1.731e-2, -1.658e+0, 1.494e-1, -1.369e-4, 5.846, 6.265e-2 +10.00, -7.507e+0, 1.522e+0, -8.393e-2, -3.290, 3.257e-1, 1.807e-1, -3.022e-3, -1.615e+0, 1.461e-1, -1.939e-4, 5.701, 6.838e-2 diff --git a/src/main/resources/gmm/coeffs/PZCT18_2ES.csv b/src/main/resources/gmm/coeffs/PZCT18_2ES.csv new file mode 100644 index 0000000000000000000000000000000000000000..4832aa30a3c54dc140affe9f85fb6b05c08fd1fd --- /dev/null +++ b/src/main/resources/gmm/coeffs/PZCT18_2ES.csv @@ -0,0 +1,24 @@ + T, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11, σReg + PGA, -7.655e-1, 8.994e-1, -7.874e-2, -3.534, 2.948e-1, -4.799e-1, 1.178e-1, -2.182e+0, 3.029e-1, -2.327e-3, 6.461, 5.836e-2 + 0.01, -1.871e-1, 8.314e-1, -7.503e-2, -3.714, 3.030e-1, -7.793e-1, 1.360e-1, -2.307e+0, 3.093e-1, -2.263e-3, 6.745, 5.637e-2 + 0.02, 8.387e-2, 7.934e-1, -7.073e-2, -3.641, 2.898e-1, -7.692e-1, 9.124e-2, -2.885e+0, 3.399e-1, -2.013e-3, 6.814, 5.774e-2 + 0.03, -1.786e-1, 8.256e-1, -7.148e-2, -3.489, 2.809e-1, -1.806e-1, 2.319e-2, -3.072e+0, 3.469e-1, -2.043e-3, 6.616, 6.483e-2 +#0.04, -4.947e-1, 8.659e-1, -7.358e-2, -3.378, 2.796e-1, 3.218e-1, -1.611e-2, -2.791e+0, 3.157e-1, -2.365e-3, 6.425, 7.330e-2 + 0.05, -7.246e-1, 8.959e-1, -7.568e-2, -3.325, 2.831e-1, 5.318e-1, -1.344e-2, -2.301e+0, 2.693e-1, -2.751e-3, 6.278, 8.156e-2 +#0.08, -1.255e+0, 1.009e+0, -8.397e-2, -3.237, 2.809e-1, 3.432e-1, 6.065e-2, -1.289e+0, 1.932e-1, -3.352e-3, 6.184, 9.085e-2 + 0.10, -1.849e+0, 1.160e+0, -9.436e-2, -3.125, 2.660e-1, 1.611e-2, 1.168e-1, -6.970e-1, 1.511e-1, -3.519e-3, -6.205, 8.908e-2 + 0.15, -2.862e+0, 1.418e+0, -1.110e-1, -2.985, 2.416e-1, -4.494e-1, 1.736e-1, -1.474e-1, 8.225e-2, -3.224e-3, -6.322, 8.524e-2 + 0.20, -3.666e+0, 1.632e+0, -1.256e-1, -2.952, 2.313e-1, -6.274e-1, 1.847e-1, -2.470e-1, 7.741e-2, -2.591e-3, 6.479, 8.309e-2 + 0.25, -4.260e+0, 1.777e+0, -1.355e-1, -2.955, 2.298e-1, -6.681e-1, 1.811e-1, -4.060e-1, 8.046e-2, -2.041e-3, 6.470, 8.142e-2 + 0.30, -4.630e+0, 1.864e+0, -1.416e-1, -2.995, 2.347e-1, -6.775e-1, 1.760e-1, -5.479e-1, 8.473e-2, -1.614e-3, 6.475, 7.692e-2 + 0.40, -5.407e+0, 2.029e+0, -1.506e-1, -2.950, 2.231e-1, -5.911e-1, 1.492e-1, -8.400e-1, 9.921e-2, -9.231e-4, 6.667, 6.876e-2 + 0.50, -5.843e+0, 2.108e+0, -1.545e-1, -2.956, 2.222e-1, -5.227e-1, 1.298e-1, -1.086e+0, 1.194e-1, -5.303e-4, 6.733, 6.338e-2 + 0.75, -6.499e+0, 2.185e+0, -1.570e-1, -2.948, 2.259e-1, -4.594e-1, 1.120e-1, -1.252e+0, 1.284e-1, -2.316e-4, 6.466, 6.348e-2 + 1.00, -6.869e+0, 2.204e+0, -1.557e-1, -2.933, 2.261e-1, -4.044e-1, 9.553e-2, -1.297e+0, 1.260e-1, -1.516e-4, 6.356, 6.575e-2 + 1.50, -7.172e+0, 2.143e+0, -1.484e-1, -2.940, 2.387e-1, -3.053e-1, 7.187e-2, -1.386e+0, 1.337e-1, -4.130e-5, 5.986, 6.896e-2 + 2.00, -7.287e+0, 2.085e+0, -1.431e-1, -2.994, 2.533e-1, -2.546e-1, 5.605e-2, -1.519e+0, 1.494e-1, -5.530e-6, -5.806, 6.909e-2 + 3.00, -7.345e+0, 1.957e+0, -1.300e-1, -3.039, 2.649e-1, -1.506e-1, 3.184e-2, -1.668e+0, 1.648e-1, 7.650e-6, -5.528, 7.060e-2 + 4.00, -7.387e+0, 1.860e+0, -1.191e-1, -3.023, 2.637e-1, -1.514e-1, 3.288e-2, -1.603e+0, 1.473e-1, -4.230e-5, -5.596, 7.011e-2 + 5.00, -7.607e+0, 1.839e+0, -1.135e-1, -2.980, 2.545e-1, -1.600e-1, 3.506e-2, -1.583e+0, 1.374e-1, -9.030e-5, -5.663, 7.060e-2 + 7.50, -8.251e+0, 1.851e+0, -1.072e-1, -2.879, 2.403e-1, -1.687e-1, 5.017e-2, -1.643e+0, 1.455e-1, -1.330e-4, -5.721, 6.987e-2 +10.00, -8.581e+0, 1.850e+0, -1.039e-1, -2.876, 2.362e-1, -8.466e-2, 5.500e-2, -1.599e+0, 1.450e-1, -2.047e-4, 5.604, 7.538e-2 diff --git a/src/main/resources/gmm/coeffs/PZCT18_sigmas.csv b/src/main/resources/gmm/coeffs/PZCT18_sigmas.csv new file mode 100644 index 0000000000000000000000000000000000000000..774afda39ba38f6652ff0f08206579a280620597 --- /dev/null +++ b/src/main/resources/gmm/coeffs/PZCT18_sigmas.csv @@ -0,0 +1,24 @@ + T, c12, c13, c14, c15, c16, c17, c18, c19, c20, c21, c22, c23, c24, c25 + PGA, 4.191e-1, 7.699e-1, -7.798e-2, 5.518e-1, -3.435e-2, 3.596e-1, -4.792e-3, 8.376e-1, -2.941e-2, 1.880e+0, -2.610e-1, 7.848e-1, -4.203e-2, 5.116e-1 + 0.01, 4.188e-1, 7.505e-1, -7.373e-2, 5.599e-1, -3.560e-2, 3.596e-1, -4.792e-3, 8.379e-1, -2.941e-2, 1.879e+0, -2.607e-1, 7.849e-1, -4.194e-2, 5.122e-1 + 0.02, 4.245e-1, 8.034e-1, -8.422e-2, 5.569e-1, -3.493e-2, 3.610e-1, -4.790e-3, 8.437e-1, -3.016e-2, 1.886e+0, -2.617e-1, 7.889e-1, -4.235e-2, 5.135e-1 + 0.03, 4.416e-1, 8.784e-1, -9.709e-2, 5.701e-1, -3.543e-2, 3.716e-1, -4.885e-3, 8.685e-1, -3.285e-2, 1.925e+0, -2.677e-1, 8.097e-1, -4.458e-2, 5.198e-1 +#0.04, 4.571e-1, 9.426e-1, -1.079e-1, 5.736e-1, -3.412e-2, 3.853e-1, -5.163e-3, 8.865e-1, -3.425e-2, 1.974e+0, -2.759e-1, 8.224e-1, -4.559e-2, 5.259e-1 + 0.05, 4.725e-1, 1.007e+0, -1.187e-1, 5.772e-1, -3.286e-2, 3.990e-1, -5.444e-3, 9.042e-1, -3.566e-2, 2.022e+0, -2.841e-1, 8.346e-1, -4.660e-2, 5.314e-1 +#0.08, 4.653e-1, 7.718e-1, -6.819e-2, 5.843e-1, -3.068e-2, 4.219e-1, -5.694e-3, 9.021e-1, -3.464e-2, 1.967e+0, -2.713e-1, 8.363e-1, -4.514e-2, 5.425e-1 + 0.10, 4.369e-1, 4.594e-1, -5.078e-3, 5.920e-1, -3.159e-2, 4.231e-1, -5.605e-3, 8.940e-1, -3.336e-2, 1.870e+0, -2.503e-1, 8.416e-1, -4.458e-2, 5.515e-1 + 0.15, 4.080e-1, 3.551e-1, 1.173e-2, 5.802e-1, -3.329e-2, 3.962e-1, -4.991e-3, 8.847e-1, -3.049e-2, 1.855e+0, -2.460e-1, 8.345e-1, -4.200e-2, 5.613e-1 + 0.20, 3.959e-1, 4.398e-1, -9.744e-3, 5.797e-1, -3.773e-2, 3.650e-1, -4.710e-3, 8.709e-1, -2.816e-2, 1.770e+0, -2.280e-1, 8.313e-1, -4.022e-2, 5.698e-1 + 0.25, 3.984e-1, 5.748e-1, -3.920e-2, 5.773e-1, -3.970e-2, 3.497e-1, -4.690e-3, 8.477e-1, -2.513e-2, 1.689e+0, -2.122e-1, 8.126e-1, -3.683e-2, 5.732e-1 + 0.30, 4.023e-1, 6.887e-1, -6.367e-2, 5.822e-1, -4.236e-2, 3.372e-1, -4.674e-3, 8.101e-1, -2.010e-2, 1.475e+0, -1.679e-1, 8.021e-1, -3.327e-2, 5.858e-1 + 0.40, 4.116e-1, 9.150e-1, -1.119e-1, 5.601e-1, -4.089e-2, 3.246e-1, -4.652e-3, 7.563e-1, -1.369e-2, 1.155e+0, -1.022e-1, 7.894e-1, -2.917e-2, 5.998e-1 + 0.50, 4.252e-1, 1.026e+0, -1.336e-1, 5.582e-1, -3.995e-2, 3.286e-1, -4.637e-3, 7.192e-1, -9.542e-3, 9.129e-1, -5.258e-2, 7.771e-1, -2.542e-2, 6.118e-1 + 0.75, 4.507e-1, 1.017e+0, -1.258e-1, 5.933e-1, -4.109e-2, 3.562e-1, -4.613e-3, 6.487e-1, -2.016e-3, 5.115e-1, 2.847e-2, 7.394e-1, -1.711e-2, 6.281e-1 + 1.00, 4.716e-1, 1.092e+0, -1.378e-1, 6.055e-1, -4.057e-2, 3.718e-1, -4.602e-3, 6.028e-1, 2.617e-3, 3.564e-1, 5.736e-2, 6.813e-1, -7.617e-3, 6.318e-1 + 1.50, 4.859e-1, 1.053e+0, -1.260e-1, 6.367e-1, -4.277e-2, 3.886e-1, -4.601e-3, 5.479e-1, 7.745e-3, 2.408e-1, 7.599e-2, 5.958e-1, 4.986e-3, 6.282e-1 + 2.00, 4.886e-1, 1.051e+0, -1.250e-1, 6.359e-1, -4.197e-2, 3.932e-1, -4.625e-3, 5.176e-1, 1.187e-2, 2.150e-1, 7.910e-2, 5.574e-1, 1.063e-2, 6.264e-1 + 3.00, 4.907e-1, 9.880e-1, -1.105e-1, 6.469e-1, -4.229e-2, 4.021e-1, -4.627e-3, 5.156e-1, 1.149e-2, 2.479e-1, 7.098e-2, 5.505e-1, 1.046e-2, 6.185e-1 + 4.00, 5.033e-1, 1.228e+0, -1.610e-1, 6.062e-1, -3.667e-2, 3.977e-1, -4.590e-3, 5.005e-1, 1.287e-2, 2.822e-1, 6.137e-2, 5.287e-1, 1.208e-2, 6.072e-1 + 5.00, 4.986e-1, 1.102e+0, -1.340e-1, 6.311e-1, -3.990e-2, 4.015e-1, -4.583e-3, 4.748e-1, 1.563e-2, 2.037e-1, 7.588e-2, 5.122e-1, 1.417e-2, 6.043e-1 + 7.50, 4.910e-1, 1.049e+0, -1.240e-1, 6.579e-1, -4.578e-2, 3.901e-1, -4.583e-3, 4.211e-1, 2.163e-2, -1.653e-2, 1.189e-1, 4.920e-1, 1.717e-2, 6.036e-1 +10.00, 4.711e-1, 8.445e-1, -8.300e-2, 6.848e-1, -5.105e-2, 3.828e-1, -4.583e-3, 3.946e-1, 2.512e-2, -5.540e-2, 1.251e-1, 4.898e-1, 1.609e-2, 5.944e-1 diff --git a/src/test/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018Tests.java b/src/test/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018Tests.java new file mode 100644 index 0000000000000000000000000000000000000000..8476d32fb23ec3ce8f695c97d1fd34d4231ac25f --- /dev/null +++ b/src/test/java/gov/usgs/earthquake/nshmp/gmm/PezeshkEtAl_2018Tests.java @@ -0,0 +1,43 @@ +package gov.usgs.earthquake.nshmp.gmm; + +import static gov.usgs.earthquake.nshmp.gmm.Gmm.PZCT18_1SS; +import static gov.usgs.earthquake.nshmp.gmm.Gmm.PZCT18_2ES; + +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; + +public class PezeshkEtAl_2018Tests implements ArgumentsProvider { + + private static String GMM_INPUTS = "/gmm/pzct18_inputs.csv"; + private static String GMM_RESULTS = "/gmm/pzct18_results.csv"; + + @ParameterizedTest(name = "({index}) {0} {2} {1}") + @ArgumentsSource(PezeshkEtAl_2018Tests.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( + PZCT18_1SS, + PZCT18_2ES); + + private static Set<Imt> imts = PezeshkEtAl_2018.COEFFS_1SS.imts(); + + public static void main(String[] args) throws IOException { + GmmTest.generateResults(gmms, imts, GMM_INPUTS, GMM_RESULTS); + } +} diff --git a/src/test/resources/gmm/pzct18_inputs.csv b/src/test/resources/gmm/pzct18_inputs.csv new file mode 100644 index 0000000000000000000000000000000000000000..b046d6b8e1b3d0ec5f0e1c3b05687e4a1cddb7df --- /dev/null +++ b/src/test/resources/gmm/pzct18_inputs.csv @@ -0,0 +1,17 @@ +#Mw,rjB,rRup,rX,dip,width,zTor,zHyp,rake,vs30,vsInf,z1p0,z2p5 +5.0,0,5.0,0,90,0,0,0,0,3000,0,NaN,NaN +6.0,0,5.0,0,90,0,0,0,0,3000,0,NaN,NaN +7.0,0,5.0,0,90,0,0,0,0,3000,0,NaN,NaN +8.0,0,5.0,0,90,0,0,0,0,3000,0,NaN,NaN +5.0,0,10.0,0,90,0,0,0,0,3000,0,NaN,NaN +6.0,0,10.0,0,90,0,0,0,0,3000,0,NaN,NaN +7.0,0,10.0,0,90,0,0,0,0,3000,0,NaN,NaN +8.0,0,10.0,0,90,0,0,0,0,3000,0,NaN,NaN +5.0,0,30.0,0,90,0,0,0,0,3000,0,NaN,NaN +6.0,0,30.0,0,90,0,0,0,0,3000,0,NaN,NaN +7.0,0,30.0,0,90,0,0,0,0,3000,0,NaN,NaN +8.0,0,30.0,0,90,0,0,0,0,3000,0,NaN,NaN +5.0,0,70.0,0,90,0,0,0,0,3000,0,NaN,NaN +6.0,0,70.0,0,90,0,0,0,0,3000,0,NaN,NaN +7.0,0,70.0,0,90,0,0,0,0,3000,0,NaN,NaN +8.0,0,70.0,0,90,0,0,0,0,3000,0,NaN,NaN diff --git a/src/test/resources/gmm/pzct18_results.csv b/src/test/resources/gmm/pzct18_results.csv new file mode 100644 index 0000000000000000000000000000000000000000..978fafeca5a1cc80cb37262a6b5e9cbc6cf2892d --- /dev/null +++ b/src/test/resources/gmm/pzct18_results.csv @@ -0,0 +1,672 @@ +0-PZCT18_1SS-PGA,0.6948163938,0.6892205743 +1-PZCT18_1SS-PGA,1.2296015484,0.6349744518 +2-PZCT18_1SS-PGA,1.5010623981,0.6066688348 +3-PZCT18_1SS-PGA,1.2640757226,0.6041068761 +4-PZCT18_1SS-PGA,0.3168717254,0.6892205743 +5-PZCT18_1SS-PGA,0.6507119077,0.6349744518 +6-PZCT18_1SS-PGA,0.9217943030,0.6066688348 +7-PZCT18_1SS-PGA,0.9007811413,0.6041068761 +8-PZCT18_1SS-PGA,0.0425512079,0.6892205743 +9-PZCT18_1SS-PGA,0.1266018809,0.6349744518 +10-PZCT18_1SS-PGA,0.2598414454,0.6066688348 +11-PZCT18_1SS-PGA,0.3678889037,0.6041068761 +12-PZCT18_1SS-PGA,0.0090641644,0.6892205743 +13-PZCT18_1SS-PGA,0.0354993598,0.6349744518 +14-PZCT18_1SS-PGA,0.0959076547,0.6066688348 +15-PZCT18_1SS-PGA,0.1787417560,0.6041068761 +0-PZCT18_1SS-SA0P01,1.0447753885,0.6906588684 +1-PZCT18_1SS-SA0P01,1.7437198685,0.6358379649 +2-PZCT18_1SS-SA0P01,2.0290394858,0.6071748967 +3-PZCT18_1SS-SA0P01,1.6461305738,0.6046150823 +4-PZCT18_1SS-SA0P01,0.4675935147,0.6906588684 +5-PZCT18_1SS-SA0P01,0.9067827756,0.6358379649 +6-PZCT18_1SS-SA0P01,1.2260213839,0.6071748967 +7-PZCT18_1SS-SA0P01,1.1557205300,0.6046150823 +8-PZCT18_1SS-SA0P01,0.0568376314,0.6906588684 +9-PZCT18_1SS-SA0P01,0.1618717953,0.6358379649 +10-PZCT18_1SS-SA0P01,0.3214151885,0.6071748967 +11-PZCT18_1SS-SA0P01,0.4449606516,0.6046150823 +12-PZCT18_1SS-SA0P01,0.0108229485,0.6906588684 +13-PZCT18_1SS-SA0P01,0.0412092612,0.6358379649 +14-PZCT18_1SS-SA0P01,0.1093966601,0.6071748967 +15-PZCT18_1SS-SA0P01,0.2024758812,0.6046150823 +0-PZCT18_1SS-SA0P02,1.6259017789,0.6925745736 +1-PZCT18_1SS-SA0P02,2.7128098767,0.6376850495 +2-PZCT18_1SS-SA0P02,3.2126573637,0.6090310755 +3-PZCT18_1SS-SA0P02,2.7004084691,0.6064689872 +4-PZCT18_1SS-SA0P02,0.7360108851,0.6925745736 +5-PZCT18_1SS-SA0P02,1.4154673935,0.6376850495 +6-PZCT18_1SS-SA0P02,1.9321262532,0.6090310755 +7-PZCT18_1SS-SA0P02,1.8719365918,0.6064689872 +8-PZCT18_1SS-SA0P02,0.0913669336,0.6925745736 +9-PZCT18_1SS-SA0P02,0.2534049214,0.6376850495 +10-PZCT18_1SS-SA0P02,0.4988397358,0.6090310755 +11-PZCT18_1SS-SA0P02,0.6969907848,0.6064689872 +12-PZCT18_1SS-SA0P02,0.0172634567,0.6925745736 +13-PZCT18_1SS-SA0P02,0.0627704531,0.6376850495 +14-PZCT18_1SS-SA0P02,0.1619954312,0.6090310755 +15-PZCT18_1SS-SA0P02,0.2967360212,0.6064689872 +0-PZCT18_1SS-SA0P03,1.6211176610,0.7059688042 +1-PZCT18_1SS-SA0P03,2.8403876763,0.6494790827 +2-PZCT18_1SS-SA0P03,3.5334642544,0.6197049088 +3-PZCT18_1SS-SA0P03,3.1209282527,0.6170588225 +4-PZCT18_1SS-SA0P03,0.7458792934,0.7059688042 +5-PZCT18_1SS-SA0P03,1.5021229161,0.6494790827 +6-PZCT18_1SS-SA0P03,2.1478423069,0.6197049088 +7-PZCT18_1SS-SA0P03,2.1805161796,0.6170588225 +8-PZCT18_1SS-SA0P03,0.1001407693,0.7059688042 +9-PZCT18_1SS-SA0P03,0.2868137303,0.6494790827 +10-PZCT18_1SS-SA0P03,0.5832422175,0.6197049088 +11-PZCT18_1SS-SA0P03,0.8420889470,0.6170588225 +12-PZCT18_1SS-SA0P03,0.0209518373,0.7059688042 +13-PZCT18_1SS-SA0P03,0.0768012320,0.6494790827 +14-PZCT18_1SS-SA0P03,0.1998822650,0.6197049088 +15-PZCT18_1SS-SA0P03,0.3693518997,0.6170588225 +0-PZCT18_1SS-SA0P05,1.2467544221,0.7299208861 +1-PZCT18_1SS-SA0P05,2.3609479793,0.6726480518 +2-PZCT18_1SS-SA0P05,3.1364034202,0.6423620441 +3-PZCT18_1SS-SA0P05,2.9229233018,0.6393193574 +4-PZCT18_1SS-SA0P05,0.5876719318,0.7299208861 +5-PZCT18_1SS-SA0P05,1.2828820918,0.6726480518 +6-PZCT18_1SS-SA0P05,1.9646198912,0.6423620441 +7-PZCT18_1SS-SA0P05,2.1106213696,0.6393193574 +8-PZCT18_1SS-SA0P05,0.0876128298,0.7299208861 +9-PZCT18_1SS-SA0P05,0.2709609079,0.6726480518 +10-PZCT18_1SS-SA0P05,0.5878757826,0.6423620441 +11-PZCT18_1SS-SA0P05,0.8947557563,0.6393193574 +12-PZCT18_1SS-SA0P05,0.0208632003,0.7299208861 +13-PZCT18_1SS-SA0P05,0.0814856377,0.6726480518 +14-PZCT18_1SS-SA0P05,0.2232652899,0.6423620441 +15-PZCT18_1SS-SA0P05,0.4291423272,0.6393193574 +0-PZCT18_1SS-SA0P1,0.8573554757,0.7555838339 +1-PZCT18_1SS-SA0P1,1.8138205621,0.7011332441 +2-PZCT18_1SS-SA0P1,2.4791780374,0.6719409113 +3-PZCT18_1SS-SA0P1,2.1892789984,0.6687547216 +4-PZCT18_1SS-SA0P1,0.4164490873,0.7555838339 +5-PZCT18_1SS-SA0P1,1.0088458611,0.7011332441 +6-PZCT18_1SS-SA0P1,1.5789476043,0.6719409113 +7-PZCT18_1SS-SA0P1,1.5965797649,0.6687547216 +8-PZCT18_1SS-SA0P1,0.0664099497,0.7555838339 +9-PZCT18_1SS-SA0P1,0.2237740578,0.7011332441 +10-PZCT18_1SS-SA0P1,0.4871540567,0.6719409113 +11-PZCT18_1SS-SA0P1,0.6851772107,0.6687547216 +12-PZCT18_1SS-SA0P1,0.0162611442,0.7555838339 +13-PZCT18_1SS-SA0P1,0.0698081951,0.7011332441 +14-PZCT18_1SS-SA0P1,0.1936162118,0.6719409113 +15-PZCT18_1SS-SA0P1,0.3469421200,0.6687547216 +0-PZCT18_1SS-SA0P15,0.6345724562,0.7495425688 +1-PZCT18_1SS-SA0P15,1.5222201242,0.6957413755 +2-PZCT18_1SS-SA0P15,2.1831003578,0.6675092847 +3-PZCT18_1SS-SA0P15,1.8718452534,0.6648213504 +4-PZCT18_1SS-SA0P15,0.3145530357,0.7495425688 +5-PZCT18_1SS-SA0P15,0.8537145175,0.6957413755 +6-PZCT18_1SS-SA0P15,1.3852605060,0.6675092847 +7-PZCT18_1SS-SA0P15,1.3438482759,0.6648213504 +8-PZCT18_1SS-SA0P15,0.0521776297,0.7495425688 +9-PZCT18_1SS-SA0P15,0.1918569836,0.6957413755 +10-PZCT18_1SS-SA0P15,0.4217653316,0.6675092847 +11-PZCT18_1SS-SA0P15,0.5543245906,0.6648213504 +12-PZCT18_1SS-SA0P15,0.0129330426,0.7495425688 +13-PZCT18_1SS-SA0P15,0.0602386161,0.6957413755 +14-PZCT18_1SS-SA0P15,0.1677448764,0.6675092847 +15-PZCT18_1SS-SA0P15,0.2792695222,0.6648213504 +0-PZCT18_1SS-SA0P2,0.4775389663,0.7415143737 +1-PZCT18_1SS-SA0P2,1.2761499309,0.6876855552 +2-PZCT18_1SS-SA0P2,1.9071927521,0.6594815850 +3-PZCT18_1SS-SA0P2,1.5939969054,0.6571228366 +4-PZCT18_1SS-SA0P2,0.2396971341,0.7415143737 +5-PZCT18_1SS-SA0P2,0.7199202160,0.6876855552 +6-PZCT18_1SS-SA0P2,1.2092216012,0.6594815850 +7-PZCT18_1SS-SA0P2,1.1358668482,0.6571228366 +8-PZCT18_1SS-SA0P2,0.0406310875,0.7415143737 +9-PZCT18_1SS-SA0P2,0.1632606472,0.6876855552 +10-PZCT18_1SS-SA0P2,0.3668635641,0.6594815850 +11-PZCT18_1SS-SA0P2,0.4610283669,0.6571228366 +12-PZCT18_1SS-SA0P2,0.0103373056,0.7415143737 +13-PZCT18_1SS-SA0P2,0.0523946168,0.6876855552 +14-PZCT18_1SS-SA0P2,0.1485135162,0.6594815850 +15-PZCT18_1SS-SA0P2,0.2354210501,0.6571228366 +0-PZCT18_1SS-SA0P25,0.3612137271,0.7333985547 +1-PZCT18_1SS-SA0P25,1.0467888548,0.6819113098 +2-PZCT18_1SS-SA0P25,1.6208914588,0.6549540724 +3-PZCT18_1SS-SA0P25,1.3410617556,0.6526979335 +4-PZCT18_1SS-SA0P25,0.1810941362,0.7333985547 +5-PZCT18_1SS-SA0P25,0.5893010179,0.6819113098 +6-PZCT18_1SS-SA0P25,1.0246369056,0.6549540724 +7-PZCT18_1SS-SA0P25,0.9519251419,0.6526979335 +8-PZCT18_1SS-SA0P25,0.0311170478,0.7333985547 +9-PZCT18_1SS-SA0P25,0.1350333420,0.6819113098 +10-PZCT18_1SS-SA0P25,0.3131004162,0.6549540724 +11-PZCT18_1SS-SA0P25,0.3879058005,0.6526979335 +12-PZCT18_1SS-SA0P25,0.0081712835,0.7333985547 +13-PZCT18_1SS-SA0P25,0.0446553514,0.6819113098 +14-PZCT18_1SS-SA0P25,0.1303937405,0.6549540724 +15-PZCT18_1SS-SA0P25,0.2034416997,0.6526979335 +0-PZCT18_1SS-SA0P3,0.2880676051,0.7355401909 +1-PZCT18_1SS-SA0P3,0.8856274414,0.6859973703 +2-PZCT18_1SS-SA0P3,1.4164710656,0.6602052168 +3-PZCT18_1SS-SA0P3,1.1785943377,0.6580626694 +4-PZCT18_1SS-SA0P3,0.1441040702,0.7355401909 +5-PZCT18_1SS-SA0P3,0.4981308463,0.6859973703 +6-PZCT18_1SS-SA0P3,0.8957991049,0.6602052168 +7-PZCT18_1SS-SA0P3,0.8380652098,0.6580626694 +8-PZCT18_1SS-SA0P3,0.0248899523,0.7355401909 +9-PZCT18_1SS-SA0P3,0.1150710790,0.6859973703 +10-PZCT18_1SS-SA0P3,0.2767631768,0.6602052168 +11-PZCT18_1SS-SA0P3,0.3462983550,0.6580626694 +12-PZCT18_1SS-SA0P3,0.0066810771,0.7355401909 +13-PZCT18_1SS-SA0P3,0.0389712080,0.6859973703 +14-PZCT18_1SS-SA0P3,0.1182609203,0.6602052168 +15-PZCT18_1SS-SA0P3,0.1866975594,0.6580626694 +0-PZCT18_1SS-SA0P4,0.1835282146,0.7356060427 +1-PZCT18_1SS-SA0P4,0.6475431998,0.6903163347 +2-PZCT18_1SS-SA0P4,1.1445496957,0.6671169802 +3-PZCT18_1SS-SA0P4,1.0134458801,0.6650936802 +4-PZCT18_1SS-SA0P4,0.0934529301,0.7356060427 +5-PZCT18_1SS-SA0P4,0.3678666035,0.6903163347 +6-PZCT18_1SS-SA0P4,0.7254171223,0.6671169802 +7-PZCT18_1SS-SA0P4,0.7166139231,0.6650936802 +8-PZCT18_1SS-SA0P4,0.0166462571,0.7356060427 +9-PZCT18_1SS-SA0P4,0.0863555053,0.6903163347 +10-PZCT18_1SS-SA0P4,0.2244209828,0.6671169802 +11-PZCT18_1SS-SA0P4,0.2921708543,0.6650936802 +12-PZCT18_1SS-SA0P4,0.0046886095,0.7356060427 +13-PZCT18_1SS-SA0P4,0.0302814553,0.6903163347 +14-PZCT18_1SS-SA0P4,0.0979736810,0.6671169802 +15-PZCT18_1SS-SA0P4,0.1587969322,0.6650936802 +0-PZCT18_1SS-SA0P5,0.1288895057,0.7420673824 +1-PZCT18_1SS-SA0P5,0.4950381124,0.7011008675 +2-PZCT18_1SS-SA0P5,0.9389872905,0.6797048859 +3-PZCT18_1SS-SA0P5,0.8795910304,0.6776974414 +4-PZCT18_1SS-SA0P5,0.0659426973,0.7420673824 +5-PZCT18_1SS-SA0P5,0.2819684136,0.7011008675 +6-PZCT18_1SS-SA0P5,0.5954348012,0.6797048859 +7-PZCT18_1SS-SA0P5,0.6209662081,0.6776974414 +8-PZCT18_1SS-SA0P5,0.0118682580,0.7420673824 +9-PZCT18_1SS-SA0P5,0.0666356989,0.7011008675 +10-PZCT18_1SS-SA0P5,0.1847680754,0.6797048859 +11-PZCT18_1SS-SA0P5,0.2530153753,0.6776974414 +12-PZCT18_1SS-SA0P5,0.0034277627,0.7420673824 +13-PZCT18_1SS-SA0P5,0.0238247851,0.7011008675 +14-PZCT18_1SS-SA0P5,0.0817800241,0.6797048859 +15-PZCT18_1SS-SA0P5,0.1386326343,0.6776974414 +0-PZCT18_1SS-SA0P75,0.0672498827,0.7603050851 +1-PZCT18_1SS-SA0P75,0.2955624068,0.7250381543 +2-PZCT18_1SS-SA0P75,0.6412197016,0.7067012454 +3-PZCT18_1SS-SA0P75,0.6866960052,0.7045988544 +4-PZCT18_1SS-SA0P75,0.0339521688,0.7603050851 +5-PZCT18_1SS-SA0P75,0.1665960733,0.7250381543 +6-PZCT18_1SS-SA0P75,0.4035171530,0.7067012454 +7-PZCT18_1SS-SA0P75,0.4824576984,0.7045988544 +8-PZCT18_1SS-SA0P75,0.0062060790,0.7603050851 +9-PZCT18_1SS-SA0P75,0.0400039195,0.7250381543 +10-PZCT18_1SS-SA0P75,0.1272880965,0.7067012454 +11-PZCT18_1SS-SA0P75,0.1999277321,0.7045988544 +12-PZCT18_1SS-SA0P75,0.0018630254,0.7603050851 +13-PZCT18_1SS-SA0P75,0.0147900513,0.7250381543 +14-PZCT18_1SS-SA0P75,0.0579589608,0.7067012454 +15-PZCT18_1SS-SA0P75,0.1121170079,0.7045988544 +0-PZCT18_1SS-SA1P0,0.0406645325,0.7590225557 +1-PZCT18_1SS-SA1P0,0.1958181671,0.7314962365 +2-PZCT18_1SS-SA1P0,0.4730316261,0.7172795072 +3-PZCT18_1SS-SA1P0,0.5732279721,0.7151122431 +4-PZCT18_1SS-SA1P0,0.0204216898,0.7590225557 +5-PZCT18_1SS-SA1P0,0.1097731337,0.7314962365 +6-PZCT18_1SS-SA1P0,0.2960059563,0.7172795072 +7-PZCT18_1SS-SA1P0,0.4004099289,0.7151122431 +8-PZCT18_1SS-SA1P0,0.0037638469,0.7590225557 +9-PZCT18_1SS-SA1P0,0.0264873834,0.7314962365 +10-PZCT18_1SS-SA1P0,0.0935074401,0.7172795072 +11-PZCT18_1SS-SA1P0,0.1655973100,0.7151122431 +12-PZCT18_1SS-SA1P0,0.0011446616,0.7590225557 +13-PZCT18_1SS-SA1P0,0.0098614525,0.7314962365 +14-PZCT18_1SS-SA1P0,0.0426191398,0.7172795072 +15-PZCT18_1SS-SA1P0,0.0923992554,0.7151122431 +0-PZCT18_1SS-SA1P5,0.0185683780,0.7511721257 +1-PZCT18_1SS-SA1P5,0.0985564256,0.7321074505 +2-PZCT18_1SS-SA1P5,0.2769465253,0.7222542561 +3-PZCT18_1SS-SA1P5,0.4120091842,0.7199950356 +4-PZCT18_1SS-SA1P5,0.0092005807,0.7511721257 +5-PZCT18_1SS-SA1P5,0.0548965802,0.7321074505 +6-PZCT18_1SS-SA1P5,0.1734104809,0.7222542561 +7-PZCT18_1SS-SA1P5,0.2900049391,0.7199950356 +8-PZCT18_1SS-SA1P5,0.0017475011,0.7511721257 +9-PZCT18_1SS-SA1P5,0.0137494356,0.7321074505 +10-PZCT18_1SS-SA1P5,0.0572733090,0.7222542561 +11-PZCT18_1SS-SA1P5,0.1263047488,0.7199950356 +12-PZCT18_1SS-SA1P5,0.0005528499,0.7511721257 +13-PZCT18_1SS-SA1P5,0.0053229073,0.7321074505 +14-PZCT18_1SS-SA1P5,0.0271325490,0.7222542561 +15-PZCT18_1SS-SA1P5,0.0732204330,0.7199950356 +0-PZCT18_1SS-SA2P0,0.0105667928,0.7444368677 +1-PZCT18_1SS-SA2P0,0.0589785876,0.7303300889 +2-PZCT18_1SS-SA2P0,0.1818216298,0.7228911679 +3-PZCT18_1SS-SA2P0,0.3095969863,0.7205937829 +4-PZCT18_1SS-SA2P0,0.0051863368,0.7444368677 +5-PZCT18_1SS-SA2P0,0.0327449862,0.7303300889 +6-PZCT18_1SS-SA2P0,0.1141902194,0.7228911679 +7-PZCT18_1SS-SA2P0,0.2199444481,0.7205937829 +8-PZCT18_1SS-SA2P0,0.0009898295,0.7444368677 +9-PZCT18_1SS-SA2P0,0.0083256995,0.7303300889 +10-PZCT18_1SS-SA2P0,0.0386795206,0.7228911679 +11-PZCT18_1SS-SA2P0,0.0992524951,0.7205937829 +12-PZCT18_1SS-SA2P0,0.0003155772,0.7444368677 +13-PZCT18_1SS-SA2P0,0.0032616481,0.7303300889 +14-PZCT18_1SS-SA2P0,0.0186195264,0.7228911679 +15-PZCT18_1SS-SA2P0,0.0587083888,0.7205937829 +0-PZCT18_1SS-SA3P0,0.0047548250,0.7436585843 +1-PZCT18_1SS-SA3P0,0.0284477244,0.7284659314 +2-PZCT18_1SS-SA3P0,0.1023671295,0.7205750992 +3-PZCT18_1SS-SA3P0,0.2215509001,0.7182120697 +4-PZCT18_1SS-SA3P0,0.0022926007,0.7436585843 +5-PZCT18_1SS-SA3P0,0.0156264904,0.7284659314 +6-PZCT18_1SS-SA3P0,0.0640611127,0.7205750992 +7-PZCT18_1SS-SA3P0,0.1579528643,0.7182120697 +8-PZCT18_1SS-SA3P0,0.0004344136,0.7436585843 +9-PZCT18_1SS-SA3P0,0.0039863245,0.7284659314 +10-PZCT18_1SS-SA3P0,0.0220009721,0.7205750992 +11-PZCT18_1SS-SA3P0,0.0730316465,0.7182120697 +12-PZCT18_1SS-SA3P0,0.0001386414,0.7436585843 +13-PZCT18_1SS-SA3P0,0.0015679016,0.7284659314 +14-PZCT18_1SS-SA3P0,0.0106646021,0.7205750992 +15-PZCT18_1SS-SA3P0,0.0436285237,0.7182120697 +0-PZCT18_1SS-SA4P0,0.0026428575,0.7251957684 +1-PZCT18_1SS-SA4P0,0.0165547211,0.7145294849 +2-PZCT18_1SS-SA4P0,0.0665012762,0.7087547283 +3-PZCT18_1SS-SA4P0,0.1713160656,0.7063981883 +4-PZCT18_1SS-SA4P0,0.0012951826,0.7251957684 +5-PZCT18_1SS-SA4P0,0.0092410918,0.7145294849 +6-PZCT18_1SS-SA4P0,0.0422839507,0.7087547283 +7-PZCT18_1SS-SA4P0,0.1240760126,0.7063981883 +8-PZCT18_1SS-SA4P0,0.0002485938,0.7251957684 +9-PZCT18_1SS-SA4P0,0.0023969961,0.7145294849 +10-PZCT18_1SS-SA4P0,0.0148219139,0.7087547283 +11-PZCT18_1SS-SA4P0,0.0587761644,0.7063981883 +12-PZCT18_1SS-SA4P0,0.0000794738,0.7251957684 +13-PZCT18_1SS-SA4P0,0.0009473539,0.7145294849 +14-PZCT18_1SS-SA4P0,0.0072420309,0.7087547283 +15-PZCT18_1SS-SA4P0,0.0355032814,0.7063981883 +0-PZCT18_1SS-SA5P0,0.0016716206,0.7256925038 +1-PZCT18_1SS-SA5P0,0.0110578953,0.7142132864 +2-PZCT18_1SS-SA5P0,0.0483779239,0.7082717611 +3-PZCT18_1SS-SA5P0,0.1399787084,0.7058921992 +4-PZCT18_1SS-SA5P0,0.0008281137,0.7256925038 +5-PZCT18_1SS-SA5P0,0.0062291875,0.7142132864 +6-PZCT18_1SS-SA5P0,0.0309893749,0.7082717611 +7-PZCT18_1SS-SA5P0,0.1019610338,0.7058921992 +8-PZCT18_1SS-SA5P0,0.0001593766,0.7256925038 +9-PZCT18_1SS-SA5P0,0.0016200460,0.7142132864 +10-PZCT18_1SS-SA5P0,0.0108910624,0.7082717611 +11-PZCT18_1SS-SA5P0,0.0484231632,0.7058921992 +12-PZCT18_1SS-SA5P0,0.0000506840,0.7256925038 +13-PZCT18_1SS-SA5P0,0.0006370179,0.7142132864 +14-PZCT18_1SS-SA5P0,0.0052950746,0.7082717611 +15-PZCT18_1SS-SA5P0,0.0291093495,0.7058921992 +0-PZCT18_1SS-SA7P5,0.0006335668,0.7197849129 +1-PZCT18_1SS-SA7P5,0.0046576533,0.7077473905 +2-PZCT18_1SS-SA7P5,0.0232124994,0.7017909691 +3-PZCT18_1SS-SA7P5,0.0784253840,0.6994640563 +4-PZCT18_1SS-SA7P5,0.0003219359,0.7197849129 +5-PZCT18_1SS-SA7P5,0.0026946088,0.7077473905 +6-PZCT18_1SS-SA7P5,0.0152898083,0.7017909691 +7-PZCT18_1SS-SA7P5,0.0588150346,0.6994640563 +8-PZCT18_1SS-SA7P5,0.0000644877,0.7197849129 +9-PZCT18_1SS-SA7P5,0.0007341417,0.7077473905 +10-PZCT18_1SS-SA7P5,0.0056658206,0.7017909691 +11-PZCT18_1SS-SA7P5,0.0296432453,0.6994640563 +12-PZCT18_1SS-SA7P5,0.0000211408,0.7197849129 +13-PZCT18_1SS-SA7P5,0.0002988621,0.7077473905 +14-PZCT18_1SS-SA7P5,0.0028641760,0.7017909691 +15-PZCT18_1SS-SA7P5,0.0186083825,0.6994640563 +0-PZCT18_1SS-SA10P0,0.0003478779,0.7137816613 +1-PZCT18_1SS-SA10P0,0.0026715324,0.6978945806 +2-PZCT18_1SS-SA10P0,0.0139390788,0.6901559077 +3-PZCT18_1SS-SA10P0,0.0494137376,0.6878382735 +4-PZCT18_1SS-SA10P0,0.0001735735,0.7137816613 +5-PZCT18_1SS-SA10P0,0.0015270631,0.6978945806 +6-PZCT18_1SS-SA10P0,0.0091278939,0.6901559077 +7-PZCT18_1SS-SA10P0,0.0370701366,0.6878382735 +8-PZCT18_1SS-SA10P0,0.0000340246,0.7137816613 +9-PZCT18_1SS-SA10P0,0.0004113107,0.6978945806 +10-PZCT18_1SS-SA10P0,0.0033782161,0.6901559077 +11-PZCT18_1SS-SA10P0,0.0188514535,0.6878382735 +12-PZCT18_1SS-SA10P0,0.0000111701,0.7137816613 +13-PZCT18_1SS-SA10P0,0.0001681761,0.6978945806 +14-PZCT18_1SS-SA10P0,0.0017203294,0.6901559077 +15-PZCT18_1SS-SA10P0,0.0119563556,0.6878382735 +0-PZCT18_2ES-PGA,0.7325561933,0.6892205743 +1-PZCT18_2ES-PGA,1.4689932724,0.6349744518 +2-PZCT18_2ES-PGA,2.0498338913,0.6066688348 +3-PZCT18_2ES-PGA,1.9903870245,0.6041068761 +4-PZCT18_2ES-PGA,0.3305552670,0.6892205743 +5-PZCT18_2ES-PGA,0.7406893621,0.6349744518 +6-PZCT18_2ES-PGA,1.1549099768,0.6066688348 +7-PZCT18_2ES-PGA,1.2530839252,0.6041068761 +8-PZCT18_2ES-PGA,0.0425048280,0.6892205743 +9-PZCT18_2ES-PGA,0.1259087761,0.6349744518 +10-PZCT18_2ES-PGA,0.2595336902,0.6066688348 +11-PZCT18_2ES-PGA,0.3722643784,0.6041068761 +12-PZCT18_2ES-PGA,0.0087888269,0.6892205743 +13-PZCT18_2ES-PGA,0.0323216283,0.6349744518 +14-PZCT18_2ES-PGA,0.0827133387,0.6066688348 +15-PZCT18_2ES-PGA,0.1472915633,0.6041068761 +0-PZCT18_2ES-SA0P01,1.1043311379,0.6906588684 +1-PZCT18_2ES-SA0P01,2.1338570683,0.6358379649 +2-PZCT18_2ES-SA0P01,2.9185781973,0.6071748967 +3-PZCT18_2ES-SA0P01,2.8256433313,0.6046150823 +4-PZCT18_2ES-SA0P01,0.4884294723,0.6906588684 +5-PZCT18_2ES-SA0P01,1.0532788306,0.6358379649 +6-PZCT18_2ES-SA0P01,1.6077734377,0.6071748967 +7-PZCT18_2ES-SA0P01,1.7371860893,0.6046150823 +8-PZCT18_2ES-SA0P01,0.0566014934,0.6906588684 +9-PZCT18_2ES-SA0P01,0.1620728488,0.6358379649 +10-PZCT18_2ES-SA0P01,0.3284979335,0.6071748967 +11-PZCT18_2ES-SA0P01,0.4712973718,0.6046150823 +12-PZCT18_2ES-SA0P01,0.0104232849,0.6906588684 +13-PZCT18_2ES-SA0P01,0.0373449715,0.6358379649 +14-PZCT18_2ES-SA0P01,0.0947108328,0.6071748967 +15-PZCT18_2ES-SA0P01,0.1700228056,0.6046150823 +0-PZCT18_2ES-SA0P02,1.7129695371,0.6925745736 +1-PZCT18_2ES-SA0P02,3.2941810985,0.6376850495 +2-PZCT18_2ES-SA0P02,4.5738856441,0.6090310755 +3-PZCT18_2ES-SA0P02,4.5852512031,0.6064689872 +4-PZCT18_2ES-SA0P02,0.7668819238,0.6925745736 +5-PZCT18_2ES-SA0P02,1.6364342227,0.6376850495 +6-PZCT18_2ES-SA0P02,2.5212075876,0.6090310755 +7-PZCT18_2ES-SA0P02,2.8045207503,0.6064689872 +8-PZCT18_2ES-SA0P02,0.0909700741,0.6925745736 +9-PZCT18_2ES-SA0P02,0.2543922120,0.6376850495 +10-PZCT18_2ES-SA0P02,0.5136284538,0.6090310755 +11-PZCT18_2ES-SA0P02,0.7487456936,0.6064689872 +12-PZCT18_2ES-SA0P02,0.0166625353,0.6925745736 +13-PZCT18_2ES-SA0P02,0.0573739028,0.6376850495 +14-PZCT18_2ES-SA0P02,0.1426355243,0.6090310755 +15-PZCT18_2ES-SA0P02,0.2560241969,0.6064689872 +0-PZCT18_2ES-SA0P03,1.6998827992,0.7059688042 +1-PZCT18_2ES-SA0P03,3.3713431856,0.6494790827 +2-PZCT18_2ES-SA0P03,4.8109068507,0.6197049088 +3-PZCT18_2ES-SA0P03,4.9395903384,0.6170588225 +4-PZCT18_2ES-SA0P03,0.7745850691,0.7059688042 +5-PZCT18_2ES-SA0P03,1.7038613354,0.6494790827 +6-PZCT18_2ES-SA0P03,2.6967418476,0.6197049088 +7-PZCT18_2ES-SA0P03,3.0710331268,0.6170588225 +8-PZCT18_2ES-SA0P03,0.0997889180,0.7059688042 +9-PZCT18_2ES-SA0P03,0.2859054020,0.6494790827 +10-PZCT18_2ES-SA0P03,0.5893895190,0.6197049088 +11-PZCT18_2ES-SA0P03,0.8742234586,0.6170588225 +12-PZCT18_2ES-SA0P03,0.0203115670,0.7059688042 +13-PZCT18_2ES-SA0P03,0.0704925782,0.6494790827 +14-PZCT18_2ES-SA0P03,0.1760286509,0.6197049088 +15-PZCT18_2ES-SA0P03,0.3162738799,0.6170588225 +0-PZCT18_2ES-SA0P05,1.2992770060,0.7299208861 +1-PZCT18_2ES-SA0P05,2.7113290442,0.6726480518 +2-PZCT18_2ES-SA0P05,3.9930311673,0.6423620441 +3-PZCT18_2ES-SA0P05,4.1501440751,0.6393193574 +4-PZCT18_2ES-SA0P05,0.6069369984,0.7299208861 +5-PZCT18_2ES-SA0P05,1.4128298632,0.6726480518 +6-PZCT18_2ES-SA0P05,2.3210052326,0.6423620441 +7-PZCT18_2ES-SA0P05,2.6909298301,0.6393193574 +8-PZCT18_2ES-SA0P05,0.0871459748,0.7299208861 +9-PZCT18_2ES-SA0P05,0.2657511775,0.6726480518 +10-PZCT18_2ES-SA0P05,0.5719302441,0.6423620441 +11-PZCT18_2ES-SA0P05,0.8686622098,0.6393193574 +12-PZCT18_2ES-SA0P05,0.0202343372,0.7299208861 +13-PZCT18_2ES-SA0P05,0.0744698463,0.6726480518 +14-PZCT18_2ES-SA0P05,0.1934246994,0.6423620441 +15-PZCT18_2ES-SA0P05,0.3545548283,0.6393193574 +0-PZCT18_2ES-SA0P1,0.8829902112,0.7555838339 +1-PZCT18_2ES-SA0P1,2.0312652388,0.7011332441 +2-PZCT18_2ES-SA0P1,3.0259195961,0.6719409113 +3-PZCT18_2ES-SA0P1,2.9189599169,0.6687547216 +4-PZCT18_2ES-SA0P1,0.4252351099,0.7555838339 +5-PZCT18_2ES-SA0P1,1.0851325218,0.7011332441 +6-PZCT18_2ES-SA0P1,1.7931492711,0.6719409113 +7-PZCT18_2ES-SA0P1,1.9188017610,0.6687547216 +8-PZCT18_2ES-SA0P1,0.0655292050,0.7555838339 +9-PZCT18_2ES-SA0P1,0.2156794577,0.7011332441 +10-PZCT18_2ES-SA0P1,0.4596874702,0.6719409113 +11-PZCT18_2ES-SA0P1,0.6344487553,0.6687547216 +12-PZCT18_2ES-SA0P1,0.0156357240,0.7555838339 +13-PZCT18_2ES-SA0P1,0.0626850712,0.7011332441 +14-PZCT18_2ES-SA0P1,0.1627384787,0.6719409113 +15-PZCT18_2ES-SA0P1,0.2735875814,0.6687547216 +0-PZCT18_2ES-SA0P15,0.6552878515,0.7495425688 +1-PZCT18_2ES-SA0P15,1.7076779105,0.6957413755 +2-PZCT18_2ES-SA0P15,2.6691919402,0.6675092847 +3-PZCT18_2ES-SA0P15,2.5023822389,0.6648213504 +4-PZCT18_2ES-SA0P15,0.3221898399,0.7495425688 +5-PZCT18_2ES-SA0P15,0.9211959002,0.6957413755 +6-PZCT18_2ES-SA0P15,1.5797637783,0.6675092847 +7-PZCT18_2ES-SA0P15,1.6249212090,0.6648213504 +8-PZCT18_2ES-SA0P15,0.0515875495,0.7495425688 +9-PZCT18_2ES-SA0P15,0.1856508491,0.6957413755 +10-PZCT18_2ES-SA0P15,0.4007273469,0.6675092847 +11-PZCT18_2ES-SA0P15,0.5188012600,0.6648213504 +12-PZCT18_2ES-SA0P15,0.0124562995,0.7495425688 +13-PZCT18_2ES-SA0P15,0.0541899286,0.6957413755 +14-PZCT18_2ES-SA0P15,0.1413995789,0.6675092847 +15-PZCT18_2ES-SA0P15,0.2212980840,0.6648213504 +0-PZCT18_2ES-SA0P2,0.4938191931,0.7415143737 +1-PZCT18_2ES-SA0P2,1.4293434715,0.6876855552 +2-PZCT18_2ES-SA0P2,2.3200921745,0.6594815850 +3-PZCT18_2ES-SA0P2,2.1119025574,0.6571228366 +4-PZCT18_2ES-SA0P2,0.2460210162,0.7415143737 +5-PZCT18_2ES-SA0P2,0.7767412599,0.6876855552 +6-PZCT18_2ES-SA0P2,1.3752464809,0.6594815850 +7-PZCT18_2ES-SA0P2,1.3654780066,0.6571228366 +8-PZCT18_2ES-SA0P2,0.0402284289,0.7415143737 +9-PZCT18_2ES-SA0P2,0.1580814299,0.6876855552 +10-PZCT18_2ES-SA0P2,0.3483602869,0.6594815850 +11-PZCT18_2ES-SA0P2,0.4305032393,0.6571228366 +12-PZCT18_2ES-SA0P2,0.0099897731,0.7415143737 +13-PZCT18_2ES-SA0P2,0.0472007837,0.6876855552 +14-PZCT18_2ES-SA0P2,0.1250670174,0.6594815850 +15-PZCT18_2ES-SA0P2,0.1858387931,0.6571228366 +0-PZCT18_2ES-SA0P25,0.3738354938,0.7333985547 +1-PZCT18_2ES-SA0P25,1.1719061165,0.6819113098 +2-PZCT18_2ES-SA0P25,1.9683626192,0.6549540724 +3-PZCT18_2ES-SA0P25,1.7714030850,0.6526979335 +4-PZCT18_2ES-SA0P25,0.1862338260,0.7333985547 +5-PZCT18_2ES-SA0P25,0.6365134586,0.6819113098 +6-PZCT18_2ES-SA0P25,1.1656189957,0.6549540724 +7-PZCT18_2ES-SA0P25,1.1436829372,0.6526979335 +8-PZCT18_2ES-SA0P25,0.0308568908,0.7333985547 +9-PZCT18_2ES-SA0P25,0.1310880923,0.6819113098 +10-PZCT18_2ES-SA0P25,0.2983831697,0.6549540724 +11-PZCT18_2ES-SA0P25,0.3639028046,0.6526979335 +12-PZCT18_2ES-SA0P25,0.0079252243,0.7333985547 +13-PZCT18_2ES-SA0P25,0.0404193493,0.6819113098 +14-PZCT18_2ES-SA0P25,0.1104503432,0.6549540724 +15-PZCT18_2ES-SA0P25,0.1617129582,0.6526979335 +0-PZCT18_2ES-SA0P3,0.2979462095,0.7355401909 +1-PZCT18_2ES-SA0P3,0.9880322874,0.6859973703 +2-PZCT18_2ES-SA0P3,1.7068855792,0.6602052168 +3-PZCT18_2ES-SA0P3,1.5361643479,0.6580626694 +4-PZCT18_2ES-SA0P3,0.1481776218,0.7355401909 +5-PZCT18_2ES-SA0P3,0.5366942678,0.6859973703 +6-PZCT18_2ES-SA0P3,1.0126778095,0.6602052168 +7-PZCT18_2ES-SA0P3,0.9954412425,0.6580626694 +8-PZCT18_2ES-SA0P3,0.0246535468,0.7355401909 +9-PZCT18_2ES-SA0P3,0.1115013839,0.6859973703 +10-PZCT18_2ES-SA0P3,0.2627127238,0.6602052168 +11-PZCT18_2ES-SA0P3,0.3224645344,0.6580626694 +12-PZCT18_2ES-SA0P3,0.0064798738,0.7355401909 +13-PZCT18_2ES-SA0P3,0.0352701163,0.6859973703 +14-PZCT18_2ES-SA0P3,0.1000108770,0.6602052168 +15-PZCT18_2ES-SA0P3,0.1477364232,0.6580626694 +0-PZCT18_2ES-SA0P4,0.1888164183,0.7356060427 +1-PZCT18_2ES-SA0P4,0.7142948511,0.6903163347 +2-PZCT18_2ES-SA0P4,1.3505643573,0.6671169802 +3-PZCT18_2ES-SA0P4,1.2763008209,0.6650936802 +4-PZCT18_2ES-SA0P4,0.0956990646,0.7356060427 +5-PZCT18_2ES-SA0P4,0.3928473359,0.6903163347 +6-PZCT18_2ES-SA0P4,0.8060090757,0.6671169802 +7-PZCT18_2ES-SA0P4,0.8265251093,0.6650936802 +8-PZCT18_2ES-SA0P4,0.0164304309,0.7356060427 +9-PZCT18_2ES-SA0P4,0.0831629927,0.6903163347 +10-PZCT18_2ES-SA0P4,0.2103833019,0.6671169802 +11-PZCT18_2ES-SA0P4,0.2660066073,0.6650936802 +12-PZCT18_2ES-SA0P4,0.0045367612,0.7356060427 +13-PZCT18_2ES-SA0P4,0.0272981044,0.6903163347 +14-PZCT18_2ES-SA0P4,0.0820954452,0.6671169802 +15-PZCT18_2ES-SA0P4,0.1233972935,0.6650936802 +0-PZCT18_2ES-SA0P5,0.1336739209,0.7420673824 +1-PZCT18_2ES-SA0P5,0.5492451947,0.7011008675 +2-PZCT18_2ES-SA0P5,1.1078624560,0.6797048859 +3-PZCT18_2ES-SA0P5,1.0969969596,0.6776974414 +4-PZCT18_2ES-SA0P5,0.0681279312,0.7420673824 +5-PZCT18_2ES-SA0P5,0.3034338056,0.7011008675 +6-PZCT18_2ES-SA0P5,0.6634416413,0.6797048859 +7-PZCT18_2ES-SA0P5,0.7121008261,0.6776974414 +8-PZCT18_2ES-SA0P5,0.0118363591,0.7420673824 +9-PZCT18_2ES-SA0P5,0.0649091052,0.7011008675 +10-PZCT18_2ES-SA0P5,0.1747403133,0.6797048859 +11-PZCT18_2ES-SA0P5,0.2309301290,0.6776974414 +12-PZCT18_2ES-SA0P5,0.0033514023,0.7420673824 +13-PZCT18_2ES-SA0P5,0.0217661209,0.7011008675 +14-PZCT18_2ES-SA0P5,0.0693961704,0.6797048859 +15-PZCT18_2ES-SA0P5,0.1086150416,0.6776974414 +0-PZCT18_2ES-SA0P75,0.0691631643,0.7603050851 +1-PZCT18_2ES-SA0P75,0.3191464949,0.7250381543 +2-PZCT18_2ES-SA0P75,0.7146695913,0.7067012454 +3-PZCT18_2ES-SA0P75,0.7766413831,0.7045988544 +4-PZCT18_2ES-SA0P75,0.0348180876,0.7603050851 +5-PZCT18_2ES-SA0P75,0.1749201780,0.7250381543 +6-PZCT18_2ES-SA0P75,0.4264566669,0.7067012454 +7-PZCT18_2ES-SA0P75,0.5045566166,0.7045988544 +8-PZCT18_2ES-SA0P75,0.0061633478,0.7603050851 +9-PZCT18_2ES-SA0P75,0.0383465048,0.7250381543 +10-PZCT18_2ES-SA0P75,0.1157803599,0.7067012454 +11-PZCT18_2ES-SA0P75,0.1696461491,0.7045988544 +12-PZCT18_2ES-SA0P75,0.0018116378,0.7603050851 +13-PZCT18_2ES-SA0P75,0.0133493594,0.7250381543 +14-PZCT18_2ES-SA0P75,0.0477363782,0.7067012454 +15-PZCT18_2ES-SA0P75,0.0828396981,0.7045988544 +0-PZCT18_2ES-SA1P0,0.0417822803,0.7590225557 +1-PZCT18_2ES-SA1P0,0.2077505840,0.7314962365 +2-PZCT18_2ES-SA1P0,0.5043039489,0.7172795072 +3-PZCT18_2ES-SA1P0,0.5976438908,0.7151122431 +4-PZCT18_2ES-SA1P0,0.0209600643,0.7590225557 +5-PZCT18_2ES-SA1P0,0.1136192765,0.7314962365 +6-PZCT18_2ES-SA1P0,0.3006847685,0.7172795072 +7-PZCT18_2ES-SA1P0,0.3884819442,0.7151122431 +8-PZCT18_2ES-SA1P0,0.0037510562,0.7590225557 +9-PZCT18_2ES-SA1P0,0.0252107926,0.7314962365 +10-PZCT18_2ES-SA1P0,0.0827216950,0.7172795072 +11-PZCT18_2ES-SA1P0,0.1325111279,0.7151122431 +12-PZCT18_2ES-SA1P0,0.0011162492,0.7590225557 +13-PZCT18_2ES-SA1P0,0.0088647587,0.7314962365 +14-PZCT18_2ES-SA1P0,0.0343694570,0.7172795072 +15-PZCT18_2ES-SA1P0,0.0650546615,0.7151122431 +0-PZCT18_2ES-SA1P5,0.0188227601,0.7511721257 +1-PZCT18_2ES-SA1P5,0.0995908959,0.7321074505 +2-PZCT18_2ES-SA1P5,0.2660455332,0.7222542561 +3-PZCT18_2ES-SA1P5,0.3588329408,0.7199950356 +4-PZCT18_2ES-SA1P5,0.0093298365,0.7511721257 +5-PZCT18_2ES-SA1P5,0.0543310263,0.7321074505 +6-PZCT18_2ES-SA1P5,0.1597429805,0.7222542561 +7-PZCT18_2ES-SA1P5,0.2371349998,0.7199950356 +8-PZCT18_2ES-SA1P5,0.0017263861,0.7511721257 +9-PZCT18_2ES-SA1P5,0.0126575723,0.7321074505 +10-PZCT18_2ES-SA1P5,0.0468557503,0.7222542561 +11-PZCT18_2ES-SA1P5,0.0875740381,0.7199950356 +12-PZCT18_2ES-SA1P5,0.0005348861,0.7511721257 +13-PZCT18_2ES-SA1P5,0.0046583397,0.7321074505 +14-PZCT18_2ES-SA1P5,0.0204833524,0.7222542561 +15-PZCT18_2ES-SA1P5,0.0454748123,0.7199950356 +0-PZCT18_2ES-SA2P0,0.0107828830,0.7444368677 +1-PZCT18_2ES-SA2P0,0.0585655112,0.7303300889 +2-PZCT18_2ES-SA2P0,0.1645693607,0.7228911679 +3-PZCT18_2ES-SA2P0,0.2392522475,0.7205937829 +4-PZCT18_2ES-SA2P0,0.0052961834,0.7444368677 +5-PZCT18_2ES-SA2P0,0.0319257805,0.7303300889 +6-PZCT18_2ES-SA2P0,0.0995680567,0.7228911679 +7-PZCT18_2ES-SA2P0,0.1606565874,0.7205937829 +8-PZCT18_2ES-SA2P0,0.0009881219,0.7444368677 +9-PZCT18_2ES-SA2P0,0.0076188228,0.7303300889 +10-PZCT18_2ES-SA2P0,0.0303924145,0.7228911679 +11-PZCT18_2ES-SA2P0,0.0627252691,0.7205937829 +12-PZCT18_2ES-SA2P0,0.0003091083,0.7444368677 +13-PZCT18_2ES-SA2P0,0.0028526742,0.7303300889 +14-PZCT18_2ES-SA2P0,0.0136205205,0.7228911679 +15-PZCT18_2ES-SA2P0,0.0336461386,0.7205937829 +0-PZCT18_2ES-SA3P0,0.0049472943,0.7436585843 +1-PZCT18_2ES-SA3P0,0.0283439414,0.7284659314 +2-PZCT18_2ES-SA3P0,0.0892385984,0.7205750992 +3-PZCT18_2ES-SA3P0,0.1543992691,0.7182120697 +4-PZCT18_2ES-SA3P0,0.0023785519,0.7436585843 +5-PZCT18_2ES-SA3P0,0.0152598812,0.7284659314 +6-PZCT18_2ES-SA3P0,0.0538009159,0.7205750992 +7-PZCT18_2ES-SA3P0,0.1042385081,0.7182120697 +8-PZCT18_2ES-SA3P0,0.0004418529,0.7436585843 +9-PZCT18_2ES-SA3P0,0.0036769510,0.7284659314 +10-PZCT18_2ES-SA3P0,0.0168150424,0.7205750992 +11-PZCT18_2ES-SA3P0,0.0422579168,0.7182120697 +12-PZCT18_2ES-SA3P0,0.0001388303,0.7436585843 +13-PZCT18_2ES-SA3P0,0.0013889642,0.7284659314 +14-PZCT18_2ES-SA3P0,0.0076365619,0.7205750992 +15-PZCT18_2ES-SA3P0,0.0230730343,0.7182120697 +0-PZCT18_2ES-SA4P0,0.0027746727,0.7251957684 +1-PZCT18_2ES-SA4P0,0.0167467976,0.7145294849 +2-PZCT18_2ES-SA4P0,0.0584052354,0.7087547283 +3-PZCT18_2ES-SA4P0,0.1176987320,0.7063981883 +4-PZCT18_2ES-SA4P0,0.0013479732,0.7251957684 +5-PZCT18_2ES-SA4P0,0.0090966421,0.7145294849 +6-PZCT18_2ES-SA4P0,0.0354716139,0.7087547283 +7-PZCT18_2ES-SA4P0,0.0799246582,0.7063981883 +8-PZCT18_2ES-SA4P0,0.0002533941,0.7251957684 +9-PZCT18_2ES-SA4P0,0.0022139706,0.7145294849 +10-PZCT18_2ES-SA4P0,0.0111775662,0.7087547283 +11-PZCT18_2ES-SA4P0,0.0326078864,0.7063981883 +12-PZCT18_2ES-SA4P0,0.0000799023,0.7251957684 +13-PZCT18_2ES-SA4P0,0.0008386953,0.7145294849 +14-PZCT18_2ES-SA4P0,0.0050868555,0.7087547283 +15-PZCT18_2ES-SA4P0,0.0178276696,0.7063981883 +0-PZCT18_2ES-SA5P0,0.0017794567,0.7256925038 +1-PZCT18_2ES-SA5P0,0.0115955111,0.7142132864 +2-PZCT18_2ES-SA5P0,0.0448014822,0.7082717611 +3-PZCT18_2ES-SA5P0,0.1026348534,0.7058921992 +4-PZCT18_2ES-SA5P0,0.0008686142,0.7256925038 +5-PZCT18_2ES-SA5P0,0.0062979477,0.7142132864 +6-PZCT18_2ES-SA5P0,0.0270751708,0.7082717611 +7-PZCT18_2ES-SA5P0,0.0690149783,0.7058921992 +8-PZCT18_2ES-SA5P0,0.0001631461,0.7256925038 +9-PZCT18_2ES-SA5P0,0.0015168371,0.7142132864 +10-PZCT18_2ES-SA5P0,0.0083618283,0.7082717611 +11-PZCT18_2ES-SA5P0,0.0273315050,0.7058921992 +12-PZCT18_2ES-SA5P0,0.0000511695,0.7256925038 +13-PZCT18_2ES-SA5P0,0.0005681295,0.7142132864 +14-PZCT18_2ES-SA5P0,0.0037401039,0.7082717611 +15-PZCT18_2ES-SA5P0,0.0145988950,0.7058921992 +0-PZCT18_2ES-SA7P5,0.0007008336,0.7197849129 +1-PZCT18_2ES-SA7P5,0.0053585485,0.7077473905 +2-PZCT18_2ES-SA7P5,0.0250080263,0.7017909691 +3-PZCT18_2ES-SA7P5,0.0712379923,0.6994640563 +4-PZCT18_2ES-SA7P5,0.0003481993,0.7197849129 +5-PZCT18_2ES-SA7P5,0.0029424054,0.7077473905 +6-PZCT18_2ES-SA7P5,0.0151766942,0.7017909691 +7-PZCT18_2ES-SA7P5,0.0477806286,0.6994640563 +8-PZCT18_2ES-SA7P5,0.0000674614,0.7197849129 +9-PZCT18_2ES-SA7P5,0.0007205623,0.7077473905 +10-PZCT18_2ES-SA7P5,0.0046977261,0.7017909691 +11-PZCT18_2ES-SA7P5,0.0186940700,0.6994640563 +12-PZCT18_2ES-SA7P5,0.0000217485,0.7197849129 +13-PZCT18_2ES-SA7P5,0.0002753910,0.7077473905 +14-PZCT18_2ES-SA7P5,0.0021284853,0.7017909691 +15-PZCT18_2ES-SA7P5,0.0100413395,0.6994640563 +0-PZCT18_2ES-SA10P0,0.0003852318,0.7137816613 +1-PZCT18_2ES-SA10P0,0.0031597492,0.6978945806 +2-PZCT18_2ES-SA10P0,0.0160613873,0.6901559077 +3-PZCT18_2ES-SA10P0,0.0505956868,0.6878382735 +4-PZCT18_2ES-SA10P0,0.0001877718,0.7137816613 +5-PZCT18_2ES-SA10P0,0.0017019172,0.6978945806 +6-PZCT18_2ES-SA10P0,0.0095597472,0.6901559077 +7-PZCT18_2ES-SA10P0,0.0332777841,0.6878382735 +8-PZCT18_2ES-SA10P0,0.0000353922,0.7137816613 +9-PZCT18_2ES-SA10P0,0.0004042622,0.6978945806 +10-PZCT18_2ES-SA10P0,0.0028616626,0.6901559077 +11-PZCT18_2ES-SA10P0,0.0125537558,0.6878382735 +12-PZCT18_2ES-SA10P0,0.0000113806,0.7137816613 +13-PZCT18_2ES-SA10P0,0.0001538237,0.6978945806 +14-PZCT18_2ES-SA10P0,0.0012884895,0.6901559077 +15-PZCT18_2ES-SA10P0,0.0066886519,0.6878382735