00001 subroutine calc_physq_center 00002 use def_matter, only : emd 00003 use def_matter_parameter, only : pinx 00004 use def_quantities, only : rho_max, pre_max, epsi_max, q_max 00005 implicit none 00006 ! 00007 rho_max = emd(0,0,0)**pinx 00008 pre_max = emd(0,0,0)**(pinx+1.0d0) 00009 epsi_max = emd(0,0,0)**pinx*(1.0d0+pinx*emd(0,0,0)) 00010 q_max = emd(0,0,0) 00011 ! 00012 end subroutine calc_physq_center