code: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          
  
          COEFFS = self.COEFFS[imt]
          R = self._compute_term_r(COEFFS, rup.mag, dists.rrup)
  
          mean = 10 ** (self._compute_mean(COEFFS, rup.mag, R))
  
          # Convert units to g,
          # but only for PGA and SA (not PGV):
          if imt.name in "SA PGA":
              mean = np.log(mean / (g*100.))
          else:
              # PGV:
              mean = np.log(mean)
  
          c1_rrup = _compute_C1_term(COEFFS, dists.rrup)
          log_phi_ss = 1.00
  
          stddevs = self._get_stddevs(
              COEFFS, stddev_types, sites.vs30.shape[0], rup.mag, c1_rrup,
              log_phi_ss, COEFFS['mean_phi_ss']
          )
  
          return mean, stddevs
docs: See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for spec of input and result values.