Cocal: /home/uryu/Cocal_svn_checkinout/Cocal/code/Main_code/Main_poisson_1bh_2pot

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/weight_midpoint_binary_excision.f90'
00006 include '../Module/radial_green_fn_grav_bhex_nb.f90'
00007 include '../Module/radial_green_fn_grav_bhex_di.f90'
00008 include '../Module/radial_green_fn_grav_bhex_dd.f90'
00009 include '../Module/radial_green_fn_grav_bhex_nd.f90'
00010 include '../Module/radial_green_fn_grav_bhex_dh.f90'
00011 include '../Module/radial_green_fn_grav_bhex_nh.f90'
00012 include '../Module/copy_array_3d.f90'
00013 !
00014 include '../Include_file/include_interface_modulefiles_peos.f90'
00015 include '../Module_interface/interface_error_metric_type0.f90'
00016 include '../Module_interface/interface_error_metric_type1.f90'
00017 include '../Module_interface/interface_grdr_gridpoint_type0.f90'
00018 include '../Module_interface/interface_grdr_gridpoint_type0_nosym.f90'
00019 include '../Module_interface/interface_grdr_gridpoint_type0_3rd_nosym.f90'
00020 include '../Module_interface/interface_grgrad_midpoint_r3rd_type0.f90'
00021 include '../Module_interface/interface_grgrad_midpoint_r4th_type0.f90'
00022 include '../Module_interface/interface_sourceterm_poisson_solver_test.f90'
00023 include '../Module_interface/interface_sourceterm_exsurf_eqm_binary.f90'
00024 include '../Module_interface/interface_sourceterm_surface_int.f90'
00025 include '../Module_interface/interface_sourceterm_volume_int_bbh_2pot_test.f90'
00026 include '../Module_interface/interface_sourceterm_surface_int_homosol.f90'
00027 include '../Module_interface/interface_sourceterm_1bh_2pot_test.f90'
00028 include '../Module_interface/interface_poisson_solver_binary_bhex_homosol.f90'
00029 include '../Module_interface/interface_poisson_solver_binary_bhex.f90'
00030 include '../Module_interface/interface_poisson_solver_binary_surf_int.f90'
00031 include '../Module_interface/interface_poisson_solver_binary_vol_int.f90'
00032 include '../Module_interface/interface_poisson_solver_bhex_surf_int.f90'
00033 include '../Module_interface/interface_poisson_solver_bhex_surf_int_all.f90'
00034 include '../Module_interface/interface_poisson_solver_1bh_homosol.f90'
00035 include '../Module_interface/interface_copy_to_hgfn_and_gfnsf.f90'
00036 !include '../Module_interface/interface_bh_boundary_test.f90'
00037 include '../Module_interface/interface_bh_boundary_nh_psi_test.f90'
00038 include '../Module_interface/interface_bh_boundary_nh_alph_test.f90'
00039 include '../Module_interface/interface_bh_boundary_dh_psi_test.f90'
00040 include '../Module_interface/interface_bh_boundary_dh_alph_test.f90'
00041 include '../Module_interface/interface_bh_boundary_1bh_nh_psi_test.f90'
00042 include '../Module_interface/interface_bh_boundary_1bh_nh_alph_test.f90'
00043 include '../Module_interface/interface_interpolation_fillup_binary.f90'
00044 include '../Module_interface/interface_interpo_gr2gr_4th.f90'
00045 !
00046 include '../Include_file/include_subroutines_peos.f90'
00047 include '../Subroutine/IO_output_poisson_test_3D.f90'
00048 include '../Subroutine/iteration_poisson_1bh_2pot_test.f90'
00049 include '../Subroutine/error_metric_type0.f90'
00050 include '../Subroutine/error_metric_type1.f90'
00051 include '../Subroutine/printout_error_metric.f90'
00052 include '../Subroutine/printout_error_all_metric.f90'
00053 include '../Subroutine/coordinate_patch_kit_bhex.f90'
00054 include '../Subroutine/copy_hgfn_nb_to_hgfn.f90'
00055 include '../Subroutine/copy_hgfn_di_to_hgfn.f90'
00056 include '../Subroutine/copy_hgfn_dd_to_hgfn.f90'
00057 include '../Subroutine/copy_hgfn_dh_to_hgfn.f90'
00058 include '../Subroutine/copy_to_hgfn_and_gfnsf.f90'
00059 include '../Subroutine/sourceterm_poisson_solver_test.f90'
00060 include '../Subroutine/sourceterm_exsurf_eqm_binary.f90'
00061 include '../Subroutine/sourceterm_surface_int.f90'
00062 include '../Subroutine/sourceterm_volume_int_bbh_2pot_test.f90'
00063 include '../Subroutine/sourceterm_surface_int_homosol.f90'
00064 include '../Subroutine/sourceterm_1bh_2pot_test.f90'
00065 include '../Subroutine/poisson_solver_binary_bhex.f90'
00066 include '../Subroutine/poisson_solver_binary_bhex_homosol.f90'
00067 include '../Subroutine/poisson_solver_binary_vol_int.f90'
00068 include '../Subroutine/poisson_solver_binary_surf_int.f90'
00069 include '../Subroutine/poisson_solver_bhex_surf_int.f90'
00070 include '../Subroutine/poisson_solver_bhex_surf_int_all.f90'
00071 include '../Subroutine/poisson_solver_1bh_homosol.f90'
00072 include '../Subroutine/grdr_gridpoint_type0.f90'
00073 include '../Subroutine/grdr_gridpoint_type0_nosym.f90'
00074 include '../Subroutine/grdr_gridpoint_type0_3rd_nosym.f90'
00075 include '../Subroutine/grgrad_midpoint_r3rd_type0.f90'
00076 include '../Subroutine/grgrad_midpoint_r4th_type0.f90'
00077 include '../Subroutine/reset_bh_boundary.f90'
00078 !include '../Subroutine/bh_boundary_test.f90'
00079 include '../Subroutine/bh_boundary_nh_psi_test.f90'
00080 include '../Subroutine/bh_boundary_nh_alph_test.f90'
00081 include '../Subroutine/bh_boundary_dh_psi_test.f90'
00082 include '../Subroutine/bh_boundary_dh_alph_test.f90'
00083 include '../Subroutine/bh_boundary_1bh_nh_psi_test.f90'
00084 include '../Subroutine/bh_boundary_1bh_nh_alph_test.f90'
00085 include '../Subroutine/allocate_poisson_bbh_test.f90'
00086 include '../Subroutine/test_analytic_solution_bhex_psialph.f90'
00087 include '../Subroutine/test_analytic_solution_1bh_psialph.f90'
00088 include '../Subroutine/interpolation_fillup_binary.f90'
00089 include '../Subroutine/interpo_gr2gr_4th.f90'
00090 !______________________________________________
00091 !
00092 !              Main Program
00093 !______________________________________________
00094 PROGRAM Main_poisson_1bh_2pot_test
00095 !
00096   use grid_parameter, only : outdata_type, iter_max
00097   use radial_green_fn_grav
00098 !  use radial_green_fn_grav_bhex_di
00099 !  use radial_green_fn_grav_bhex_dd
00100   use radial_green_fn_grav_bhex_nb
00101   use radial_green_fn_grav_bhex_dh
00102   use radial_green_fn_grav_bhex_nh
00103   use radial_green_fn_grav_bhex_dd
00104   use radial_green_fn_grav_bhex_nd
00105   implicit none
00106   integer :: iseq, iter_count, total_iteration
00107 !
00108   call coordinate_patch_kit_grav
00109   call allocate_hgfn_bhex
00110 !  call allocate_hgfn_bhex_di
00111 !  call calc_hgfn_bhex_di
00112 !  call copy_hgfn_di_to_hgfn
00113   call allocate_hgfn_bhex_dd
00114   call calc_hgfn_bhex_dd
00115   call allocate_hgfn_bhex_nd
00116   call calc_hgfn_bhex_nd
00117 !  call copy_hgfn_dd_to_hgfn
00118 !
00119 ! -- No boundary Green's fn
00120   call allocate_hgfn_bhex_nb
00121   call calc_hgfn_bhex_nb
00122 !  call copy_to_hgfn_and_gfnsf(hgfn_nb,gfnsf_nb)
00123 !
00124 ! -- Dirichlet with homogeneous sol
00125   call allocate_hgfn_bhex_dh
00126   call calc_hgfn_bhex_dh
00127 !  call copy_to_hgfn_and_gfnsf(hgfn_dh,gfnsf_dh)
00128 ! -- Neumann with homogeneous sol
00129   call allocate_hgfn_bhex_nh
00130   call calc_hgfn_bhex_nh
00131 !n  call copy_to_hgfn_and_gfnsf(hgfn_nh,gfnsf_nh)
00132 ! --
00133   call allocate_poisson_bbh_test
00134   call test_analytic_solution_1bh_psialph
00135 ! 
00136   call iteration_poisson_1bh_2pot_test(iter_count)
00137   if (total_iteration.ge.iter_max) then
00138     write(6,*)' ** Solution did not converge **'
00139   end if
00140 !  call printout_NS_shape_seq(iseq)
00141   if (outdata_type.eq.'3D') call IO_output_poisson_test_3D
00142 !!  call printout_debug
00143 !  call printout_NS_shape
00144 !
00145 END PROGRAM Main_poisson_1bh_2pot_test
/home/uryu/Cocal_svn_checkinout/Cocal/code/Main_code/Main_poisson_1bh_2pot_test.f90
