Subdomain scalar realloc

- 1st patch is there for RELAP-7 (I need a programmatic way of turning that flag on/off)
- 2nd patch fixes the assumption that coupling matrix is symmetrical (it could be part of the third patch if needed)
- 3rd patch is the ham ;-)