anchor.diff:                    ! apply vibrational population
                     read(quantum_numbers(QN%dens_col-4,ileveli),*) ndensity
                     abscoef=cmcoef*acoef*gtot(ilevelf)*exp(-c2/temp*ener_rot)*ndensity*&
  -                          (1.0_rk-exp(-c2/ener_rot*tranfreq))/(tranfreq**2*partfunc)
  +                          (1.0_rk-exp(-c2/temp*ener_rot*tranfreq))/(tranfreq**2*partfunc)
                     ! split into a product of vib and rot parts 
                  else
                     abscoef=cmcoef*acoef*gtot(ilevelf)*exp(-c2/temp*ener_rot)*exp(-c2/temp_vib*ener_vib)*&
  
positive.code:   type(HitranErrorT),target,save :: HITRAN_E(HITRAN_max_ierr),HITRAN_S(HITRAN_max_ierr)
    type(HitranErrorT),target,save :: HITRAN_Air,HITRAN_Self,HITRAN_N,HITRAN_Delta
    !
    integer(ik),allocatable,save  :: gamma_idx(:,:) !Used for indexing the gamma for the fast_voigt
    real(rk),allocatable,save  :: gamma_comb(:,:) !Used for indexing the gamma for the fast_voigt
    logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran_do = .false.,  histogramJ = .false., &
  
negative.code:               call readi(QN%dens_col)
                !
                vibpopulation_do = .true.
                vibtemperature_do  = .true.
                !
              case ("VIB","VIBRATIONAL")
                !
  