diff options
author | knarf <knarf@a192deed-9f5d-41fb-bb60-7f59df5b1b8d> | 2009-12-02 22:22:53 +0000 |
---|---|---|
committer | knarf <knarf@a192deed-9f5d-41fb-bb60-7f59df5b1b8d> | 2009-12-02 22:22:53 +0000 |
commit | 9bdbc9209f0004429393aad47d0973046a4ee0dc (patch) | |
tree | 1468bb1aa622e40126a1f25dd42b7ffe6b23ccdc /src/EOS_GP_Setup.F90 | |
parent | 794e79d409732b82bb8871444d6e17c80b0e879c (diff) |
thorns for general eos interface (moved from Whisky_Dev repositories)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOSG_Polytrope/trunk@2 a192deed-9f5d-41fb-bb60-7f59df5b1b8d
Diffstat (limited to 'src/EOS_GP_Setup.F90')
-rwxr-xr-x | src/EOS_GP_Setup.F90 | 46 |
1 files changed, 46 insertions, 0 deletions
diff --git a/src/EOS_GP_Setup.F90 b/src/EOS_GP_Setup.F90 new file mode 100755 index 0000000..4601e92 --- /dev/null +++ b/src/EOS_GP_Setup.F90 @@ -0,0 +1,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 |