Oops, I spoke too fast. I can make a slice but if I try a projection, the same kind of error occurs: http://paste.yt-project.org/show/7QiwjVoYXKYxCn3Iu2u9/ This is the script I am using: ------------------------------------------------------------------------- def pf_compare(pf1, pf2): def _NewField(field, data): data2 = pf2.h.grids[data.id - data._id_offset] return abs(data2["TotalEnergy"]-data["TotalEnergy"])/data["TotalEnergy"] pf1.field_info.add_field("CompareEnergy", function = _NewField, validators = [ValidateGridType()]) pf1 = load('helix/DD0000/CE0000') pf2 = load('parallel/DD0000/CE0000') # Comparison pf_compare(pf1, pf2) region = pf1.h.region([0.5, 0.5, 0.5], [0.0, 0.0, 0.0], [1.0, 1.0, 1.0]) compare_energy = region["CompareEnergy"] pc = PlotCollection(pf1) pc.add_slice("CompareEnergy",2,center=[0.5, 0.5, 0.5]) pc.add_projection("CompareEnergy",2) pc.set_lim((0.4, 0.6, 0.4, 0.6)) pc.save() ------------------------------------------------------------------------- Any idea? Cheers, JC On 5/1/12 3:35 PM, Jean-Claude Passy wrote:
Thanks Britton, it works perfectly.
Cheers,
JC