Hi Jean-Claude, There are a couple aspects to this. The first is that Enzo doesn't output the ghost zones -- so any ghost zones handled inside yt are generated by yt. Were Enzo to output the ghost zones, we would probably be able to handle this, but it does not. Derived fields can depend on the generation of ghost zones, but keep in mind that these are ghost zones generated by yt. These ghost zones are constructed in a similar manner to how Enzo generates them, but there may be minor differences. You can manually inspect ghost zones on a *grid* by calling retrieve_ghost_zones on that grid. If you could tell us a bit more about your goal, maybe we could help out a bit more? -Matt On Wed, Jun 2, 2010 at 9:32 AM, Jean-Claude Passy <jcpassy@gmail.com> wrote:
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
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org