/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/IO_input_physq_plot.f90
Go to the documentation of this file.00001 subroutine IO_input_physq_plot(iseq,flag_restmass) 00002 use def_quantities 00003 use def_quantities_derived 00004 implicit none 00005 integer, intent(out) :: iseq, flag_restmass 00006 ! 00007 if (iseq.eq.1) then 00008 open(210,file='rnsphyplot_all_input.dat',status='old') 00009 end if 00010 ! 00011 ! Reference for the data 00012 ! Slot # 00013 ! 1 '== Sequence number == ', iseq 00014 ! 2 '#### Solution did not converge #### ', flag_restmass 00015 ! '## Coordinate and Proper Radii in K = 1 unit ##' 00016 ! 3, 4 ' NS radius along x = ', coord_radius_x, proper_radius_x 00017 ! 5, 6 ' NS radius along y = ', coord_radius_y, proper_radius_y 00018 ! 7, 8 ' NS radius along z = ', coord_radius_z, proper_radius_z 00019 ! 9 10 ' Axis ratio y/x = ', coord_axis_ratio_yx, proper_axis_ratio_yx 00020 ! 11 12 ' Axis ratio z/x = ', coord_axis_ratio_zx, proper_axis_ratio_zx 00021 ! '## Radii in Req = 1 unit ##' 00022 ! 13 ' NS radius along x = ', coord_radius_normalized_x 00023 ! 14 ' NS radius along y = ', coord_radius_normalized_y 00024 ! 15 ' NS radius along z = ', coord_radius_normalized_z 00025 ! '## Eccentricity in coordinate and Proper Radii ##' 00026 ! 16 17 ' sqrt(1-(Rz/Rx)^2) = ', coord_eccentricity_zx,proper_eccentricity_zx 00027 ! 00028 ! '## M = spherical gravitational mass (K = 1 unit) ##' 00029 ! 18 19 ' Omega M and Omega = ', omega_M, omega 00030 ! 20 ' M_ADM = ', admmass 00031 ! 21 ' M_K Komar mass = ', komarmass 00032 ! 22 ' Rest mass M_0 1star= ', restmass 00033 ! 23 ' Proper mass M_p = ', propermass 00034 ! 24 ' Angular momentum J = ', angmom 00035 ! 25 ' Spherical rest mass= ', restmass_sph 00036 ! 26 ' Spherical grav mass= ', gravmass_sph 00037 ! 27 ' Spherical M/R = ', MoverR_sph 00038 ! 28 ' Schwarzschildradius= ', schwarz_radi_sph 00039 ! 00040 ! 29 ' E/M = (M_ADM-M)/M = ', E_bind_over_M 00041 ! 30 31 ' J/M^2 and J/M_ADM^2= ', J_over_M2, J_over_Madm2 00042 ! 00043 ! 32 ' T kinetic energy = ', T_kinene 00044 ! 33 ' W grav energy = ', W_gravene 00045 ! 34 ' T/|W| = ', ToverW 00046 ! 35 ' Moment of inertia = ', I_inertia 00047 ! 00048 ! '## Virial and Rest mass accuracy ##' 00049 ! 36 ' 1 - M_K/M_ADM = ', Madm_MK_virial 00050 ! 37 ' M_0/M_0sph - 1 = ', M0_minus_M0sph_normalized 00051 ! 00052 ! '## Maximums of the density, pressure and p/rho ##' 00053 ! 38 ' rho_max = ', rho_c_normalized 00054 ! 39 ' pre_max = ', pre_c_normalized 00055 ! 40 ' epsilon_max = ', epsilon_c_normalized 00056 ! 41 ' (p/rho)_max = ', p_over_rho_c_normalized 00057 ! 00058 ! '## Red and blue shift ##' 00059 ! 42 43 ' surface on x axis = ', zrb_xp_plus, zrb_xp_minus 00060 ! 44, 45 ' surface on y axis = ', zrb_yp_plus, zrb_yp_minus 00061 ! 46, 47 ' surface on z axis = ', zrb_zp_plus, zrb_zp_minus 00062 ! 00063 ! 48 ' (Omega M)**2 ' omega2_M2 00064 ! 00065 ! 00066 read(210,'(2i5,1p,48e23.15)',ERR=9999,END=9999) & 00067 iseq, flag_restmass, & 00068 coord_radius_x, proper_radius_x, & 00069 coord_radius_y, proper_radius_y, & 00070 coord_radius_z, proper_radius_z, & 00071 coord_axis_ratio_yx, proper_axis_ratio_yx, & 00072 coord_axis_ratio_zx, proper_axis_ratio_zx, & 00073 coord_radius_normalized_x, & 00074 coord_radius_normalized_y, & 00075 coord_radius_normalized_z, & 00076 coord_eccentricity_zx, proper_eccentricity_zx, & 00077 ! 00078 omega_M, omega, & 00079 admmass, komarmass, restmass, propermass, angmom, restmass_sph, & 00080 gravmass_sph, MoverR_sph, schwarz_radi_sph, & 00081 ! 00082 E_bind_over_M, J_over_M2, J_over_Madm2, & 00083 T_kinene, W_gravene, ToverW, I_inertia, & 00084 Madm_MK_virial, M0_minus_M0sph_normalized, & 00085 rho_c_normalized, pre_c_normalized, & 00086 epsilon_c_normalized, p_over_rho_c_normalized, & 00087 zrb_xp_plus,zrb_xp_minus,zrb_yp_plus,zrb_yp_minus,zrb_zp_plus,zrb_zp_minus, & 00088 omega2_M2 00089 ! 00090 return 00091 9999 continue 00092 flag_restmass = 9999 00093 return 00094 ! 00095 end subroutine IO_input_physq_plot
