After a short motivation for non-minimal SUSY models I’ll introduce a
framework to study SUSY models beyond the MSSM using the Mathematica
package SARAH. The main part of the talk will focus on two aspects.
First, the automatized calculation of two-loop Higgs masses using
SARAH/SPheno for a large variety of SUSY models with the same
precision public codes provide today for the MSSM is explained. I will
also discuss results for important corrections  in the NMSSM in the
large $\lambda$-limit which have been calculated for the first time
with this setup. Second, the automatized calculation of quark and
lepton flavour violating observables at one-loop using the FlavorKit
interface between FeynArts/FormCalc, SARAH, and SPheno is presented.
I’ll also discuss at short examples how the user can implement easily
new observables and -- if there is time -- also new Wilson
coefficients.