anchor.diff:         case ("ENERMAX")
             !
             call readf(enermax)
  -          !
  -          ! energy threshold based on the uncertainty 
  -          !
  -       case('UNCERTAINTY_MAX','UNC_MAX','UNC')
  -          !
  -          uncertainty_filter = .true.
  -          !
  -          call readf(unc_threshold)
             !
          case('GAUSSIAN','GAUSS','DOPPL','DOPPLER','RECT','BOX','BIN','STICKS','STICK','GAUS0','DOPP0',&
               'LOREN','LORENTZIAN','LORENTZ','MAX','VOIGT','PSEUDO','PSE-ROCCO','PSE-LIU','VOI-QUAD','PHOENIX',&
  
positive.code:   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., &
               vibpopulation_do = .false., super_energies_do = .false., super_Einstein_do = .false., &
               uncertainty_filter = .false.
    logical :: lineprofile_do = .false., use_width_offset = .false.
    logical :: microns = .false.
    logical :: use_resolving_power = .false.  ! using resolving for creating the grid
  
negative.code:    real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,jf,ji,acoef,j0rk,gfcoef
     real(rk)    :: acoef_,cutoff,ndensity,abscoef_ref
     integer(ik) :: Jmax,Jp,Jpp,Noffset,Nspecies_,Nvib_states,ivib1,ivib2,ivib,JmaxAll,imin
     real(rk)    :: gamma_,n_,gamma_s,ener_vib,ener_rot,J_,pf_1,pf_2,t_1,t_2,unc_i,unc_f
     character(len=cl) :: ioname,intname
     !
     real(rk),allocatable :: freq(:),intens(:),jrot(:),pf(:,:),energies(:),Asum(:),weight(:),abciss(:),bnormq(:)
  