docs: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          """
          See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for spec of input and result values.
          """
          C = self.COEFFS[imt]
          mean = (
              self._get_magnitude_term(C, rup.mag) +
              self._get_distance_term(C, rup.mag, dists.rrup) +
              self._get_style_of_faulting_term(C, rup.rake) +
              self._get_site_response_term(C, sites.vs30))
  
          stddevs = self._get_stddevs(C, stddev_types, len(sites.vs30))
          return mean, stddevs