anchor.diff:         !
     end subroutine voi_quad 
   
  -
  -
  +  subroutine get_Voigt_gamma_n(Nspecies,Jf,Ji)
  +     !
  +     implicit none 
  +     !
  +     integer(ik),intent(in) :: Nspecies
  +     real(rk),intent(in) :: Jf,Ji
  +     real(rk) :: halfwidth,gamma_,n_
  +     integer(ik) :: ispecies,Jpp,Jp
  +     
  +     !
  +     halfwidth = 0
  +     !
  +     do ispecies =1,Nspecies
  +       !
  +       if ( trim(species(ispecies)%filename)/="" ) then
  +         !
  +         Jpp = nint(Ji)
  +         Jp  = nint(Jf)
  +         !
  +         gamma_ = species(ispecies)%gammaQN(Jpp,Jp-Jpp)
  +         n_     = species(ispecies)%nQN(Jpp,Jp-Jpp)
  +         !
  +       else
  +         !
  +         gamma_ = species(ispecies)%gamma
  +         n_     = species(ispecies)%n
  +         !
  +       endif
  +       !
  +       halfwidth =  halfwidth + species(ispecies)%ratio*gamma_*(species(ispecies)%T0/Temp)**N_*pressure/species(ispecies)%P0
  +       !
  +     enddo
  +
  +  end subroutine get_Voigt_gamma_n
   
     !
     function voigt_faddeeva(nu,nu0,alpha,gamma) result(f)
  +   !
  +   implicit none 
      !
      real(rk),intent(in) :: nu,nu0,alpha,gamma
      real(rk) :: f,sigma,x,y
  
positive.code: module spectrum
    !
    use accuracy
    use timer
  
negative.code:     integer(ik) :: ierr = 0           ! Number of QN-entries
      !
    end type HitranErrorT
    !
    type(selectT) :: upper(filtermax),lower(filtermax)
    !
  