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.
          """
  
          sites.vs30 = 700 * np.ones(len(sites.vs30))
  
          mean, stddevs = super().get_mean_and_stddevs(
              sites, rup, dists, imt, stddev_types)
  
          C = CauzziFaccioli2008SWISS01.COEFFS
  
          tau_ss = 'tau'
          log_phi_ss = np.log(10)
          mean, stddevs = _apply_adjustments(
              C, self.COEFFS_FS_ROCK[imt], tau_ss,
              mean, stddevs, sites, rup, dists.rhypo, imt, stddev_types,
              log_phi_ss)
  
          return mean, np.log(10 ** np.array(stddevs))