














































 














































































 


































































































































































 
























































































































































































 























































 










 








































































































































































































































 










                                
                                

                                
                                

                                
                                

                                















































 

      integer function EOS_Hybrid_Startup()

      USE EOS_Polytrope_Scalars
      USE EOS_Hybrid_Scalars

      implicit none

      REAL*8 eos_gamma_supernuclear
      integer, parameter :: cctki_use_eos_gamma_supernuclear = kind(eos_
     &gamma_supernuclear)
      REAL*8 eos_gamma_th
      integer, parameter :: cctki_use_eos_gamma_th = kind(eos_gamma_th)
      REAL*8 rho_nuc
      integer, parameter :: cctki_use_rho_nuc = kind(rho_nuc)
      COMMON /EOS_Hybridrest/eos_gamma_supernuclear, eos_gamma_th, rho_n
     &uc
      REAL*8 eos_gamma
      integer, parameter :: cctki_use_eos_gamma = kind(eos_gamma)
      REAL*8 eos_k
      integer, parameter :: cctki_use_eos_k = kind(eos_k)
      REAL*8 gamma_ini
      integer, parameter :: cctki_use_gamma_ini = kind(gamma_ini)
      INTEGER*4 CCTKH0
      integer, parameter :: cctki_use_CCTKH0 = kind(CCTKH0)
      COMMON /EOS_2D_POLYTROPErest/eos_gamma, eos_k, gamma_ini, CCTKH0

      integer handle, ierr
      
      external EOS_Hybrid_Pressure
      external EOS_Hybrid_SpecificIE
      external EOS_Hybrid_RestMassDens
      external EOS_Hybrid_DPressByDRho
      external EOS_Hybrid_DPressByDEps

      call EOS_RegisterMethod(handle,"EOS_Hybrid")

      if (handle .ge. 0) then

        call EOS_RegisterPressure(ierr, handle, EOS_Hybrid_Pressure)
        call EOS_RegisterSpecificIntEnergy(ierr, handle, 
     .                                     EOS_Hybrid_SpecificIE)
        call EOS_RegisterRestMassDens(ierr, handle, 
     .                                EOS_Hybrid_RestMassDens)
        call EOS_RegisterDPressByDRho(ierr, handle, 
     .                                EOS_Hybrid_DPressByDRho)
        call EOS_RegisterDPressByDEps(ierr, handle, 
     .                                EOS_Hybrid_DPressByDEps)

      else

        call CCTK_Warn(0,60,"EOS_Hybrid_Startup.F","EOS_Hybrid", "Unable
     & to register the EOS method!")

      end if

c

c

c
c
c
c
c
c
c
c
c
c

      eos_k_supernuclear_cgs = eos_k_cgs * (rho_nuc * rho_geom_factor_in
     &v)**
     .     (eos_gamma - eos_gamma_supernuclear)

      EOS_Hybrid_Startup = 0
        
      end function EOS_Hybrid_Startup
