<div dir="ltr">Hello,<div><br></div><div>I&#39;d like to point out also my python library/script called &quot;rugutils&quot;. This one started off inspired by David&#39;s scidata (and pygraph), but has developed into a set of tools to look in a quick and hassle-free way at Carpet data. It&#39;s written in python 3.5</div><div><br></div><div><a href="https://bitbucket.org/fguercilena/rugutils">https://bitbucket.org/fguercilena/rugutils</a><br></div><div><a href="https://docs.einsteintoolkit.org/et-docs/Analysis_and_post-processing">https://docs.einsteintoolkit.org/et-docs/Analysis_and_post-processing</a><br></div><div><br></div><div>Federico</div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-05-17 19:44 GMT+02:00 Roland Haas <span dir="ltr">&lt;<a href="mailto:rhaas@aei.mpg.de" target="_blank">rhaas@aei.mpg.de</a>&gt;</span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Frank, Michael,<br>
<br>
@Michael: both of what Frank or David describe should work. If you know<br>
how to compile the VisIt reader yourselv, you can also try the copy of<br>
CarpetHDF5 at <a href="https://bitbucket.org/rhaas80/carpethdf5" rel="noreferrer" target="_blank">https://bitbucket.org/rhaas80/carpethdf5</a> either the<br>
&quot;multifile&quot; branch (which is the default one that has all features) or<br>
the &quot;for_VisIt&quot; branch which has fewer features and is what I wanted at<br>
one point propose for inclusion in VisIt. Both will let you<br>
create .visit files listing your hdf5 files and VisIt will treat them a<br>
a single database. You will have to create .visit files that look like<br>
this:<br>
<br>
!NBLOCKS 3<br>
output-0000/rho.file_0.h5<br>
output-0000/rho.file_1.h5<br>
output-0000/rho.file_2.h5<br>
output-0001/rho.file_0.h5<br>
output-0001/rho.file_1.h5<br>
output-0001/rho.file_2.h5<br>
<br>
for a 3 process run and similar for runs with more processes. You can probably (I have not tested this) also list both rho and eps:<br>
<br>
!NBLOCKS 3<br>
output-0000/rho.file_0.h5<br>
output-0000/rho.file_1.h5<br>
output-0000/rho.file_2.h5<br>
output-0000/eps.file_0.h5<br>
output-0000/eps.file_1.h5<br>
output-0000/eps.file_2.h5<br>
output-0001/rho.file_0.h5<br>
output-0001/rho.file_1.h5<br>
output-0001/rho.file_2.h5<br>
output-0001/eps.file_0.h5<br>
output-0001/eps.file_1.h5<br>
output-0001/eps.file_2.h5<br>
<br>
possibly you have to set NBLOCKS to 6 (2*3) for this, though I have not tested<br>
this either.<br>
<br>
@Frank:<br>
&gt; Also, in case you have an old version of VisIt, you might need to enable<br>
&gt; the use of soft links in the VisIt Carpet reader. I am currently not<br>
&gt; sure if that change made it back upstream yet. The corresponding patch<br>
&gt; is:<br>
&gt;<br>
&gt; Index: avtCarpetHDF5FileFormat.C<br>
&gt; ===================================================================<br>
&gt; --- avtCarpetHDF5FileFormat.C   (revision 57)<br>
&gt; +++ avtCarpetHDF5FileFormat.C   (working copy)<br>
&gt; @@ -1044,7 +1044,7 @@<br>
&gt;     sprintf(fullname, &quot;%s%s%s&quot;, rootname, rootname[strlen(rootname)-1]==&#39;/&#39; ? &quot;&quot; : &quot;/&quot;,member_name);<br>
&gt;<br>
&gt;     // we are interested in datasets only - skip anything else<br>
&gt; -   H5Gget_objinfo (group_id, member_name, 0, &amp;object_info);<br>
&gt; +   H5Gget_objinfo (group_id, member_name, 1, &amp;object_info);<br>
&gt;     if (object_info.type != H5G_DATASET)<br>
&gt;     {<br>
&gt;        if (object_info.type == H5G_GROUP)<br>
As far as I can tell this code change has not even made it into the ET<br>
maintained copy of the reader:<br>
<a href="https://svn.cactuscode.org/VizTools/CarpetHDF5/" rel="noreferrer" target="_blank">https://svn.cactuscode.org/VizTools/CarpetHDF5/</a><br>
<br>
where if this does not cause any other problems I think it would be<br>
welcome as well.<br>
<br>
Yours,<br>
Roland<br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
<br>
--<br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://keys.gnupg.net" rel="noreferrer" target="_blank">http://keys.gnupg.net</a>.<br>
</font></span><br>_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>Federico Guercilena</div><div>Institut für Theoretische Physik</div><div>Johann Wolfgang Goethe-Universität</div><div>Max-von-Laue-Str. 1</div><div>60438 Frankfurt am Main, Germany</div><div>Telephone: +49 69 798 47887</div><div>Email: guercilena[at]<a href="http://th.physik.uni-frankfurt.de" target="_blank">th.physik.uni-frankfurt.de</a></div></div></div>
</div>