/home/uryu/Cocal_svn_checkinout/Cocal/code/Main_utility/interpolation_contour_wave_binary_test_mpt.f90
Go to the documentation of this file.00001 include '../Include_file/include_modulefiles_plot_helm.f90' 00002 include '../Include_file/include_modulefiles_analysis_plot_helm.f90' 00003 include '../Include_file/include_interface_modulefiles_plot_helm.f90' 00004 include '../Include_file/include_interface_modulefiles_analysis_plot_helm.f90' 00005 include '../Include_file/include_subroutines_plot_helm.f90' 00006 include '../Include_file/include_subroutines_analysis_plot_helm.f90' 00007 include '../Include_file/include_functions.f90' 00008 !______________________________________________ 00009 ! 00010 ! Interpolation Program 00011 !______________________________________________ 00012 PROGRAM interpolation_contour_wave_binary_test_mpt 00013 ! 00014 use grid_parameter_binary_excision 00015 use phys_constant 00016 use def_metric, only : psi 00017 use def_metric_cartesian, only : psica 00018 use interface_modules_cartesian 00019 use grid_points_binary_excision 00020 use grid_points_asymptotic_patch 00021 use def_vector_x 00022 implicit none 00023 integer :: impt 00024 ! 00025 ! -- Read parameters 00026 call allocate_grid_parameter_mpt 00027 call allocate_grid_parameter_binary_excision_mpt 00028 call allocate_def_bh_parameter_mpt 00029 do impt = 1, nmpt 00030 call read_parameter_mpt(impt) 00031 call copy_grid_parameter_to_mpt(impt) 00032 call read_parameter_binary_excision_mpt(impt) 00033 call copy_grid_parameter_binary_excision_to_mpt(impt) 00034 call read_parameter_bh_mpt(impt) 00035 call copy_def_bh_parameter_to_mpt(impt) 00036 end do 00037 ! 00038 ! -- Allocate arrays 00039 call set_allocate_size_mpt 00040 ! 00041 call allocate_grid_points_binary_excision 00042 call allocate_grid_points_asymptotic_patch 00043 call allocate_metric_and_matter_BHNS_test 00044 call allocate_grid_points_binary_excision_mpt 00045 call allocate_grid_points_asymptotic_patch_mpt 00046 call allocate_metric_and_matter_BHNS_test_mpt 00047 call coordinate_patch_cartesian 00048 call allocate_metric_and_matter_cartesian 00049 call allocate_mpatch_all_BBH_CF 00050 call allocate_coordinate_patch_kit_grav_mpt 00051 call allocate_vector_x 00052 ! 00053 do impt = 1, nmpt 00054 call copy_grid_parameter_from_mpt(impt) 00055 call copy_grid_parameter_binary_excision_from_mpt(impt) 00056 call copy_def_bh_parameter_from_mpt(impt) 00057 call coordinate_patch_kit_grav_noGreen_mpt 00058 call calc_parameter_binary_excision 00059 call calc_grid_points_binary_excision 00060 call copy_to_mpatch_interpolation_utility(impt) 00061 ! 00062 call IO_input_potential_test_3D_mpt(impt) 00063 call copy_def_metric_to_mpt(impt) 00064 call copy_def_binary_parameter_to_mpt(impt) 00065 end do 00066 ! 00067 call copy_from_mpatch_interpolation_utility(nmpt) 00068 ! -- coordinates of asymptotic patch in central patches 00069 do impt = 1, nmpt 00070 call calc_grid_points_asymptotic_patch(impt,nmpt) 00071 call copy_grid_points_asymptotic_patch_to_mpt(impt) 00072 end do 00073 ! 00074 do impt = 1, nmpt 00075 call copy_from_mpatch_interpolation_utility(impt) 00076 call copy_def_metric_from_mpt(impt) 00077 ! call IO_output_plot_xyz_mpt(impt) 00078 call IO_output_plot_xyz_BBH_CF_mpt(impt) 00079 call calc_potential_minus_rinv_excision 00080 call interpolation_fillup_cartesian(psi,psica) 00081 call IO_output_cartesian_contour_potential_test_mpt(impt) 00082 end do 00083 ! 00084 END PROGRAM interpolation_contour_wave_binary_test_mpt
