00001 subroutine IO_output_cartesian_contour_FLU
00002 use phys_constant, only : long
00003 use def_matter_cartesian, only : emdca, vxca, vyca, vzca, utca
00004 use grid_parameter_cartesian, only : nx, ny, nz, nx_mid
00005 use coordinate_grav_xyz, only : x, y, z
00006 implicit none
00007 real(long) :: small = 1.0d-20
00008 real(long) :: emdc, vxc, vyc, vzc, utc
00009 integer :: ix, iy, iz
00010
00011 open(14,file='rns_contour_xy_FLU.dat',status='unknown')
00012 iz = nx_mid
00013 do iy = 1, ny
00014 write(14,*) ' '
00015 do ix = 1, nx
00016 emdc = emdca(ix,iy,iz)
00017 vxc = vxca(ix,iy,iz)
00018 vyc = vyca(ix,iy,iz)
00019 vzc = vzca(ix,iy,iz)
00020 utc = utca(ix,iy,iz)
00021 if (dabs(emdc).le.small) emdc = 0.0d0
00022 if (dabs(vxc).le.small) vxc = 0.0d0
00023 if (dabs(vyc).le.small) vyc = 0.0d0
00024 if (dabs(vzc).le.small) vzc = 0.0d0
00025 if (dabs(utc).le.small) utc = 0.0d0
00026
00027 write(14,'(20es14.6)') x(ix), y(iy), emdc, vxc, vyc, vzc, utc
00028
00029 end do
00030 end do
00031 close(14)
00032
00033 open(14,file='rns_contour_xz_FLU.dat',status='unknown')
00034 iy = nx_mid
00035 do iz = 1, nz
00036 write(14,*) ' '
00037 do ix = 1, nx
00038 emdc = emdca(ix,iy,iz)
00039 vxc = vxca(ix,iy,iz)
00040 vyc = vyca(ix,iy,iz)
00041 vzc = vzca(ix,iy,iz)
00042 utc = utca(ix,iy,iz)
00043 if (dabs(emdc).le.small) emdc = 0.0d0
00044 if (dabs(vxc).le.small) vxc = 0.0d0
00045 if (dabs(vyc).le.small) vyc = 0.0d0
00046 if (dabs(vzc).le.small) vzc = 0.0d0
00047 if (dabs(utc).le.small) utc = 0.0d0
00048
00049 write(14,'(20es14.6)') x(ix), z(iz), emdc, vxc, vyc, vzc, utc
00050
00051 end do
00052 end do
00053 close(14)
00054
00055 open(14,file='rns_contour_yz_FLU.dat',status='unknown')
00056 ix = nx_mid
00057 do iz = 1, nz
00058 write(14,*) ' '
00059 do iy = 1, ny
00060 emdc = emdca(ix,iy,iz)
00061 vxc = vxca(ix,iy,iz)
00062 vyc = vyca(ix,iy,iz)
00063 vzc = vzca(ix,iy,iz)
00064 utc = utca(ix,iy,iz)
00065 if (dabs(emdc).le.small) emdc = 0.0d0
00066 if (dabs(vxc).le.small) vxc = 0.0d0
00067 if (dabs(vyc).le.small) vyc = 0.0d0
00068 if (dabs(vzc).le.small) vzc = 0.0d0
00069 if (dabs(utc).le.small) utc = 0.0d0
00070
00071 write(14,'(20es14.6)') y(iy), z(iz), emdc, vxc, vyc, vzc, utc
00072
00073 end do
00074 end do
00075 close(14)
00076
00077 end subroutine IO_output_cartesian_contour_FLU
