[yt-users] fields in ghost zones

Dear all, in order to set up my simulation, I need to access the values of certain fields in the ghost zones. If I do something like: pf = load(data) region = pf.h.region(...) x = region[field] x contains the values of field for the physical grid only. Is there a way I can get the same 1D-array but with the ghost zones included as well ? Thanks for your help, Jean-Claude