I'm dealing with some athena data right now, and running into some issues with defining regions.
In : center=[0.5,0.5,0.5]
In : lc = [0,0,0.4]
In : rc=[1.0,1.0,0.6]
In : region = pf.h.region(center, lc, rc)
In : region.volume()
I'm able to get at region.quantities, but trying to make a projection plot fails.
Interestingly, or maybe not, if I try to do something similar with data that has its domain centered on the origin, the volume is finite and everything I've tried with plotting works fine. I guess athena support is still a bit preliminary, but I thought I'd throw this out there.