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]
          # Get the mean
          mean = self.get_mean(C, rup, sites, dists)
          if imt.name == "PGV":
              # Convert from log10 m/s to ln cm/s
              mean = np.log((10.0 ** mean) * 100.)
          else:
              # convert from log10 m/s/s to ln g
              mean = np.log((10.0 ** mean) / g)
  
          # Get the standard deviations
          stddevs = self.get_stddevs(C, mean.shape, stddev_types)
          return mean, stddevs