anchor.diff:    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., &
  -             stick_hitran = .false.,stick_oxford = .false.,vibtemperature_do = .false., spectra_do = .false.
  +             stick_hitran = .false.,stick_oxford = .false.,vibtemperature_do = .false., spectra_do = .false., vibpopulation_do = .false.
     logical :: lineprofile_do = .false., use_width_offset = .false.
     logical :: microns = .false.
     logical :: use_resolving_power = .false.  ! using resolving for creating the grid
  
positive.code:                !
                 if (.not.vibtemperature_do) then
                    abscoef=cmcoef*acoef*gtot(ilevelf)*exp(-beta*energyi)*(1.0_rk-exp(-beta*tranfreq))/(tranfreq**2*partfunc)
                 else 
                    ! split into a product of vib and rot parts 
                    abscoef=cmcoef*acoef*gtot(ilevelf)*exp(-c2/temp*ener_rot)*exp(-c2/temp_vib*ener_vib)*&
                            (1.0_rk-exp(-c2/temp_vib*tranfreq))/(tranfreq**2*partfunc)
                 endif 
  
negative.code:                !
                 if (.not.vibtemperature_do) then
                    abscoef=emcoef*acoef*gtot(ilevelf)*exp(-beta*energyf)*tranfreq/(partfunc)
  
                 elseif(vibpopulation_do) then
                    ! apply vibrational population
                    read(quantum_numbers(QN%dens_col-4,ileveli),*) ndensity
                    !
                    abscoef=emcoef*acoef*gtot(ilevelf)*exp(-c2/temp*ener_rot)*ndensity*&
                            tranfreq/(partfunc)
                    !
                 else 
                    ! split into a product of vib and rot parts 
                    abscoef=emcoef*acoef*gtot(ilevelf)*exp(-c2/temp*ener_rot)*&
                            exp(-c2/temp_vib*ener_vib)*tranfreq/(partfunc)
  
                    !
                 endif 
                 !
  