anchor.diff:    real(rk),allocatable,save  :: gamma_comb(:,:) !Used for indexing the gamma for the fast_voigt
     logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran_do = .false.,  histogramJ = .false., &
                stick_hitran = .false.,stick_oxford = .false.,vibtemperature_do = .false., spectra_do = .false., &
  -             vibpopulation_do = .false.
  +             vibpopulation_do = .false., super_energies_do = .false.
     logical :: lineprofile_do = .false., use_width_offset = .false.
     logical :: microns = .false.
     logical :: use_resolving_power = .false.  ! using resolving for creating the grid
  
positive.code:      stop 'Undefined PF'
     endif
     !
     if (spectra_do.and.(partfunc<0.0.or.partfunc_ref<0.0)) then
       write(out,"('For SPECTRA pf and ref-pf must be defined using PF XXXX REF YYYY')")
       stop 'Undefined PF'
     endif
     !
     if (spectra_do.and.(temp_ref<0.0)) then
       write(out,"('For SPECTRA reference temperatures must be defined using temperature XXXX REF YYYY')")
       stop 'Undefined reference Temperature'
     endif
  
negative.code:         case ("SPECTRA")
            !
            spectra_do = .true.
            !
          case ("SUPER-ENERGIES")
            !
            super_energies_do = .true.
            !
          case ("HITRAN","OXFORD")
            !
  