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

Go to the documentation of this file.

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

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


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

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

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

  f4(1:4) = enei(i0:i0+3)
!  write(6,*) "ene: x4:", (x4(i),i=1,4)
!  write(6,*) "ene: f4:", (f4(i),i=1,4)
  ened    = fn_lagint(x4,f4,q)
!
end subroutine teos_q2hprho

---