/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/calc_soundspeed_peos.f90

Go to the documentation of this file.

subroutine calc_soundspeed_peos(c_s)
  use phys_constant, only : long
  use def_peos_parameter     ! abc,abi,rhoi,qi,hi,nphase
  use grid_parameter, only : nrf, ntf, npf, ntfeq, ntfpolp, npfxzp, npfyzp
  use def_matter, only : emd
  use make_array_3d
  implicit none
  real(long), pointer  :: c_s(:,:,:)
  real(long) :: qq, hh, pre, rho, epsi, abin
  integer :: irf, itf, ipf, iphase
!
  do ipf = 0, npf
    do itf = 0, ntf
      do irf = 0, nrf
        qq = emd(irf,itf,ipf)
        call peos_q2hprho(qq, hh, pre, rho, epsi)
        call peos_lookup(qq, qi, iphase)
        abin = abi(iphase)
        c_s(irf,itf,ipf) = sqrt(abin*pre/rho/hh)
      end do
    end do
  end do
!
!
end subroutine calc_soundspeed_peos

---