#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
  character(len=512) eosfilename
  CCTK_INT fslen

  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)

  if(nuceos_read_table.ne.0) then
     ! call EOS table reader
     call CCTK_FortranString(fslen,nuceos_table_name,eosfilename)
     call CCTK_INFO("##################################################")
     call CCTK_INFO("EOS_Omni: Reading finite-T nuclear EOS table")
     call CCTK_INFO(eosfilename)
     call CCTK_INFO("##################################################")
     call nuc_eos_readtable(eosfilename)
  endif
  

end subroutine EOS_Omni_Startup