docs: def JMS_to_array(C, sectors=None):
      """For a dictionary with JMS Wilson coefficients, return a dictionary
      of arrays."""
      if sectors is None:
          wc_keys = wcxf.Basis['WET', 'JMS'].all_wcs
      else:
          try:
              wc_keys = [k for s in sectors for k in wcxf.Basis['WET', 'JMS'].sectors[s]]
          except KeyError:
              print(sectors)
      # fill in zeros for missing coefficients
      C_complete = {k: C.get(k, 0) for k in wc_keys}
      Ca = _scalar2array(C_complete)
      for k in Ca:
          if k in C_symm_keys[5]:
              Ca[k] = _symm_herm(Ca[k])
          if k in C_symm_keys[41]:
              Ca[k] = _symm_current(Ca[k])
          if k in C_symm_keys[4]:
              Ca[k] = _symm_herm(_symm_current(Ca[k]))
          if k in C_symm_keys[9]:
              Ca[k] = _antisymm_12(Ca[k])
      return Ca