/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine_mpatch/printout_physq_BNS_all_mpt.f90

Go to the documentation of this file.

subroutine printout_physq_BNS_all_mpt(filename)
  use phys_constant,  only : nmpt, g, c, solmas
  use grid_parameter, only : ntfeq, npfxzp, npfyzp, ntfpolp, sw_mass_iter
  use def_matter, only : emd, rs
  use def_matter_parameter, only : ome, radi, pinx
  use def_quantities
  use def_quantities_mpt
  use def_matter_parameter_mpt
  use grid_parameter_mpt
  use def_binary_parameter
  use def_binary_parameter_mpt
  implicit none
  real(8) :: MM = solmas, LL = g*solmas/c**2, TT = g*solmas/c**3
  real(8) :: Msph_tot
  real(8) :: fixeddlm, fixedvir
  integer :: i,ia,ib
  character(30) :: filename
!
  Msph_tot = def_quantities_real_(28,1) + def_quantities_real_(28,2)

  open(100,file=filename,status='unknown')

  write(100,'(a62)') '## If units are not mentioned, G=c=Msol=1 is implied        ##'
  write(100,'(a62)') '## M is total spherical gravitational mass (G=c=Msol=1)     ##'
  write(100,'(a62)') '## R0 is the rescaling factor of all coordinate patches     ##'

  write(100,'(a72)') '========================================================================'

  write(100,'(a24,1p,1e23.15)') 'Omega R0              = ', def_matter_param_real_(3,1)
  write(100,'(a24,1p,1e23.15)') 'R0                    = ', def_matter_param_real_(5,1)
  write(100,'(a24,1p,1e23.15)') 'Omega                 = ',  &
  &  def_matter_param_real_(3,1)/def_matter_param_real_(5,1)

  write(100,'(a24,1p,1e23.15)') 'Omega [rad/sec]       = ', def_quantities_real_(122,1)

  write(100,'(a24,1p,1e23.15)') 'Omega M               = ',  &
  &  def_matter_param_real_(3,1)/def_matter_param_real_(5,1)*Msph_tot 

  write(100,'(a24,1p,2e23.15)') 'Distance d_s =sepa*R0 = ',  &
  &  sepa_(1)*def_matter_param_real_(5,1), sepa_(2)*def_matter_param_real_(5,1)
!  &  2.0d0*grid_param_real_(2,1)*def_matter_param_real_(5,1)

  write(100,'(a24,1p,2e23.15)') 'Distance d_s [km]     = ',  &
  &  sepa_(1)*def_matter_param_real_(5,1)*LL*1.0d-5, sepa_(2)*def_matter_param_real_(5,1)*LL*1.0d-5

  write(100,'(a24,1p,2e23.15)') 'Distance d =dis*R0    = ',  &
  &  dis_(1)*def_matter_param_real_(5,1), dis_(2)*def_matter_param_real_(5,1)
!  &  2.0d0*grid_param_real_(2,1)*def_matter_param_real_(5,1)

  write(100,'(a24,1p,2e23.15)') 'Distance d [km]       = ',  &
  &  dis_(1)*def_matter_param_real_(5,1)*LL*1.0d-5, dis_(2)*def_matter_param_real_(5,1)*LL*1.0d-5

  write(100,'(a24,1p,2e23.15)') 'Radius [km]           = ',  &
  &  surf_param_real_(1,1)*def_matter_param_real_(5,1)*LL*1.0d-5,  &
  &  surf_param_real_(1,2)*def_matter_param_real_(5,1)*LL*1.0d-5

  write(100,'(a72)') '========================================================================'

  i=1     !1
  write(100,'(a24,1p,2e23.15,a30)') 'admmass               = ', & 
  &  def_quantities_real_(i,1), def_quantities_real_(i,2), '     nc: Not valid for COCP1,2'

  i=i+1   !2
  write(100,'(a24,1p,2e23.15)') 'komarmass             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !3
  write(100,'(a24,1p,2e23.15,a30)') 'komarmass_nc          = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2), '     nc: Not valid for COCP1,2'

  i=i+1   !4
  write(100,'(a24,1p,2e23.15)') 'restmass              = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !5
  write(100,'(a24,1p,2e23.15)') 'propermass            = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !6
  write(100,'(a24,1p,2e23.15)') 'angmom                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a24,1p,2e23.15)') 'M_0/M_0sph - 1        = ', &
  & (def_quantities_real_(4,1) - def_quantities_real_(29,1))/def_quantities_real_(29,1),  &
  & (def_quantities_real_(4,2) - def_quantities_real_(29,2))/def_quantities_real_(29,2)

  write(100,'(a72)') '========================================================================'

  i=i+1   !7
  write(100,'(a24,1p,2e23.15)') 'admmass_asymp         = ', def_quantities_real_(i,nmpt)

  i=i+1   !8
  write(100,'(a24,1p,2e23.15)') 'komarmass_asymp       = ', def_quantities_real_(i,nmpt)

  i=i+1   !9
  write(100,'(a24,1p,2e23.15)') 'angmom_asymp          = ', def_quantities_real_(i,nmpt)

  i=i+1   !10
  write(100,'(a24,1p,2e23.15)') 'admmom_asymp(1)       = ', def_quantities_real_(i,nmpt)

  i=i+1   !11
  write(100,'(a24,1p,2e23.15)') 'admmom_asymp(2)       = ', def_quantities_real_(i,nmpt)

  i=i+1   !12
  write(100,'(a24,1p,2e23.15)') 'admmom_asymp(3)       = ', def_quantities_real_(i,nmpt)

  write(100,'(a24,1p,2e23.15)') 'E/M = (M_ADM-M)/M     = ', &
  &  def_quantities_real_(7,nmpt)/Msph_tot - 1.0d0

  write(100,'(a24,1p,2e23.15)') 'J/M^2                 = ', &
  &  def_quantities_real_(9,nmpt)/Msph_tot**2.0d0

  write(100,'(a24,1p,2e23.15)') 'J/M_ADM^2             = ', &
  &  def_quantities_real_(9,nmpt)/def_quantities_real_(7,nmpt)**2.0d0

  write(100,'(a24,1p,2e23.15)') '1 - M_K/M_ADM         = ', &
  &  (def_quantities_real_(7,nmpt) - def_quantities_real_(8,nmpt))/def_quantities_real_(7,nmpt) 

  write(100,'(a72)') '========================================================================'

  i=i+1   !13
  write(100,'(a24,1p,2e23.15)') 'T_kinene              = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !14
  write(100,'(a24,1p,2e23.15)') 'W_gravene             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !15
  write(100,'(a24,1p,2e23.15)') 'P_intene              = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !16
  write(100,'(a24,1p,2e23.15)') 'M_emfene              = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !17
  write(100,'(a24,1p,2e23.15)') 'M_torBene             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !18
  write(100,'(a24,1p,2e23.15)') 'M_polBene             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !19
  write(100,'(a24,1p,2e23.15)') 'M_eleEene             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !20
  write(100,'(a24,1p,2e23.15)') 'Virial                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !21
  write(100,'(a24,1p,2e23.15)') 'ToverW                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !22
  write(100,'(a24,1p,2e23.15)') 'PoverW                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !23
  write(100,'(a24,1p,2e23.15)') 'MoverW                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !24
  write(100,'(a24,1p,2e23.15)') 'MtorBoverW            = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !25
  write(100,'(a24,1p,2e23.15)') 'MpolBoverW            = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !26
  write(100,'(a24,1p,2e23.15)') 'MeleEoverW            = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !27
  write(100,'(a24,1p,2e23.15)') 'I_inertia             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a72)') '========================================================================'

  i=i+1   !28
  write(100,'(a24,1p,2e23.15)') 'gravmass_sph          = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !29
  write(100,'(a24,1p,2e23.15)') 'restmass_sph          = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !30
  write(100,'(a24,1p,2e23.15)') 'propermass_sph        = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !31
  write(100,'(a24,1p,2e23.15)') 'MoverR_sph            = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !32
  write(100,'(a24,1p,2e23.15)') 'schwarz_radi_sph      = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !33
  write(100,'(a24,1p,2e23.15)') 'schwarz_radi_sph_km   = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a72)') '========================================================================'

  i=i+1   !34
  write(100,'(a24,1p,2e23.15)') 'coord_radius_x        = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !35
  write(100,'(a24,1p,2e23.15)') 'coord_radius_y        = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !36
  write(100,'(a24,1p,2e23.15)') 'coord_radius_z        = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a24,1p,2e23.15)') 'Coord axis ratio y/x  = ', &
  &  def_quantities_real_(35,1)/def_quantities_real_(34,1),  &
  &  def_quantities_real_(35,2)/def_quantities_real_(34,2)      

  write(100,'(a24,1p,2e23.15)') 'Coord axis ratio z/x  = ', &
  &  def_quantities_real_(36,1)/def_quantities_real_(34,1),  &
  &  def_quantities_real_(36,2)/def_quantities_real_(34,2)

  write(100,'(a24,1p,2e23.15)') 'sqrt(1-(Rz/Rx)^2)     = ', &
  &           sqrt(1.0d0-(def_quantities_real_(36,1)/def_quantities_real_(34,1))**2), &
  &           sqrt(1.0d0-(def_quantities_real_(36,2)/def_quantities_real_(34,2))**2)

  i=i+1   !37
  write(100,'(a24,1p,2e23.15)') 'proper_radius_x       = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !38
  write(100,'(a24,1p,2e23.15)') 'proper_radius_y       = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !39
  write(100,'(a24,1p,2e23.15)') 'proper_radius_z       = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a24,1p,2e23.15)') 'Proper axis ratio y/x = ', &
  &  def_quantities_real_(38,1)/def_quantities_real_(37,1),  &
  &  def_quantities_real_(38,2)/def_quantities_real_(37,2)      

  write(100,'(a24,1p,2e23.15)') 'Proper axis ratio z/x = ', &
  &  def_quantities_real_(39,1)/def_quantities_real_(37,1),  &
  &  def_quantities_real_(39,2)/def_quantities_real_(37,2)

  write(100,'(a24,1p,2e23.15)') 'sqrt(1-(Rz/Rx)^2)     = ', &
  &           sqrt(1.0d0-(def_quantities_real_(39,1)/def_quantities_real_(37,1))**2), &
  &           sqrt(1.0d0-(def_quantities_real_(39,2)/def_quantities_real_(37,2))**2)

  write(100,'(a72)') '========================================================================'

  i=i+1   !40
  write(100,'(a24,1p,2e23.15)') 'rho_c                 = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !41
  write(100,'(a24,1p,2e23.15)') 'pre_c                 = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !42
  write(100,'(a24,1p,2e23.15)') 'q_c                   = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !43
  write(100,'(a24,1p,2e23.15)') 'rho_max               = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !44
  write(100,'(a24,1p,2e23.15)') 'pre_max               = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !45
  write(100,'(a24,1p,2e23.15)') 'epsi_max              = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !46
  write(100,'(a24,1p,2e23.15)') 'q_max                 = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a72)') '========================================================================'

  i=i+1   !47
  write(100,'(a24,1p,2e23.15)') 'coord_radius_x_km     = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !48
  write(100,'(a24,1p,2e23.15)') 'coord_radius_y_km     = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !49
  write(100,'(a24,1p,2e23.15)') 'coord_radius_z_km     = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !50
  write(100,'(a24,1p,2e23.15)') 'proper_radius_x_km    = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !51
  write(100,'(a24,1p,2e23.15)') 'proper_radius_y_km    = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !52
  write(100,'(a24,1p,2e23.15)') 'proper_radius_z_km    = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a72)') '========================================================================'

  i=i+1   !53
  write(100,'(a24,1p,2e23.15)') 'rho_c_cgs             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !54
  write(100,'(a24,1p,2e23.15)') 'pre_c_cgs             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !55
  write(100,'(a24,1p,2e23.15)') 'epsi_c_cgs            = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !56
  write(100,'(a24,1p,2e23.15)') 'q_c_cgs               = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !57
  write(100,'(a24,1p,2e23.15)') 'rho_max_cgs           = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !58
  write(100,'(a24,1p,2e23.15)') 'pre_max_cgs           = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !59
  write(100,'(a24,1p,2e23.15)') 'epsi_max_cgs          = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !60
  write(100,'(a24,1p,2e23.15)') 'q_max_cgs             = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a72)') '========================================================================'

  i=i+1   !61
  write(100,'(a24,1p,2e23.15)') 'zrb_xp_plus           = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !62
  write(100,'(a24,1p,2e23.15)') 'zrb_xp_minus          = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !63
  write(100,'(a24,1p,2e23.15)') 'zrb_yp_plus           = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !64
  write(100,'(a24,1p,2e23.15)') 'zrb_yp_minus          = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !65
  write(100,'(a24,1p,2e23.15)') 'zrb_zp_plus           = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !66
  write(100,'(a24,1p,2e23.15)') 'zrb_zp_minus          = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !67
  write(100,'(a24,1p,2e23.15)') 'dhdr_x                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !68
  write(100,'(a24,1p,2e23.15)') 'dhdr_y                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !69
  write(100,'(a24,1p,2e23.15)') 'dhdr_z                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  write(100,'(a72)') '========================================================================'

  do ib  = 1, 3
    do ia  = 1, 3
      i=i+1   !70,75,80,85,90,95,100,105,110
      write(100,'(a4,i1,a1,i1,a17,1p,2e23.15)') 'Iij(',ia,',',ib,')             = ', &
      &  def_quantities_real_(i,1), def_quantities_real_(i,2)
    
      i=i+1   !71,76,81,86,91,96,101,106,111
      write(100,'(a4,i1,a1,i1,a17,1p,2e23.15)') 'Itf(',ia,',',ib,')             = ', &
      &  def_quantities_real_(i,1), def_quantities_real_(i,2)

      i=i+1   !72,77,82,87,92,97,102,107,112
      write(100,'(a7,i1,a1,i1,a14,1p,2e23.15)') 'dt1Itf(',ia,',',ib,')          = ', &
      &  def_quantities_real_(i,1), def_quantities_real_(i,2)

      i=i+1   !73,78,83,88,93,98,103,108,113
      write(100,'(a7,i1,a1,i1,a14,1p,2e23.15)') 'dt2Itf(',ia,',',ib,')          = ', &
      &  def_quantities_real_(i,1), def_quantities_real_(i,2)

      i=i+1   !74,79,84,89,94,99,104,109,114
      write(100,'(a7,i1,a1,i1,a14,1p,2e23.15)') 'dt3Itf(',ia,',',ib,')          = ', &
      &  def_quantities_real_(i,1), def_quantities_real_(i,2)
    end do
  end do

  i=i+1   !115
  write(100,'(a24,1p,2e23.15)')   'LGW                   = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  do ia  = 1, 3
    i=i+1   !116,117,118
    write(100,'(a5,i1,a18,1p,2e23.15)') 'dJdt(',ia,')               = ', &
    &   def_quantities_real_(i,1), def_quantities_real_(i,2)
  end do

  i=i+1   !119
  write(100,'(a24,1p,2e23.15)')   'hplus                 = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !120
  write(100,'(a24,1p,2e23.15)')   'hcross                = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  i=i+1   !121
  write(100,'(a24,1p,2e23.15)')   'chi_cusp              = ', &
  &  def_quantities_real_(i,1), def_quantities_real_(i,2)

  close(100)
!
end subroutine printout_physq_BNS_all_mpt

---