Added null() static construction method.
ImpesTPFAAD now takes a linear solver parameter.

Also work in progress on adding well bhp to the system.
Work in progress on adding well eqns.