! File: submodule_diffstar_base_io.f90 ! Authors: Francesco Torsello (FT) !************************************************************************ ! Copyright (C) 2020-2023 Francesco Torsello * ! * ! This file is part of SPHINCS_ID * ! * ! SPHINCS_ID is free software: you can redistribute it and/or modify * ! it under the terms of the GNU General Public License as published by * ! the Free Software Foundation, either version 3 of the License, or * ! (at your option) any later version. * ! * ! SPHINCS_ID is distributed in the hope that it will be useful, * ! but WITHOUT ANY WARRANTY; without even the implied warranty of * ! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * ! GNU General Public License for more details. * ! * ! You should have received a copy of the GNU General Public License * ! along with SPHINCS_ID. If not, see <https://www.gnu.org/licenses/>. * ! The copy of the GNU General Public License should be in the file * ! 'COPYING'. * !************************************************************************ SUBMODULE (diffstar_base) io !*************************************************** ! !# This submodule contains the implementation of the ! methods of TYPE diffstarbase that handle I/O (input/output) ! ! FT 5.11.2021 ! !*************************************************** IMPLICIT NONE CONTAINS !------------------------------! !-- OVERRIDING SUBROUTINES --! !------------------------------! MODULE PROCEDURE print_summary_drs !************************************************ ! !# Prints a summary of the physical properties the system ! to the standard output and, optionally, to a formatted ! file whose name is given as the optional argument `filename` ! ! FT 5.11.2021 ! !************************************************ USE utility, ONLY: density_si2cu USE constants, ONLY: kg2g, m2cm IMPLICIT NONE PRINT *, " * Differentially rotating star (DRS):" PRINT * PRINT *, " Baryon mass of the DRS=", this% mass, "Msun" PRINT *, " Gravitational mass of the DRS=", this% mass_grav, "Msun" PRINT * PRINT *, " Equatorial (not areal) radius of the DRS at phi=0 = ", & this% r_eq, "Msun_geo" PRINT *, " Equatorial (not areal) radius of the DRS at phi=pi/2 = ", & this% r_eq_pi2, "Msun_geo" PRINT *, " Equatorial (not areal) radius of the DRS at phi=pi = ", & this% r_eq_pi, "Msun_geo" PRINT *, " Equatorial (not areal) radius of the DRS at phi=3pi/2 = ", & this% r_eq_3pi2, "Msun_geo" PRINT *, " Polar radius of the DRS= ", this% r_pole, "Msun_geo" PRINT *, " Ratio between polar radius at equatiorial radius at phi=0= ", & this% r_ratio, "Msun_geo" PRINT * PRINT *, " EOS for the DRS= ", this% eos PRINT * PRINT *, " Central baryon mass density for the DRS= ", & this% rho_center, "Msun/Msun_geo**3= ", & this% rho_center/density_si2cu*kg2g/(m2cm**3), "g cm^{-3}" PRINT * PRINT *, " Ratio T/|W| between the rotaional kinetic energy and ", & "the gravitational binding energy: ", this% tsw PRINT *, " For axisymmetric configurations as this one, the ", & "threshold for dynamical bar-mode instability is T/|W|~0.25 ", & " [Masaru Shibata et al 2000 ApJ 542 453, ", & "https://arxiv.org/pdf/astro-ph/0005378.pdf]. See also ", & "[Manca et al., Classical and Quantum Gravity, 24, 171, ", & "https://arxiv.org/abs/0705.1826], ", & "Sec.3.3 in [Galeazzi et al., Astron Astrophys 541:A156, ", & "arXiv:1101.2664], and Sec.5.1.3 in ", & "[Paschalidis, V., Stergioulas, N., Rotating stars in ", & "relativity. Living Rev Relativ 20, 7 (2017), ", & "https://link.springer.com/article/10.1007%2Fs41114-017-0008-x]." PRINT * END PROCEDURE print_summary_drs END SUBMODULE io