We present an extension to the Mathematica package SARAH which allows for Higgs mass calculations at the two-loop level in a wide range of supersymmetric (SUSY) models beyond the MSSM. These calculations are based on the effective potential approach and include all two-loop corrections which are independent of electroweak gauge couplings. For the numerical evaluation Fortran code for SPheno is generated by SARAH. This allows the prediction of the Higgs mass in more complicated SUSY models with the same precision that most state-of-the-art spectrum generators provide for the MSSM.