code: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          
          C = self.COEFFS[imt]
  
          mean = (self._compute_magnitude_term(C, rup.mag) +
                  self._compute_distance_term(C, dists.rhypo, rup.mag))
          stddevs = self._get_stddevs(C, dists.rhypo, stddev_types)
          return mean, stddevs
docs: See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for spec of input and result values.