/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/printout_physq_plot.f90
Go to the documentation of this file.00001 subroutine printout_physq_plot(iseq,flag_restmass) 00002 use grid_parameter, only : ntfeq, npfxzp, npfyzp, ntfpolp, sw_mass_iter 00003 use def_matter, only : emd, rs 00004 use def_matter_parameter, only : ome, radi, pinx 00005 use def_quantities 00006 implicit none 00007 real(8) :: fixeddlm, fixedvir 00008 integer, intent(in) :: iseq, flag_restmass 00009 ! 00010 if (iseq.eq.1) then 00011 open(200,file='rnsphyplot_all.dat',status='unknown') 00012 ! open(210,file='rnsphyplot.dat',status='unknown') 00013 else 00014 open(200,file='rnsphyplot_all.dat',status='old', position="append") 00015 end if 00016 ! 00017 ! Reference for the data 00018 ! Slot # 00019 ! 1 '== Sequence number == ', iseq 00020 ! 2 '#### Solution did not converge #### ', flag_restmass 00021 ! '## Coordinate and Proper Radii in K = 1 unit ##' 00022 ! 3, 4 ' NS radius along x = ', coord_radius_x, proper_radius_x 00023 ! 5, 6 ' NS radius along y = ', coord_radius_y, proper_radius_y 00024 ! 7, 8 ' NS radius along z = ', coord_radius_z, proper_radius_z 00025 ! 9 ' Axis ratio y/x = ', coord_radius_y/coord_radius_x 00026 ! 10 proper_radius_y/proper_radius_x 00027 ! 11 ' Axis ratio z/x = ', coord_radius_z/coord_radius_x 00028 ! 12 proper_radius_z/proper_radius_x 00029 ! '## Radii in Req = 1 unit ##' 00030 ! 13 ' NS radius along x = ', rs(ntfeq,npfxzp) 00031 ! 14 ' NS radius along y = ', rs(ntfeq,npfyzp) 00032 ! 15 ' NS radius along z = ', rs(ntfpolp,npfxzp) 00033 ! '## Eccentricity in coordinate and Proper Radii ##' 00034 ! 16 ' sqrt(1-(Rz/Rx)^2) = ', sqrt(1.0d0-(coord_radius_z/coord_radius_x)**2) 00035 ! 17 sqrt(1.0d0-(proper_radius_z/proper_radius_x)**2) 00036 ! 00037 ! '## M = spherical gravitational mass (K = 1 unit) ##' 00038 ! 18 19 ' Omega M and Omega = ', ome/radi*gravmass_sph, ome/radi 00039 ! 20 ' M_ADM = ', admmass 00040 ! 21 ' M_K Komar mass = ', komarmass 00041 ! 22 ' Rest mass M_0 1star= ', restmass 00042 ! 23 ' Proper mass M_p = ', propermass 00043 ! 24 ' Angular momentum J = ', angmom 00044 ! 25 ' Spherical rest mass= ', restmass_sph 00045 ! 26 ' Spherical grav mass= ', gravmass_sph 00046 ! 27 ' Spherical M/R = ', MoverR_sph 00047 ! 28 ' Schwarzschildradius= ', schwarz_radi_sph 00048 ! 00049 ! 29 ' E/M = (M_ADM-M)/M = ', admmass/gravmass_sph - 1.0d0 00050 ! 30 ' J/M^2 and J/M_ADM^2= ', angmom/gravmass_sph**2.0d0 00051 ! 31 angmom/admmass**2.0d0 00052 ! 00053 ! 32 ' T kinetic energy = ', T_kinene 00054 ! 33 ' W grav energy = ', W_gravene 00055 ! 34 ' T/|W| = ', ToverW 00056 ! 35 ' Moment of inertia = ', I_inertia 00057 ! 00058 ! '## Virial and Rest mass accuracy ##' 00059 ! 36 ' 1 - M_K/M_ADM = ', fixedvir = (admmass - komarmass)/admmass 00060 ! 37 ' M_0/M_0sph - 1 = ', fixeddlm = (restmass - restmass_sph) 00061 ! /restmass_sphfixeddlm 00062 ! '## Maximums of the density, pressure and p/rho ##' 00063 ! 38 ' rho_max = ', rho_max 00064 ! 39 ' pre_max = ', pre_max 00065 ! 40 ' epsilon_max = ', epsi_max 00066 ! 41 ' (p/rho)_max = ', q_max 00067 ! 00068 ! '## Red and blue shift ##' 00069 ! 42 43 ' surface on x axis = ', zrb_xp_plus, zrb_xp_minus 00070 ! 44, 45 ' surface on y axis = ', zrb_yp_plus, zrb_yp_minus 00071 ! 46, 47 ' surface on z axis = ', zrb_zp_plus, zrb_zp_minus 00072 ! 00073 ! 48 ' (Omega M)**2 ' (ome/radi*gravmass_sph)**2 00074 ! 00075 ! 49 ' Omega[rad/s] = ', ome_cgs 00076 ! 50 ' Quasi local spin S = ', qua_loc_spin 00077 ! 51, 52 ' S/M^2 and S/M_ADM^2= ', qua_loc_spin/gravmass_sph**2.0d0, qua_loc_spin/admmass**2.0d0 00078 00079 fixedvir = (admmass - komarmass)/admmass 00080 fixeddlm = (restmass - restmass_sph)/restmass_sph 00081 ! 00082 write(200,'(2i5,1p,52e23.15)') & 00083 iseq, flag_restmass, & 00084 coord_radius_x, proper_radius_x, & 00085 coord_radius_y, proper_radius_y, & 00086 coord_radius_z, proper_radius_z, & 00087 coord_radius_y/coord_radius_x, proper_radius_y/proper_radius_x, & 00088 coord_radius_z/coord_radius_x, proper_radius_z/proper_radius_x, & 00089 rs(ntfeq,npfxzp), rs(ntfeq,npfyzp), rs(ntfpolp,npfxzp), & 00090 sqrt(1.0d0-(coord_radius_z/coord_radius_x)**2), & 00091 sqrt(1.0d0-(proper_radius_z/proper_radius_x)**2), & 00092 ! 00093 ome/radi*gravmass_sph, ome/radi, & 00094 admmass, komarmass, restmass, propermass, angmom, restmass_sph, & 00095 gravmass_sph, MoverR_sph, schwarz_radi_sph, & 00096 ! 00097 admmass/gravmass_sph - 1.0d0, & 00098 angmom/gravmass_sph**2.0d0, angmom/admmass**2.0d0, & 00099 T_kinene, W_gravene, ToverW, I_inertia, fixedvir, fixeddlm, & 00100 rho_max, pre_max, epsi_max, q_max, & 00101 zrb_xp_plus,zrb_xp_minus,zrb_yp_plus,zrb_yp_minus,zrb_zp_plus,zrb_zp_minus, & 00102 (ome/radi*gravmass_sph)**2, & 00103 ome_cgs, & 00104 qua_loc_spin, qua_loc_spin/gravmass_sph**2.0d0, qua_loc_spin/admmass**2.0d0 00105 00106 ! 00107 close(200) 00108 ! 00109 end subroutine printout_physq_plot
