anchor.diff:                 species(2)%gamma = gamma_s
                  !species(2)%n = n_
                  !
  -               gamma_RAM(iswap_) = gamma_ * pressure * (296.0/temp)**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()
                  !
                endif
  
positive.code:    call ArrayStart('swap:ilevelf_RAM',info,size(ilevelf_RAM),kind(ilevelf_RAM))
     !
     if (lineprofile_do) then
       allocate(gamma_RAM(N_to_RAM),stat=info)
       allocate(gamma_idx_RAM(N_to_RAM),stat=info)
       call ArrayStart('swap:gamma_RAM',info,size(gamma_RAM),kind(gamma_RAM))
     endif
     !
     if (verbose>=4) call MemoryReport
  
negative.code:                 abscoef = abscoef_ram(iswap)
                  tranfreq = nu_ram(iswap)
                  !
                  call do_Voigt_Fast(tranfreq,abscoef,use_resolving_power, &
                       freq,gamma_idx_RAM(iswap),dpwcoef,offset,freql,intens_omp(:,iomp))
                  !
                enddo
                !
  