/home/uryu/Cocal_svn_checkinout/Cocal/code/EOS/Subroutine/teos_h2qprho.f90

Go to the documentation of this file.

subroutine teos_h2qprho(h,q,pre,rho,ened)
!
  use def_teos_parameter   !rhoi,qi,hi,nphase
  implicit none
!
  real(8), intent(inout)  :: h
  real(8), intent(out) :: q, pre, rho, ened
  real(8), external :: fn_lagint, fn_lagint_2nd
  real(8) :: x4(4), f4(4), x2(2), f2(2)
  integer              :: i, i0, ii
  integer, save        :: iphase
!
  if (h < hi(0)) then
    q=0.0d0;  h=1.0d0;  pre=0.0d0;  rho=0.0d0;  ened=0.0d0
    return
  end if

  ii=iphase
  if(h.lt.hi(iphase)) then
    do i = ii, 1, -1
      if( (h.ge.hi(i-1)) .and. (h.le.hi(i))) then
        iphase = i
        exit
      end if
    end do
  else
    do i = ii, nphase
      if( (h.ge.hi(i-1)) .and. (h.le.hi(i))) then
        iphase = i
        exit
      end if
    end do
  end if

!  do i = 1, nphase
!    if( (h.ge.hi(i-1)) .and. (h.le.hi(i))) then
!      iphase = i
!      exit
!    end if
!  end do
!   write(6,*) "h2qprho: iphase=", iphase
!
  i0 = min0(max0(iphase-2,0),nphase-3)
  x4(1:4) = hi(i0:i0+3)
  f4(1:4) = qi(i0:i0+3)
  q       = fn_lagint(x4,f4,h)

  f4(1:4) = prei(i0:i0+3)
  pre     = fn_lagint(x4,f4,h)

  f4(1:4) = rhoi(i0:i0+3)
  rho     = fn_lagint(x4,f4,h)

  f4(1:4) = enei(i0:i0+3)
  ened    = fn_lagint(x4,f4,h)

!
end subroutine teos_h2qprho

---