blob: d080e3ef0fae757b0f61469a0c25d7d3b616446b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"
#include "cctk_Functions.h"
subroutine EOS_Omni_Startup(CCTK_ARGUMENTS)
use EOS_Omni_Module
implicit none
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
gl_k_cgs = gl_k * rho_gf**poly_gamma_ini / press_gf
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)
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
|
