blob: 672b22e84c99f09c72d10bbfde9e8bab60160863 (
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
37
38
39
40
41
42
|
#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
if(poly_gamma_initial .gt. 0d0) then
poly_gamma_ini = poly_gamma_initial
else
poly_gamma_ini = poly_gamma
end if
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)
end subroutine EOS_Omni_Startup
subroutine EOS_Omni_Get_Energy_Shift(CCTK_ARGUMENTS)
use EOS_Omni_Module
implicit none
DECLARE_CCTK_PARAMETERS
DECLARE_CCTK_ARGUMENTS
call nuc_eos_c_get_energy_shift(energy_shift,eos_tempmin,eos_tempmax,&
eos_yemin,eos_yemax)
end subroutine EOS_Omni_Get_Energy_Shift
|