anchor.diff:      real(rk)       :: ratio = 1.0_rk  ! Ratio
       real(rk)       :: T0    = 296_rk  ! Reference T, K
       real(rk)       :: P0    = 1.0_rk  ! Reference P, bar
  -    character(len=cl) :: name = "NA"         ! Broadener number of quadrature points
  +    character(len=cl) :: name = "NA"         ! Broadener name
       character(len=cl) :: filename=""  ! File name with QN-dependent broadening parameters
       character(len=cl) :: model="const"  ! Broadening model
       real(rk),pointer  :: gammaQN(:,:) ! Voigt parameter gamma as a function of QN
  
positive.code:             !
          case ('PHOEN')
              !
              call do_gf_oscillator_strength_Phoenix(i,ichunk,iso,nswap,nlines,energies,Jrot,gtot,ilevelf_ram,ileveli_ram,acoef_ram,abscoef_ram,gfunit)
              !
              call TimerStop('Calc')
              cycle loop_tran
  
negative.code:    use  input
     use VoigtKampff_module
     !
     integer(ik) :: info,ipoint,nlevels,i,itemp,enunit,tunit,sunit,bunit,j,j0,ilevelf,ileveli,indexi,indexf,iline,maxitems,gfunit
     integer(ik) :: indexf_,indexi_,kitem,nlines,ifilter
     real(rk)    :: beta,ln2,ln22,dtemp,dfreq,temp0,beta0,intband,dpwcoef,tranfreq,abscoef,halfwidth0,tranfreq0
     real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,jf,ji,acoef,j0rk,gfcoef
     real(rk)    :: acoef_,cutoff
     integer(ik) :: Jmax,Jp,Jpp,J_,Noffset,Nspecies_,Nvib_states,ivib1,ivib2,ivib,JmaxAll
     real(rk)    :: gamma_,n_,gamma_s,ener_vib,ener_rot
     character(len=cl) :: ioname
     !
  