diff options
Diffstat (limited to 'src')
| -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 |