00001 subroutine read_parameter_binary_excision_mpt(impt)
00002 use phys_constant, only : long
00003 use grid_parameter, only : nrf
00004 use coordinate_grav_r, only : rg
00005 use grid_parameter_binary_excision, only : ex_nrg, ex_ndis
00006 implicit none
00007 integer,intent(in) :: impt
00008 character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
00009 open(1,file='bin_ex_par_mpt'//np(impt)//'.dat',status='old')
00010 read(1,'(2i5)') ex_nrg, ex_ndis
00011 close(1)
00012
00013 if (ex_nrg.ne.0) then
00014 if (nrf.ge.ex_nrg) write(6,*) '** Warning ** nrf > ex_nrg '
00015 if (nrf.ge.ex_nrg) write(6,*) 'nrf = ', nrf, ' ex_nrg = ', ex_nrg
00016 end if
00017 if (ex_nrg.eq.0) write(6,*) 'No binary excision for patch =',impt
00018 end subroutine read_parameter_binary_excision_mpt
00019
00020 subroutine read_parameter_binary_excision_mpt_cactus(impt, dir_path)
00021 use phys_constant, only : long
00022 use grid_parameter, only : nrf
00023 use coordinate_grav_r, only : rg
00024 use grid_parameter_binary_excision, only : ex_nrg, ex_ndis
00025 implicit none
00026 integer,intent(in) :: impt
00027 character*400, intent(in) :: dir_path
00028 character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
00029 open(1,file=trim(dir_path)//'/'//'bin_ex_par_mpt'//np(impt)//'.dat',status='old')
00030 read(1,'(2i5)') ex_nrg, ex_ndis
00031 close(1)
00032
00033 if (ex_nrg.ne.0) then
00034 if (nrf.ge.ex_nrg) write(6,*) '** Warning ** nrf > ex_nrg '
00035 if (nrf.ge.ex_nrg) write(6,*) 'nrf = ', nrf, ' ex_nrg = ', ex_nrg
00036 end if
00037 if (ex_nrg.eq.0) write(6,*) 'No binary excision for patch =',impt
00038 end subroutine read_parameter_binary_excision_mpt_cactus
