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

Go to the documentation of this file.

subroutine copy_msec_BNS_iterqt_to_mpt(x_vector,sw_master)
  use phys_constant, only : long
  use grid_parameter_mpt
  use def_quantities_mpt
  use def_matter_parameter_mpt
  use def_binary_parameter_mpt 
  implicit none
  real(long), pointer  :: x_vector(:)
  character(3), intent(in) :: sw_master
!
  select case (sw_master)
    case ("010", "030", "050")
      def_quantities_real_(49,1)  = x_vector(1)  !emdc1
      def_matter_param_real_(2,1) = x_vector(1)  !emdc1
!     Second star has same q
      def_quantities_real_(49,2)  = x_vector(1)  !emdc1
      def_matter_param_real_(2,2) = x_vector(1)  !emdc1
    case ("020", "040", "060")
      write(6,*) "***** Inside copy_msec_BNS_iterqt_to_mpt: dis1, dis2, emdc1, emdc2, rsf2:", &
      & dis_(1), dis_(2), def_quantities_real_(49,1), def_quantities_real_(49,2), surf_param_real_(1,2)

      dis_(1) = x_vector(1)                      !dis1
      def_quantities_real_(49,1)  = x_vector(2)  !emdc1
      def_matter_param_real_(2,1) = x_vector(2)  !emdc1

      dis_(2) = sepa_(2) - x_vector(1)           !dis2
      surf_param_real_(1,2)       = x_vector(3)  !rsf2 
      def_quantities_real_(49,2)  = x_vector(4)  !emdc2
      def_matter_param_real_(2,2) = x_vector(4)  !emdc2
    case ("110", "130", "150", "210", "230", "250")
      write(6,*) "***** Inside copy_msec_BNS_iterqt_to_mpt: rsurf, emdc:", &
      & surf_param_real_(1,1), def_quantities_real_(49,1)

      surf_param_real_(1,1)       = x_vector(1)  !rsf1 
      def_quantities_real_(49,1)  = x_vector(2)  !emdc1       
      def_matter_param_real_(2,1) = x_vector(2)  !emdc1
!     Second star has same rs, q
      surf_param_real_(1,2)       = x_vector(1)  !rsf1 
      def_quantities_real_(49,2)  = x_vector(2)  !emdc1       
      def_matter_param_real_(2,2) = x_vector(2)  !emdc1
    case ("120", "140", "160", "220", "240", "260")
      surf_param_real_(1,1)       = x_vector(1)  !rsf1 
      dis_(1)                     = x_vector(2)  !dis1 
      def_quantities_real_(49,1)  = x_vector(3)  !emdc1  
      def_matter_param_real_(2,1) = x_vector(3)  !emdc1
      surf_param_real_(1,2)       = x_vector(4)  !rsf2
      def_quantities_real_(49,2)  = x_vector(5)  !emdc2  
      def_matter_param_real_(2,2) = x_vector(5)  !emdc2
  end select


!  do ii = 1, 2
!    call copy_def_matter_parameter_from_mpt(ii)
!    call copy_def_binary_parameter_from_mpt(ii)
!    if (ii.eq.1) then
!      emdc = x_vector(1)
!      dis  = x_vector(3)
!    end if
!    if (ii.eq.2) then
!      emdc = x_vector(2)
!      dis  = sepa - x_vector(3)
!    end if
!    call copy_def_matter_parameter_to_mpt(ii)
!    call copy_def_binary_parameter_to_mpt(ii)
!  end do

!
end subroutine copy_msec_BNS_iterqt_to_mpt

---