diff options
author | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2013-03-23 20:08:12 +0000 |
---|---|---|
committer | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2013-03-23 20:08:12 +0000 |
commit | 0e045b8b3c7caeb62c491a4b3a6534890d4276cf (patch) | |
tree | 89c8398b7f17c9417b697b6152ed4b9b5e8bbb67 | |
parent | d9199f1361e1a22b63e59f11dcf575fc4b6fcfc4 (diff) |
* Low-density polytrope for tabulated cold EOS. This
part of the commit was forgotten.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@77 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
-rw-r--r-- | src/EOS_Omni_ColdEOSReadTable.F90 | 20 |
1 files changed, 16 insertions, 4 deletions
diff --git a/src/EOS_Omni_ColdEOSReadTable.F90 b/src/EOS_Omni_ColdEOSReadTable.F90 index b73bd4d..faee59d 100644 --- a/src/EOS_Omni_ColdEOSReadTable.F90 +++ b/src/EOS_Omni_ColdEOSReadTable.F90 @@ -19,7 +19,8 @@ subroutine EOS_Omni_ReadColdTable(CCTK_ARGUMENTS) LOGICAL :: tablethere character(len=256) :: buffer1,buffer2,buffer3 CCTK_INT :: lnrho,lnye,lntemp - integer :: i + CCTK_REAL :: temp_press + integer :: i ! convert file name to fortran string call CCTK_FortranString ( slength, coldeos_table_name, & @@ -47,6 +48,7 @@ subroutine EOS_Omni_ReadColdTable(CCTK_ARGUMENTS) endif coldeos_nrho = lnrho + ! allocate vars allocate(coldeos_logrho(coldeos_nrho)) allocate(coldeos_eps(coldeos_nrho)) @@ -55,8 +57,8 @@ subroutine EOS_Omni_ReadColdTable(CCTK_ARGUMENTS) ! read rho min and max read(667,"(A16,A16)") buffer1,buffer2 - read(buffer1(9:),*) coldeos_rhomin - read(buffer2(9:),*) coldeos_rhomax + read(buffer1(10:),*) coldeos_rhomin + read(buffer2(10:),*) coldeos_rhomax ! read heat capacity? Not sure what exactly that is used for, just ignore it ! for now @@ -74,6 +76,7 @@ subroutine EOS_Omni_ReadColdTable(CCTK_ARGUMENTS) read(667,"(A30)") buffer1 read(buffer1(9:),*) coldeos_kappa + ! make sure density is in log spacing read(667,"(A30)") buffer1 if(.not.(trim(adjustl(buffer1)).eq."RhoSpacing = Log")) then @@ -94,7 +97,15 @@ subroutine EOS_Omni_ReadColdTable(CCTK_ARGUMENTS) coldeos_logrho(i) = log10(coldeos_rhomin) + (i-1)*coldeos_dlrho enddo -#if 0 + ! set up low_kappa to be able to extrapolate to very low densities + ! we are using gamma=2 + coldeos_low_gamma = 2.0d0 + temp_press = (10.0d0**coldeos_logrho(1))**coldeos_gamma(1) + coldeos_low_kappa = temp_press / & + ((10.0d0**coldeos_logrho(1))**coldeos_low_gamma) + + +#if 0 ! debug output do i=1,coldeos_nrho write(6,"(i5,1P10E15.6)") i, coldeos_logrho(i), & @@ -102,4 +113,5 @@ subroutine EOS_Omni_ReadColdTable(CCTK_ARGUMENTS) enddo #endif + end subroutine EOS_Omni_ReadColdTable |