anchor.diff:               !
                select case (trim(cutoff_model))
                   !
  +             case ("EXP")
  +                !
  +                continue
  +                !
                case ("EXP-WEAK")
                   !
  +                abscoef_ref =cmcoef*acoef*gtot(ilevelf)*exp(-beta_ref*energyi)*(1.0_rk-exp(-beta_ref*tranfreq))/&
  +                             (tranfreq**2*partfunc_ref)
  +                !
                   cutoff = S_crit*exp(-tranfreq/nu_crit)
                   !
  -                if (trim(proftype(1:5))/='TRANS'.and.(abscoef>=cutoff.or.energyi<=ener_cut)) then
  +                if ((abscoef>=cutoff.or.abscoef_ref>=cutoff).or.energyi<=ener_cut) then
                     cycle loop_swap
                   endif
                   !
                case ("EXP-STRONG")
                   !
  +                abscoef_ref =cmcoef*acoef*gtot(ilevelf)*exp(-beta_ref*energyi)*(1.0_rk-exp(-beta_ref*tranfreq))/&
  +                             (tranfreq**2*partfunc_ref)
  +                !
                   cutoff = S_crit*exp(-tranfreq/nu_crit)
                   !
  -                if (trim(proftype(1:5))/='TRANS'.and.(abscoef<cutoff.and.energyi>ener_cut)) then
  +                if ((abscoef<cutoff.and.abscoef_ref<cutoff).and.energyi>ener_cut) then
                     cycle loop_swap
                   endif
                   !
  
positive.code:           case("HITRAN")
               write(out,"(10x,/'Stick pectra of type ',a,' with the HITRAN cut-off model, Scrit = ',e18.5,' nu_crit = ',f16.4)") &
               trim(proftype),S_crit,nu_crit
            case("EXP","EXP-WEAK","EXP-STRONG")
               write(out,"(10x,/'Strong/weak lines ',a,' with the EXP cut-off model, Scrit = ',e18.5,' nu_crit = ',f16.4)") &
               trim(proftype),S_crit,nu_crit
            case default
               write(out,"(10x,/'Stick spectra of type ',a,' stronger than ',e18.5)") trim(proftype),thresh
            end select 
  
negative.code:               exit loop_tran
             endif
             !
          elseif (super_Einstein_do) then ! super_energies_do
             !
             iswap_ = 0
             !
  