c: example read for a Meraxes grid

example_grid_read.c

#include <stdlib.h>
#include <hdf5.h>
#include <hdf5_hl.h>
#include <math.h>

int main(int argc, char *argv[])
{
  hsize_t dims[1];
  size_t type_size;
  H5T_class_t class_id;
  hid_t fd;
  float *grid;
  int dim;
  int snapshot = 51;
  char ds_name[] = {"delta_T"};
  char ds_path[128];

  printf("Reading grid...\n");
  fd = H5Fopen("/Users/smutch/data_temp/meraxes.hdf5", H5F_ACC_RDONLY, H5P_DEFAULT);
  sprintf(ds_path, "Snap%03d/Grids/%s", snapshot, ds_name);
  H5LTget_dataset_info(fd, ds_path, dims, &class_id, &type_size);
  grid = malloc(sizeof(float) * dims[0]);
  dim = (int)(cbrt((double)(dims[0])));
  H5LTread_dataset_float(fd, ds_path, grid);
  H5Fclose(fd);
  printf("...done\n");
  
  return 0;
}