anchor.diff:        !
         halfwidth = 0 
         !
  +      lineprofile_do = .true.
  +      !
         do i=1,Nspecies
           !
           halfwidth =  halfwidth + species(i)%ratio*species(i)%gamma*(species(i)%T0/Temp)**species(i)%N*pressure/species(i)%P0
  
positive.code:             !
          case ('VOI-Q') ! VOIGT-QUADRATURES
              !
              halfwidth=get_Voigt_gamma_n(Nspecies,Jf,Ji)
              !
              intens_omp = 0
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
                !
                do iswap = iomp,nswap,N_omp_procs
                  !
                  abscoef = abscoef_ram(iswap)
                  tranfreq = nu_ram(iswap)
                  !
                  call do_Voi_Q(tranfreq,abscoef,dfreq,freq,abciss,weight,halfwidth,offset,freql,dpwcoef,intens_omp(:,iomp))
                  !
  
negative.code:      real(rk),intent(in) :: Jf,Ji
       real(rk) :: halfwidth,gamma_,n_,f
       integer(ik) :: ispecies,Jpp,Jp
       !
       halfwidth = 0
       !
  