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.
          """
          # Check that the GSIM supports the standard deviations requested
          assert all(stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
                     for stddev_type in stddev_types)
          coe = self.COEFFS[imt]
          mean = (self._get_magnitude_scaling(coe, rup.mag) +
                  self._get_distance_scaling(coe, dists.rhypo) +
                  self._get_azimuth_correction(coe, dists.azimuth))
          # Convert to g
          if imt.name in "SA PGA":
              mean = np.log(np.exp(mean) / (100.0 * g))
          # Compute std
          stddevs = self._compute_std(coe, stddev_types, dists.azimuth.shape)
          return mean, stddevs