anchor.diff:                   !
                  enddo
                  !
  +               ! filter based on the uncertainty 
  +               if (uncertainty_filter_do) then
  +                  !
  +                  ! reconstruct the uncertainty from 
  +                  read(quantum_numbers(QN%unc_col-4,ilevelf),*) unc_f
  +                  read(quantum_numbers(QN%unc_col-4,ileveli),*) unc_i
  +                  !
  +                  if (unc_f>unc_threshold.or.unc_i>unc_threshold) cycle loop_swap
  +                  !
  +               endif
  +               !
                endif
                !
                if (vibtemperature_do) then
  
positive.code:     !
    end type CutoffT
  
    type QNT ! broadener
      !
      integer(ik) :: Kcol=1                ! Column with K 
      integer(ik) :: statecol=7            ! State column
  
negative.code:                 call readi(upper(ifilter)%i) ; upper(ifilter)%i = upper(ifilter)%i - 4
                  call reada(upper(ifilter)%mask)
                  !
                case('UNC','UNCERTAINTY')
                  !
                  ! no need for the counter to change here
                  ifilter = ifilter - 1
                  !
                  uncertainty_filter_do = .true.
                  !
                  call readf(unc_threshold)
                  !
                case default
                  !
                  call report ("Unrecognized unit name "//trim(w),.true.)
  