anchor.diff:               !
                iswap_ = iswap_ + 1
                !
  -             abscoef_ram(iswap_)=abscoef*exp(-beta    *energyi)*(1.0_rk-exp(-beta    *tranfreq))*partfunc_ref/&
  -                                        (exp(-beta_ref*energyi)*(1.0_rk-exp(-beta_ref*tranfreq))*partfunc   ) 
  +             !abscoef_ram(iswap_)=abscoef*exp(-beta    *energyi)*(1.0_rk-exp(-beta    *tranfreq))*partfunc_ref/&
  +             !                           (exp(-beta_ref*energyi)*(1.0_rk-exp(-beta_ref*tranfreq))*partfunc   ) 
  +             !
  +             ! gtot is set to 1
  +             !
  +             abscoef_ram(iswap_)=cmcoef*abscoef*exp(-beta*energyi)*(1.0_rk-exp(-beta*tranfreq))/(tranfreq**2*partfunc)
                !
                nu_ram(iswap_) = tranfreq
                !
  
positive.code:         case ("SPECTRA")
            !
            spectra_do = .true.
            !
          case ("SUPER-ENERGIES")
            !
  
negative.code:               !
                S_crit = thresh
                !
                !write(out,"('EXP intensity cut-off model I0*exp(-nu/nucrit)')")
                !
                call readu(w)
                !
  