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