`SubDomainSet` should be overridable at runtime

[EdCaunt]

It would be good if a SubDomainSet could be passed when running an operator, along the lines of op.apply(sds=new_sds).