00001 !______________________________________________ 00002 include '../Include_file/include_modulefiles_peos.f90' 00003 include '../Include_file/include_interface_utilities.f90' 00004 include '../Include_file/include_subroutines_utilities.f90' 00005 ! 00006 include '../Module/def_matter_velocity.f90' 00007 ! 00008 include '../Subroutine/allocate_matter_emdrsrhof.f90' 00009 include '../Subroutine/allocate_matter_3velocity.f90' 00010 include '../Subroutine/allocate_matter_4velocity.f90' 00011 include '../Subroutine/calc_3velocity_corot.f90' 00012 include '../Subroutine/calc_4velocity_corot.f90' 00013 include '../Subroutine/calc_matter_rhof.f90' 00014 ! 00015 include '../Subroutine/IO_input_converged_solution_fluid_3D.f90' 00016 include '../Subroutine/IO_output_solution_fluid_3D.f90' 00017 ! 00018 include '../Analysis/Subroutine/coordinate_patch_kit_grav_noGreen.f90' 00019 !______________________________________________ 00020 ! 00021 ! Interpolation Program 00022 !______________________________________________ 00023 PROGRAM output_fluid_velocity 00024 ! 00025 use def_vector_phi 00026 implicit none 00027 ! 00028 call coordinate_patch_kit_grav_noGreen 00029 call peos_initialize 00030 call allocate_matter_emdrsrhof 00031 call allocate_matter_3velocity 00032 call allocate_matter_4velocity 00033 call allocate_vector_phi 00034 call IO_input_converged_solution_fluid_3D 00035 ! 00036 call calc_vector_phi_matter(1) 00037 call calc_3velocity_corot 00038 call calc_4velocity_corot 00039 call calc_matter_rhof 00040 ! 00041 call IO_output_solution_fluid_3D 00042 ! 00043 END PROGRAM output_fluid_velocity