anchor.diff:             !
              do while(iswap_<N_to_RAM)
                !
  -             !read(tunit,*,end=119) indexf_RAM(iswap),indexi_RAM(iswap),acoef_RAM(iswap)
  -             read(tunit,*,end=120) &
  -                  imol,isotope,tranfreq,abscoef,energyi,gamma_,n_
  +             !
  +             read(tunit,*,end=120) imol,isotope,tranfreq,abscoef,energyi,gamma_,n_
  +             !
  +             if (imol/=imolecule.or.iso/=isotope) cycle
                !
                energyf = tranfreq + energyi
  -             !
  -             if (imol/=imolecule.or.iso/=isotope) cycle
                !
                if (tranfreq<small_) cycle
                !
  
positive.code:         case ("SPECTRA")
            !
            spectra_do = .true.
            !
          case ("HITRAN","OXFORD")
            !
  
negative.code:   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., 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
  