00001 function quark_rho2h_dot(rho)
00002 use def_qeos_parameter
00003 implicit none
00004 real(8) :: rho, quark_rho2h_dot
00005 integer :: ii
00006 quark_rho2h_dot=0.0d0
00007 do ii=1,2*nphase
00008 quark_rho2h_dot=abchdot(ii)*rho**abihdot(ii)+quark_rho2h_dot
00009 end do
00010 return
00011 END function quark_rho2h_dot