00001 subroutine printout_physq_BNS_all_mpt(filename)
00002 use phys_constant, only : nmpt, g, c, solmas
00003 use grid_parameter, only : ntfeq, npfxzp, npfyzp, ntfpolp, sw_mass_iter
00004 use def_matter, only : emd, rs
00005 use def_matter_parameter, only : ome, radi, pinx
00006 use def_quantities
00007 use def_quantities_mpt
00008 use def_matter_parameter_mpt
00009 use grid_parameter_mpt
00010 use def_binary_parameter
00011 use def_binary_parameter_mpt
00012 implicit none
00013 real(8) :: MM = solmas, LL = g*solmas/c**2, TT = g*solmas/c**3
00014 real(8) :: Msph_tot
00015 real(8) :: fixeddlm, fixedvir
00016 integer :: i,ia,ib
00017 character(30) :: filename
00018
00019 Msph_tot = def_quantities_real_(28,1) + def_quantities_real_(28,2)
00020
00021 open(100,file=filename,status='unknown')
00022
00023 write(100,'(a62)') '## If units are not mentioned, G=c=Msol=1 is implied ##'
00024 write(100,'(a62)') '## M is total spherical gravitational mass (G=c=Msol=1) ##'
00025 write(100,'(a62)') '## R0 is the rescaling factor of all coordinate patches ##'
00026
00027 write(100,'(a72)') '========================================================================'
00028
00029 write(100,'(a24,1p,1e23.15)') 'Omega*R0 = ', def_matter_param_real_(3,1)
00030
00031 write(100,'(a24,1p,1e23.15)') 'Omega = ', &
00032 & def_matter_param_real_(3,1)/def_matter_param_real_(5,1)
00033
00034 write(100,'(a24,1p,1e23.15)') 'Omega M_ADM = ', &
00035 & def_matter_param_real_(3,1)/def_matter_param_real_(5,1)*def_quantities_real_(7,nmpt)
00036
00037 write(100,'(a24,1p,1e23.15)') 'Omega [rad/sec] = ', def_quantities_real_(78,1)
00038
00039 write(100,'(a24,1p,1e23.15)') 'R0 = ', def_matter_param_real_(5,1)
00040
00041 write(100,'(a24,1p,2e23.15)') 'r_surf = ', surf_param_real_(1,1), surf_param_real_(1,2)
00042
00043 write(100,'(a24,1p,2e23.15)') 'Radius Rx [km] = ', &
00044 & surf_param_real_(1,1)*def_matter_param_real_(5,1)*LL*1.0d-5, &
00045 & surf_param_real_(1,2)*def_matter_param_real_(5,2)*LL*1.0d-5
00046
00047 write(100,'(a24,1p,2e23.15)') 'Separation = ', sepa_(1)
00048
00049 write(100,'(a24,1p,2e23.15)') 'Separation*R0 = ', sepa_(1)*def_matter_param_real_(5,1)
00050
00051 write(100,'(a24,1p,2e23.15)') 'Separation*R0 [km] = ', sepa_(1)*def_matter_param_real_(5,1)*LL*1.0d-5
00052
00053 write(100,'(a24,1p,2e23.15)') 'Distance to CM (d) = ', dis_(1), dis_(2)
00054
00055 write(100,'(a24,1p,2e23.15)') 'd*R0 = ', dis_(1)*def_matter_param_real_(5,1), &
00056 & dis_(2)*def_matter_param_real_(5,2)
00057
00058 write(100,'(a24,1p,2e23.15)') 'd*R0 [km] = ', dis_(1)*def_matter_param_real_(5,1)*LL*1.0d-5, &
00059 & dis_(2)*def_matter_param_real_(5,2)*LL*1.0d-5
00060
00061 write(100,'(a72)') '========================================================================'
00062 write(100,'(a24,1p,2e23.15,a30)') 'admmass = ', def_quantities_real_(1,1), &
00063 & def_quantities_real_(1,2), ' nc: Not valid for COCP1,2'
00064 write(100,'(a24,1p,2e23.15)') 'komarmass = ', def_quantities_real_(2,1), &
00065 & def_quantities_real_(2,2)
00066 write(100,'(a24,1p,2e23.15,a30)') 'komarmass_nc = ', def_quantities_real_(3,1), &
00067 & def_quantities_real_(3,2), ' nc: Not valid for COCP1,2'
00068 write(100,'(a24,1p,2e23.15)') 'restmass = ', def_quantities_real_(4,1), def_quantities_real_(4,2)
00069 write(100,'(a24,1p,2e23.15)') 'propermass = ', def_quantities_real_(5,1), def_quantities_real_(5,2)
00070 write(100,'(a24,1p,2e23.15)') 'angmom = ', def_quantities_real_(6,1), def_quantities_real_(6,2)
00071
00072 write(100,'(a72)') '========================================================================'
00073 write(100,'(a24,1p,2e23.15)') 'admmass_asymp = ', def_quantities_real_(7,nmpt)
00074 write(100,'(a24,1p,2e23.15)') 'komarmass_asymp = ', def_quantities_real_(8,nmpt)
00075 write(100,'(a24,1p,2e23.15)') 'angmom_asymp = ', def_quantities_real_(9,nmpt)
00076 write(100,'(a24,1p,2e23.15)') 'admmom_asymp(1) = ', def_quantities_real_(10,nmpt)
00077 write(100,'(a24,1p,2e23.15)') 'admmom_asymp(2) = ', def_quantities_real_(11,nmpt)
00078 write(100,'(a24,1p,2e23.15)') 'admmom_asymp(3) = ', def_quantities_real_(12,nmpt)
00079 write(100,'(a24,1p,2e23.15)') 'E/M = (M_ADM-M)/M = ', def_quantities_real_(7,nmpt)/Msph_tot - 1.0d0
00080 write(100,'(a24,1p,2e23.15)') 'J/M^2 = ', def_quantities_real_(9,nmpt)/Msph_tot**2.0d0
00081 write(100,'(a24,1p,2e23.15)') 'J/M_ADM^2 = ', &
00082 & def_quantities_real_(9,nmpt)/def_quantities_real_(7,nmpt)**2.0d0
00083 write(100,'(a24,1p,2e23.15)') '1 - M_K/M_ADM = ', &
00084 & 1.0d0 - def_quantities_real_(8,nmpt)/def_quantities_real_(7,nmpt)
00085
00086 write(100,'(a72)') '========================================================================'
00087 write(100,'(a24,1p,2e23.15)') 'T_kinene = ', def_quantities_real_(13,1), def_quantities_real_(13,2)
00088 write(100,'(a24,1p,2e23.15)') 'W_gravene = ', def_quantities_real_(14,1), def_quantities_real_(14,2)
00089 write(100,'(a24,1p,2e23.15)') 'P_intene = ', def_quantities_real_(15,1), def_quantities_real_(15,2)
00090 write(100,'(a24,1p,2e23.15)') 'M_emfene = ', def_quantities_real_(16,1), def_quantities_real_(16,2)
00091 write(100,'(a24,1p,2e23.15)') 'M_torBene = ', def_quantities_real_(17,1), def_quantities_real_(17,2)
00092 write(100,'(a24,1p,2e23.15)') 'M_polBene = ', def_quantities_real_(18,1), def_quantities_real_(18,2)
00093 write(100,'(a24,1p,2e23.15)') 'M_eleEene = ', def_quantities_real_(19,1), def_quantities_real_(19,2)
00094
00095 write(100,'(a24,1p,2e23.15)') 'Virial = ', def_quantities_real_(20,1), def_quantities_real_(20,2)
00096
00097 write(100,'(a24,1p,2e23.15)') 'ToverW = ', def_quantities_real_(21,1), def_quantities_real_(21,2)
00098 write(100,'(a24,1p,2e23.15)') 'PoverW = ', def_quantities_real_(22,1), def_quantities_real_(22,2)
00099 write(100,'(a24,1p,2e23.15)') 'MoverW = ', def_quantities_real_(23,1), def_quantities_real_(23,2)
00100 write(100,'(a24,1p,2e23.15)') 'MtorBoverW = ', def_quantities_real_(24,1), def_quantities_real_(24,2)
00101 write(100,'(a24,1p,2e23.15)') 'MpolBoverW = ', def_quantities_real_(25,1), def_quantities_real_(25,2)
00102 write(100,'(a24,1p,2e23.15)') 'MeleEoverW = ', def_quantities_real_(26,1), def_quantities_real_(26,2)
00103
00104 write(100,'(a24,1p,2e23.15)') 'I_inertia = ', def_quantities_real_(27,1), def_quantities_real_(27,2)
00105
00106 write(100,'(a72)') '========================================================================'
00107 write(100,'(a24,1p,2e23.15)') 'gravmass_sph = ', def_quantities_real_(28,1), def_quantities_real_(28,2)
00108 write(100,'(a24,1p,2e23.15)') 'restmass_sph = ', def_quantities_real_(29,1), def_quantities_real_(29,2)
00109 write(100,'(a24,1p,2e23.15)') 'propermass_sph = ', def_quantities_real_(30,1), def_quantities_real_(30,2)
00110 write(100,'(a24,1p,2e23.15)') 'MoverR_sph = ', def_quantities_real_(31,1), def_quantities_real_(31,2)
00111 write(100,'(a24,1p,2e23.15)') 'M_0/M_0sph - 1 = ', def_quantities_real_(4,1)/def_quantities_real_(29,1) - 1.0d0, &
00112 & def_quantities_real_(4,2)/def_quantities_real_(29,2) - 1.0d0
00113
00114 write(100,'(a72)') '========================================================================'
00115 write(100,'(a24,1p,2e23.15)') 'schwarz_radi_sph = ', def_quantities_real_(32,1), def_quantities_real_(32,2)
00116 write(100,'(a24,1p,2e23.15)') 'coord_radius_x = ', def_quantities_real_(33,1), def_quantities_real_(33,2)
00117 write(100,'(a24,1p,2e23.15)') 'coord_radius_y = ', def_quantities_real_(34,1), def_quantities_real_(34,2)
00118 write(100,'(a24,1p,2e23.15)') 'coord_radius_z = ', def_quantities_real_(35,1), def_quantities_real_(35,2)
00119 write(100,'(a24,1p,2e23.15)') 'Coord axis ratio y/x = ', def_quantities_real_(34,1)/def_quantities_real_(33,1), &
00120 & def_quantities_real_(34,2)/def_quantities_real_(33,2)
00121 write(100,'(a24,1p,2e23.15)') 'Coord axis ratio z/x = ', def_quantities_real_(35,1)/def_quantities_real_(33,1), &
00122 & def_quantities_real_(35,2)/def_quantities_real_(33,2)
00123 write(100,'(a24,1p,2e23.15)') 'sqrt(1-(Ry/Rx)^2) = ', &
00124 & sqrt(1.0d0-(def_quantities_real_(34,1)/def_quantities_real_(33,1))**2), &
00125 & sqrt(1.0d0-(def_quantities_real_(34,2)/def_quantities_real_(33,2))**2)
00126 write(100,'(a24,1p,2e23.15)') 'sqrt(1-(Rz/Rx)^2) = ', &
00127 & sqrt(1.0d0-(def_quantities_real_(35,1)/def_quantities_real_(33,1))**2), &
00128 & sqrt(1.0d0-(def_quantities_real_(35,2)/def_quantities_real_(33,2))**2)
00129
00130 write(100,'(a24,1p,2e23.15)') 'proper_radius_x = ', def_quantities_real_(36,1), def_quantities_real_(36,2)
00131 write(100,'(a24,1p,2e23.15)') 'proper_radius_y = ', def_quantities_real_(37,1), def_quantities_real_(37,2)
00132 write(100,'(a24,1p,2e23.15)') 'proper_radius_z = ', def_quantities_real_(38,1), def_quantities_real_(38,2)
00133 write(100,'(a24,1p,2e23.15)') 'Proper axis ratio y/x = ', def_quantities_real_(37,1)/def_quantities_real_(36,1), &
00134 & def_quantities_real_(37,2)/def_quantities_real_(36,2)
00135 write(100,'(a24,1p,2e23.15)') 'Proper axis ratio z/x = ', def_quantities_real_(38,1)/def_quantities_real_(36,1), &
00136 & def_quantities_real_(38,2)/def_quantities_real_(36,2)
00137 write(100,'(a24,1p,2e23.15)') 'sqrt(1-(pRy/pRx)^2) = ', &
00138 & sqrt(1.0d0-(def_quantities_real_(37,1)/def_quantities_real_(36,1))**2), &
00139 & sqrt(1.0d0-(def_quantities_real_(37,2)/def_quantities_real_(36,2))**2)
00140 write(100,'(a24,1p,2e23.15)') 'sqrt(1-(pRz/pRx)^2) = ', &
00141 & sqrt(1.0d0-(def_quantities_real_(38,1)/def_quantities_real_(36,1))**2), &
00142 & sqrt(1.0d0-(def_quantities_real_(38,2)/def_quantities_real_(36,2))**2)
00143
00144 write(100,'(a72)') '========================================================================'
00145 write(100,'(a24,1p,2e23.15)') 'schwarz_radi_sph_km = ', def_quantities_real_(39,1), def_quantities_real_(39,2)
00146 write(100,'(a24,1p,2e23.15)') 'coord_radius_x_km = ', def_quantities_real_(40,1), def_quantities_real_(40,2)
00147 write(100,'(a24,1p,2e23.15)') 'coord_radius_y_km = ', def_quantities_real_(41,1), def_quantities_real_(41,2)
00148 write(100,'(a24,1p,2e23.15)') 'coord_radius_z_km = ', def_quantities_real_(42,1), def_quantities_real_(42,2)
00149 write(100,'(a24,1p,2e23.15)') 'proper_radius_x_km = ', def_quantities_real_(43,1), def_quantities_real_(43,2)
00150 write(100,'(a24,1p,2e23.15)') 'proper_radius_y_km = ', def_quantities_real_(44,1), def_quantities_real_(44,2)
00151 write(100,'(a24,1p,2e23.15)') 'proper_radius_z_km = ', def_quantities_real_(45,1), def_quantities_real_(45,2)
00152
00153 write(100,'(a72)') '========================================================================'
00154 write(100,'(a24,1p,2e23.15)') 'rho_c = ', def_quantities_real_(46,1), def_quantities_real_(46,2)
00155 write(100,'(a24,1p,2e23.15)') 'pre_c = ', def_quantities_real_(47,1), def_quantities_real_(47,2)
00156 write(100,'(a24,1p,2e23.15)') 'epsi_c = ', def_quantities_real_(48,1), def_quantities_real_(48,2)
00157 write(100,'(a24,1p,2e23.15)') 'q_c = ', def_quantities_real_(49,1), def_quantities_real_(49,2)
00158 write(100,'(a24,1p,2e23.15)') 'rho_max = ', def_quantities_real_(50,1), def_quantities_real_(50,2)
00159 write(100,'(a24,1p,2e23.15)') 'pre_max = ', def_quantities_real_(51,1), def_quantities_real_(51,2)
00160 write(100,'(a24,1p,2e23.15)') 'epsi_max = ', def_quantities_real_(52,1), def_quantities_real_(52,2)
00161 write(100,'(a24,1p,2e23.15)') 'q_max = ', def_quantities_real_(53,1), def_quantities_real_(53,2)
00162
00163 write(100,'(a72)') '========================================================================'
00164 write(100,'(a24,1p,2e23.15)') 'rho_c_cgs = ', def_quantities_real_(54,1), def_quantities_real_(54,2)
00165 write(100,'(a24,1p,2e23.15)') 'pre_c_cgs = ', def_quantities_real_(55,1), def_quantities_real_(55,2)
00166 write(100,'(a24,1p,2e23.15)') 'epsi_c_cgs = ', def_quantities_real_(56,1), def_quantities_real_(56,2)
00167 write(100,'(a24,1p,2e23.15)') 'q_c_cgs = ', def_quantities_real_(57,1), def_quantities_real_(57,2)
00168 write(100,'(a24,1p,2e23.15)') 'rho_max_cgs = ', def_quantities_real_(58,1), def_quantities_real_(58,2)
00169 write(100,'(a24,1p,2e23.15)') 'pre_max_cgs = ', def_quantities_real_(59,1), def_quantities_real_(59,2)
00170 write(100,'(a24,1p,2e23.15)') 'epsi_max_cgs = ', def_quantities_real_(60,1), def_quantities_real_(60,2)
00171 write(100,'(a24,1p,2e23.15)') 'q_max_cgs = ', def_quantities_real_(61,1), def_quantities_real_(61,2)
00172
00173 write(100,'(a72)') '========================================================================'
00174 write(100,'(a24,1p,2e23.15)') 'zrb_xp_plus = ', def_quantities_real_(62,1), def_quantities_real_(62,2)
00175 write(100,'(a24,1p,2e23.15)') 'zrb_xp_minus = ', def_quantities_real_(63,1), def_quantities_real_(63,2)
00176 write(100,'(a24,1p,2e23.15)') 'zrb_yp_plus = ', def_quantities_real_(64,1), def_quantities_real_(64,2)
00177 write(100,'(a24,1p,2e23.15)') 'zrb_yp_minus = ', def_quantities_real_(65,1), def_quantities_real_(65,2)
00178 write(100,'(a24,1p,2e23.15)') 'zrb_zp_plus = ', def_quantities_real_(66,1), def_quantities_real_(66,2)
00179 write(100,'(a24,1p,2e23.15)') 'zrb_zp_minus = ', def_quantities_real_(67,1), def_quantities_real_(67,2)
00180
00181 write(100,'(a72)') '========================================================================'
00182 write(100,'(a24,1p,2e23.15)') 'dhdr_x = ', def_quantities_real_(68,1), def_quantities_real_(68,2)
00183 write(100,'(a24,1p,2e23.15)') 'dhdr_y = ', def_quantities_real_(69,1), def_quantities_real_(69,2)
00184 write(100,'(a24,1p,2e23.15)') 'dhdr_z = ', def_quantities_real_(70,1), def_quantities_real_(70,2)
00185 write(100,'(a24,1p,2e23.15)') 'chi_cusp = ', def_quantities_real_(71,1), def_quantities_real_(71,2)
00186
00187 write(100,'(a72)') '========================================================================'
00188 write(100,'(a24,1p,2e23.15)') 'circ_line_xy = ', def_quantities_real_(72,1), def_quantities_real_(72,2)
00189 write(100,'(a24,1p,2e23.15)') 'circ_line_yz = ', def_quantities_real_(73,1), def_quantities_real_(73,2)
00190 write(100,'(a24,1p,2e23.15)') 'circ_line_zx = ', def_quantities_real_(74,1), def_quantities_real_(74,2)
00191 write(100,'(a24,1p,2e23.15)') 'circ_surf_xy = ', def_quantities_real_(75,1), def_quantities_real_(75,2)
00192 write(100,'(a24,1p,2e23.15)') 'circ_surf_yz = ', def_quantities_real_(76,1), def_quantities_real_(76,2)
00193 write(100,'(a24,1p,2e23.15)') 'circ_surf_zx = ', def_quantities_real_(77,1), def_quantities_real_(77,2)
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237 close(100)
00238
00239 end subroutine printout_physq_BNS_all_mpt
