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