From 1b6828ab44209fd2084765bfdb1605f6e553547c Mon Sep 17 00:00:00 2001 From: cott Date: Sat, 5 Feb 2011 17:07:51 +0000 Subject: * 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 --- src/nuc_eos/readtable.F90 | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) 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 -- cgit v1.2.3