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

00001 !______________________________________________
00002 include '../Include_file/include_modulefiles_BBH_CF.f90'
00003 include '../Include_file/include_interface_modulefiles_BBH_CF.f90'
00004 include '../Include_file/include_subroutines_BBH_CF.f90'
00005 include '../Include_file/include_functions.f90'
00006 !______________________________________________
00007 !
00008 !              Main Program
00009 !______________________________________________
00010 PROGRAM Main_BBH_CF
00011 !
00012   use grid_parameter, only : indata_type, outdata_type, iter_max
00013   use radial_green_fn_grav
00014   use radial_green_fn_grav_bhex_nb
00015   use radial_green_fn_grav_bhex_dd
00016   use radial_green_fn_grav_bhex_nd
00017   use grid_parameter_binary_excision
00018   use grid_points_binary_excision
00019   use weight_midpoint_binary_excision
00020   use interface_copy_to_hgfn_and_gfnsf
00021   implicit none
00022   integer :: iseq, iter_count, total_iteration
00023 !
00024   call coordinate_patch_kit_bhex
00025   call allocate_hgfn_bhex
00026   call allocate_hgfn_bhex_dd
00027   call calc_hgfn_bhex_dd
00028   call allocate_hgfn_bhex_nd
00029   call calc_hgfn_bhex_nd
00030 ! -- No boundary Green's fn
00031   call allocate_hgfn_bhex_nb
00032   call calc_hgfn_bhex_nb
00033 !  call copy_to_hgfn_and_gfnsf(hgfn_nb,gfnsf_nb)
00034 ! --
00035   call read_parameter_bh
00036   call read_parameter_binary_excision
00037   call calc_parameter_binary_excision
00038   call allocate_grid_points_binary_excision
00039   call calc_grid_points_binary_excision
00040   call allocate_weight_midpoint_binary_excision
00041   call calc_weight_midpoint_binary_excision
00042   call allocate_BBH_CF
00043   if (indata_type.eq.'IN') call initial_metric_CF
00044   if (indata_type.eq.'3D') call IO_input_initial_3D
00045 !  call test_analytic_solution_bhex_psialph
00046 ! 
00047   call iteration_BBH_CF(iter_count)
00048   if (total_iteration.ge.iter_max) then
00049     write(6,*)' ** Solution did not converge **'
00050   end if
00051 !  call printout_NS_shape_seq(iseq)
00052   if (outdata_type.eq.'3D') call IO_output_solution_3D
00053 !!  call printout_debug
00054 !  call printout_NS_shape
00055 !
00056 END PROGRAM Main_BBH_CF
/home/uryu/Cocal_svn_checkinout/Cocal/code/Main_code/Main_BBH_CF.f90
