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:        !
         if (verbose>=2) then
            !
            write(out,"(10x,/'Stick pectra of type ',a,' stronger than ',e18.5)") trim(proftype),thresh
            write(out,"(10x,'Range = ',f18.7,'-',f18.7)") freql,freqr
            write(out,"(10x,'Temperature = ',f18.7)") temp
            write(out,"(10x,'Partition function = ',f17.4)") partfunc
  
negative.code:    integer(ik) :: indexf_,indexi_
     real(rk)    :: beta,ln2,ln22,dtemp,dfreq,temp0,beta0,intband,dpwcoef,x0,tranfreq,tranfreq_i,abscoef,dfreq_,xp,xm,de,lor,b,lor2,dfreq_2,halfwidth0,dnu_half,dxp
     real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,elow,jf,ji,acoef,j0rk,bnorm,f,eta1,eta2,intens1,intens2,Va0,gammaV,a,wg,d
     real(rk)    :: sigma,alpha,gamma,y,x1,x2,voigt_,dx2,xi,L1,L2,acoef_,cutoff
     integer(ik) :: Jmax,Jp,Jpp,Kp,Kpp,J_,ispecies,alloc_p
     real(rk)    :: gamma_,n_
     character(len=cl) :: ioname
  