code: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          
  
          # Extracting dictionary of coefficients specific to required
          # intensity measure type.
          C = self.COEFFS[imt]
  
          imean = (self._compute_magnitude(rup, C) +
                   self._compute_attenuation(rup, dists, imt, C) +
                   self._compute_distance(rup, dists, imt, C))
  
          mean = np.log(10.0 ** (imean))
  
          istddevs = self._get_stddevs(C, stddev_types, rup, imt,
                                       num_sites=len(dists.rrup))
  
          stddevs = np.log(10.0 ** np.array(istddevs))
  
          return mean, stddevs
docs: See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for spec of input and result values.