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]
          C_SIG = self.SIGMA_COEFFS[imt]
  
          mean = (self.get_magnitude_scaling_term(C, rup.mag) +
                  self.get_distance_scaling_term(C, dists.rhypo))
  
          std_devs = self.get_stddevs(C_SIG, stddev_types, len(dists.rhypo))
  
          #: Mean ground motions initially returned in cm/s/s (for PGA, SA)
          #: and cm/s for PGV
          if not imt.name == "PGV":
              # Convert mean from log(cm/s/s) to g
              mean = np.log(np.exp(mean) / (100. * g))
  
          return mean, std_devs