/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/IO_output_1D_general.f90

Go to the documentation of this file.

subroutine IO_output_1D_general(filename,fg,mg,pot,jr,jt,jp)
  use phys_constant, only : long
  use grid_parameter, only  :   nrg, ntg, npg, nrf, ntf, npf
  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, ii, ir0,it0,ip0
  real(long), pointer  :: pot(:,:,:)
  integer, intent(in) :: jr, jt, jp
  character(len=1), intent(in) :: fg,mg
  character(len=30), intent(in) :: filename
!
  if(fg == 'f') then
    nrtmp=nrf;  nttmp=ntf;  nptmp=npf
  else
    nrtmp=nrg;  nttmp=ntg;  nptmp=npg
  end if

  if (mg=='m') then
    ir0=1;  it0=1;  ip0=1
  else
    ir0=0;  it0=0;  ip0=0
  end if

  open(12,file=filename, status='unknown')
  if (jr==-1)  then
    do ii=ir0,nrtmp
      write(12,'(1p,2e20.12)') rg(ii), pot(ii,jt,jp)
    end do
  end if
  if (jt==-1)  then
    do ii=it0,nttmp
      write(12,'(1p,2e20.12)') thg(ii), pot(jr,ii,jp)
    end do
  end if
  if (jp==-1)  then
    do ii=ip0,nptmp
      write(12,'(1p,2e20.12)') phig(ii), pot(jr,jt,ii)
    end do
  end if
  close(12)
!
end subroutine IO_output_1D_general

---