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

Go to the documentation of this file.

subroutine IO_output_cartesian_contour_potential_2pot_test
  use phys_constant, only : long
  use def_metric_cartesian, only  : psica, alphca
  use grid_parameter_cartesian, only  : nx, ny, nz, nx_mid
  use coordinate_grav_xyz, only  : x, y, z
  implicit none
  real(long) :: small = 1.0d-20
  real(long) :: psic, alphc
  integer :: ix, iy, iz
!
  open(14,file='rns_contour_xy.dat',status='unknown')
  iz = nx_mid
  do iy = 1, ny
    write(14,*) ' '
    do ix = 1, nx
      psic = psica(ix,iy,iz)
      alphc = alphca(ix,iy,iz)
      if (dabs(psic).le.small) psic = 0.0d0
      if (dabs(alphc).le.small) alphc = 0.0d0
      write(14,'(20es14.6)') x(ix), y(iy), psic, alphc
!
    end do
  end do
  close(14)
!
  open(14,file='rns_contour_xz.dat',status='unknown')
  iy = nx_mid
  do iz = 1, nz
    write(14,*) ' '
    do ix = 1, nx
      psic = psica(ix,iy,iz)
      alphc = alphca(ix,iy,iz)
      if (dabs(psic).le.small) psic = 0.0d0
      if (dabs(alphc).le.small) alphc = 0.0d0
      write(14,'(20es14.6)') x(ix), z(iz), psic, alphc
!
    end do
  end do
  close(14)
!
  open(14,file='rns_contour_yz.dat',status='unknown')
  ix = nx_mid
  do iz = 1, nz
    write(14,*) ' '
    do iy = 1, ny
      psic = psica(ix,iy,iz)
      alphc = alphca(ix,iy,iz)
      if (dabs(psic).le.small) psic = 0.0d0
      if (dabs(alphc).le.small) alphc = 0.0d0
      write(14,'(20es14.6)') y(iy), z(iz), psic, alphc
!
    end do
  end do
  close(14)
!
end subroutine IO_output_cartesian_contour_potential_2pot_test

---