/home/uryu/Cocal_svn_checkinout/Cocal/code/Main_code/Main_AHfinder.f90
Go to the documentation of this file.00001 include '../Include_file/include_modulefiles_AHfinder.f90' 00002 include '../Include_file/include_interface_modulefiles_AHfinder.f90' 00003 include '../Include_file/include_subroutines_AHfinder.f90' 00004 include '../Include_file/include_functions.f90' 00005 ! 00006 program Main_AHfinder 00007 use grid_parameter, only : indata_type, outdata_type, iter_max 00008 use grid_parameter_binary_excision 00009 use grid_points_binary_excision 00010 use weight_midpoint_binary_excision 00011 implicit none 00012 integer :: iseq, iter_count, total_iteration 00013 ! 00014 call coordinate_patch_kit_bhex 00015 call read_parameter_bh 00016 call read_parameter_binary_excision 00017 call calc_parameter_binary_excision 00018 call allocate_grid_points_binary_excision 00019 call calc_grid_points_binary_excision 00020 call allocate_weight_midpoint_binary_excision 00021 call calc_weight_midpoint_binary_excision 00022 call allocate_BBH_CF 00023 call allocate_horizon 00024 call initial_AHfinder 00025 if (indata_type.eq.'IN') call initial_metric_CF 00026 if (indata_type.eq.'3D') call IO_input_initial_3D 00027 ! 00028 call calc_vector_bh(2) 00029 call excurve 00030 call iteration_AHfinder(iter_count) 00031 if (total_iteration.ge.iter_max) then 00032 write(6,*)' ** Solution did not converge **' 00033 end if 00034 ! 00035 call calc_AHarea_AHfinder 00036 call IO_output_AHfinder 00037 call IO_output_AHfinder_gnuplot 00038 ! 00039 end program Main_AHfinder
