anchor.diff:                !
                 S_crit = thresh
                 !
  -              write(out,"('EXP intensity cut-off model I0*exp(-nu/nucrit)')")
  +              !write(out,"('EXP intensity cut-off model I0*exp(-nu/nucrit)')")
                 !
                 call readu(w)
                 !
  
positive.code:              !
               select case (trim(cutoff_model))
                  !
               case ("EXP-WEAK")
                  !
                  cutoff = S_crit*exp(-tranfreq/nu_crit)
                  !
                  if (trim(proftype(1:5))/='TRANS'.and.(abscoef>=cutoff.or.energyi<=ener_cut)) then
                    cycle loop_swap
                  endif
                  !
               case ("EXP-STRONG")
                  !
                  cutoff = S_crit*exp(-tranfreq/nu_crit)
                  !
                  if (trim(proftype(1:5))/='TRANS'.and.(abscoef<cutoff.and.energyi>ener_cut)) then
                    cycle loop_swap
                  endif
                  !
  
negative.code:         case ("SPECTRA")
            !
            spectra_do = .true.
            !
          case ("SUPER-EINSTEIN")
            !
            super_Einstein_do = .true.
            !
          case ("SUPER-ENERGIES")
            !
  