blob: 4601e9225c448ebd5981c8fb7c2c0610eb2b638d (
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
43
44
45
46
|
/*@@
@file EOS_GP_Setup.F90
@date Mon Mar 14 17:41:20 2005
@author Ian Hawke
@desc
Setup the scalar variables
@enddesc
@@*/
#include "cctk.h"
#include "cctk_Arguments.h"
#include "cctk_Parameters.h"
subroutine EOS_GP_Setup()
USE EOS_GP_Scalars
implicit none
DECLARE_CCTK_PARAMETERS
m_solar_cgs = 1.987d33
c_cgs = 2.9979d10
G_cgs = 6.6732d-8
m_solar_geom = G_cgs / c_cgs**2 * m_solar_cgs
eos_gamma_local = eos_gamma
if (use_cgs .ne. 0) then
p_geom_factor = G_cgs / c_cgs**4 * m_solar_geom**2
rho_geom_factor = p_geom_factor * c_cgs**2
rho_geom_factor_inv = 1.d0 / rho_geom_factor
eos_k_cgs = eos_k * rho_geom_factor**gamma_ini / p_geom_factor
else
p_geom_factor = 1.d0
rho_geom_factor = 1.d0
rho_geom_factor_inv = 1.d0
eos_k_cgs = eos_k
end if
end subroutine EOS_GP_Setup
|