Cocal: /home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/calc_physq_center

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, emdmax, xsol
00006   integer    :: ir, iremax
00007 !
00008   call search_emdmax_xaxis_grid(iremax)
00009   if (iremax.eq.0) then 
00010     q_c = emd(0,0,0)
00011   else
00012     call search_emdmax_xaxis(iremax,emdmax,xsol)
00013     q_c = emdmax
00014   end if
00015   call peos_q2hprho(q_c, hh, pre_c, rho_c, epsi_c)
00016 !
00017 end subroutine calc_physq_center_peos

/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/calc_physq_center_peos.f90