/*@@
   @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