Rob Low? Really--Probably and I mean probably the spatio-temporal
part of your spatio can be assumed to go to zero. Then use the same time constant on each equation reducing it to a 1st ODE and then use MathCads solver. I make this assumption because the advection is usually uniform with the solution. If this is not the case you have a non-linear dynamic model and you must isolate the dampening terms.
