[Commits] [svn:einsteintoolkit] EOS_Omni/trunk/src/nuc_eos/ (Rev. 36)

[cott at tapir.caltech.edu]

User: cott
Date: 2011/02/05 11:07 AM

Modified:
 /trunk/src/nuc_eos/
  readtable.F90

Log:
 * 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).

File Changes:

Directory: /trunk/src/nuc_eos/
==============================

File [modified]: readtable.F90
Delta lines: +20 -0
===================================================================
--- trunk/src/nuc_eos/readtable.F90	2011-01-26 16:15:44 UTC (rev 35)
+++ trunk/src/nuc_eos/readtable.F90	2011-02-05 17:07:51 UTC (rev 36)
@@ -238,7 +238,27 @@
 
 !  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