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

Go to the documentation of this file.

subroutine IO_output_2D_general(filename,fg,mg,pot,pl)
  use phys_constant, only : long
  use grid_parameter, only  :   nrg, ntg, npg, nrf, ntf, npf
  use trigonometry_grav_theta, only : sinthg, costhg
  use trigonometry_grav_phi,   only : sinphig, cosphig
  use coordinate_grav_r, only : rg
  implicit none
  integer :: ir, it, ip, nrtmp, nttmp, nptmp, ii, ir0,it0,ip0
  real(long), pointer  :: pot(:,:,:)
  real(long) :: xx, yy, zz
  character(len=1), intent(in) :: fg,mg
  character(len=2), intent(in) :: pl
  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 (pl=='xy' .or. pl=='yx')  then
    it = nttmp/2
    do ip = ip0, nptmp
      do ir = ir0, nrtmp
        xx = rg(ir)*sinthg(it)*cosphig(ip)
        yy = rg(ir)*sinthg(it)*sinphig(ip)
        write(12,'(1p,3e23.15)') xx, yy, pot(ir,it,ip)
      end do
    end do
  end if
!
  if (pl=='yz' .or. pl=='zy')  then
    ip = nptmp/4
    do it = it0, nttmp
      do ir = ir0, nrtmp
        yy = rg(ir)*sinthg(it)*sinphig(ip)
        zz = rg(ir)*costhg(it)
        write(12,'(1p,3e23.15)') yy, zz, pot(ir,it,ip)
      end do
    end do
!
    ip = 3*nptmp/4
    do it = it0, nttmp
      do ir = ir0, nrtmp
        yy = rg(ir)*sinthg(it)*sinphig(ip)
        zz = rg(ir)*costhg(it)
        write(12,'(1p,3e23.15)') yy, zz, pot(ir,it,ip)
      end do
    end do
  end if
!
  if (pl=='zx' .or. pl=='xz')  then
    ip = 0
    do it = it0, nttmp
      do ir = ir0, nrtmp
        xx = rg(ir)*sinthg(it)*cosphig(ip)
        zz = rg(ir)*costhg(it)
        write(12,'(1p,3e23.15)') xx, zz, pot(ir,it,ip)
      end do
    end do
!
    ip = nptmp/2
    do it = it0, nttmp
      do ir = ir0, nrtmp
        xx = rg(ir)*sinthg(it)*cosphig(ip)
        zz = rg(ir)*costhg(it)
        write(12,'(1p,3e23.15)') xx, zz, pot(ir,it,ip)
      end do
    end do
  end if
!
  close(12)
!
end subroutine IO_output_2D_general

---