diff options
author | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2010-10-24 14:06:52 +0000 |
---|---|---|
committer | cott <cott@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2010-10-24 14:06:52 +0000 |
commit | 5621733be96b97235891c5da5e1aaa52967058c3 (patch) | |
tree | 11f840e137111193f641e553600c86e7cebc310c /src/EOS_Omni_Startup.F90 | |
parent | 7449c12a9f5976d5d9eff74ea4a43741a64e8cb3 (diff) |
* add capability allowing to read in finite-T EOS tables.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/EOS_Omni@11 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
Diffstat (limited to 'src/EOS_Omni_Startup.F90')
-rw-r--r-- | src/EOS_Omni_Startup.F90 | 17 |
1 files changed, 15 insertions, 2 deletions
diff --git a/src/EOS_Omni_Startup.F90 b/src/EOS_Omni_Startup.F90 index d13a0f6..d080e3e 100644 --- a/src/EOS_Omni_Startup.F90 +++ b/src/EOS_Omni_Startup.F90 @@ -10,6 +10,8 @@ subroutine EOS_Omni_Startup(CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS DECLARE_CCTK_ARGUMENTS + character(len=512) eosfilename + CCTK_INT fslen poly_k_cgs = poly_k * rho_gf**poly_gamma_ini / press_gf @@ -17,7 +19,18 @@ subroutine EOS_Omni_Startup(CCTK_ARGUMENTS) hybrid_k1_cgs = hybrid_k1 * rho_gf**poly_gamma_ini / press_gf - hybrid_k2_cgs = hybrid_k1_cgs * (hybrid_rho_nuc * inv_rho_gf)**(hybrid_gamma1-hybrid_gamma2) - + hybrid_k2_cgs = hybrid_k1_cgs * & + (hybrid_rho_nuc * inv_rho_gf)**(hybrid_gamma1-hybrid_gamma2) + + if(nuceos_read_table.ne.0) then + ! call EOS table reader + call CCTK_FortranString(fslen,nuceos_table_name,eosfilename) + call CCTK_INFO("##################################################") + call CCTK_INFO("EOS_Omni: Reading finite-T nuclear EOS table") + call CCTK_INFO(eosfilename) + call CCTK_INFO("##################################################") + call nuc_eos_readtable(eosfilename) + endif + end subroutine EOS_Omni_Startup |