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_term(C, rup.mag) +
                  self._get_distance_term(C, dists.rjb, rup.mag) +
                  self._get_site_term(C, sites.vs30))
  
          # Units of GMPE are in terms of m/s (corrected in an Erratum)
          # Convert to g
          if imt.name in "SA PGA":
              mean = np.log(np.exp(mean) / g)
          else:
              # For PGV convert from m/s to cm/s/s
              mean = np.log(np.exp(mean) * 100.)
  
          # Get standard deviations
          stddevs = self._get_stddevs(C, stddev_types, dists.rjb.shape)
          return mean, stddevs