00001 subroutine adjust_calc_fncval_Virial_Py_Mratio_mpt(niq,msec_f_oold)
00002 use phys_constant, only : long, nmpt
00003 use def_quantities, only : admmass_asymp, komarmass_asymp, admmom_asymp
00004 use def_quantities_bh, only : AHmass
00005 use def_binary_parameter, only : mass_ratio
00006 implicit none
00007 integer :: niq
00008 real(long) :: msec_f_oold(niq), Mirr1, Mirr2
00009
00010 call copy_def_quantities_from_mpt(nmpt)
00011 msec_f_oold(1) = (admmass_asymp - komarmass_asymp)/admmass_asymp
00012 msec_f_oold(2) = admmom_asymp(2)
00013 call copy_def_quantities_bh_from_mpt(1)
00014 Mirr1 = AHmass
00015 call copy_def_quantities_bh_from_mpt(2)
00016 Mirr2 = AHmass
00017 msec_f_oold(3) = Mirr2/(Mirr1*mass_ratio) - 1.0d0
00018
00019 end subroutine adjust_calc_fncval_Virial_Py_Mratio_mpt
