code: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          
          assert all(stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
                     for stddev_type in stddev_types)
  
          C = self.COEFFS[imt]
  
          imean = (self._compute_magnitude_scaling(C, rup.mag) +
                   self._compute_distance_scaling(C, dists.rrup, rup.mag))
          # Original GMPE returns log10 acceleration in cm/s/s
          # Converts to natural logarithm of g
          mean = np.log((10.0 ** (imean - 2.0)) / g)
          mean = self._compute_site_scaling(sites.vs30, mean)
          istddevs = self._compute_stddevs(
              C, dists.rrup.shape, stddev_types
          )
          # Convert from common logarithm to natural logarithm
          stddevs = np.log(10 ** np.array(istddevs))
          return mean, stddevs
docs: See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for spec of input and result values.