Extend the spatial database implementation to allow specification of time dependence. Starting case could be x,y,z,t. * Need to adjust interface to allow time * PyLith - Would we need an 4D auxiliary field?