00001 subroutine calc_physq_center_peos
00002   use def_matter, only : emd
00003   use def_quantities, only : rho_c, pre_c, epsi_c, q_c 
00004   implicit none
00005   real(8) :: hh
00006 
00007   q_c = emd(0,0,0)
00008   call peos_q2hprho(q_c, hh, pre_c, rho_c, epsi_c)
00009 
00010 end subroutine calc_physq_center_peos