00001 subroutine printout_physq_BH(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 write(110,'(a22,1p,2e23.15)') ' Radius & dis 1st BH= ', rgin, dis
00019 write(110,'(a22,1p,2e23.15)') ' Separation = ', sepa
00020
00021 write(110,*) '## BH orbital and spin angular velocity parameters ##'
00022 write(110,'(a22,1p,2e23.15)') ' Omega = ', ome_bh
00023 write(110,'(a22,1p,2e23.15)') ' Spin 1st BH = ', spin_bh
00024 write(110,'(a22,1p,2e23.15)') ' Spin axis 1st BH = ', th_spin_bh_deg, &
00025 & phi_spin_bh_deg
00026 write(110,'(a22,1p,2e23.15)') ' Puncture spin Sx = ', spin_pBH(1)
00027 write(110,'(a22,1p,2e23.15)') ' Puncture spin Sy = ', spin_pBH(2)
00028 write(110,'(a22,1p,2e23.15)') ' Puncture spin Sz = ', spin_pBH(3)
00029 write(110,*) '## BH linear momentum parameters ##'
00030 write(110,'(a22,1p,2e23.15)') ' Px asympto & punct = ', admmom_asymp(1), &
00031 & mom_pBH(1)
00032 write(110,'(a22,1p,2e23.15)') ' Py asympto & punct = ', admmom_asymp(2), &
00033 & mom_pBH(2)
00034 write(110,'(a22,1p,2e23.15)') ' Pz asympto & punct = ', admmom_asymp(3), &
00035 & mom_pBH(3)
00036
00037 write(110,*) '## Mass and angular momentum (G = c = 1 unit) ##'
00038 write(110,'(a22,1p,2e23.15)') ' Puncture mass para = ', mass_pBH
00039 write(110,'(a22,1p,2e23.15)') ' M_ADM (asymp & vol)= ', admmass_asymp, &
00040 & admmass
00041 write(110,'(a22,1p,2e23.15)') ' M_K (asymp & vol)= ', komarmass_asymp, &
00042 & komarmass
00043 write(110,'(a22,1p,2e23.15)') ' J (asymptotic) = ', angmom_asymp
00044 write(110,'(a22,1p,2e23.15)') ' AH mass 1st BH = ', AHmass
00045 write(110,'(a22,1p,2e23.15)') ' AH area 1st BH = ', AHarea
00046 write(110,'(a22,1p,2e23.15)') ' AH spin 1st BH = ', AHspin
00047
00048 fixedvir = (admmass_asymp - komarmass_asymp)/admmass_asymp
00049 write(110,'(a22,1p,2e23.15)') ' 1 - M_K/M_ADM = ', fixedvir
00050
00051 end subroutine printout_physq_BH
