00001
00002 include '../Include_file/include_modulefiles_peos.f90'
00003 include '../Module/grid_parameter_binary_excision.f90'
00004 include '../Module/grid_points_binary_excision.f90'
00005 include '../Module/grid_points_asymptotic_patch.f90'
00006 include '../Module/make_char1_array_2d.f90'
00007 include '../Module/make_char2_array_2d.f90'
00008 include '../Module/make_int_array_3d.f90'
00009 include '../Module/weight_midpoint_binary_excision.f90'
00010 include '../Module_mpatch/grid_points_asymptotic_patch_mpt.f90'
00011 include '../Module_interface/interface_grgrad_midpoint_r3rd_type0.f90'
00012
00013 include '../Include_file/include_interface_modulefiles_peos.f90'
00014 include '../Include_file/include_subroutines_peos.f90'
00015 include '../Include_file/include_modulefiles_mpatch.f90'
00016
00017 include '../Analysis/Module/grid_parameter_cartesian.f90'
00018 include '../Analysis/Module/coordinate_grav_xyz.f90'
00019 include '../Analysis/Module/def_metric_cartesian.f90'
00020 include '../Analysis/Module/def_matter_cartesian.f90'
00021 include '../Analysis/Module/interface_modules_cartesian.f90'
00022
00023 include '../Analysis/Subroutine/coordinate_patch_kit_grav_noGreen.f90'
00024 include '../Analysis/Subroutine/coordinate_patch_kit_grav_noGreen_mpt.f90'
00025 include '../Analysis/Subroutine/coordinate_patch_cartesian.f90'
00026 include '../Analysis/Subroutine/allocate_metric_and_matter_cartesian.f90'
00027 include '../Analysis/Subroutine/IO_output_cartesian_contour_potential_test.f90'
00028 include '../Analysis/Subroutine/IO_output_cartesian_contour_potential_test_mpt.f90'
00029 include '../Analysis/Subroutine/interpo_gr2cgr_4th.f90'
00030 include '../Analysis/Subroutine/interpo_fl2cgr_4th.f90'
00031 include '../Analysis/Subroutine/interpolation_fillup_cartesian.f90'
00032 include '../Analysis/Subroutine/interpolation_fillup_cartesian_mpt.f90'
00033 include '../Analysis/Subroutine/interpolation_matter.f90'
00034 include '../Analysis/Subroutine/IO_output_plot_xyz.f90'
00035 include '../Analysis/Subroutine/IO_output_plot_xyz_mpt.f90'
00036 include '../Analysis/Subroutine/IO_output_plot_averaged_error_mpt.f90'
00037 include '../Analysis/Subroutine/copy_to_mpatch_interpolation_utility.f90'
00038 include '../Analysis/Subroutine/copy_from_mpatch_interpolation_utility.f90'
00039
00040 include '../Subroutine/allocate_metric_and_matter_BHNS_test.f90'
00041 include '../Subroutine/allocate_poisson_solver_test.f90'
00042 include '../Subroutine/interpo_lag4th_2Dsurf.f90'
00043 include '../Subroutine/IO_input_potential_test_3D.f90'
00044 include '../Subroutine/grgrad_midpoint_r3rd_type0.f90'
00045 include '../Subroutine_mpatch/IO_input_potential_test_3D_mpt.f90'
00046 include '../Subroutine_mpatch/test_analytic_BHNS_solution_mpt.f90'
00047 include '../Subroutine_mpatch/test_analytic_bns_solution_mpt.f90'
00048 include '../Subroutine_mpatch/test_analytic_solution_bhex_mpt.f90'
00049 include '../Subroutine_mpatch/test2_analytic_solution_bhex_mpt.f90'
00050 include '../Subroutine_mpatch/allocate_grid_points_asymptotic_patch_mpt.f90'
00051 include '../Subroutine_mpatch/allocate_metric_and_matter_BHNS_test_mpt.f90'
00052 include '../Subroutine_mpatch/copy_grid_points_asymptotic_patch_to_mpt.f90'
00053 include '../Subroutine_mpatch/copy_grid_points_asymptotic_patch_from_mpt.f90'
00054 include '../Include_file/include_subroutines_mpatch.f90'
00055
00056
00057
00058
00059 PROGRAM interpolation_contour_potential_binary_test2_mpt
00060
00061 use grid_parameter_binary_excision
00062 use grid_points_asymptotic_patch
00063 use phys_constant
00064 use def_metric, only : psi
00065 use def_metric_cartesian, only : psica
00066 use interface_modules_cartesian
00067 use grid_points_binary_excision
00068 use trigonometry_grav_phi
00069 implicit none
00070 integer :: impt, impt_ex
00071
00072
00073 call allocate_grid_parameter_mpt
00074 call allocate_grid_parameter_binary_excision_mpt
00075 do impt = 1, nmpt
00076 call read_parameter_mpt(impt)
00077 call copy_grid_parameter_to_mpt(impt)
00078 call read_parameter_binary_excision_mpt(impt)
00079 call copy_grid_parameter_binary_excision_to_mpt(impt)
00080 end do
00081
00082
00083 call set_allocate_size_mpt
00084
00085 call allocate_grid_points_binary_excision
00086 call allocate_grid_points_asymptotic_patch
00087 call allocate_metric_and_matter_BHNS_test
00088 call allocate_grid_points_binary_excision_mpt
00089 call allocate_grid_points_asymptotic_patch_mpt
00090 call allocate_metric_and_matter_BHNS_test_mpt
00091 call allocate_trig_grav_mphi
00092 call coordinate_patch_cartesian
00093 call allocate_metric_and_matter_cartesian
00094 call allocate_mpatch_all_test
00095 call allocate_coordinate_patch_kit_grav_mpt
00096
00097 do impt = 1, nmpt
00098 call copy_grid_parameter_from_mpt(impt)
00099 call copy_grid_parameter_binary_excision_from_mpt(impt)
00100 call coordinate_patch_kit_grav_noGreen_mpt
00101 call calc_parameter_binary_excision
00102 call calc_grid_points_binary_excision
00103 call copy_to_mpatch_interpolation_utility(impt)
00104 call IO_input_potential_test_3D_mpt(impt)
00105 call copy_poisson_solver_test_to_mpt(impt)
00106 end do
00107 if (nmpt.gt.2) then
00108 call copy_from_mpatch_interpolation_utility(nmpt)
00109 do impt = 1, nmpt
00110
00111 call calc_grid_points_asymptotic_patch(impt,nmpt)
00112 call copy_grid_points_asymptotic_patch_to_mpt(impt)
00113 end do
00114 end if
00115
00116 do impt = 1, nmpt
00117 call copy_from_mpatch_interpolation_utility(impt)
00118 call copy_poisson_solver_test_from_mpt(impt)
00119
00120 call test2_analytic_solution_bhex_mpt(impt)
00121 call IO_output_plot_xyz_mpt(impt)
00122 call IO_output_plot_averaged_error_mpt(impt)
00123 if(impt.eq.1) impt_ex = 2
00124 if(impt.eq.2) impt_ex = 1
00125 call interpolation_fillup_cartesian_mpt(psi, psica, impt, impt_ex)
00126 call IO_output_cartesian_contour_potential_test_mpt(impt)
00127 end do
00128
00129 END PROGRAM interpolation_contour_potential_binary_test2_mpt
