anchor.diff:                  ilevelf = ilevelf_ram(iswap)
                   ileveli = ileveli_ram(iswap)
                   !
  -                cutoff = S_crit*exp(-tranfreq/nu_crit)
  +                !
  +                select case (trim(cutoff_model))
  +                   !
  +                case ("EXP","EXP-WEAK","EXP-STRONG")
  +                  cutoff = S_crit*exp(-tranfreq/nu_crit)
  +                case ("GRID")
  +                  inu  = min(max(1,nint(tranfreq/ ( (freqr-freql)/cutoffs%npoints ) )),cutoffs%npoints)
  +                  cutoff = cutoffs%intensity(inu)*cutoffs%factor
  +                end select 
                   !
                   ileveli = ileveli_ram(iswap)
                   energyi = energies(ileveli)
                   !
  -                if (abscoef<cutoff.and.energyi>ener_cut) then
  +                if (abscoef<cutoff.and.energyi>cutoffs%enercut) then
                       write(wunit,my_fmt) ilevelf,ileveli,acoef !,tranfreq
                   else
                       write(sunit,my_fmt) ilevelf,ileveli,acoef,tranfreq
  
positive.code:    !
     real(rk)    :: hitran_Tref = 296_rk
     integer(ik) :: info,ipoint,ipoint_,nlevels,i,itemp,enunit,tunit,sunit,wunit,bunit,pfunit,j,j0,ilevelf,ileveli,indexi,indexf
     integer(ik) :: indexf_,indexi_,kitem,nlines,ifilter,k,igrid,maxitems,iline
     real(rk)    :: beta,ln2,ln22,dtemp,dfreq,temp0,beta0,intband,dpwcoef,tranfreq,abscoef,halfwidth0,tranfreq0,delta_air,beta_ref
     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
     character(len=cl) :: ioname
     !
     real(rk),allocatable :: freq(:),intens(:),jrot(:),pf(:,:),energies(:),Asum(:),weight(:),abciss(:),bnormq(:)
     integer(ik),allocatable :: gtot(:),indices(:)
  
negative.code:   logical :: accepted_work = .false.
    logical :: completed_work = .true.
    logical :: all_done = .false.
    !
    type(VoigtKampffCollection),save :: fast_voigt
    
    
    !
  