docs: def _compute_stress_drop_adjustment(self, SC, mag, scale_fac):
          """
          Compute equation (6) p. 2200
          """
          return scale_fac * np.minimum(
              SC['delta'] + 0.05,
              0.05 + SC['delta'] * (
                  np.maximum(mag - SC['M1'], 0) / (SC['Mh'] - SC['M1'])
              )
          )