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, omeca
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, omec, 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 omec = omeca(ix,iy,iz)
00019 vxc = vxca(ix,iy,iz)
00020 vyc = vyca(ix,iy,iz)
00021 vzc = vzca(ix,iy,iz)
00022 psic = psica(ix,iy,iz)
00023 alphc=alphca(ix,iy,iz)
00024 bvxc = bvxdca(ix,iy,iz)
00025 bvyc = bvydca(ix,iy,iz)
00026 bvzc = bvzdca(ix,iy,iz)
00027 if (dabs(emdc).le.small) emdc = 0.0d0
00028 if (dabs(omec).le.small) omec = 0.0d0
00029 if (dabs(vxc).le.small) vxc = 0.0d0
00030 if (dabs(vyc).le.small) vyc = 0.0d0
00031 if (dabs(vzc).le.small) vzc = 0.0d0
00032 if (dabs(psic).le.small) psic = 0.0d0
00033 if (dabs(alphc).le.small) alphc = 0.0d0
00034 if (dabs(bvxc).le.small) bvxc = 0.0d0
00035 if (dabs(bvyc).le.small) bvyc = 0.0d0
00036 if (dabs(bvzc).le.small) bvzc = 0.0d0
00037
00038 write(14,'(20es14.6)') x(ix), y(iy), &
00039 & emdc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc, omec
00040
00041 end do
00042 end do
00043 close(14)
00044
00045 open(14,file='rns_contour_xz.dat',status='unknown')
00046 iy = nx_mid
00047 do iz = 1, nz
00048 write(14,*) ' '
00049 do ix = 1, nx
00050 emdc = emdca(ix,iy,iz)
00051 omec = omeca(ix,iy,iz)
00052 vxc = vxca(ix,iy,iz)
00053 vyc = vyca(ix,iy,iz)
00054 vzc = vzca(ix,iy,iz)
00055 psic = psica(ix,iy,iz)
00056 alphc=alphca(ix,iy,iz)
00057 bvxc = bvxdca(ix,iy,iz)
00058 bvyc = bvydca(ix,iy,iz)
00059 bvzc = bvzdca(ix,iy,iz)
00060 if (dabs(emdc).le.small) emdc = 0.0d0
00061 if (dabs(omec).le.small) omec = 0.0d0
00062 if (dabs(vxc).le.small) vxc = 0.0d0
00063 if (dabs(vyc).le.small) vyc = 0.0d0
00064 if (dabs(vzc).le.small) vzc = 0.0d0
00065 if (dabs(psic).le.small) psic = 0.0d0
00066 if (dabs(alphc).le.small) alphc = 0.0d0
00067 if (dabs(bvxc).le.small) bvxc = 0.0d0
00068 if (dabs(bvyc).le.small) bvyc = 0.0d0
00069 if (dabs(bvzc).le.small) bvzc = 0.0d0
00070
00071 write(14,'(20es14.6)') x(ix), z(iz), &
00072 & emdc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc, omec
00073
00074 end do
00075 end do
00076 close(14)
00077
00078 open(14,file='rns_contour_yz.dat',status='unknown')
00079 ix = nx_mid
00080 do iz = 1, nz
00081 write(14,*) ' '
00082 do iy = 1, ny
00083 emdc = emdca(ix,iy,iz)
00084 omec = omeca(ix,iy,iz)
00085 vxc = vxca(ix,iy,iz)
00086 vyc = vyca(ix,iy,iz)
00087 vzc = vzca(ix,iy,iz)
00088 psic = psica(ix,iy,iz)
00089 alphc=alphca(ix,iy,iz)
00090 bvxc = bvxdca(ix,iy,iz)
00091 bvyc = bvydca(ix,iy,iz)
00092 bvzc = bvzdca(ix,iy,iz)
00093 if (dabs(emdc).le.small) emdc = 0.0d0
00094 if (dabs(omec).le.small) omec = 0.0d0
00095 if (dabs(vxc).le.small) vxc = 0.0d0
00096 if (dabs(vyc).le.small) vyc = 0.0d0
00097 if (dabs(vzc).le.small) vzc = 0.0d0
00098 if (dabs(psic).le.small) psic = 0.0d0
00099 if (dabs(alphc).le.small) alphc = 0.0d0
00100 if (dabs(bvxc).le.small) bvxc = 0.0d0
00101 if (dabs(bvyc).le.small) bvyc = 0.0d0
00102 if (dabs(bvzc).le.small) bvzc = 0.0d0
00103
00104 write(14,'(20es14.6)') y(iy), z(iz), &
00105 & emdc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc, omec
00106
00107 end do
00108 end do
00109 close(14)
00110
00111 end subroutine IO_output_cartesian_contour
