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:                !
                 Jf = 0.0_rk
                 Ji = 0.0_rk
                 species(1)%gamma = gamma_
                 species(1)%n = n_
                 !
                 gamma_ram(iswap_) = gamma_ * pressure * (temp_ref/temp)**n_
  
                 !Below: probably not correct, but gamma_idx_RAM definitely needs
                 !to be initialized somehow
  
negative.code:            !
             do while(iswap_<N_to_RAM)
               !
               !
               read(tunit,*,end=120) imol,isotope,tranfreq,abscoef,energyi,gamma_,n_
               !
               if (imol/=imolecule.or.iso/=isotope) cycle
               !
               energyf = tranfreq + energyi
               !
               if (tranfreq<small_) cycle
               !
  