anchor.diff:                 species(2)%gamma = gamma_s
                  !species(2)%n = n_
                  !
  -               gamma_ram(iswap_) = gamma_ * pressure * (hitran_Tref/temp)**n_
  -
  +               gamma_ram(iswap_) = gamma_ * pressure * (hitran_Tref/temp)**n_*species(1)%ratio
  +               if (Nspecies>1) gamma_ram(iswap_) = gamma_ram(iswap_) + gamma_s * pressure * species(2)%ratio
  +               !
                  !Below: probably not correct, but gamma_idx_RAM definitely needs
                  !to be initialized somehow
                  gamma_idx_RAM(iswap_) = fast_voigt%get_size()
  
positive.code:              ji = jrot(ileveli)
               !
               tranfreq = energyf-energyi
               !
               tranfreq0 = tranfreq
               !
  
negative.code:   type speciesT ! broadener
      !
      real(rk)       :: N = 0.5_rk      ! Voigt parameter N
      real(rk)       :: gamma = 0.0_rk  ! Voigt parameter gamma
      real(rk)       :: ratio = 1.0_rk  ! Ratio
      real(rk)       :: T0    = 296_rk  ! Reference T, K
      real(rk)       :: P0    = 1.0_rk  ! Reference P, bar
  