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