<div dir="ltr">Hi there,<div><br></div><div style>I'm dealing with some athena data right now, and running into some issues with defining regions. </div><div style><br></div><div style>an example, with domain_left_edge = [0, 0, 0], and right_edge = [1, 1, 1]</div>

<div style><br></div><div style><div>In [18]: center=[0.5,0.5,0.5]</div><div>In [19]: lc = [0,0,0.4]</div><div>In [20]: rc=[1.0,1.0,0.6]</div><div>In [21]: region = pf.h.region(center, lc, rc)</div><div>In [22]: region.volume()</div>

<div>Out[22]: -inf</div><div><br></div><div style>I'm able to get at region.quantities, but trying to make a projection plot fails.  </div><div><br></div><div style>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.</div>

<div style><br></div><div style>best,</div><div style>Nick</div></div></div>