/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine_mpatch/IO_input_potential_test_3D_mpt.f90

Go to the documentation of this file.

subroutine IO_input_potential_test_3D_mpt(impt)
  use phys_constant, only : long, nnrg, nntg, nnpg
  use def_metric, only : alph, psi, bvxd, bvyd, bvzd
  use def_matter, only : emd, rs
  use def_matter_parameter, only : ome, ber, radi
  use def_binary_parameter, only : dis
  use coordinate_grav_r, only : rg
  use coordinate_grav_theta, only : thg
  use coordinate_grav_phi, only : phig 
  implicit none
  integer :: ir, it, ip, nrtmp, nttmp, nptmp
  integer :: impt
  real(long) :: rgtmp(0:nnrg), thgtmp(0:nntg), phigtmp(0:nnpg)
  character(len=1) :: np(5) = (/'1', '2','3', '4', '5'/)
!
!
! --- Metric potentials.
  open(13,file='rnsgra_3D_mpt'//np(impt)//'.las',status='old')
  read(13,'(5i5)') nrtmp, nttmp, nptmp
  do ip = 0, nptmp
    do it = 0, nttmp
      do ir = 0, nrtmp
        read(13,'(1p,6e20.12)')  psi(ir,it,ip), &
    &                           alph(ir,it,ip)
      end do
    end do
  end do
  read(13,'(1p,6e20.12)') ome, ber, radi, dis
  close(13)
!
! --- Coordinate grids.
  open(14,file='rnsgrids_3D_mpt'//np(impt)//'.las',status='old')
  read(14,'(5i5)') nrtmp, nttmp, nptmp
  do ir = 0, nrtmp
    read(14,'(1p,6e20.12)') rgtmp(ir)
  end do
  do it = 0, nttmp
    read(14,'(1p,6e20.12)') thgtmp(it)
  end do
  do ip = 0, nptmp
    read(14,'(1p,6e20.12)') phigtmp(ip)
  end do
  close(14)
!
end subroutine IO_input_potential_test_3D_mpt

---