Wow, thanks a million, this worked just fine and I think I understood at least parts of the method :)<br><br>Cheers<div><br>Thomas</div><div><br></div><div><br></div><div><span class="Apple-style-span" style="font-family: &#39;Times New Roman&#39;; font-size: medium; "><i><span class="Apple-style-span" style="font-style: normal; "><pre>

On Tue, 2010-10-19 at 14:20 +0200, Nadav Horesh wrote:</pre><pre><span class="Apple-style-span" style="font-family: &#39;Times New Roman&#39;; white-space: normal; "><i><span class="Apple-style-span" style="font-style: normal; "><pre style="display: inline !important; ">

&gt; Of course there is an (at least one) error:</pre></span></i></span></pre></span></i></span></div><div><span class="Apple-style-span" style="font-family: &#39;Times New Roman&#39;; font-size: medium; "><pre>&gt; the line should be:</pre>

</span><span class="Apple-style-span" style="font-family: &#39;Times New Roman&#39;; font-size: medium; "><pre>&gt; 
&gt; XYZ =
&gt; np.mgrid[lwrbnd[0]:uprbnd[0]:shape[0]*1j,lwrbnd[1]:uprbnd[1]:shape[1]*1j, lwrbnd[2]:uprbnd[2]:shape[2]*1j]


&gt; On Tue, 2010-10-19 at 14:10 +0200, Nadav Horesh wrote:
&gt; &gt;<i> You can aid mgrid, riughy as the follows (I may have mistakes, but the
</i>&gt; &gt;<i> direction should be clear):
</i>&gt; &gt;<i> 
</i>&gt; &gt;<i> def transform_3d_data_(field,lwrbnd,uprbnd): 
</i>&gt; &gt;<i>   shape = field.shape 
</i>&gt; &gt;<i>   XYZ = np.mgrid[lwrbnd[0]:uprbnd[0]:shape[0],
</i>&gt; &gt;<i> lwrbnd[1]:uprbnd[1]:shape[1], lwrbnd[2]:uprbnd[2]:shape[2]] 
</i>&gt;<i> </i><i><pre style="display: inline !important; ">&gt; </pre></i><i>  vectors = fields.reshape(-1,3) </i></pre></span><span class="Apple-style-span" style="font-family: &#39;Times New Roman&#39;; font-size: medium; "><pre>

&gt; &gt;<i>   np.savetxt(np.hstack((XYZ.reshape(3,-1).T, vectors)))
</i>&gt; &gt;<i> 
</i>&gt; &gt;<i> 
</i>&gt; &gt;<i>   Nadav </i></pre></span></div>