Skip to content

Misplaced brace in Kuehn 2020 basin term gives nonzero values when basin == false

A misplaced brace at KuenEtAl_2020.java L481 results in the cap for the Seattle basin being applied outside Seattle (seattle == false) and for basin == true and basin == false. I believe the correct code is

    if (seattle) { // basin is implicitely true for seattle() and m9()
      fb = c.θ11S; // period-dependent seattle scale factor
      if (m9 && !slab && z2p5 > 6.0 && (c.imt.isSA() && c.imt.period() > 1.9)) {
        fb = log(2.0); // M9 scaling
      } else {
	  fb = Math.min(c.θ11 + c.θ12 * δlnZ, c.θ11S); // Seattle cap
	  /*
	   * NOTE: If we stop using the USGS basin model in NSHMs, then at shorter
	   * periods c.θ11S can often impose significant deamplification, more than
	   * c.θ11 + c.θ12 * δlnZ and this was not likely the intent of the
	   * developers.
	   */
      }
    }
    if (basin) {
      fb *= GmmUtils.deltaZ25scale(c.imt, z2p5);
    }

If the cap is meant to be applied for all cases when seattle == true (case with M9 scaling), then the else should be removed.