Cocal: /home/uryu/Cocal_svn_checkinout/Cocal/code/Analysis/Subroutine/coordinate_patch_kit_grav_grid_coc2cac

00001 subroutine coordinate_patch_kit_grav_grid_coc2cac_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   implicit none
00012   integer :: igrid
00013 ! call subroutines. the order is important.
00014   if (igrid==1) then
00015     write(6,*) "Default grid..."
00016     call grid_r
00017   else if (igrid==2) then
00018     write(6,*) "Black hole grid..."
00019     call grid_r_bhex('eBH')
00020   else if (igrid==3) then
00021     write(6,*) "Neutron star grid with rg(nrf)<1..."
00022     call grid_r_bns
00023   else if (igrid==4) then
00024     write(6,*) "Neutron star grid with rg(nrf)=1 and constant dr until r=7..."
00025     call grid_r_bns_const
00026   else
00027     write(6,*) "Choose a grid: 1 or 2 or 3 or 4...exiting"
00028     stop
00029   end if
00030 
00031   call grid_theta
00032   call trig_grav_theta
00033   call grid_phi
00034   call trig_grav_phi
00035 !  call weight_calc_midpoint_grav
00036 !  call weight_calc_midpoint_grav_th4th
00037 !  call weight_calc_midpoint_fluid
00038   call grid_extended
00039 end subroutine coordinate_patch_kit_grav_grid_coc2cac_mpt
00040

/home/uryu/Cocal_svn_checkinout/Cocal/code/Analysis/Subroutine/coordinate_patch_kit_grav_grid_coc2cac_mpt.f90