anchor.diff:               ji = jrot(ileveli)
                !
                tranfreq = energyf-energyi
  +             !
  +             tranfreq = tranfreq + lineshift
                !
                tranfreq0 = tranfreq
                !
  
positive.code:                species(2)%gamma = gamma_s
                 !species(2)%n = n_
                 !
                 gamma_ram(iswap_) = gamma_ * pressure * (hitran_Tref/temp)**n_
  
                 !Below: probably not correct, but gamma_idx_RAM definitely needs
                 !to be initialized somehow
                 gamma_idx_RAM(iswap_) = fast_voigt%get_size()
  
negative.code:                 tranfreq = nu_ram(iswap)
                  halfwidth = gamma_ram(iswap)
                  !
                  call do_Voi_Q(tranfreq,abscoef,dfreq,freq,abciss,weight,halfwidth,offset,freql,dpwcoef,intens_omp(:,iomp))
                  !
                enddo
                !
              enddo
              !$omp end parallel do
              !
          case ('VOI-9') ! VOIGT-Alg-916
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq,halfwidth) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
  