GitLab

Summary

The asciih5block script returns error messages and an empty (header only) file

Steps to reproduce

module use unstable; module load gcc/7.3.0 mpich/3.3 hdf5/1.10.4 H5hut/2.0.0rc5

I believe this will happen with any input file, an example is given here.

Example input file: NewChimney_R000_D000E.txt

Modify ascii2h5block_asgic.cpp for the right dimensions: ascii2h5block_asgic.cpp

$ export LDLIBS=-lH5hut
$ make ascii2h5block_asgic
$ ./ascii2h5block_asgic NewChimney_R000_D000E.txt out
--------------------------------------------------------
Using NewChimney_R000_D000E.txt to create test_CYC.h5part
Lines in finE: 67370
Hardcoded limits: x(-90/-10) cm
Hardcoded limits: y(-30/30) cm
Hardcoded limits: z(-1/1) cm
Hardcoded spacing: 0.5 cm
Hardcoded nlines: 67367
File nlines: 67365
Grid dimensions: Px = 23 , Py = 29 , Pz = 101
E_max = (29.0302, -5.6267, -23.2274) V/cm at index 20641.
Converting from V/cm and cm to kV/mm and mm before saving h5part
[proc 0] E: H5Block3dSetView: Iteration is invalid! Have you set the time step?
[proc 0] E: H5Block3dWriteVector3dFieldFloat64: No view has been defined!
[proc 0] E: H5Block3dWriteVector3dFieldFloat64: No view has been defined!
Done, bye ...

What is the current bug behavior?

Header only H5-file

What is the expected correct behavior?

H5-file with correct field

Possible fixes

The H5SetStep should be set before the View:

        H5SetStep(file, 0);
        h5_err_t h5err = H5Block3dSetView(file,
                                          0, gridPx - 1,
                                          0, gridPy - 1,
                                          0, gridPz - 1);

This modified ascii2h5block_asgic.cpp runs successfully.