00001 Subroutine quark_rho2phenedpdrho(rho,pre,hh,ene,dpdrho)
00002 use def_qeos_parameter
00003 implicit none
00004 integer :: ii
00005 real(8),intent(inout) :: rho, pre, hh, ene, dpdrho
00006 pre=0.0d0
00007 hh =0.0d0
00008 ene=0.0d0
00009 dpdrho=0.0d0
00010 if (rho.le.0.0d0) rho=0.0d0
00011 do ii=1,nphase
00012 pre=abc(ii)*rho**abi(ii)+pre
00013 hh =abch(ii)*rho**abih(ii)+hh
00014 ene=abcene(ii)*rho**abiene(ii)+ene
00015 dpdrho=abi(ii)*abc(ii)*rho**(abi(ii)-1.0d0)+dpdrho
00016 end do
00017 do ii=nphase+1,2*nphase
00018 hh=abch(ii)*rho**abih(ii)+hh
00019 end do
00020 ene=ene+rho+rho*eneconst_gcm1
00021 hh=hh+1.0d0+eneconst_gcm1
00022 END Subroutine quark_rho2phenedpdrho
