Roy,

What's up with making this const?  We have some code that depends on elem's being non-const... like changing their subdomain ids after reading the mesh.

So what gives?  I see your message that says that this is "correct"... but I'm not convinced...

Derek