Improved mixing_zone functions and added plot of h and Xi.

mixing_zone was over-built, spitting out boundary regions and
plotting instead of just the boundary width.  Now it leaves
visualization to the driver and just outputs h and Xi (from Cabot 2006).

Added plot of h, Xi for series input in the driver.