/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine_mpatch/read_parameter_binary_excision_mpt.f90

Go to the documentation of this file.

subroutine read_parameter_binary_excision_mpt(impt)
  use phys_constant, only : long
  use grid_parameter, only : nrf
  use coordinate_grav_r, only : rg
  use grid_parameter_binary_excision, only : ex_nrg, ex_ndis
  implicit none
  integer,intent(in)  :: impt
  character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
  open(1,file='bin_ex_par_mpt'//np(impt)//'.dat',status='old')
  read(1,'(2i5)') ex_nrg, ex_ndis
  close(1)
!  if (nrf.ge.ex_nrg) stop ' nrf > ex_nrg '
  if (ex_nrg.ne.0) then 
    if (nrf.ge.ex_nrg) write(6,*) '** Warning ** nrf > ex_nrg '
    if (nrf.ge.ex_nrg) write(6,*) 'nrf = ', nrf, '  ex_nrg = ', ex_nrg
  end if
  if (ex_nrg.eq.0) write(6,*) 'No binary excision for patch =',impt 
end subroutine read_parameter_binary_excision_mpt

---