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
00007 include '../Include_file/include_interface_modulefiles_peos.f90'
00008 include '../Module_interface/interface_error_metric.f90'
00009 include '../Module_interface/interface_grdr_gridpoint_type0.f90'
00010 include '../Module_interface/interface_sourceterm_exsurf_eqm_binary.f90'
00011 include '../Module_interface/interface_sourceterm_poisson_solver_test.f90'
00012 include '../Module_interface/interface_poisson_solver_binary.f90'
00013 include '../Module_interface/interface_poisson_solver_binary_surf_int.f90'
00014 include '../Module_interface/interface_poisson_solver_binary_vol_int.f90'
00015
00016 include '../Include_file/include_subroutines_peos.f90'
00017 include '../Subroutine/IO_output_poisson_test_3D.f90'
00018 include '../Subroutine/iteration_poisson_solver_test.f90'
00019 include '../Subroutine/error_metric.f90'
00020 include '../Subroutine/printout_error_metric.f90'
00021 include '../Subroutine/sourceterm_poisson_solver_test.f90'
00022 include '../Subroutine/sourceterm_exsurf_eqm_binary.f90'
00023 include '../Subroutine/poisson_solver_binary.f90'
00024 include '../Subroutine/poisson_solver_binary_vol_int.f90'
00025 include '../Subroutine/poisson_solver_binary_surf_int.f90'
00026 include '../Subroutine/grdr_gridpoint_type0.f90'
00027 include '../Subroutine/test_source.f90'
00028 include '../Subroutine/allocate_poisson_solver_test.f90'
00029 include '../Subroutine/test_analytic_solution.f90'
00030
00031
00032
00033
00034 PROGRAM Main_poisson_binary_test
00035
00036 use grid_parameter, only : indata_type, outdata_type, &
00037 & iter_max
00038 use grid_parameter_binary_excision
00039 use grid_points_binary_excision
00040 use weight_midpoint_binary_excision
00041 implicit none
00042 integer :: iseq, iter_count, total_iteration
00043
00044 call coordinate_patch_kit_grav
00045 call read_parameter_binary_excision
00046 call calc_parameter_binary_excision
00047 call allocate_grid_points_binary_excision
00048 call calc_grid_points_binary_excision
00049 call allocate_weight_midpoint_binary_excision
00050 call calc_weight_midpoint_binary_excision
00051 call allocate_poisson_solver_test
00052 call test_source
00053 call test_analytic_solution
00054
00055 call iteration_poisson_solver_test(iter_count)
00056 if (total_iteration.ge.iter_max) then
00057 write(6,*)' ** Solution did not converge **'
00058 end if
00059
00060 if (outdata_type.eq.'3D') call IO_output_poisson_test_3D
00061
00062
00063
00064 END PROGRAM Main_poisson_binary_test