00001 subroutine coordinate_patch_kit_grav_grid_mpt(igrid)
00002 use grid_parameter
00003 use coordinate_grav_r
00004 use coordinate_grav_phi
00005 use coordinate_grav_theta
00006 use coordinate_grav_extended
00007 use weight_midpoint_grav
00008 use weight_midpoint_fluid
00009 use trigonometry_grav_theta
00010 use trigonometry_grav_phi
00011 use legendre_fn_grav
00012 implicit none
00013 integer :: igrid
00014
00015 if (igrid==1) then
00016 write(6,*) "***Default grid... ****"
00017 call grid_r
00018 else if (igrid==2) then
00019 write(6,*) "***Black hole grid... ****"
00020 call grid_r_bhex('eBH')
00021 else if (igrid==3) then
00022 write(6,*) "***Neutron star grid with rg(nrf)<1... ****"
00023 call grid_r_bns
00024 else if (igrid==4) then
00025 write(6,*) "***Neutron star grid with rg(nrf)=1 and constant dr until r=7... ****"
00026 call grid_r_bns_const
00027 else
00028 write(6,*) "***Choose a grid: 1 or 2 or 3 or 4...exiting****"
00029 stop
00030 end if
00031
00032 call grid_theta
00033 call trig_grav_theta
00034 call legendre
00035 call grid_phi
00036 call trig_grav_phi
00037 call weight_calc_midpoint_grav
00038 call weight_calc_midpoint_grav_th4th
00039 call weight_calc_midpoint_fluid
00040 call grid_extended
00041 end subroutine coordinate_patch_kit_grav_grid_mpt
00042
