[Users] Merging carpet hdf5 data without creating additional files

[Roland Haas]

Hello Frank, Michael,

@Michael: both of what Frank or David describe should work. If you know
how to compile the VisIt reader yourselv, you can also try the copy of
CarpetHDF5 at https://bitbucket.org/rhaas80/carpethdf5 either the
"multifile" branch (which is the default one that has all features) or
the "for_VisIt" branch which has fewer features and is what I wanted at
one point propose for inclusion in VisIt. Both will let you
create .visit files listing your hdf5 files and VisIt will treat them a
a single database. You will have to create .visit files that look like
this:

!NBLOCKS 3
output-0000/rho.file_0.h5
output-0000/rho.file_1.h5
output-0000/rho.file_2.h5
output-0001/rho.file_0.h5
output-0001/rho.file_1.h5
output-0001/rho.file_2.h5

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:

!NBLOCKS 3
output-0000/rho.file_0.h5
output-0000/rho.file_1.h5
output-0000/rho.file_2.h5
output-0000/eps.file_0.h5
output-0000/eps.file_1.h5
output-0000/eps.file_2.h5
output-0001/rho.file_0.h5
output-0001/rho.file_1.h5
output-0001/rho.file_2.h5
output-0001/eps.file_0.h5
output-0001/eps.file_1.h5
output-0001/eps.file_2.h5

possibly you have to set NBLOCKS to 6 (2*3) for this, though I have not tested
this either.

@Frank:
> Also, in case you have an old version of VisIt, you might need to enable
> the use of soft links in the VisIt Carpet reader. I am currently not
> sure if that change made it back upstream yet. The corresponding patch
> is:
> 
> Index: avtCarpetHDF5FileFormat.C
> ===================================================================
> --- avtCarpetHDF5FileFormat.C   (revision 57)
> +++ avtCarpetHDF5FileFormat.C   (working copy)
> @@ -1044,7 +1044,7 @@
>     sprintf(fullname, "%s%s%s", rootname, rootname[strlen(rootname)-1]=='/' ? "" : "/",member_name);
>     
>     // we are interested in datasets only - skip anything else
> -   H5Gget_objinfo (group_id, member_name, 0, &object_info);
> +   H5Gget_objinfo (group_id, member_name, 1, &object_info);
>     if (object_info.type != H5G_DATASET)
>     {
>        if (object_info.type == H5G_GROUP)
As far as I can tell this code change has not even made it into the ET
maintained copy of the reader:
https://svn.cactuscode.org/VizTools/CarpetHDF5/

where if this does not cause any other problems I think it would be
welcome as well.

Yours,
Roland



-- 
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://keys.gnupg.net.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 5517 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20160517/c8d3df1f/attachment.bin