!====================================================================
!   SUBROUTINE   ---- sub_eos_saumon_read_table ----              ===
!                                                                 ===
!   Subroutine for reading table of Saumon EOS                    ===
!                                                      2012/09/12 ===
!====================================================================

subroutine sub_eos_saumon_read_table
  use mod_parameter_saumon_eos
  implicit none

!---- Variables Used in This Subroutine -----------------------------
  integer(4) :: i,j
!--------------------------------------------------------------------

  open(50,file=TABLE_SAUMON_EOS)

  read(50,'(d15.5)') d_log_min_Saumon
  read(50,'(d15.5)') d_log_max_Saumon
  read(50,'(d15.5)') dd_log_Saumon
  read(50,'(i7)')    id_Saumon
  read(50,'(d15.5)') e_log_min_Saumon
  read(50,'(d15.5)') e_log_max_Saumon
  read(50,'(d15.5)') de_log_Saumon
  read(50,'(i7)')    ie_Saumon
  read(50,'(d15.5)') p_log_out_Saumon

  !---- Allocate ----------------------------------
  allocate(d_log_Saumon(0:id_Saumon))
  allocate(e_log_Saumon(0:ie_Saumon))
  allocate(p_log_Saumon(0:id_Saumon,0:ie_Saumon))
  allocate(t_log_Saumon(0:id_Saumon,0:ie_Saumon))
  !------------------------------------------------

  do i=0,id_Saumon
     read(50,'(d15.5)') d_log_Saumon(i)
     do j=0,ie_Saumon
        read(50,'(3d15.5)') e_log_Saumon(j),    &
        &                   p_log_Saumon(i,j),  &
        &                   t_log_Saumon(i,j)
     end do
  end do

  close(50)

  return
end subroutine sub_eos_saumon_read_table