The 13th edition's treatment of WF torsion cannot be reasonably implemented in any software program. The problem as I see it is an inconsistency with the rest of the code. As written, the code would have you move from a moment capacity of Mp to a moment capacity of My when you have even the smallest iota of torque in the member. You can't have that type of large discontinuity in your combined stress equations.
In RISA, we consider torsion in wide flanges for the 13th edition. But, we use the recomendations of the AISC Design Guide on Torsion. Those recomendations seem more rational to me. The recommendations consist of calculating torsional shear stresses and torsional warping normal stresses and adding them to the existing flexural shear stresses and weak axis bending stresses. Done in this manner this is not any discontinuity between the combined stress equations with or without torsion.
Also, FWIW, I think the reason why the code is so conservative with torsion of wide flanges is because it is assuming that whatever analysis was done is ignoring the torsional warping effects on the wide flange members. Since these can be significant, it makes sense to be conservative.
Therefore, ewhen a program like STAAD (which I don't believe copnsiders torsional warping effects), I would always recommend running a side calculation related to warping effects.... at least until you are comfortable with what your program does related to warping.