/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine_mpatch/adjust_copy_trpPunc_to_mpt.f90
Go to the documentation of this file.00001 subroutine adjust_copy_trpPunc_to_mpt(niq,msec_x) 00002 use phys_constant, only : long, nmpt 00003 use def_bh_parameter, only : mom_pBH, mass_pBH 00004 use def_binary_parameter, only : dis, sepa 00005 implicit none 00006 integer :: niq, ii 00007 real(long) :: msec_x(niq) 00008 ! 00009 do ii = 1, 2 00010 call copy_def_bh_parameter_from_mpt(ii) 00011 call copy_def_binary_parameter_from_mpt(ii) 00012 if (ii.eq.1) then 00013 mom_pBH(2) = msec_x(1) 00014 dis = msec_x(2) 00015 end if 00016 if (ii.eq.2) then 00017 mom_pBH(2) = msec_x(1) 00018 dis = sepa - msec_x(2) 00019 mass_pBH = msec_x(3) 00020 end if 00021 call copy_def_bh_parameter_to_mpt(ii) 00022 call copy_def_binary_parameter_to_mpt(ii) 00023 end do 00024 ! 00025 end subroutine adjust_copy_trpPunc_to_mpt
