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]
          C_SITE = self.SITE_COEFFS[imt]
          s_c, idx = self._get_site_classification(sites.vs30)
          sa_rock = (self.get_magnitude_scaling_term(C, rup) +
                     self.get_sof_term(C, rup) +
                     self.get_depth_term(C, rup) +
                     self.get_distance_term(C, dists, rup))
  
          sa_soil = self.add_site_amplification(C, C_SITE, sites,
                                                sa_rock, idx, rup)
  
          stddevs = self.get_stddevs(C, sites.vs30.shape, idx, stddev_types)
          return sa_soil, stddevs
docs: See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for spec of input and result values.