00001 subroutine update_coordinates_mpt 00002 use phys_constant, only : nmpt 00003 use grid_parameter, only : rgin 00004 use def_bh_parameter, only : ome_bh 00005 use def_binary_parameter, only : dis 00006 use grid_parameter_binary_excision 00007 use grid_points_binary_excision 00008 use grid_points_asymptotic_patch 00009 use grid_points_binary_in_asympto 00010 use weight_midpoint_binary_excision 00011 implicit none 00012 integer :: impt 00013 ! 00014 do impt = 1, nmpt 00015 call copy_grid_parameter_from_mpt(impt) 00016 call copy_def_binary_parameter_from_mpt(impt) 00017 call copy_grid_parameter_binary_excision_from_mpt(impt) 00018 call copy_def_bh_parameter_from_mpt(impt) 00019 call coordinate_patch_kit_grav_mpt 00020 !!! call calc_parameter_binary_excision 00021 if (impt.ne.nmpt) then 00022 call calc_grid_points_binary_excision 00023 call calc_grid_points_binary_in_asympto(impt,nmpt) 00024 call copy_grid_points_binary_in_asympto_to_mpt(impt) 00025 end if 00026 call calc_weight_midpoint_binary_excision 00027 call calc_vector_x_grav(0) 00028 call calc_vector_phi_grav(0) 00029 call calc_vector_bh(0) 00030 ! 00031 call copy_to_mpatch_all_BBH_CF(impt) 00032 end do 00033 call copy_from_mpatch_all_BBH_CF(nmpt) 00034 ! -- coordinates of asymptotic patch in central patches 00035 do impt = 1, 2 00036 call calc_grid_points_asymptotic_patch(impt,nmpt) 00037 call copy_grid_points_asymptotic_patch_to_mpt(impt) 00038 end do 00039 ! 00040 end subroutine update_coordinates_mpt