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.