diff --git a/bfps/cpp/full_code/NSVE.cpp b/bfps/cpp/full_code/NSVE.cpp
index 292d05298a236fa17ff8f30e1fbfb26eb96f3dc2..2a6ae9379020192895ccf8f92921f3a08bb8b333 100644
--- a/bfps/cpp/full_code/NSVE.cpp
+++ b/bfps/cpp/full_code/NSVE.cpp
@@ -246,97 +246,6 @@ int NSVE<rnumber>::finalize(void)
     return EXIT_SUCCESS;
 }
 
-//template <typename rnumber>
-//int NSVE<rnumber>::read_parameters(void)
-//{
-//    hid_t parameter_file;
-//    hid_t dset, memtype, space;
-//    char *string_data;
-//    parameter_file = H5Fopen(
-//            (this->simname + std::string(".h5")).c_str(),
-//            H5F_ACC_RDONLY,
-//            H5P_DEFAULT);
-//    dset = H5Dopen(parameter_file, "/parameters/checkpoints_per_file", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &checkpoints_per_file);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/dealias_type", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &dealias_type);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/dkx", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &dkx);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/dky", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &dky);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/dkz", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &dkz);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/dt", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &dt);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/famplitude", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &famplitude);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/fk0", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &fk0);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/fk1", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &fk1);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/fmode", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &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(forcing_type, "%s", string_data);
-//    free(string_data);
-//    H5Sclose(space);
-//    H5Tclose(memtype);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/histogram_bins", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &histogram_bins);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/max_velocity_estimate", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &max_velocity_estimate);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/max_vorticity_estimate", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &max_vorticity_estimate);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/niter_out", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &niter_out);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/niter_part", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &niter_part);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/niter_stat", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &niter_stat);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/niter_todo", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &niter_todo);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/nparticles", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &nparticles);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/nu", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, &nu);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/nx", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &nx);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/ny", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &ny);
-//    H5Dclose(dset);
-//    dset = H5Dopen(parameter_file, "/parameters/nz", H5P_DEFAULT);
-//    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &nz);
-//    H5Dclose(dset);
-//    H5Fclose(parameter_file);
-//    return EXIT_SUCCESS;
-//}
-
-
 template <typename rnumber>
 int NSVE<rnumber>::grow_file_datasets()
 {