  We describe an extension of the SOFTSUSY spectrum calculator to include
two-loop supersymmetric QCD (SUSYQCD) corrections of order
$\mathcal{O}(\alpha_s^2)$ to gluino and squark pole masses, either in the
minimal supersymmetric standard model (MSSM) or the next-to-minimal
supersymmetric standard model (NMSSM). This document provides an overview of
the program and acts as a manual for the new version of SOFTSUSY, which
includes the increase in accuracy in squark and gluino pole mass predictions.
