anchor.diff:                 cycle loop_swap
                endif
                !
  -             ilevelf_ram(iswap_) = ilevelf
  -             ileveli_ram(iswap_) = ileveli
  -             abscoef_ram(iswap_) = abscoef
  -             acoef_ram(iswap_) = acoef
  -             !
  -             nu_ram(iswap_) = tranfreq
  +             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)
  +               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
  -           !
  +           !$omp end parallel do
  +           !
  +           iswap_ = 1
  +           !
  +           do iswap = 1,nswap_
  +             !
  +             if (ileveli_ram(iswap)>0) then 
  +                ilevelf_ram(iswap_) = ilevelf_ram(iswap)
  +                ileveli_ram(iswap_) = ileveli_ram(iswap)
  +                abscoef_ram(iswap_) = abscoef_ram(iswap)
  +                acoef_ram(iswap_)   = acoef_ram(iswap)
  +                !
  +                nu_ram(iswap_) = nu_ram(iswap)
  +                !
  +                if (lineprofile_do) then
  +                  gamma_RAM(iswap_)=gamma_RAM(iswap)
  +                endif
  +                !
  +                iswap_ = iswap_ + 1
  +                !
  +             endif
  +             !
  +           enddo
              nswap = iswap_-1
  +           !
  +           intband = intband + sum(abscoef_ram(1:nswap))
              !
           endif
           !
  
positive.code:            !
             !read(tunit,*,end=20) indexf,indexi,acoef
             iswap_ = 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)
             loop_swap : do iswap = 1,nswap_
               !
               indexf = indexf_RAM(iswap)
  
negative.code:                !
               case ('LIFETIME')
                 !
                 !$omp critical
                 !
                 if (Asum(ilevelf)<0) Asum(ilevelf) = 0
                 !
                 Asum(ilevelf) = Asum(ilevelf) + acoef
                 !$omp end critical
                 !
                 !print*,ilevelf,indexf,indexi,acoef
                 !
  