docs: def get_mean_values(self, C, sites, rup, dists, a1100):
          """
          Returns the mean values for a specific IMT
          """
          if isinstance(a1100, np.ndarray):
              # Site model defined
              temp_vs30 = sites.vs30
              temp_z2pt5 = sites.z2pt5
          else:
              # Default site and basin model
              temp_vs30 = 1100.0 * np.ones(len(sites.vs30))
              temp_z2pt5 = self._select_basin_model(1100.0) *\
                  np.ones_like(temp_vs30)
  
          return (self._get_magnitude_term(C, rup.mag) +
                  self._get_geometric_attenuation_term(C, rup.mag, dists.rrup) +
                  self._get_style_of_faulting_term(C, rup) +
                  self._get_hanging_wall_term(C, rup, dists) +
                  self._get_shallow_site_response_term(C, temp_vs30, a1100) +
                  self._get_basin_response_term(C, temp_z2pt5) +
                  self._get_hypocentral_depth_term(C, rup) +
                  self._get_fault_dip_term(C, rup) +
                  self._get_anelastic_attenuation_term(C, dists.rrup))