diff options
Diffstat (limited to 'src/nuc_eos/readtable.c')
| -rw-r--r-- | src/nuc_eos/readtable.c | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/src/nuc_eos/readtable.c b/src/nuc_eos/readtable.c index 72077e1..786467b 100644 --- a/src/nuc_eos/readtable.c +++ b/src/nuc_eos/readtable.c @@ -43,6 +43,12 @@ void EOS_OMNI_ReadTable(CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS DECLARE_CCTK_ARGUMENTS + CCTK_Info(CCTK_THORNSTRING,"*******************************"); + CCTK_Info(CCTK_THORNSTRING,"Reading nuc_eos table file:"); + CCTK_Info(CCTK_THORNSTRING,nuceos_table_name); + CCTK_Info(CCTK_THORNSTRING,"*******************************"); + + hid_t file; if (!file_is_readable(nuceos_table_name)) CCTK_VWarn(CCTK_WARN_ABORT, __LINE__, __FILE__, CCTK_THORNSTRING, @@ -74,6 +80,7 @@ void EOS_OMNI_ReadTable(CCTK_ARGUMENTS) READ_EOS_HDF5("pointstemp", &ntemp, H5T_NATIVE_INT, H5S_ALL); READ_EOS_HDF5("pointsye", &nye, H5T_NATIVE_INT, H5S_ALL); + // Allocate memory for tables CCTK_REAL *alltables, *logrho, *logtemp, *ye, energy_shift; @@ -116,6 +123,7 @@ void EOS_OMNI_ReadTable(CCTK_ARGUMENTS) // Gamma READ_EOSTABLE_HDF5("gamma", 18); + // Read additional tables and variables READ_EOS_HDF5("logrho", logrho, H5T_NATIVE_DOUBLE, H5S_ALL); READ_EOS_HDF5("logtemp", logtemp, H5T_NATIVE_DOUBLE, H5S_ALL); @@ -129,6 +137,7 @@ void EOS_OMNI_ReadTable(CCTK_ARGUMENTS) CCTK_FNAME(allocate_eosmodule) (&nrho, &ntemp, &nye, alltables, logrho, logtemp, ye, &energy_shift); + // Free the memory again because fortran copied the whole thing now free(ye); free(logtemp); |