00001 subroutine read_parameter_bh
00002 use phys_constant, only : long, pi
00003 use def_bh_parameter
00004 implicit none
00005 open(1,file='bbhpar.dat',status='old')
00006 read(1,'(1p,2e14.6)') ome_bh, spin_bh
00007 read(1,'(1p,2e14.6)') alph_bh, psi_bh
00008 read(1,'(1p,2e14.6)') th_spin_bh_deg, phi_spin_bh_deg
00009 read(1,'(2(2x,a2))') bh_bctype, bh_soltype
00010 read(1,'(1p,2e14.6)') mass_pBH
00011 close(1)
00012 th_spin_bh = pi*th_spin_bh_deg/180.0d0
00013 phi_spin_bh = pi*phi_spin_bh_deg/180.0d0
00014 bh_sptype = 'SP'
00015 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq. 0.0) bh_sptype = 'Xp'
00016 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq. 90.0) bh_sptype = 'Yp'
00017 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq.180.0) bh_sptype = 'Xm'
00018 if (th_spin_bh_deg.eq.90.0.or.phi_spin_bh_deg.eq.270.0) bh_sptype = 'Ym'
00019 if (th_spin_bh_deg.eq. 0.0) bh_sptype = 'Zp'
00020 if (th_spin_bh_deg.eq.180.0) bh_sptype = 'Zm'
00021 end subroutine read_parameter_bh
