00001 function quark_h2rho_guo(hh)
00002 use def_qeos_parameter
00003 implicit none
00004 real(8) :: hh, quark_h2rho_guo
00005 integer :: ii
00006 quark_h2rho_guo=((hh-1.0d0-eneconst_gcm1)/(abch(1)+abch(1+nphase)))**&
00007 &(1.0d0/abih(1))
00008 return
00009 END function quark_h2rho_guo