00001 subroutine read_parameter_bh_mpt(impt)
00002 use phys_constant, only : long, pi
00003 use grid_parameter, only : num_sol_seq
00004 use def_bh_parameter
00005 use def_binary_parameter, only : mass_ratio
00006 implicit none
00007 integer,intent(in) :: impt
00008 character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
00009 open(1,file='bbhpar_mpt'//np(impt)//'.dat',status='old')
00010 read(1,'(1p,2e14.6)') ome_bh, spin_bh
00011 read(1,'(1p,2e14.6)') alph_bh, psi_bh
00012 read(1,'(1p,2e14.6)') th_spin_bh_deg, phi_spin_bh_deg
00013 read(1,'(2(2x,a2),6x,1p,e14.6)') bh_bctype, bh_soltype, mass_ratio
00014 read(1,'(1p,e14.6,i5)') rgin_deform, num_sol_seq
00015 close(1)
00016 th_spin_bh = pi*th_spin_bh_deg/180.0d0
00017 phi_spin_bh = pi*phi_spin_bh_deg/180.0d0
00018 bh_sptype = 'SP'
00019 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq. 0.0) bh_sptype = 'Xp'
00020 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq. 90.0) bh_sptype = 'Yp'
00021 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq.180.0) bh_sptype = 'Xm'
00022 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq.270.0) bh_sptype = 'Ym'
00023 if (th_spin_bh_deg.eq. 0.0) bh_sptype = 'Zp'
00024 if (th_spin_bh_deg.eq.180.0) bh_sptype = 'Zm'
00025 end subroutine read_parameter_bh_mpt
