use read_value in postprocess

bfps/cpp/full_code/postprocess.cpp

@@ -39,58 +39,36 @@ int postprocess::read_parameters()
     char *string_data;
     sprintf(fname, "%s.h5", this->simname.c_str());
     parameter_file = H5Fopen(fname, H5F_ACC_RDONLY, H5P_DEFAULT);
-    dset = H5Dopen(parameter_file, "/parameters/dealias_type", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->dealias_type);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/dkx", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->dkx);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/dky", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->dky);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/dkz", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->dkz);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/dt", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->dt);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/famplitude", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->famplitude);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/friction_coefficient", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->friction_coefficient);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/fk0", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->fk0);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/fk1", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->fk1);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/fmode", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->fmode);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/forcing_type", H5P_DEFAULT);
-    space = H5Dget_space(dset);
-    memtype = H5Dget_type(dset);
-    string_data = (char*)malloc(256);
-    H5Dread(dset, memtype, H5S_ALL, H5S_ALL, H5P_DEFAULT, &string_data);
-    sprintf(this->forcing_type, "%s", string_data);
-    free(string_data);
-    H5Sclose(space);
-    H5Tclose(memtype);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/nu", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->nu);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/nx", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->nx);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/ny", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->ny);
-    H5Dclose(dset);
-    dset = H5Dopen(parameter_file, "/parameters/nz", H5P_DEFAULT);
-    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &this->nz);
-    H5Dclose(dset);
+    this->dealias_type = hdf5_tools::read_value<int>(
+            parameter_file, "/parameters/dealias_type");
+    this->dkx = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/dkx");
+    this->dky = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/dky");
+    this->dkz = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/dkz");
+    this->dt = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/dt");
+    this->famplitude = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/famplitude");
+    this->friction_coefficient = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/friction_coefficient");
+    this->fk0 = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/fk0");
+    this->fk1 = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/fk1");
+    this->fmode = hdf5_tools::read_value<int>(
+            parameter_file, "/parameters/fmode");
+    sprintf(this->forcing_type, "%s",
+            hdf5_tools::read_string(parameter_file, "/parameters/forcing_type").c_str());
+    this->nu = hdf5_tools::read_value<double>(
+            parameter_file, "/parameters/nu");
+    this->nx = hdf5_tools::read_value<int>(
+            parameter_file, "/parameters/nx");
+    this->ny = hdf5_tools::read_value<int>(
+            parameter_file, "/parameters/ny");
+    this->nz = hdf5_tools::read_value<int>(
+            parameter_file, "/parameters/nz");
     H5Fclose(parameter_file);
     return 0;
 }