00001 subroutine IO_output_cartesian_contour
00002 use phys_constant, only : long
00003 use def_metric_cartesian, only : alphca, psica, bvxdca, bvydca, bvzdca
00004 use def_matter_cartesian, only : emdca, vxca, vyca, vzca
00005 use grid_parameter_cartesian, only : nx, ny, nz, nx_mid
00006 use coordinate_grav_xyz, only : x, y, z
00007 implicit none
00008 real(long) :: small = 1.0d-20
00009 real(long) :: emdc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
00010 integer :: ix, iy, iz
00011
00012 open(14,file='rns_contour_xy.dat',status='unknown')
00013 iz = nx_mid
00014 do iy = 1, ny
00015 write(14,*) ' '
00016 do ix = 1, nx
00017 emdc = emdca(ix,iy,iz)
00018 vxc = vxca(ix,iy,iz)
00019 vyc = vyca(ix,iy,iz)
00020 vzc = vzca(ix,iy,iz)
00021 psic = psica(ix,iy,iz)
00022 alphc=alphca(ix,iy,iz)
00023 bvxc = bvxdca(ix,iy,iz)
00024 bvyc = bvydca(ix,iy,iz)
00025 bvzc = bvzdca(ix,iy,iz)
00026 if (dabs(emdc).le.small) emdc = 0.0d0
00027 if (dabs(vxc).le.small) vxc = 0.0d0
00028 if (dabs(vyc).le.small) vyc = 0.0d0
00029 if (dabs(vzc).le.small) vzc = 0.0d0
00030 if (dabs(psic).le.small) psic = 0.0d0
00031 if (dabs(alphc).le.small) alphc = 0.0d0
00032 if (dabs(bvxc).le.small) bvxc = 0.0d0
00033 if (dabs(bvyc).le.small) bvyc = 0.0d0
00034 if (dabs(bvzc).le.small) bvzc = 0.0d0
00035
00036 write(14,'(20es14.6)') x(ix), y(iy), &
00037 & emdc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
00038
00039 end do
00040 end do
00041 close(14)
00042
00043 open(14,file='rns_contour_xz.dat',status='unknown')
00044 iy = nx_mid
00045 do iz = 1, nz
00046 write(14,*) ' '
00047 do ix = 1, nx
00048 emdc = emdca(ix,iy,iz)
00049 vxc = vxca(ix,iy,iz)
00050 vyc = vyca(ix,iy,iz)
00051 vzc = vzca(ix,iy,iz)
00052 psic = psica(ix,iy,iz)
00053 alphc=alphca(ix,iy,iz)
00054 bvxc = bvxdca(ix,iy,iz)
00055 bvyc = bvydca(ix,iy,iz)
00056 bvzc = bvzdca(ix,iy,iz)
00057 if (dabs(emdc).le.small) emdc = 0.0d0
00058 if (dabs(vxc).le.small) vxc = 0.0d0
00059 if (dabs(vyc).le.small) vyc = 0.0d0
00060 if (dabs(vzc).le.small) vzc = 0.0d0
00061 if (dabs(psic).le.small) psic = 0.0d0
00062 if (dabs(alphc).le.small) alphc = 0.0d0
00063 if (dabs(bvxc).le.small) bvxc = 0.0d0
00064 if (dabs(bvyc).le.small) bvyc = 0.0d0
00065 if (dabs(bvzc).le.small) bvzc = 0.0d0
00066
00067 write(14,'(20es14.6)') x(ix), z(iz), &
00068 & emdc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
00069
00070 end do
00071 end do
00072 close(14)
00073
00074 open(14,file='rns_contour_yz.dat',status='unknown')
00075 ix = nx_mid
00076 do iz = 1, nz
00077 write(14,*) ' '
00078 do iy = 1, ny
00079 emdc = emdca(ix,iy,iz)
00080 vxc = vxca(ix,iy,iz)
00081 vyc = vyca(ix,iy,iz)
00082 vzc = vzca(ix,iy,iz)
00083 psic = psica(ix,iy,iz)
00084 alphc=alphca(ix,iy,iz)
00085 bvxc = bvxdca(ix,iy,iz)
00086 bvyc = bvydca(ix,iy,iz)
00087 bvzc = bvzdca(ix,iy,iz)
00088 if (dabs(emdc).le.small) emdc = 0.0d0
00089 if (dabs(vxc).le.small) vxc = 0.0d0
00090 if (dabs(vyc).le.small) vyc = 0.0d0
00091 if (dabs(vzc).le.small) vzc = 0.0d0
00092 if (dabs(psic).le.small) psic = 0.0d0
00093 if (dabs(alphc).le.small) alphc = 0.0d0
00094 if (dabs(bvxc).le.small) bvxc = 0.0d0
00095 if (dabs(bvyc).le.small) bvyc = 0.0d0
00096 if (dabs(bvzc).le.small) bvzc = 0.0d0
00097
00098 write(14,'(20es14.6)') y(iy), z(iz), &
00099 & emdc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
00100
00101 end do
00102 end do
00103 close(14)
00104
00105 end subroutine IO_output_cartesian_contour
