diff options
Diffstat (limited to 'src/nuc_eos/readtable.F90')
-rw-r--r-- | src/nuc_eos/readtable.F90 | 12 |
1 files changed, 11 insertions, 1 deletions
diff --git a/src/nuc_eos/readtable.F90 b/src/nuc_eos/readtable.F90 index 1eee600..a06cff8 100644 --- a/src/nuc_eos/readtable.F90 +++ b/src/nuc_eos/readtable.F90 @@ -1,3 +1,7 @@ +#include "cctk.h" +#include "cctk_Parameters.h" + + subroutine nuc_eos_readtable(eos_filename) ! This routine reads the table and initializes ! all variables in the module. @@ -7,6 +11,8 @@ subroutine nuc_eos_readtable(eos_filename) implicit none + DECLARE_CCTK_PARAMETERS; + character(*) eos_filename character(len=100) message @@ -205,13 +211,17 @@ subroutine nuc_eos_readtable(eos_filename) call h5dread_f(dset_id, H5T_NATIVE_DOUBLE, energy_shift, dims1, error) call h5dclose_f(dset_id,error) accerr=accerr+error - + if(accerr.ne.0) then stop "Problem reading EOS table file" endif call h5fclose_f (file_id,error) + if(do_energy_shift.ne.1) then + energy_shift = 0.0d0 + endif + ! not a good idea to call this function as it will screw recovery ! call h5close_f (error) |