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.
          """
          # get the necessary set of coefficients
          C = self.COEFFS[imt]
          # compute median sa on rock (vs30=1180m/s). Used for site response
          # term calculation
          sa1180 = np.exp(self._get_sa_at_1180(C, imt, sites, rup, dists))
  
          # get the mean value
          mean = (self._get_basic_term(C, rup, dists) +
                  self._get_faulting_style_term(C, rup) +
                  self._get_site_response_term(C, imt, sites.vs30, sa1180) +
                  self._get_hanging_wall_term(C, dists, rup) +
                  self._get_top_of_rupture_depth_term(C, imt, rup) +
                  self._get_soil_depth_term(C, sites.z1pt0 / METRES_PER_KM,
                                            sites.vs30)
                  )
          mean += self._get_regional_term(C, imt, sites.vs30, dists.rrup)
          # get standard deviations
          stddevs = self._get_stddevs(C, imt, rup, sites, stddev_types, sa1180,
                                      dists)
          return mean, stddevs