00001 subroutine printout_physq_BBH_mpt(iseq)
00002 use phys_constant, only : long, nmpt
00003 use grid_parameter, only : rgin
00004 use def_quantities
00005 use def_bh_parameter
00006 use def_binary_parameter, only : sepa, dis
00007 use def_quantities_bh
00008 implicit none
00009 real(long) :: fixedvir
00010 integer, intent(in) :: iseq
00011 integer :: impt
00012
00013 if (iseq.eq.1) then
00014 open(110,file='bbhphyseq.dat',status='unknown')
00015 end if
00016 write(110,*) '== Sequence number == ', iseq
00017 write(110,*) '## BH Radii in G = c = 1 unit ##'
00018 call copy_grid_parameter_from_mpt(1)
00019 call copy_def_binary_parameter_from_mpt(1)
00020 write(110,'(a22,1p,2e23.15)') ' Radius & dis 1st BH= ', rgin, dis
00021 call copy_grid_parameter_from_mpt(2)
00022 call copy_def_binary_parameter_from_mpt(2)
00023 write(110,'(a22,1p,2e23.15)') ' Radius & dis 2nd BH= ', rgin, dis
00024 call copy_def_binary_parameter_from_mpt(1)
00025 write(110,'(a22,1p,2e23.15)') ' Separation = ', sepa
00026
00027 write(110,*) '## BBH orbital and spin angular velocity parameters ##'
00028 call copy_def_bh_parameter_from_mpt(1)
00029 write(110,'(a22,1p,2e23.15)') ' Omega = ', ome_bh
00030 write(110,'(a22,1p,2e23.15)') ' Spin 1st BH = ', spin_bh
00031 write(110,'(a22,1p,2e23.15)') ' Spin axis 1st BH = ', th_spin_bh_deg, &
00032 & phi_spin_bh_deg
00033 call copy_def_bh_parameter_from_mpt(2)
00034 write(110,'(a22,1p,2e23.15)') ' Spin 2nd BH = ', spin_bh
00035 write(110,'(a22,1p,2e23.15)') ' Spin axis 2nd BH = ', th_spin_bh_deg, &
00036 & phi_spin_bh_deg
00037 write(110,*) '## Mass and angular momentum (G = c = 1 unit) ##'
00038 call copy_def_quantities_from_mpt(nmpt)
00039 write(110,'(a22,1p,2e23.15)') ' M_ADM (asymptotic) = ', admmass_asymp
00040 write(110,'(a22,1p,2e23.15)') ' M_K Komar mass = ', komarmass_asymp
00041 write(110,'(a22,1p,2e23.15)') ' J (asymptotic) = ', angmom_asymp
00042 call copy_def_quantities_bh_from_mpt(1)
00043 write(110,'(a22,1p,2e23.15)') ' AH mass 1st BH = ', AHmass
00044 write(110,'(a22,1p,2e23.15)') ' AH area 1st BH = ', AHarea
00045 write(110,'(a22,1p,2e23.15)') ' AH spin 1st BH = ', AHspin
00046 call copy_def_quantities_bh_from_mpt(2)
00047 write(110,'(a22,1p,2e23.15)') ' AH mass 2nd BH = ', AHmass
00048 write(110,'(a22,1p,2e23.15)') ' AH area 2nd BH = ', AHarea
00049 write(110,'(a22,1p,2e23.15)') ' AH spin 2nd BH = ', AHspin
00050
00051 fixedvir = (admmass_asymp - komarmass_asymp)/admmass_asymp
00052 write(110,'(a22,1p,2e23.15)') ' 1 - M_K/M_ADM = ', fixedvir
00053
00054 end subroutine printout_physq_BBH_mpt
