[Paraview] problems with Xdmf-file

[Natalie Happenhofer]

Hi!
I´m trying to write my data in the xdmf-format and read it in paraview. I have a 3D - curvilinear grid, i.e. I have to specify every point with its coordinates and then I have scalar attribute data, to be assigned to each node. 
The .xmf meta-data file is recognised by paraview and works fine, but when I click the "Apply"- button, paraview crashes. I believe that I get the .h5 files wrong, when I write them (as I don´t get the Xdmf-API to work, I writing the .xmf-file "manually" with a program and the .h5-files with the HDF5-API. 
Does anyone know how to set the parameters for writing the correct .h5-file (Rank, and the shape of the array which is written to the file)? 
This is as I have it now:

The coordinate file, which has (in this case) the dimensions 350 140 1. 
I write the coordinates in the following array:
float* Coordinates;
Coordinates = new float [3*xdims];
int j=0; //counter

//write xCoord to array
for(i=0;i<xdims;i=i+3) 
    {Coordinates[i] = xPoints[j];
     j++;
    }
//write yCoord to array
j=0;
for(i=1;i<xdims;i=i+3) 
    {Coordinates[i] = yPoints[j];
     j++;
    }
//write zCoord to array
for(i=2;i<xdims;i=i+3) 
    {Coordinates[i] = 0.0;
    }

This array obviously has rank 1, so I write the .h5 file with the following parameters:

file = H5Fcreate(H5FileName.c_str(), H5F_ACC_TRUNC , H5P_DEFAULT, H5P_DEFAULT);
dataspace = H5Screate_simple(1, dims, NULL);
datatype = H5Tcopy(H5T_NATIVE_FLOAT);
status = H5Tset_order(datatype,H5T_ORDER_LE);
dataset = H5Dcreate2(file,DATASETNAME1,datatype,dataspace,H5P_DEFAULT,H5P_DEFAULT,H5P_DEFAULT);
status = H5Dwrite(dataset, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, Coordinates);


The attribute data files I write analogue. 

My .xmf file looks like this:
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []>

<Xdmf>
<Domain>
<Grid Name="Curvilinear" GridType="Uniform">
<Topology TopologyType="3DSMesh" Dimensions ="350 140 1">
</Topology><Geometry GeometryType="XYZ">
<DataItem DataType="Float" Precision="4" Dimensions="350 140 1"  Format="HDF">
C:\Dokumente und Einstellungen\User1\Eigene Dateien\xdmf_files\c1d-T132-30SEPT-deep_coord_00000.h5/Curvilinear/XYZ
</DataItem>
</Geometry>
<Attribute Name="Density" Active="1" Type="Scalar" Center="Node">
<DataItem DataType="Float" Precision="4" Dimensions="350 140 1" Format="HDF">
C:\Dokumente und Einstellungen\User1\Eigene Dateien\xdmf_files\c1d-T132-30SEPT-deep.p_00000.h5:/Curvilinear/Density
</DataItem>
</Attribute>
</Grid>
</Domain>
</Xdmf>

Any help would be appreciated!
Thx,
NH




_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20081013/290c88e8/attachment-0001.htm>