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 original values
          hslab = 50  # See info in GMPEt_Inslab_med.dat
          rjb, rrup = utils.get_equivalent_distance_inslab(rup.mag, dists.repi,
                                                           hslab)
          dists.rjb = rjb
          dists.rrup = rrup
          mean, stddevs = super().get_mean_and_stddevs(sites, rup, dists, imt,
                                                       stddev_types)
          cff = self.SITE_COEFFS[imt]
          mean_adj = np.log(np.exp(mean) * 10**cff['mf'])
          stddevs = [np.ones(len(dists.rrup))*get_sigma(imt)]
          return mean_adj, stddevs