/orange/adamginsburg/galactic_plane_surveys/uchii_2026/pg2026_model/mass_resolved_boost.py:56: IntegrationWarning: The maximum number of subdivisions (80) has been achieved.
  If increasing the limit yields no improvement it is advised to analyze 
  the integrand in order to determine the difficulties.  If the position of a 
  local difficulty can be determined (singularity, discontinuity) one will 
  probably gain from splitting up the interval and calling the integrator 
  on the subranges.  Perhaps a special-purpose integrator should be used.
  val, _ = quad(integrand, lo, hi, limit=80)