anchor.diff:              !
            endif
             !
  -       case ("THRESHOLD")
  -          !
  -          call readf(thresh)
  +       case ("THRESHOLD","CUTOFF")
  +          !
  +          if (Nitems>2) then 
  +            !
  +            call readu(w)
  +            !
  +            if (trim(w)/="HITRAN") call report ("Expected: either HITRAN with the a cutoff Intensity or just one number = cutoff "//trim(w),.true.)
  +            cutoff_model = "HITRAN"
  +            call readf(thresh)
  +            !
  +            S_crit = thresh
  +            !
  +            write(out,"('HITRAN intensity cut-off model is used, see HITRAN 2012 paper')")
  +            !
  +            call readu(w)
  +            !
  +            if (trim(w)/="NU_CRIT".and.trim(w)/="NUCRIT") call report ("Expected: keyword = NU_CRIT with the switching frquency"//trim(w),.true.)
  +            !
  +            call readf(nu_crit)
  +            !
  +          else
  +            !
  +            call readf(thresh)
  +            !
  +          endif
             !
          case ("ENERMAX")
             !
  
positive.code:   integer(ik)   :: GNS=1,npoints=1001,nchar=1,nfiles=1,ipartf=0,verbose=3,ioffset = 10,iso=1
    real(rk)      :: temp=298.0,partfunc=-1.0,freql=-small_,freqr= 20000.0,thresh=1.0d-90,halfwidth=1e-2,meanmass=1.0,maxtemp=10000.0
    real(rk)      :: voigt_gamma = 0.05, voigt_n = 0.44, offset = -25.0, pressure = 1.0_rk 
    real(rk)      :: enermax = 1e6, abscoef_thresh = 1.0d-50
    integer(ik)   :: nquad = 20      ! Number of quadrature points 
    !
    character(len=cl) :: specttype="absorption",enrfilename="none",intfilename(nfiles_max),proftype="DOPPL",output="output"
  
negative.code:         !
          Nintens(ilog) = Nintens(ilog)+1
          !
          !
          ! HITAN cutoff model as an intensity threshold
          if (trim(cutoff_model)=="HITRAN") then
             if (tranfreq<=nu_crit) then
               cutoff = S_crit*tranfreq/nu_crit*tanh(0.5_rk*c2/temp*tranfreq)
             else
               cutoff = S_crit
             endif
          else
             cutoff = thresh
          endif
          !
          if (abscoef<cutoff) cycle
          !
          ! if only stick spectrum needed
          if (proftype(1:5)=='STICK') then
  