asciih5block not working with recent H5Hut library?
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.