There is a generic mechanism for surface loads in form of the
NeumannEnergy class (to be replaced soon(TM) by something
even more generic).  This is what should be used, instead
of ad hoc implementations such as the one in SimoFoxEnergy.