Shifted van-Genuchten curve added.  Fixes #3350

This is for 2-phase flow.  It may make the Jacobian nonsingular in
2phase models when one phase disappears and residual saturation is
zero.  Time will tell whether this is useful.

Included are doco and tests
Shifted van-Genuchten curve added.  Fixes #3350

This is for 2-phase flow.  It makes the Jacobian nonsingular
in 2phase models when one phase disappears and the residual
saturation is zero.  In limited testing this appears to give
good convergence characteristics.

Included are doco and tests.

Also, I have changed RichardsSeff::dseff and d2seff so they
accept references to the answers, rather than providing the
answers as return values.  According to jwpeterson this might
speed up computations by not having to do extra copying.
Since I have so many RichardsSeff type of classes, this required
a lot of minor changes
Cleaning up the warehouses, making most methods const, avoiding vector copies
closes #3358, #3368
Merge pull request #3377 from permcody/warehouse_cleanup

Cleaning up the warehouses
Removed vector copy
(refs #3378)
Merge pull request #3386 from aeslaughter/output_warehouse-3378

Removed a vector copy from OutputWarehouse
How 'bout this fix, eh?
closes #3388
Merge pull request #3390 from permcody/max_parallel

How 'bout this fix, eh?
Merge pull request #3353 from WilkAndy/shiftedVG_3350

Shifted van-Genuchten curve added.  Fixes #3350
Doco for VGshifted added into user.tex.

Fixes #3393

This document had fallen a few weeks behind the code
so was substantially revised to with my current thoughts.
Merge pull request #3395 from WilkAndy/vgshifted_doco_3393

Doco for VGshifted added into user.tex.
Re-factored output to the screen, including MultiApp indenting
(refs #3321)
Allow printing in const functions
Get rid of one more ostringstream
InitialSetup messages now printing
(refs #3321)
Added test for const method output (#3321)