anchor.diff:          x0 = sqrt(ln2)/halfwidth*dfreq*0.5_rk
           !
           !   half width for Doppler profiling
  -        if (proftype(1:5)=='VOIGT'.or.proftype(1:3)=='PSE') then
  +        if (proftype(1:3)=='VOI'.or.proftype(1:3)=='PSE') then
             if (tranfreq<small_) cycle
             halfwidth0=dpwcoef*tranfreq
             x0 = sqrt(ln2)/halfwidth0*dfreq*0.5_rk
  
positive.code:                !
                 tranfreq_i = freq(ipoint)
                 !
                 intens(ipoint)=intens(ipoint)+voigt_humlicek(tranfreq_i,tranfreq,halfwidth0,halfwidth)*abscoef
                 !
              enddo
  
negative.code:     real(rk)       :: ratio = 0.9_rk  ! Voigt parameter gamma
      real(rk)       :: T0    = 298_rk  ! Reference T, K
      real(rk)       :: P0    = 1.0_rk  ! Reference P, bar
      character(len=cl) :: name         ! Broadener number of quadrature points 
    end type speciesT
    !
    type(selectT) :: upper,lower
  