diff options
author | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2011-02-05 17:07:51 +0000 |
---|---|---|
committer | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2011-02-05 17:07:51 +0000 |
commit | 1b6828ab44209fd2084765bfdb1605f6e553547c (patch) | |
tree | 6d37351c13fc924bebe652947e3dabe23e10b9d6 | |
parent | c8f3fb496c9afbd68a17f447daf9f5ab2ac256a5 (diff) |
* update hdf5 table reader routine so that
2006 results from Ott et al. 2007 can be reproduced
by using the Shen table of the Garching group.
Note that the code is currently commented (with #if 0 [...] #endif).
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@36 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
-rw-r--r-- | src/nuc_eos/readtable.F90 | 20 |
1 files changed, 20 insertions, 0 deletions
diff --git a/src/nuc_eos/readtable.F90 b/src/nuc_eos/readtable.F90 index b865b39..0df833f 100644 --- a/src/nuc_eos/readtable.F90 +++ b/src/nuc_eos/readtable.F90 @@ -238,6 +238,26 @@ subroutine nuc_eos_readtable(eos_filename) ! write(6,*) "Done reading eos tables" +#if 0 +! if using the actual garching table; +! keep this around for compatiblity and to +! reproduce 2006 results using the Garching Shen EOS table. + amu_cgs_andi = 1.66d-24 + do i=1,nrho + do j=1,ntemp + do k=1,nye + buffer1 = 10.0d0**alltables(i,j,k,2) + buffer1 = buffer1*amu_cgs_andi/mev_to_erg + buffer2 = (buffer1 + 939.5731d0 - 10.8d0 - amu_mev) + buffer3 = buffer2 - ye(k)*0.511d0 + buffer4 = buffer3/amu_cgs * mev_to_erg + alltables(i,j,k,2) = log10(buffer4) + enddo + enddo + enddo +#endif + + end subroutine nuc_eos_readtable |