00001 subroutine write_last_physq
00002 use def_quantities
00003 use def_bh_parameter
00004 use grid_parameter, only : rgin
00005 implicit none
00006
00007 open(15,file='last_physq.txt' ,status='unknown')
00008 write(15,'(1p,13e16.8)') rgin, spin_bh, ome_bh, admmass, komarmass, admmass_thr, &
00009 & app_hor_area_bh, irredmass, bindingene, &
00010 & angmom, angmom_thr, angmom_smarr, qua_loc_spin
00011
00012 write (15,'(a30,1p,e16.8)') ' Radius of bh = ', rgin
00013 write (15,'(a30,1p,e16.8)') ' Spin of = ', spin_bh
00014 write (15,'(a30,1p,e16.8)') ' Orbital angular velocity = ', ome_bh
00015 write (15,'(a30,1p,e16.8)') ' ADM mass = ', admmass
00016 write (15,'(a30,1p,e16.8)') ' Komar mass = ', komarmass
00017 write (15,'(a30,1p,e16.8)') ' ADM mass (throats) = ', admmass_thr
00018 write (15,'(a30,1p,e16.8)') ' Apparent horizon area = ', app_hor_area_bh
00019 write (15,'(a30,1p,e16.8)') ' Irreducible mass = ', irredmass
00020 write (15,'(a30,1p,e16.8)') ' Binding energy = ', bindingene
00021 write (15,'(a30,1p,e16.8)') ' Angular momentum = ', angmom
00022 write (15,'(a30,1p,e16.8)') ' Angular momentum (throats) = ', angmom_thr
00023 write (15,'(a30,1p,e16.8)') ' Angular momentum (Smarr) = ', angmom_smarr
00024 write (15,'(a30,1p,e16.8)') ' Quasi local spin BH = ', qua_loc_spin
00025 write (15,*) '#'
00026 close(15)
00027
00028 end subroutine write_last_physq
