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

Go to the documentation of this file.

subroutine store_msec_BNS_fiterqt_mpt(f_vector,sw_master)
  use phys_constant, only : long, nmpt
  use def_quantities
  use grid_parameter_mpt
  use def_quantities_mpt
  use def_binary_parameter_mpt
  use def_matter_parameter_mpt
  implicit none
  real(long), pointer :: f_vector(:)
  character(3), intent(in) :: sw_master
  real(long) :: omega, ome, R0
  real(long) :: Mo1, Mo2, Madm, MoR1, MoR2, Rx, Sz, Jasym, Pasym
  real(long) :: tsepa, binsepa, tqt1, tqt2, tqt, mr
!
  Madm  = def_quantities_real_( 7,nmpt)   ! admmass_asymp  
  Jasym = def_quantities_real_( 9,nmpt)   ! angmom_asymp 
  Pasym = def_quantities_real_(11,nmpt)   ! admmom_asymp(2) 

  Mo1   = def_quantities_real_( 4,1)  
  Mo2   = def_quantities_real_( 4,2)    

  R0      = def_matter_param_real_(5,1)    !  radi
  ome     = def_matter_param_real_(3,1)    !  ome
  omega   = ome/R0                         !  omega
  binsepa = sepa_(1)*R0

  write(6,*) 'KKKKK sepa, R0, ome=', sepa_(1), R0, ome

  tsepa = surf_param_real_(2,1)     !target_sepa
  tqt1  = surf_param_real_(3,1)     !target_qt
  tqt2  = surf_param_real_(3,2)     !target_qt
  tqt   = tqt1 + tqt2
  mr    = tqt2/tqt1

  select case (sw_master)
    case("010")   ! niq=1   TESTED
      f_vector(1) = (Mo1 - tqt1)/tqt1
      write(6,*) "***** Inside store_msec_BNS_fiterqt_mpt, case 010:", Mo1, tqt1
    case("030")   ! niq=1   TESTED
      f_vector(1) = (Madm - tqt)/tqt
      write(6,*) "***** Inside store_msec_BNS_fiterqt_mpt, case 030:", Madm, tqt
    case("050")   ! niq=1   NOT DONE YET
!      f_vector(1) = (MoR1 - tqt1)/tqt1  
    case("020")   ! niq=4
      f_vector(1) = Pasym
      f_vector(2) = (Mo1 - tqt1)/tqt1
      f_vector(3) = mr*Mo1/Mo2 -1.0d0
      f_vector(4) = (Mo2 - tqt2)/tqt2
      write(6,*) "***** Inside store_msec_BNS_fiterqt_mpt, case 020:", Mo1, tqt1, &
             &   Mo2, tqt2, Pasym, mr
    case("040")   ! niq=4
      f_vector(1) = Pasym
      f_vector(2) = (Mo1 - tqt1)/tqt1    ! NOT CORRECT
      f_vector(3) = (Madm - tqt)/tqt
      f_vector(4) = (Mo2 - tqt2)/tqt2
    case("060")   ! niq=3
!      f_vector(1) = Pasym
!      f_vector(2) = (MoR1 - tqt1)/tqt1
!      f_vector(3) = (MoR2 - tqt2)/tqt2
    case("110")   ! niq=2   TESTED 
      f_vector(1) = (binsepa - tsepa)/tsepa
      f_vector(2) = (Mo1 - tqt1)/tqt1
      write(6,*) "***** Inside store_msec_BNS_fiterqt_mpt, case 110:", Mo1, tqt1, &
             &   binsepa, tsepa
    case("130")   ! niq=2   TESTED
      f_vector(1) = (binsepa - tsepa)/tsepa
      f_vector(2) = (Madm - tqt)/tqt
      write(6,*) "***** Inside store_msec_BNS_fiterqt_mpt, case 130:", Madm, tqt, &
             &   binsepa, tsepa
    case("150")   ! niq=2
!      f_vector(1) = (binsepa - tsepa_coc)/tsepa_coc
!      f_vector(2) = (MoR1 - tqt1)/tqt1
    case("210")   ! niq=2   NOT TESTED
      f_vector(1) = (omega - tsepa)/tsepa
      f_vector(2) = (Mo1 - tqt1)/tqt1
      write(6,*) "***** Inside store_msec_BNS_fiterqt_mpt, case 210:", Mo1, tqt1, &
             &   omega, tsepa
    case("230")   ! niq=2   TESTED
      f_vector(1) = (omega - tsepa)/tsepa
      f_vector(2) = (Madm - tqt)/tqt
      write(6,*) "***** Inside store_msec_BNS_fiterqt_mpt, case 230:", Madm, tqt, &
             &   omega, tsepa
    case("250")   ! niq=2  NOT DONE YET
 !     f_vector(1) = (omega - tsepa_coc)/tsepa_coc
 !     f_vector(2) = (MoR1 - tqt1)/tqt1
  end select
!
end subroutine store_msec_BNS_fiterqt_mpt

---