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.
          """
          mean = self._get_mean(sites.vs30, rup.mag, dists.rrup, imt, scale_fac=0)
          stddevs = self._get_stddevs(stddev_types, num_sites=sites.vs30.size)
  
          return mean, stddevs