docs: def random_galactic_extinction(
          log,
          sampleNumber,
          extinctionType,
          extinctionConstant,
          galacticExtinctionDistributions,
          plot=False):
      """
      *Generate a Numpy array of random galactic extinctions*
  
      **Key Arguments:**
          - ``log`` -- logger
          - ``sampleNumber`` -- the sample number, i.e. array size
          - ``extinctionType`` -- constant or random?
          - ``extinctionConstant`` -- the constant value (when extinctionType is constant)
          - ``galacticExtinctionDistributions`` -- the galactic extinction distribution (when extinctionType is random)
          - ``plot`` -- generate plot?
  
      **Return:**
          - ``galacticExtinctionArray``
      """
      ################ > IMPORTS ################
      ## STANDARD LIB ##
      ## THIRD PARTY ##
      import numpy as np
      ## LOCAL APPLICATION ##
  
      ################ > VARIABLE SETTINGS #######
  
      ################ >ACTION(S) ################
      # xxx come back here and add in the random extinctions generation -- will
      # need to account for what type of SN we have ##
      if extinctionType == "constant":
          galacticExtinctionArray = np.ones(sampleNumber) * extinctionConstant
      else:
          log.error('galactic extiction distributions not included yet')
  
      # log.debug('galacticExtinctionArray %s' % (galacticExtinctionArray,))
  
      return galacticExtinctionArray