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.
          """
          # extracting dictionary of coefficients specific to required
          # intensity measure type.
  
          C = self.COEFFS[imt]
          mean = (self._get_magnitude_scaling(C, rup.mag) +
                  self._get_distance_scaling(C, dists, rup.mag) +
                  self._get_site_term(C, sites.vs30))
  
          # Mean is returned in terms of m/s^2. Need to convert to g
          mean -= np.log(g)
          stddevs = self.get_stddevs(C, sites.vs30.shape, stddev_types)
          return mean + self.adjustment_factor, stddevs