<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Oops, I spoke too fast. <br>
I can make a slice but if I try a projection, the same kind of error
occurs: <a class="moz-txt-link-freetext" href="http://paste.yt-project.org/show/7QiwjVoYXKYxCn3Iu2u9/">http://paste.yt-project.org/show/7QiwjVoYXKYxCn3Iu2u9/</a><br>
<br>
This is the script I am using:<br>
<br>
<address>-------------------------------------------------------------------------</address>
<address>def pf_compare(pf1, pf2):</address>
<address> def _NewField(field, data):</address>
<address> data2 = pf2.h.grids[data.id - data._id_offset]</address>
<address> return
abs(data2["TotalEnergy"]-data["TotalEnergy"])/data["TotalEnergy"]</address>
<address> pf1.field_info.add_field("CompareEnergy", function =
_NewField, validators = [ValidateGridType()])</address>
<address><br>
pf1 = load('helix/DD0000/CE0000')</address>
<address>pf2 = load('parallel/DD0000/CE0000')</address>
<address><br>
#
Comparison
</address>
<address>pf_compare(pf1, pf2)</address>
<address>region = pf1.h.region([0.5, 0.5, 0.5], [0.0, 0.0, 0.0],
[1.0, 1.0, 1.0])</address>
<address>compare_energy = region["CompareEnergy"]</address>
<address><br>
pc = PlotCollection(pf1)</address>
<address>pc.add_slice("CompareEnergy",2,center=[0.5, 0.5, 0.5])</address>
<address>pc.add_projection("CompareEnergy",2)</address>
<address>pc.set_lim((0.4, 0.6, 0.4, 0.6))</address>
<address>pc.save()</address>
<address>-------------------------------------------------------------------------</address>
<br>
Any idea?<br>
Cheers,<br>
<br>
JC<br>
<br>
On 5/1/12 3:35 PM, Jean-Claude Passy wrote:
<blockquote cite="mid:4FA03B1D.2070501@uvic.ca" type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
Thanks Britton, it works perfectly.<br>
<br>
Cheers,<br>
<br>
JC</blockquote>
</body>
</html>