00001 subroutine printout_iterqt_BNS_mpt(x_vector, f_vector, im,id,is)
00002 use phys_constant, only : nmpt, g, c, solmas, long
00003 use def_quantities_mpt
00004 use def_matter_parameter_mpt
00005 use def_binary_parameter_mpt
00006 implicit none
00007 real(long), pointer :: x_vector(:), f_vector(:)
00008 real(8) :: MM = solmas, LL = g*solmas/c**2, TT = g*solmas/c**3
00009 integer :: im,id,is
00010
00011 if(id.eq.1.and.im.eq.1) then
00012 call length_conversion(2, surf_param_real_(2,1), bd)
00013 bsepa = sepa_(1)
00014 write(6,'(a41,1p,e23.15)') "Binary separation in COCAL units :", bsepa
00015 write(6,'(a41,1p,e23.15)') "Target binary separation in km :", surf_param_real_(2,1)
00016 write(6,'(a41,1p,e23.15)') "Target binary separation in COCAL units :", bd
00017 end if
00018
00019
00020 end subroutine printout_iterqt_BNS_mpt