Pulled in srb-forces_us from scf_shirley.

Tested for ZnS-f in serial and parallel.  It could use more testing,
but most of the qe-tests are force-less.

Note that this is a rather inefficient way to compute the forces.
There is a better way, based on an auxillary basis for projector
derivatives, which will be implemented later.