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

Go to the documentation of this file.

subroutine update_coordinates_BNS_mpt
  use phys_constant, only : nmpt
  use grid_parameter
  use grid_parameter_mpt
  use def_binary_parameter
  use def_quantities_mpt       
  use def_matter_parameter_mpt
  use grid_parameter_binary_excision
  use grid_points_binary_excision
  use grid_points_asymptotic_patch
  use grid_points_binary_in_asympto
  use weight_midpoint_binary_excision
  implicit none
  integer :: impt
!
  write(6,*) "----------------------- Start update coordinates --------------------------"
  do impt = 1, nmpt
    call copy_from_mpatch_all_BNS_CF(impt)
    if (impt==1 .or. impt==2)  then
      write(6,*)  "==> start impt, r_surf, qc, emdc", impt, surf_param_real_(1,impt), &
        &          def_quantities_real_(49,impt), def_matter_param_real_(2,impt)
    end if
!    call copy_grid_parameter_from_mpt(impt)          ! r_surf is here
!    call copy_def_binary_parameter_from_mpt(impt)
!    call copy_grid_parameter_binary_excision_from_mpt(impt)
!    call copy_def_peos_parameter_from_mpt(impt)
    call coordinate_patch_kit_grav_grid_mpt(3)       ! 3:r_surf is used
!    call calc_parameter_binary_excision
!    call IO_printout_grid_data_mpt(impt)
    if (impt.ne.nmpt) then
      call calc_grid_points_binary_excision
      call calc_grid_points_binary_in_asympto(impt,nmpt)
      call copy_grid_points_binary_in_asympto_to_mpt(impt)
    end if
    call calc_weight_midpoint_binary_excision
    if (impt==1 .or. impt==2) then
      call calc_vector_x_grav(2)
      call calc_vector_phi_grav(2)
    else
      call calc_vector_x_grav(1)
      call calc_vector_phi_grav(1)
    end if
    call copy_to_mpatch_all_BNS_CF(impt)
    if (impt==1 .or. impt==2)  then
      write(6,*)  "==> end impt, r_surf, qc, emdc", impt, surf_param_real_(1,impt), &
        &        def_quantities_real_(49,impt), def_matter_param_real_(2,impt)
    end if
  end do
  call copy_from_mpatch_all_BNS_CF(nmpt)
! -- coordinates of asymptotic patch in central patches
  do impt = 1, 2
    call calc_grid_points_asymptotic_patch(impt,nmpt)
    call copy_grid_points_asymptotic_patch_to_mpt(impt)
  end do
  write(6,*) "------------------------ End update coordinates --------------------------"
!
end subroutine update_coordinates_BNS_mpt

---