!#######################################################################
!     MODULE PROGRAM  --- mod_parameter_tillotson_eos ---            ###
!#######################################################################
!                                                                    ###
!     Parameter Sets for Tillotson EOS                               ###
!       (Meaning of parameters is described                          ###
!        in book 'Impact Createring' by Melosh (1989).)              ###
!                                                                    ###
!     Any Marerial  : imat = 1                                       ###
!     Murchison     : imat = 2                                       ###
!     Hard Ryugu    : imat = 3                                       ###
!     Soft Ryugu    : imat = 4                                       ###
!                                                                    ###
!                                                         2005/10/07 ###
!                                                     mod 2011/06/23 ###
!                                                     mod 2011/07/13 ###
!                                                     mod 2021/11/23 ###
!                                                     mod 2021/11/25 ###
!#######################################################################

module mod_parameter_tillotson_eos

  implicit none

!==== Reference Temperature for Tillotson EOS ==========================
  real(8), parameter :: T_REF = 220.0d0

!==== Switch for Temperature Calculation in Tillotson EOS ==============
  integer(4), parameter :: SWT_CV = 1  ! = 0 : Constant Cv (= CV_CONST)
                                       ! = 1 : T-dependent Cv
  
!=== Any Material (imat = 1) ===========================================
  real(8), parameter :: any_rho0=7800.0d0          ! [kg/m3]
  real(8), parameter :: any_a=0.5d0
  real(8), parameter :: any_b=1.5d0
  real(8), parameter :: any_ap=128.0d9             ! [Pa]
  real(8), parameter :: any_bp=105.0d9             ! [Pa]
  real(8), parameter :: any_e0=9.5d6               ! [J/kg]
  real(8), parameter :: any_alph=5.0d0
  real(8), parameter :: any_beta=5.0d0
  real(8), parameter :: any_eiv=2.4d6              ! [J/kg]
  real(8), parameter :: any_ecv=8.67d6             ! [J/kg]
  real(8), parameter :: any_debye=737.0d0          ! [K] Debye Temperature
  real(8), parameter :: any_cv_lim=1170.0d0        ! [J/K/kg] Dulong-Petit limit

!=== New Murchison created by Kurosawa (imat = 2) =====================
  real(8), parameter :: murch_rho0=2200.0d0        ! [kg/m3]
  real(8), parameter :: murch_a=0.5d0
  real(8), parameter :: murch_b=0.5d0
  real(8), parameter :: murch_ap=7.67d9            ! [Pa]
  real(8), parameter :: murch_bp=15.0d9            ! [Pa]
  real(8), parameter :: murch_e0=10.0d6            ! [J/kg]
  real(8), parameter :: murch_alph=5.0d0  
  real(8), parameter :: murch_beta=5.0d0  
  real(8), parameter :: murch_eiv=7.4d12           ! [J/kg]
  real(8), parameter :: murch_ecv=1.22d13          ! [J/kg]
  real(8), parameter :: murch_debye=688.0d0        ! [K] Debye Temperature
  real(8), parameter :: murch_cv_lim=1059.0d0      ! [J/K/kg] Dulong-Petit limit
  real(8), parameter :: murch_cv_const=930.0d0    ! [J/kg/K] Constant Cv
  
!=== New Hard Ryugu created by Kurosawa (imat = 3) =====================
  real(8), parameter :: hryugu_rho0=1820.0d0       ! [kg/m3]
  real(8), parameter :: hryugu_a=0.5d0
  real(8), parameter :: hryugu_b=-0.436d0
  real(8), parameter :: hryugu_ap=3.67d9           ! [Pa]
  real(8), parameter :: hryugu_bp=7.07d9           ! [Pa]
  real(8), parameter :: hryugu_e0=0.2d6            ! [J/kg]
  real(8), parameter :: hryugu_alph=5.0d0  
  real(8), parameter :: hryugu_beta=5.0d0  
  real(8), parameter :: hryugu_eiv=6.4d12          ! [J/kg]
  real(8), parameter :: hryugu_ecv=1.1d13          ! [J/kg]
  real(8), parameter :: hryugu_debye=737.0d0       ! [K] Debye Temperature
  real(8), parameter :: hryugu_cv_lim=1170.0d0     ! [J/K/kg] Dulong-Petit limit
  real(8), parameter :: hryugu_cv_const=880.0d0    ! [J/kg/K] Constant Cv
  
!=== New Soft Ryugu created by Kurosawa (imat = 4) =====================
  real(8), parameter :: sryugu_rho0=1820.0d0        ! [kg/m3]
  real(8), parameter :: sryugu_a=0.5d0
  real(8), parameter :: sryugu_b=-0.457d0
  real(8), parameter :: sryugu_ap=2.85d9            ! [Pa]
  real(8), parameter :: sryugu_bp=5.3d7             ! [Pa]
  real(8), parameter :: sryugu_e0=0.2d6             ! [J/kg]
  real(8), parameter :: sryugu_alph=5.0d0
  real(8), parameter :: sryugu_beta=5.0d0
  real(8), parameter :: sryugu_eiv=6.4d12           ! [J/kg]
  real(8), parameter :: sryugu_ecv=1.1d13           ! [J/kg]
  real(8), parameter :: sryugu_debye=737.0d0        ! [K] Debye Temperature
  real(8), parameter :: sryugu_cv_lim=1170.0d0      ! [J/K/kg] Dulong-Petit limit
  real(8), parameter :: sryugu_cv_const=880.0d0     ! [J/kg/K] Constant Cv
  
!======================================================================

end module mod_parameter_tillotson_eos