Initial commit from Polytechnique Montreal

1 files changed, 877 insertions, 0 deletions

diff --git a/Ganlib/src/hdf5_aux.c b/Ganlib/src/hdf5_aux.c
new file mode 100644
index 0000000..74e7565
--- /dev/null
+++ b/Ganlib/src/hdf5_aux.c
@@ -0,0 +1,877 @@
+
+/**********************************/
+/* C API for hdf5 file support    */
+/* (auxiliary functions)          */
+/* author: A. Hebert (30/11/2021) */
+/**********************************/
+
+/*
+   Copyright (C) 2021 Ecole Polytechnique de Montreal
+
+   This library is free software; you can redistribute it and/or
+   modify it under the terms of the GNU Lesser General Public
+   License as published by the Free Software Foundation; either
+   version 2.1 of the License, or (at your option) any later version.
+ */
+
+#if defined(HDF5_LIB)
+#include <stdlib.h>
+#include <string.h>
+#include "hdf5_aux.h"
+#define MAX_NAME 1024
+
+hid_t h5tools_get_native_type(hid_t type);
+herr_t print_info(hid_t ifile, const char *name, void *opdata);  /* Operator function */
+
+herr_t iretcd;
+
+static char AbortString[164];
+static char name1024[1024];
+
+void hdf5_open_file_c(const char *fname, hid_t *ifile, int_32 irdonly) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Open an HDF5 file
+ *
+ * input parameters:
+ *   fname : hdf5 file namr
+ *   ifile : hdf5 file identificator.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_open_file_c";
+  switch (irdonly) {
+  case 0:
+    *ifile = H5Fopen(fname, H5F_ACC_RDWR, H5P_DEFAULT);
+    if (*ifile < 0) {
+      sprintf(AbortString,"%s: H5Fopen failure on read-write HDF5 file '%s'.",nomsub,fname);
+      xabort_c(AbortString);
+     }
+    break;
+  case 1:
+    *ifile = H5Fopen(fname, H5F_ACC_RDONLY, H5P_DEFAULT);
+    if (*ifile < 0) {
+      sprintf(AbortString,"%s: H5Fopen failure on read-only HDF5 file '%s'.",nomsub,fname);
+      xabort_c(AbortString);
+     }
+    break;
+  default:
+    sprintf(AbortString,"%s: Invalid action on HDF5 file '%.72s'.",nomsub,fname);
+    xabort_c(AbortString);
+    break;
+  }
+}
+
+void hdf5_close_file_c(hid_t *ifile) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Close an HDF5 file
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_close_file_c";
+  iretcd = H5Fclose(*ifile);
+  if (iretcd != 0) {
+    sprintf(AbortString,"%s: HDF5 close failure. iretcd=%d\n", nomsub, iretcd);
+    xabort_c(AbortString);
+  }
+}
+
+void hdf5_list_c(hid_t *ifile, const char *namp) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * List the root table of contents
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_list_c";
+  H5G_stat_t statbuf;
+  hid_t loc_id=*ifile;
+  sprintf(name1024,"/%s",namp);
+  printf("\nTable of contents --'%s'--\n", name1024);
+  if(!hdf5_group_exists_c(ifile, namp)) {
+    printf("%s: HDF5 missing group=%s\n", nomsub, namp);
+    return;
+  }
+  H5Gget_objinfo(loc_id, name1024, 0, &statbuf);
+  switch (statbuf.type) {
+  case H5G_GROUP:
+    printf("dataset/group name........................                                  ");
+    printf("type............   bytes..    shape.................\n");
+    H5Giterate(loc_id, name1024, NULL, print_info, NULL);
+    break;
+  default:
+    printf("%s is not a H5G_GROUP.\n", namp);
+    break;
+  }
+}
+
+void hdf5_get_dimensions_c(hid_t *ifile, const char *namp, int_32 *rank) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Find dataset rank (number of dimensions)
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *
+ * output parameter:
+ *   rank  : number of dimensions.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_get_dimensions_c";
+  hid_t loc_id=*ifile;
+  
+  hid_t dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  hid_t filespace = H5Dget_space (dataset);
+  *rank = (int_32)H5Sget_simple_extent_ndims(filespace);
+  if (*rank < 0) {
+      sprintf(AbortString,"%s: H5Sget_simple_extent_ndims failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Dclose(dataset);
+}
+
+void hdf5_get_num_group_c(hid_t *ifile, const char *namp, int_32 *nbobj) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Find the number of objects (daughter datasets and daughter groups) in a group
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the group.
+ *
+ * output parameter:
+ *   nbobj : number of objects
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_get_num_group_c";
+  hid_t loc_id=*ifile;
+
+  hid_t group = H5Oopen(loc_id,namp,H5P_DEFAULT);
+  if (group < 0) {
+      sprintf(AbortString,"%s: H5Oopen failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  hsize_t nobj;
+  iretcd = H5Gget_num_objs(group, &nobj);
+  if (iretcd < 0) {
+      sprintf(AbortString,"%s: H5Gget_num_objs failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  *nbobj = (int_32)nobj;
+  H5Oclose(group);
+}
+
+void hdf5_list_datasets_c(hid_t *ifile, const char *namp, int_32 *ndsets, char *idata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Recover character daughter dataset names in a group.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the group.
+ *
+ * output parameter:
+ *  ndsets : number of daughter datasets.
+ *   idata : dataset names.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_list_datasets_c";
+  hid_t loc_id=*ifile;
+  int idx;
+
+  hid_t group = H5Oopen(loc_id,namp,H5P_DEFAULT);
+  if (group < 0) {
+      sprintf(AbortString,"%s: H5Oopen failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  hsize_t nobj;
+  int otype;
+  char memb_name[MAX_NAME];
+  iretcd = H5Gget_num_objs(group, &nobj);
+  if (iretcd < 0) {
+      sprintf(AbortString,"%s: H5Gget_num_objs failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  *ndsets = 0;
+  for (idx = 0; idx < nobj; idx++) {
+    H5Gget_objname_by_idx(group, (hsize_t)idx, memb_name, (size_t)MAX_NAME);
+    otype =  H5Gget_objtype_by_idx(group, (size_t)idx);
+    switch(otype) {
+    case H5G_DATASET:
+      strncpy (idata+MAX_NAME*(*ndsets), memb_name, MAX_NAME);
+      (*ndsets)++;
+      break;
+    }
+  }
+  H5Oclose(group);
+}
+
+void hdf5_list_groups_c(hid_t *ifile, const char *namp, int_32 *ndsets, char *idata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Recover character daughter group names in a group.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the group.
+ *
+ * output parameter:
+ *  ndsets : number of daughter groups.
+ *   idata : group names.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_list_groups_c";
+  hid_t loc_id=*ifile;
+  int idx;
+
+  hid_t group = H5Oopen(loc_id,namp,H5P_DEFAULT);
+  if (group < 0) {
+      sprintf(AbortString,"%s: H5Oopen failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  hsize_t nobj;
+  int otype;
+  char memb_name[MAX_NAME];
+  iretcd = H5Gget_num_objs(group, &nobj);
+  if (iretcd < 0) {
+      sprintf(AbortString,"%s: H5Gget_num_objs failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  *ndsets = 0;
+  for (idx = 0; idx < nobj; idx++) {
+    H5Gget_objname_by_idx(group, (hsize_t)idx, memb_name, (size_t)MAX_NAME);
+    otype =  H5Gget_objtype_by_idx(group, (size_t)idx);
+    switch(otype) {
+    case H5G_GROUP:
+      strncpy (idata+MAX_NAME*(*ndsets), memb_name, MAX_NAME);
+      (*ndsets)++;
+      break;
+    }
+  }
+  H5Oclose(group);
+}
+
+void hdf5_info_c(hid_t *ifile, const char *namp, int_32 *rank, int_32 *type, int_32 *nbyte, int_32 *dimsr) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Find dataset information.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *
+ * output parameter:
+ *   rank  : number of dimensions.
+ *   type  : type of dataset.
+ *   nbyte : number of bytes per value.
+ *   dimsr : shape.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_info_c";
+  hid_t loc_id=*ifile;
+  hsize_t dimsr_t[10];
+  int i;
+  
+  hid_t dataset = -1;
+  H5E_BEGIN_TRY {
+    dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  } H5E_END_TRY;
+  if (dataset < 0) {
+      *rank=0;
+      *type=99;
+      *nbyte=0;
+      return;
+  }
+  hid_t filespace = H5Dget_space (dataset);
+  *rank = H5Sget_simple_extent_ndims (filespace);
+  if (*rank < 0) {
+      sprintf(AbortString,"%s: H5Sget_simple_extent_ndims failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  } else if (*rank > 10) {
+      sprintf(AbortString,"%s: the object '%.72s' has rank= %d > 10.",nomsub,namp,*rank);
+      xabort_c(AbortString);
+  }
+  hid_t htype = H5Dget_type(dataset);
+  if (htype < 0) {
+      sprintf(AbortString,"%s: H5Dget_type failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  size_t precision = 0;
+  switch (H5Tget_class(htype)) {
+  case H5T_INTEGER: 
+       *type=1;
+       precision = H5Tget_size(htype);
+       break;
+  case H5T_FLOAT: 
+       /*precision = H5Tget_precision(htype);*/
+       precision = H5Tget_size(htype);
+       if (precision == 4) {
+         *type=2;
+       } else if (precision == 8) {
+         *type=4;
+       } else {
+         *type=7;
+       }
+       break;
+  case H5T_STRING: 
+       *type=3;
+       precision = H5Tget_size(htype);
+       break;
+  default:
+       *type=7;
+  }
+  *nbyte = precision;
+  iretcd = H5Sget_simple_extent_dims(filespace, dimsr_t, NULL);
+  if (iretcd == -1) {
+      sprintf(AbortString,"%s: H5Sget_simple_extent_dims failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  for(i=0; i<*rank; i++) dimsr[i] = (int_32)dimsr_t[(*rank)-i-1];
+  
+  H5Tclose(htype);
+  H5Dclose(dataset);
+}
+
+int_32 hdf5_group_exists_c(hid_t *ifile, const char *namp)
+/*
+ *-----------------------------------------------------------------------
+ *
+ * Test for existence of a group.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the group.
+ *
+ * output parameter:
+ *  hdf5_group_exists_c : =0: the group does exists; =1: does not exists.
+ *
+ *-----------------------------------------------------------------------
+ */
+{
+  hid_t loc_id=*ifile;
+  H5G_stat_t statbuf;
+H5E_BEGIN_TRY
+  if(strlen(namp)==0) {
+    iretcd = H5Gget_objinfo(loc_id, "/", 0, &statbuf);
+  } else {
+    iretcd = H5Gget_objinfo(loc_id, namp, 0, &statbuf);
+  }
+H5E_END_TRY
+  if (iretcd >= 0) {
+    if (statbuf.type == H5G_GROUP) return 1;
+  }
+  return 0;
+}
+
+void hdf5_create_group_c(hid_t *ifile, const char *namp)
+/*
+ *-----------------------------------------------------------------------
+ *
+ * Create a group.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the group to create.
+ *
+ *-----------------------------------------------------------------------
+ */
+{
+  char *nomsub="hdf5_create_group_c";
+  hid_t loc_id=*ifile;
+  H5G_stat_t statbuf;
+H5E_BEGIN_TRY
+  iretcd = H5Gget_objinfo(loc_id, namp, 0, &statbuf);
+H5E_END_TRY
+  if (iretcd >= 0) {
+    if (statbuf.type == H5G_GROUP) return;
+  }
+  hid_t groupid = H5Gcreate(loc_id, namp, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
+  if (groupid < 0)  {
+    sprintf(AbortString,"%s: H5Gcreate failure on object '%.72s'.",nomsub,namp);
+    xabort_c(AbortString);
+  }
+  H5Gclose(groupid);
+}
+
+void hdf5_delete_c(hid_t *ifile, const char *namp)
+/*
+ *-----------------------------------------------------------------------
+ *
+ * Delete a group or a dataset.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the group or dataset to delete.
+ *
+ *-----------------------------------------------------------------------
+ */
+{
+  char *nomsub="hdf5_delete_c";
+  hid_t loc_id=*ifile;
+  iretcd = H5Ldelete(loc_id,namp,H5P_DEFAULT);
+  if (iretcd < 0) {
+    sprintf(AbortString,"%s: HDF5 delete failure. iretcd=%d\n", nomsub, iretcd);
+    xabort_c(AbortString);
+  }
+}
+
+void hdf5_copy_c(hid_t *ifile_s, const char *namp_s, hid_t *ifile_d, const char *namp_d)
+/*
+ *-----------------------------------------------------------------------
+ *
+ * Copy a group or a dataset from one location to another. The source and
+ * destination need not be in the same file.
+ *
+ * input parameters:
+ *   ifile_s : hdf5 source file identificator.
+ *    namp_s : character name of the source group or dataset to copy.
+ *
+ * output parameters:
+ *   ifile_d : hdf5 destination file identificator.
+ *    namp_d : character name of the destination group or dataset.
+ *
+ *-----------------------------------------------------------------------
+ */
+{
+  char *nomsub="hdf5_copy_c";
+  if(hdf5_group_exists_c(ifile_s, namp_s)) {
+    /* namp_s is a group */
+    hid_t loc_id_s=*ifile_s;
+    hid_t loc_id_d=*ifile_d;
+    iretcd = H5Ocopy(loc_id_s,namp_s,loc_id_d,namp_d,H5P_DEFAULT,H5P_DEFAULT);
+    if (iretcd < 0) {
+      sprintf(AbortString,"%s: HDF5 copy failure. iretcd=%d\n", nomsub, iretcd);
+      xabort_c(AbortString);
+    }
+  } else {
+    /* namp_s is a dataset */
+    int i;
+    int_32 rank, type, nbyte, dimsr[10], dimsr2[10], length;
+    hdf5_info_c(ifile_s, namp_s, &rank, &type, &nbyte, dimsr);
+    length = 1;
+    for(i=0; i<rank; i++) {
+      dimsr2[i] = (int_32)dimsr[(rank)-i-1];
+      length = length*dimsr2[i];
+    }
+    if (type == 1) {
+      int_32 *idata = (int_32 *)malloc(sizeof(int_32)*length);
+      hdf5_read_data_int_c(ifile_s, namp_s, idata);
+      hdf5_write_data_int_c(ifile_d, namp_s, rank, dimsr2, idata);
+      free(idata);
+    } else if (type == 2) {
+      float *rdata = (float *)malloc(sizeof(float)*length);
+      hdf5_read_data_real4_c(ifile_s, namp_s, rdata);
+      hdf5_write_data_real4_c(ifile_d, namp_s, rank, dimsr2, rdata);
+      free(rdata);
+    } else if (type == 3) {
+      char *ihdata = (char *)malloc(nbyte*length);
+      hdf5_read_data_string_c(ifile_s, namp_s, ihdata);
+      hdf5_write_data_string_c(ifile_d, namp_s, rank, nbyte, dimsr2, ihdata);
+      free(ihdata);
+    } else if (type == 4) {
+      double *rddata = (double *)malloc(sizeof(double)*length);
+      hdf5_read_data_real8_c(ifile_s, namp_s, rddata);
+      hdf5_write_data_real8_c(ifile_d, namp_s, rank, dimsr2, rddata);
+      free(rddata);
+    } else {
+      sprintf(AbortString,"%s: dataset '%.60s' has the wrong type(2).", nomsub, namp_s);
+      xabort_c(AbortString);
+    }
+  }
+}
+
+void hdf5_read_data_int_c(hid_t *ifile, const char *namp, int_32 *idata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy an integer dataset from hdf5 into memory.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *
+ * output parameter:
+ *   idata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_read_data_int_c";
+  hid_t loc_id=*ifile;
+  hid_t dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  iretcd = H5Dread(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, idata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be readed.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Dclose(dataset);
+}
+
+void hdf5_read_data_real4_c(hid_t *ifile, const char *namp, float *rdata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy a real(4) dataset from hdf5 into memory.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *
+ * output parameter:
+ *   rdata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_read_data_real4_c";
+  hid_t loc_id=*ifile;
+  hid_t dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  iretcd = H5Dread(dataset, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be readed.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Dclose(dataset);
+}
+
+void hdf5_read_data_real8_c(hid_t *ifile, const char *namp, double *rdata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy a real(8) dataset from hdf5 into memory.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *
+ * output parameter:
+ *   rdata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_read_data_real8_c";
+  hid_t loc_id=*ifile;
+  hid_t dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  iretcd = H5Dread(dataset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be readed.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Dclose(dataset);
+}
+
+void hdf5_read_data_string_c(hid_t *ifile, const char *namp, char *idata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy a string dataset from hdf5 into memory.
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *
+ * output parameter:
+ *   idata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_read_data_string_c";
+  hid_t loc_id=*ifile;
+  
+  hid_t dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  hid_t datatype = H5Dget_type(dataset);
+  if (datatype < 0) {
+      sprintf(AbortString,"%s: H5Dget_type failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  iretcd = H5Dread(dataset, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, idata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be readed.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Tclose(datatype);
+  H5Dclose(dataset);
+}
+
+void hdf5_write_data_int_c(hid_t *ifile, const char *namp, int_32 rank, int_32 *dimsr, int_32 *idata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy an integer array from memory into a hdf5 dataset
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *    rank : number of dimensions.
+ *   dimsr : number of information along each dimension.
+ *   idata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_write_data_int_c";
+  hid_t loc_id=*ifile;
+  int i;
+  hsize_t dimsr_t[10];
+  if (rank > 10) {
+      sprintf(AbortString,"%s: rank > 10 on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  hid_t dataset = -1;
+  H5E_BEGIN_TRY {
+  dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  } H5E_END_TRY;
+  if (dataset  >= 0) {
+    iretcd = H5Ldelete(loc_id,namp,H5P_DEFAULT);
+    if (iretcd < 0) {
+      sprintf(AbortString,"%s: HDF5 delete failure. iretcd=%d\n", nomsub, iretcd);
+      xabort_c(AbortString);
+    }
+    H5Dclose(dataset);
+  }
+  for (i = 0; i < rank; ++i) dimsr_t[i] = dimsr[i];
+  hid_t dataspace = H5Screate_simple(rank, dimsr_t, NULL); 
+  hid_t datatype = H5Tcopy(H5T_NATIVE_INT);
+  dataset = H5Dcreate(loc_id, namp, datatype, dataspace,
+			H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
+  if (dataset < 0) {
+      sprintf(AbortString,"%s: H5Dcreate failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  iretcd = H5Dwrite(dataset, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, idata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be written.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Sclose(dataspace);
+  H5Tclose(datatype);
+  H5Dclose(dataset);
+}
+
+void hdf5_write_data_real4_c(hid_t *ifile, const char *namp, int_32 rank, int_32 *dimsr, float *rdata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy a real(4) array from memory into a hdf5 dataset
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *    rank : number of dimensions.
+ *   dimsr : number of information along each dimension.
+ *   rdata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_write_data_real4_c";
+  hid_t loc_id=*ifile;
+  int i;
+  hsize_t dimsr_t[10];
+  if (rank > 10) {
+      sprintf(AbortString,"%s: rank > 10 on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  hid_t dataset = -1;
+  H5E_BEGIN_TRY {
+  dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  } H5E_END_TRY;
+  if (dataset  >= 0) {
+    iretcd = H5Ldelete(loc_id,namp,H5P_DEFAULT);
+    if (iretcd < 0) {
+      sprintf(AbortString,"%s: HDF5 delete failure. iretcd=%d\n", nomsub, iretcd);
+      xabort_c(AbortString);
+    }
+    H5Dclose(dataset);
+  }
+  for (i = 0; i < rank; ++i) dimsr_t[i] = dimsr[i];
+  hid_t dataspace = H5Screate_simple(rank, dimsr_t, NULL); 
+  hid_t datatype = H5Tcopy(H5T_NATIVE_FLOAT);
+  dataset = H5Dcreate(loc_id, namp, datatype, dataspace,
+			H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
+  if (dataset < 0) {
+      sprintf(AbortString,"%s: H5Dcreate failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  iretcd = H5Dwrite(dataset, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be written.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Sclose(dataspace);
+  H5Tclose(datatype);
+  H5Dclose(dataset);
+}
+
+void hdf5_write_data_real8_c(hid_t *ifile, const char *namp, int_32 rank, int_32 *dimsr, double *rdata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy a real(8) array from memory into a hdf5 dataset
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *    rank : number of dimensions.
+ *   dimsr : number of information along each dimension.
+ *   rdata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_write_data_real8_c";
+  hid_t loc_id=*ifile;
+  int i;
+  hsize_t dimsr_t[10];
+  if (rank > 10) {
+      sprintf(AbortString,"%s: rank > 10 on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  hid_t dataset = -1;
+  H5E_BEGIN_TRY {
+  dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  } H5E_END_TRY;
+  if (dataset  >= 0) {
+    iretcd = H5Ldelete(loc_id,namp,H5P_DEFAULT);
+    if (iretcd < 0) {
+      sprintf(AbortString,"%s: HDF5 delete failure. iretcd=%d\n", nomsub, iretcd);
+      xabort_c(AbortString);
+    }
+    H5Dclose(dataset);
+  }
+  for (i = 0; i < rank; ++i) dimsr_t[i] = dimsr[i];
+  hid_t dataspace = H5Screate_simple(rank, dimsr_t, NULL); 
+  hid_t datatype = H5Tcopy(H5T_NATIVE_DOUBLE);
+  dataset = H5Dcreate(loc_id, namp, datatype, dataspace,
+			H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
+  if (dataset < 0) {
+      sprintf(AbortString,"%s: H5Dcreate failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  iretcd = H5Dwrite(dataset, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be written.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Sclose(dataspace);
+  H5Tclose(datatype);
+  H5Dclose(dataset);
+}
+
+void hdf5_write_data_string_c(hid_t *ifile, const char *namp, int_32 rank, int_32 len, int_32 *dimsr, char *idata) {
+/*
+ *----------------------------------------------------------------------
+ *
+ * Copy a character array from memory into a hdf5 dataset
+ *
+ * input parameters:
+ *   ifile : hdf5 file identificator.
+ *    namp : character name of the dataset.
+ *    rank : number of dimensions.
+ *     len : length of a string element in the array (in bytes).
+ *   dimsr : number of information along each dimension.
+ *   idata : information elements.
+ *
+ *----------------------------------------------------------------------
+ */
+  char *nomsub="hdf5_write_data_string_c";
+  hid_t loc_id=*ifile;
+  int i;
+  hsize_t dimsr_t[10];
+  if (rank > 10) {
+      sprintf(AbortString,"%s: rank > 10 on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  hid_t dataset = -1;
+  H5E_BEGIN_TRY {
+  dataset = H5Dopen(loc_id,namp,H5P_DEFAULT);
+  } H5E_END_TRY;
+  if (dataset  >= 0) {
+    iretcd = H5Ldelete(loc_id,namp,H5P_DEFAULT);
+    if (iretcd < 0) {
+      sprintf(AbortString,"%s: HDF5 delete failure. iretcd=%d\n", nomsub, iretcd);
+      xabort_c(AbortString);
+    }
+    H5Dclose(dataset);
+  }
+  for (i = 0; i < rank; ++i) dimsr_t[i] = dimsr[i];
+  hid_t dataspace = H5Screate_simple(rank, dimsr_t, NULL);
+  hid_t datatype = H5Tcopy(H5T_C_S1);
+  H5Tset_size(datatype, len * sizeof(char));
+  dataset = H5Dcreate(loc_id, namp, datatype , dataspace, H5P_DEFAULT,
+      H5P_DEFAULT, H5P_DEFAULT);
+  if (dataset < 0) {
+      sprintf(AbortString,"%s: H5Dcreate failure on object '%.72s'.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  iretcd = H5Dwrite(dataset, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, idata);
+  if (iretcd != 0) {
+      sprintf(AbortString,"%s: the object '%.72s' cannot be written.",nomsub,namp);
+      xabort_c(AbortString);
+  }
+  H5Sclose(dataspace);
+  H5Tclose(datatype);
+  H5Dclose(dataset);
+}
+
+herr_t print_info(hid_t loc_id, const char *name, void *opdata) {
+    H5G_stat_t statbuf;
+    char* ctype[]={"","INTEGER","REAL","CHARACTER","DOUBLE PRECISION",
+                   "LOGICAL","COMPLEX","UNDEFINED"};
+    /*
+     * Get type of the object and display its name and type.
+     * The name of the object is passed to this function by 
+     * the Library.
+     */
+    H5Gget_objinfo(loc_id, name, 0, &statbuf);
+    int rank, type, nbyte;
+    int_32 dimsr[6];
+    switch (statbuf.type) {
+    case H5G_GROUP: 
+         printf(" '%-72s' GROUP \n", name);
+         break;
+    case H5G_DATASET: 
+         hdf5_info_c(&loc_id, name, &rank, &type, &nbyte, dimsr);
+         if (rank == 1) {
+           printf(" '%-72s' %-16s   %-10d %d \n", name,ctype[type],nbyte,dimsr[0]);
+         } else if (rank == 2) {
+           printf(" '%-72s' %-16s   %-10d %d  %d\n", name,ctype[type],nbyte,dimsr[0],dimsr[1]);
+         } else if (rank == 3) {
+           printf(" '%-72s' %-16s   %-10d %d  %d  %d\n", name,ctype[type],nbyte,dimsr[0],dimsr[1],dimsr[2]);
+         } else if (rank == 4) {
+           printf(" '%-72s' %-16s   %-10d %d  %d  %d  %d\n", name,ctype[type],nbyte,dimsr[0],dimsr[1],
+           dimsr[2],dimsr[3]);
+         } else if (rank == 5) {
+           printf(" '%-72s' %-16s   %-10d %d  %d  %d  %d  %d\n", name,ctype[type],nbyte,dimsr[0],
+           dimsr[1],dimsr[2],dimsr[3],dimsr[4]);
+         } else if (rank == 6) {
+           printf(" '%-72s' %-16s   %-10d %d  %d  %d  %d  %d  %d\n", name,ctype[type],nbyte,dimsr[0],
+           dimsr[1],dimsr[2],dimsr[3],dimsr[4],dimsr[5]);
+         }
+         break;
+    case H5G_TYPE: 
+         printf(" '%-72s' NAMED DATATYPE \n", name);
+         break;
+    default:
+         printf(" '%-72s' UNKNOWN \n", name);
+    }
+    fflush(stdout);
+    return 0;
+ }
+#endif /* defined(HDF5_LIB) */