code: def Inorm_bol_bb(Teff=5772., logg=4.43, abun=0.0, atm='blackbody', photon_weighted=False):
      
  
      if atm != 'blackbody':
          raise ValueError('atmosphere must be set to blackbody for Inorm_bol_bb.')
  
      if photon_weighted:
          factor = 2.6814126821264836e22/Teff
      else:
          factor = 1.0
  
      # convert scalars to vectors if necessary:
      if not hasattr(Teff, '__iter__'):
          Teff = np.array((Teff,))
  
      return factor * sigma_sb.value * Teff**4 / np.pi
docs: @Teff: value or array of effective temperatures
      @logg: surface gravity; not used, for class compatibility only
      @abun: abundances; not used, for class compatibility only
      @atm: atmosphere model, must be blackbody, otherwise exception is raised
      @photon_weighted: intensity weighting scheme; must be False, otherwise exception is raised
  
      Computes normal bolometric intensity using the Stefan-Boltzmann law,
      Inorm_bol_bb = 1/\pi \sigma T^4. If photon-weighted intensity is
      requested, Inorm_bol_bb is multiplied by a conversion factor that
      comes from integrating lambda/hc P(lambda) over all lambda.
  
      Input parameters mimick the Passband class Inorm method for calling
      convenience.