docs: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          # pylint: disable=too-many-arguments
          """
          See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for specification of input and result values.
          """
  
          # extract dictionary of coefficients specific to required
          # intensity measure type
          coeffs = self.COEFFS[imt]
          coeffs.update(self.CONSTS)
  
          # equation (1) is in terms of common logarithm
          log_mean = (self._compute_magnitude(rup, coeffs) +
                      self._compute_distance(dists, coeffs) +
                      self._get_site_amplification(sites, coeffs) +
                      self._get_mechanism(rup, coeffs))
          # so convert to g and thence to the natural logarithm
          mean = log_mean*np.log(10.0) - np.log(g)
  
          # convert standard deviations from common to natural logarithm
          log_stddevs = self._get_stddevs(coeffs, stddev_types, len(sites.vs30))
          stddevs = log_stddevs*np.log(10.0)
  
          return mean, stddevs