Complete rewrite, so that stabilized and unstabilized are now in one code. All parameters are available via options, which includes viscosity.