Generalize screen output of the rod and continuum forces to more than one object

The entire primal (i.e., Dirichlet-to-Neumann) should now run for
arbitrary numbers of rods and continua.

[[Imported from SVN: r6844]]