code: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          # pylint: disable=too-many-arguments
          
          ln_mean, [ln_stddev] = super().get_mean_and_stddevs(
              sites, rup, dists, imt, stddev_types)
  
          # compute site corrections, equation (9)
          coeffs = self.COEFFS_UPPER[imt]
          ln_mean += np.log(coeffs['correction'])
  
          return ln_mean, [ln_stddev]
docs: See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for specification of input and result values.
  
          Implements the correction factor for the upper crust, equation (12) on
          p. 484:
  
          ``P' = P x Correction_factor``