














































 










































































 


































































































































































 























































































































































































 





















































 










 








































































































































































































































 










                                
                                

                                
                                

                                
                                

                                















































 

      REAL*8 function EOS_Polytrope_Pressure(rho, eps)

      USE EOS_Polytrope_Scalars

      implicit none
      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 use_cgs
      integer, parameter :: cctki_use_use_cgs = kind(use_cgs)
      COMMON /EOS_2d_Polytroperest/eos_gamma, eos_k, gamma_ini, use_cgs

      REAL*8 rho
      REAL*8  eps
      integer, parameter :: cctki_use_ eps = kind( eps)

      EOS_Polytrope_Pressure = p_geom_factor * eos_k_cgs * 
     .     (rho * rho_geom_factor_inv) ** eos_gamma
      
      end function EOS_Polytrope_Pressure

 

      REAL*8 function EOS_Polytrope_SpecificIE(rho, press)

      USE EOS_Polytrope_Scalars

      implicit none
      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 use_cgs
      integer, parameter :: cctki_use_use_cgs = kind(use_cgs)
      COMMON /EOS_2d_Polytroperest/eos_gamma, eos_k, gamma_ini, use_cgs

      REAL*8 rho, press, EOS_Polytrope_Pressure

      EOS_Polytrope_SpecificIE = EOS_Polytrope_Pressure(rho, press) 
     .     / ( rho * (eos_gamma - 1.d0) )
      
      end function EOS_Polytrope_SpecificIE

 

      REAL*8 function EOS_Polytrope_RestMassDens(eps, press)

      USE EOS_Polytrope_Scalars

      implicit none
      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 use_cgs
      integer, parameter :: cctki_use_use_cgs = kind(use_cgs)
      COMMON /EOS_2d_Polytroperest/eos_gamma, eos_k, gamma_ini, use_cgs

      REAL*8 eps, press

      EOS_Polytrope_RestMassDens = press / ((eos_gamma - 1.d0) * eps)
      
      end function EOS_Polytrope_RestMassDens

 
      
      REAL*8 function EOS_Polytrope_DPressByDRho(rho, eps)

      USE EOS_Polytrope_Scalars

      implicit none
      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 use_cgs
      integer, parameter :: cctki_use_use_cgs = kind(use_cgs)
      COMMON /EOS_2d_Polytroperest/eos_gamma, eos_k, gamma_ini, use_cgs

      REAL*8 rho
      REAL*8  eps
      integer, parameter :: cctki_use_ eps = kind( eps)

      EOS_Polytrope_DPressByDRho = p_geom_factor * eos_k_cgs * 
     .     eos_gamma * rho_geom_factor_inv *  
     .     (rho*rho_geom_factor_inv) ** (eos_gamma - 1.d0)

      end function EOS_Polytrope_DPressByDRho

 

      REAL*8 function EOS_Polytrope_DPressByDEps(rho, eps)

      USE EOS_Polytrope_Scalars

      implicit none
      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 use_cgs
      integer, parameter :: cctki_use_use_cgs = kind(use_cgs)
      COMMON /EOS_2d_Polytroperest/eos_gamma, eos_k, gamma_ini, use_cgs

      REAL*8  rho
      integer, parameter :: cctki_use_ rho = kind( rho)
      REAL*8  eps
      integer, parameter :: cctki_use_ eps = kind( eps)

      EOS_Polytrope_DPressByDEps = 0.d0

      end function EOS_Polytrope_DPressByDEps

