anchor.diff:             !
              !read(tunit,*,end=20) indexf,indexi,acoef
              iswap_ = 1
  +           !
  +           ileveli_ram = -1
              abscoef_ram = 0
              !
  -           !omp parallel do private(iswap,indexf,indexi,acoef,ilevelf,ileveli,energyf,energyi,ifilter,jf,ji,tranfreq,&
  -           !omp& tranfreq0,offset,abscoef,cutoff)&
  -           !omp&  schedule(static) reduction(+:iswap_,intband) shared(ilevelf_ram,ileveli_ram,abscoef_ram,acoef_ram,nu_ram,Asum)
  +           !$omp  parallel do private(iswap,indexf,indexi,acoef,ilevelf,ileveli,energyf,energyi,ifilter,&
  +           !$omp& ivib,ener_vib,ener_rot,jf,ji,tranfreq,tranfreq0,offset,abscoef,cutoff)&
  +           !$omp& schedule(static) shared(ilevelf_ram,ileveli_ram,abscoef_ram,acoef_ram,nu_ram,Asum,gamma_RAM)
              loop_swap : do iswap = 1,nswap_
                !
                indexf = indexf_RAM(iswap)
  
positive.code:                cycle loop_swap
               endif
               !
               ilevelf_ram(iswap_) = ilevelf
               ileveli_ram(iswap_) = ileveli
               abscoef_ram(iswap_) = abscoef
               acoef_ram(iswap_) = acoef
               !
               nu_ram(iswap_) = tranfreq
               !
               if (lineprofile_do) then
                 gamma_RAM(iswap_)=get_Voigt_gamma_n(Nspecies,Jf,Ji)
               endif
               !
               intband = intband + abscoef
               !
               !do do_log_intensity
               iswap_ = iswap_ + 1
               !
             enddo loop_swap
             !omp end parallel do
             !
             nswap = iswap_-1
             !
          endif
          !
  
negative.code:            intband = intband + sum(abscoef_RAM(1:nswap_))
             !
             if (nswap_<1) then
                if (verbose>=5) write(out,"('... Finish swap')")
                exit loop_tran
             endif
             !
  