anchor.diff:    !
     logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran_do = .false.,  histogramJ = .false., stick_hitran = .false.
     logical :: microns = .false.
  +  logical :: use_resolving_power = .true.  ! using resolving for creating the grid
     !
     contains
     !
  
positive.code:   real(rk)      :: enermax = 1e6, abscoef_thresh = 1.0d-50, abundance = 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 
    character(len=cl) :: cutoff_model = "NONE"
    integer(ik)   :: nquad = 20      ! Number of quadrature points 
    !
  
negative.code:             !
              if (ipoint<1.or.ipoint>npoints) cycle 
              !
            elseif (use_resolving_power) then 
              !
              ipoint =  nint(real((log(tranfreq)-log(freql))/resolving_f,rk))+1
              !
              !nint( ( 10000.0_rk/tranfreq-freql )/dfreq )+1
              !
              if (ipoint<1.or.ipoint>npoints) cycle 
              !
            else 
              !
              ipoint =  max(nint( ( tranfreq-freql)/dfreq )+1,1)
              !
  