Cocal: /home/uryu/Cocal_svn_checkinout/Cocal/code/Main_utility/interpolation_contour_potential_BBH_CF

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