Hack for setting heterogeneous porosity/permeability maps

This is far from being a good solution, but since I have no time to re-factor the entire hierarchy of classes, I'm adding non-const versions for allowing the user to manually set the maps as in:

```cpp
perm = props.permeability();
perm[i] = 10*milli*darcy;
```