Add support for product dims

[AHsu98]

We probably want to add support for products of dims as var_groups, so that we could, for example, have intercept vary by both year and region, where it varies smoothly by year, and has a different temporal pattern for each region.

This is somewhat lower priority right now I think, but should be implemented and used at some point.