anchor.diff:               ji = jrot(ileveli)
                !
                tranfreq = energyf-energyi
  +             !
  +             tranfreq = tranfreq + lineshift
                !
                tranfreq0 = tranfreq
                !
  
positive.code:   integer(ik)   :: GNS=1,npoints=1001,nchar=1,nfiles=1,ipartf=0,verbose=2,ioffset = 10,iso=1,imolecule =-1
    real(rk)      :: temp=298.0,partfunc=-1.0,partfunc_ref=-1.0,freql=-small_,freqr= 200000.0,thresh=1.0d-70
    real(rk)      :: voigt_gamma = 0.05, voigt_n = 0.44, offset = 25.0, pressure = 1.0_rk,halfwidth=1e-2,meanmass=1.0,maxtemp=10000.0
    real(rk)      :: enermax = 1e7, abscoef_thresh = 1.0d-50, abundance = 1.0d0, gf_factor = 1.0d0
    real(rk)      :: S_crit = 1e-29      ! cm/molecule, HITRAN cut-off paramater
    real(rk)      :: nu_crit = 2000.0d0  ! cm-1, HITRAN cut-off paramater
  
negative.code:   type speciesT ! broadener
      !
      real(rk)       :: N = 0.5_rk      ! Voigt parameter N
      real(rk)       :: gamma = 0.0_rk  ! Voigt parameter gamma
      real(rk)       :: ratio = 1.0_rk  ! Ratio
      real(rk)       :: T0    = 296_rk  ! Reference T, K
      real(rk)       :: P0    = 1.0_rk  ! Reference P, bar
  