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/def_bh_parameter.f90'
00006 include '../Module/def_metric_excurve_grid.f90'
00007 include '../Module/def_vector_bh.f90'
00008 include '../Module/def_vector_irg.f90'
00009 include '../Module/weight_midpoint_binary_excision.f90'
00010 include '../Module/radial_green_fn_grav_bhex_nb.f90'
00011 include '../Module/radial_green_fn_grav_bhex_dd.f90'
00012 include '../Module/radial_green_fn_grav_bhex_nd.f90'
00013 include '../Module/radial_green_fn_grav_bhex_nh.f90'
00014 include '../Module/radial_green_fn_grav_bhex_dh.f90'
00015 include '../Module/copy_array_3d.f90'
00016 include '../Module/gnufor2.f90'
00017
00018 include '../Include_file/include_interface_modulefiles_peos.f90'
00019 include '../Module_interface/interface_error_metric_type0.f90'
00020 include '../Module_interface/interface_error_metric_type1.f90'
00021 include '../Module_interface/interface_error_metric_type2.f90'
00022 include '../Module_interface/interface_grdr_gridpoint_type0.f90'
00023 include '../Module_interface/interface_grdr_gridpoint_type0_nosym.f90'
00024 include '../Module_interface/interface_grdr_gridpoint_type0_3rd_nosym.f90'
00025 include '../Module_interface/interface_grgrad_4th_gridpoint_bhex.f90'
00026 include '../Module_interface/interface_grgrad_midpoint_r3rd_type0.f90'
00027 include '../Module_interface/interface_grgrad_midpoint_r4th_type0.f90'
00028 include '../Module_interface/interface_grgrad_midpoint_type0.f90'
00029 include '../Module_interface/interface_sourceterm_poisson_solver_test.f90'
00030 include '../Module_interface/interface_sourceterm_surface_int_homosol.f90'
00031 include '../Module_interface/interface_poisson_solver_binary_bhex_homosol.f90'
00032 include '../Module_interface/interface_sourceterm_exsurf_eqm_binary.f90'
00033 include '../Module_interface/interface_sourceterm_surface_int.f90'
00034 include '../Module_interface/interface_poisson_solver_binary_bhex.f90'
00035 include '../Module_interface/interface_poisson_solver_binary_surf_int.f90'
00036 include '../Module_interface/interface_poisson_solver_binary_vol_int.f90'
00037 include '../Module_interface/interface_poisson_solver_bhex_surf_int.f90'
00038 include '../Module_interface/interface_poisson_solver_bhex_surf_int_all.f90'
00039 include '../Module_interface/interface_copy_to_hgfn_and_gfnsf.f90'
00040 include '../Module_interface/interface_interpo_gr2gr_4th.f90'
00041 include '../Module_interface/interface_update_parameter.f90'
00042 include '../Module_interface/interface_interpolation_fillup_binary.f90'
00043 include '../Module_interface/interface_surf_source_adm_mass.f90'
00044 include '../Module_interface/interface_surf_source_komar_mass.f90'
00045 include '../Module_interface/interface_source_ang_mom_inf.f90'
00046 include '../Module_interface/interface_source_ang_mom_thr.f90'
00047 include '../Module_interface/interface_source_ang_mom_exc.f90'
00048 include '../Module_interface/interface_source_ang_mom_smarr.f90'
00049 include '../Module_interface/interface_surf_int_grav_rg.f90'
00050 include '../Module_interface/interface_vol_int_grav_bhex.f90'
00051
00052
00053 include '../Module_interface/interface_bh_boundary_d_psi.f90'
00054 include '../Module_interface/interface_bh_boundary_nh_psi_test.f90'
00055 include '../Module_interface/interface_bh_boundary_d_alps.f90'
00056 include '../Module_interface/interface_bh_boundary_d_bvxd.f90'
00057 include '../Module_interface/interface_bh_boundary_d_bvyd.f90'
00058 include '../Module_interface/interface_bh_boundary_d_bvzd.f90'
00059 include '../Module_interface/interface_bh_boundary_d_Bfun.f90'
00060 include '../Module_interface/interface_bh_boundary_n_Bfun.f90'
00061 include '../Module_interface/interface_bh_boundary_d_potx.f90'
00062 include '../Module_interface/interface_bh_boundary_d_poty.f90'
00063 include '../Module_interface/interface_bh_boundary_d_potz.f90'
00064 include '../Module_interface/interface_bh_boundary_AH.f90'
00065 include '../Module_interface/interface_bh_boundary_CF.f90'
00066 include '../Module_interface/interface_outer_boundary_d_psi.f90'
00067 include '../Module_interface/interface_outer_boundary_d_alps.f90'
00068 include '../Module_interface/interface_outer_boundary_d_bvxd.f90'
00069 include '../Module_interface/interface_outer_boundary_d_bvyd.f90'
00070 include '../Module_interface/interface_outer_boundary_d_bvzd.f90'
00071 include '../Module_interface/interface_outer_boundary_d_Bfun.f90'
00072 include '../Module_interface/interface_outer_boundary_n_Bfun.f90'
00073 include '../Module_interface/interface_outer_boundary_d_potx.f90'
00074 include '../Module_interface/interface_outer_boundary_d_poty.f90'
00075 include '../Module_interface/interface_outer_boundary_d_potz.f90'
00076 include '../Module_interface/interface_sourceterm_HaC_CF.f90'
00077 include '../Module_interface/interface_sourceterm_trG_CF.f90'
00078 include '../Module_interface/interface_sourceterm_MoC_CF.f90'
00079 include '../Module_interface/interface_sourceterm_MoC_CF_type1_bhex.f90'
00080 include '../Module_interface/interface_sourceterm_Bfun.f90'
00081 include '../Module_interface/interface_sourceterm_exsurf_eqm_binary_parity.f90'
00082 include '../Module_interface/interface_interpolation_fillup_binary_parity.f90'
00083 include '../Module_interface/interface_compute_shift_v2.f90'
00084 include '../Module_interface/interface_compute_dBfun.f90'
00085
00086 include '../Include_file/include_subroutines_peos.f90'
00087 include '../Subroutine/read_parameter_bh.f90'
00088 include '../Subroutine/excurve_CF_gridpoint.f90'
00089 include '../Subroutine/excurve_CF_gridpoint_bhex.f90'
00090 include '../Subroutine/calc_vector_bh.f90'
00091 include '../Subroutine/calc_vector_irg.f90'
00092
00093 include '../Subroutine/IO_output_poisson_test_3D.f90'
00094 include '../Subroutine/error_metric_type0.f90'
00095 include '../Subroutine/error_metric_type1.f90'
00096 include '../Subroutine/error_metric_type2.f90'
00097 include '../Subroutine/printout_error_metric.f90'
00098 include '../Subroutine/printout_error_all_metric.f90'
00099 include '../Subroutine/coordinate_patch_kit_bhex.f90'
00100 include '../Subroutine/copy_hgfn_nb_to_hgfn.f90'
00101 include '../Subroutine/copy_hgfn_dd_to_hgfn.f90'
00102 include '../Subroutine/copy_to_hgfn_and_gfnsf.f90'
00103 include '../Subroutine/sourceterm_poisson_solver_test.f90'
00104 include '../Subroutine/sourceterm_surface_int_homosol.f90'
00105 include '../Subroutine/sourceterm_exsurf_eqm_binary.f90'
00106 include '../Subroutine/sourceterm_surface_int.f90'
00107 include '../Subroutine/poisson_solver_binary_bhex.f90'
00108 include '../Subroutine/poisson_solver_binary_bhex_homosol.f90'
00109 include '../Subroutine/poisson_solver_binary_vol_int.f90'
00110 include '../Subroutine/poisson_solver_binary_surf_int.f90'
00111 include '../Subroutine/poisson_solver_bhex_surf_int.f90'
00112 include '../Subroutine/poisson_solver_bhex_surf_int_all.f90'
00113 include '../Subroutine/grdr_gridpoint_type0.f90'
00114 include '../Subroutine/grdr_gridpoint_type0_nosym.f90'
00115 include '../Subroutine/grdr_gridpoint_type0_3rd_nosym.f90'
00116 include '../Subroutine/grgrad_4th_gridpoint_bhex.f90'
00117 include '../Subroutine/grgrad_midpoint_r3rd_type0.f90'
00118 include '../Subroutine/grgrad_midpoint_r4th_type0.f90'
00119 include '../Subroutine/grgrad_midpoint_type0.f90'
00120 include '../Subroutine/interpo_linear1p_type0_2Dsurf.f90'
00121 include '../Subroutine/interpo_gr2gr_4th.f90'
00122 include '../Subroutine/interpolation_fillup_binary.f90'
00123
00124 include '../Subroutine/allocate_BBH_CF.f90'
00125 include '../Subroutine/allocate_BBH_CF_AH.f90'
00126 include '../Subroutine/iteration_BBH_CF.f90'
00127 include '../Subroutine/iteration_BBH_CF_type1.f90'
00128 include '../Subroutine/iteration_BBH_CF_AH.f90'
00129 include '../Subroutine/bh_boundary_d_psi.f90'
00130 include '../Subroutine/bh_boundary_nh_psi_test.f90'
00131 include '../Subroutine/bh_boundary_d_alps.f90'
00132 include '../Subroutine/bh_boundary_d_bvxd.f90'
00133 include '../Subroutine/bh_boundary_d_bvyd.f90'
00134 include '../Subroutine/bh_boundary_d_bvzd.f90'
00135 include '../Subroutine/bh_boundary_d_Bfun.f90'
00136 include '../Subroutine/bh_boundary_n_Bfun.f90'
00137 include '../Subroutine/bh_boundary_d_potx.f90'
00138 include '../Subroutine/bh_boundary_d_poty.f90'
00139 include '../Subroutine/bh_boundary_d_potz.f90'
00140 include '../Subroutine/bh_boundary_AH.f90'
00141 include '../Subroutine/bh_boundary_CF.f90'
00142 include '../Subroutine/outer_boundary_d_psi.f90'
00143 include '../Subroutine/outer_boundary_d_alps.f90'
00144 include '../Subroutine/outer_boundary_d_bvxd.f90'
00145 include '../Subroutine/outer_boundary_d_bvyd.f90'
00146 include '../Subroutine/outer_boundary_d_bvzd.f90'
00147 include '../Subroutine/outer_boundary_d_Bfun.f90'
00148 include '../Subroutine/outer_boundary_n_Bfun.f90'
00149 include '../Subroutine/outer_boundary_d_potx.f90'
00150 include '../Subroutine/outer_boundary_d_poty.f90'
00151 include '../Subroutine/outer_boundary_d_potz.f90'
00152 include '../Subroutine/sourceterm_HaC_CF.f90'
00153 include '../Subroutine/sourceterm_trG_CF.f90'
00154 include '../Subroutine/sourceterm_MoC_CF.f90'
00155 include '../Subroutine/sourceterm_MoC_CF_type1_bhex.f90'
00156 include '../Subroutine/sourceterm_Bfun.f90'
00157 include '../Subroutine/interpolation_fillup_binary_parity.f90'
00158 include '../Subroutine/printout_error_metric_combined.f90'
00159 include '../Subroutine/printout_error_metric_combined_B.f90'
00160 include '../Subroutine/sourceterm_exsurf_eqm_binary_parity.f90'
00161 include '../Subroutine/compute_shift_v2.f90'
00162 include '../Subroutine/compute_dBfun.f90'
00163 include '../Subroutine/reset_bh_boundary_all.f90'
00164 include '../Subroutine/reset_metric_CF.f90'
00165 include '../Subroutine/reset_metric_CF_AH.f90'
00166 include '../Subroutine/initial_metric_CF.f90'
00167 include '../Subroutine/reset_bh_boundary_AH.f90'
00168
00169 include '../Subroutine/calc_physical_quantities_BBH_CF.f90'
00170 include '../Subroutine/calc_mass_BBH_CF.f90'
00171 include '../Subroutine/calc_mass_BBH_CF_adm_vol.f90'
00172
00173 include '../Subroutine/surf_source_adm_mass.f90'
00174 include '../Subroutine/surf_source_komar_mass.f90'
00175 include '../Subroutine/source_ang_mom_inf.f90'
00176 include '../Subroutine/source_ang_mom_thr.f90'
00177 include '../Subroutine/source_ang_mom_exc.f90'
00178 include '../Subroutine/source_ang_mom_smarr.f90'
00179 include '../Subroutine/printout_physq_console_BBH.f90'
00180 include '../Subroutine/calc_mass_ir.f90'
00181 include '../Subroutine/calc_circular_orbit.f90'
00182 include '../Subroutine/calc_ang_mom_BBH_CF_inf.f90'
00183 include '../Subroutine/calc_ang_mom_BBH_CF_thr.f90'
00184 include '../Subroutine/calc_ang_mom_BBH_CF_smarr.f90'
00185 include '../Subroutine/save_solution.f90'
00186 include '../Subroutine/write_omega_last.f90'
00187 include '../Subroutine/test_excurve.f90'
00188 include '../Subroutine/modify_r0_excurve.f90'
00189 include '../Subroutine/calc_app_hor_area_BBH_CF.f90'
00190 include '../Subroutine/surf_int_grav_rg.f90'
00191 include '../Subroutine/vol_int_grav_bhex.f90'
00192
00193 include '../Function/lagint_2nd.f90'
00194
00195 include '../Analysis/Subroutine/IO_output_BBH_CF.f90'
00196
00197
00198
00199
00200
00201 PROGRAM Main_BBH_CF_seq
00202
00203 use grid_parameter, only : indata_type, outdata_type, iter_max
00204 use radial_green_fn_grav
00205 use radial_green_fn_grav_bhex_nb
00206 use radial_green_fn_grav_bhex_dd
00207 use radial_green_fn_grav_bhex_nd
00208 use grid_parameter_binary_excision
00209 use grid_points_binary_excision
00210 use weight_midpoint_binary_excision
00211 use interface_copy_to_hgfn_and_gfnsf
00212 implicit none
00213 integer :: iseq, iter_count, total_iteration, flag_omega, count_adj
00214
00215 call coordinate_patch_kit_bhex
00216 call allocate_hgfn_bhex
00217 call allocate_hgfn_bhex_dd
00218 call calc_hgfn_bhex_dd
00219 call allocate_hgfn_bhex_nd
00220 call calc_hgfn_bhex_nd
00221
00222 call allocate_hgfn_bhex_nb
00223 call calc_hgfn_bhex_nb
00224
00225
00226 call read_parameter_bh
00227 call read_parameter_binary_excision
00228 call calc_parameter_binary_excision
00229 call IO_printout_grid_data
00230 call allocate_grid_points_binary_excision
00231 call calc_grid_points_binary_excision
00232 call allocate_weight_midpoint_binary_excision
00233 call calc_weight_midpoint_binary_excision
00234 call allocate_BBH_CF_AH
00235 if (indata_type.eq.'IN') call initial_metric_CF
00236 if (indata_type.eq.'3D') call IO_input_initial_3D
00237
00238
00239 total_iteration = 0
00240 call calc_circular_orbit(iter_count)
00241
00242 if (outdata_type.eq.'3D') call IO_output_solution_3D
00243
00244 END PROGRAM Main_BBH_CF_seq
