00001 subroutine write_last_physq_BNS_mpt
00002 use phys_constant, only : long, nmpt
00003 use def_quantities_mpt
00004 use def_matter_parameter_mpt
00005 use grid_parameter_mpt
00006 implicit none
00007 real(long) :: emdenc, omega, radius, m_rest, j_star, m_adm, m_kom, j_adm
00008 real(long) :: m_grav_sph, m_rest_sph, M_R_sph, rhoc, qc, rhomax, qmax
00009 real(long) :: chicusp, rsurf, rgmid
00010 character(40) :: char1, char2, char3, char4, char5
00011 character(100) :: dircommand
00012
00013 emdenc = def_matter_param_real_(2,1)
00014 omega = def_matter_param_real_(3,1)
00015 radius = def_matter_param_real_(5,1)
00016 m_rest = def_quantities_real_(4,1)
00017 j_star = def_quantities_real_(6,1)
00018 m_adm = def_quantities_real_(7,nmpt)
00019 m_kom = def_quantities_real_(8,nmpt)
00020 j_adm = def_quantities_real_(9,nmpt)
00021 m_grav_sph = def_quantities_real_(28,1)
00022 m_rest_sph = def_quantities_real_(29,1)
00023 M_R_sph = def_quantities_real_(31,1)
00024 rhoc = def_quantities_real_(40,1)
00025 qc = def_quantities_real_(42,1)
00026 rhomax = def_quantities_real_(43,1)
00027 qmax = def_quantities_real_(46,1)
00028 chicusp = def_quantities_real_(121,1)
00029 rsurf = surf_param_real_(1,1)
00030 rgmid = grid_param_real_(2,1)
00031
00032 open(15,file='physical_quantities.txt',status='unknown')
00033 write(15,'(1p,20e20.12)') omega, radius, m_rest, j_star, m_adm, m_kom, j_adm, &
00034 & m_grav_sph, m_rest_sph, M_R_sph, rhoc, qc, rhomax, qmax, chicusp, rsurf, rgmid
00035
00036 write (15,'(a30,1p,e20.12)') '#Omega*Radius = ', omega
00037 write (15,'(a30,1p,e20.12)') '#Radius = ', radius
00038 write (15,'(a30,1p,e20.12)') '#Rest mass = ', m_rest
00039 write (15,'(a30,1p,e20.12)') '#NS angular momentum = ', j_star
00040 write (15,'(a30,1p,e20.12)') '#ADM mass = ', m_adm
00041 write (15,'(a30,1p,e20.12)') '#Komar mass = ', m_kom
00042 write (15,'(a30,1p,e20.12)') '#Angular momentum = ', j_adm
00043 write (15,'(a30,1p,e20.12)') '#Gravitational spher. mass = ', m_grav_sph
00044 write (15,'(a30,1p,e20.12)') '#Rest spher. mass = ', m_rest_sph
00045 write (15,'(a30,1p,e20.12)') '#Spherical compactness M/R = ', M_R_sph
00046 write (15,'(a30,1p,e20.12)') '#Central restmass density = ', rhoc
00047 write (15,'(a30,1p,e20.12)') '#Central Emden function = ', qc
00048 write (15,'(a30,1p,e20.12)') '#Maximum restmass density = ', rhomax
00049 write (15,'(a30,1p,e20.12)') '#Maximum Emden function = ', qmax
00050 write (15,'(a30,1p,e20.12)') '#chi at cusp = ', chicusp
00051 write (15,'(a30,1p,e20.12)') '#r_surf = ', rsurf
00052 write (15,'(a30,1p,e20.12)') '#rgmid = ', rgmid
00053 close(15)
00054
00055 end subroutine write_last_physq_BNS_mpt
