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

Go to the documentation of this file.

subroutine calc_physq_center_peos_grid
  use grid_parameter, only  :  nrf, ntf, npf
  use def_matter, only : emd
  use def_quantities, only : rho_c, pre_c, epsi_c, q_c, &
  &                          rho_max, pre_max, epsi_max, q_max 
  implicit none
  real(8) :: hh, emdmax=-1.0d0
  integer :: irf,itf,ipf, irfmax, itfmax, ipfmax
!
  q_c = emd(0,0,0)
  call peos_q2hprho(q_c, hh, pre_c, rho_c, epsi_c)
  write(6,'(a6,1p,e23.15,a21,i3,a1,i3,a1,i3,a1)') "q_c  =", q_c,"   at (irf,itf,ipf)=(", &
  &                                                0, ",", 0, ",", 0,")"
!
  q_max = 0.0d0
  do irf=0, nrf
    do itf=0, ntf
      do ipf=0, npf
        if ( emd(irf,itf,ipf) > q_max )  then
          irfmax = irf
          itfmax = itf
          ipfmax = ipf
          q_max = emd(irf,itf,ipf)
        end if
      end do
    end  do 
  end do
!
  write(6,'(a6,1p,e23.15,a21,i3,a1,i3,a1,i3,a1)') "q_max=", q_max,"   at (irf,itf,ipf)=(", &
  &                                                irfmax, ",", itfmax, ",", ipfmax,")"
!
  call peos_q2hprho(q_max, hh, pre_max, rho_max, epsi_max)
!
end subroutine calc_physq_center_peos_grid

---