00001 subroutine IO_output_cartesian_contour_qeos
00002 use phys_constant, only : long
00003 use def_metric_cartesian, only : alphca, psica, bvxdca, bvydca, bvzdca
00004 use def_matter_parameter, only : rhos_qs
00005 use def_matter_cartesian_qeos, only : rhoca, vxca, vyca, vzca
00006 use grid_parameter_cartesian, only : nx, ny, nz, nx_mid
00007 use coordinate_grav_xyz, only : x, y, z
00008
00009 implicit none
00010 real(long) :: small = 1.0d-20
00011 real(long) :: rhoc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
00012 integer :: ix, iy, iz
00013
00014 open(14,file='rns_contour_xy.dat',status='unknown')
00015 iz = nx_mid
00016 do iy = 1, ny
00017 write(14,*) ' '
00018 do ix = 1, nx
00019 rhoc = rhoca(ix,iy,iz)
00020 vxc = vxca(ix,iy,iz)
00021 vyc = vyca(ix,iy,iz)
00022 vzc = vzca(ix,iy,iz)
00023 psic = psica(ix,iy,iz)
00024 alphc=alphca(ix,iy,iz)
00025 bvxc = bvxdca(ix,iy,iz)
00026 bvyc = bvydca(ix,iy,iz)
00027 bvzc = bvzdca(ix,iy,iz)
00028 if (dabs(rhoc).le.rhos_qs) rhoc = 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 & rhoc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
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 rhoc = rhoca(ix,iy,iz)
00051 vxc = vxca(ix,iy,iz)
00052 vyc = vyca(ix,iy,iz)
00053 vzc = vzca(ix,iy,iz)
00054 psic = psica(ix,iy,iz)
00055 alphc=alphca(ix,iy,iz)
00056 bvxc = bvxdca(ix,iy,iz)
00057 bvyc = bvydca(ix,iy,iz)
00058 bvzc = bvzdca(ix,iy,iz)
00059 if (dabs(rhoc).le.rhos_qs) rhoc = 0.0d0
00060 if (dabs(vxc).le.small) vxc = 0.0d0
00061 if (dabs(vyc).le.small) vyc = 0.0d0
00062 if (dabs(vzc).le.small) vzc = 0.0d0
00063 if (dabs(psic).le.small) psic = 0.0d0
00064 if (dabs(alphc).le.small) alphc = 0.0d0
00065 if (dabs(bvxc).le.small) bvxc = 0.0d0
00066 if (dabs(bvyc).le.small) bvyc = 0.0d0
00067 if (dabs(bvzc).le.small) bvzc = 0.0d0
00068
00069 write(14,'(20es14.6)') x(ix), z(iz), &
00070 & rhoc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
00071
00072 end do
00073 end do
00074 close(14)
00075
00076 open(14,file='rns_contour_yz.dat',status='unknown')
00077 ix = nx_mid
00078 do iz = 1, nz
00079 write(14,*) ' '
00080 do iy = 1, ny
00081 rhoc = rhoca(ix,iy,iz)
00082 vxc = vxca(ix,iy,iz)
00083 vyc = vyca(ix,iy,iz)
00084 vzc = vzca(ix,iy,iz)
00085 psic = psica(ix,iy,iz)
00086 alphc=alphca(ix,iy,iz)
00087 bvxc = bvxdca(ix,iy,iz)
00088 bvyc = bvydca(ix,iy,iz)
00089 bvzc = bvzdca(ix,iy,iz)
00090 if (dabs(rhoc).le.rhos_qs) rhoc = 0.0d0
00091 if (dabs(vxc).le.small) vxc = 0.0d0
00092 if (dabs(vyc).le.small) vyc = 0.0d0
00093 if (dabs(vzc).le.small) vzc = 0.0d0
00094 if (dabs(psic).le.small) psic = 0.0d0
00095 if (dabs(alphc).le.small) alphc = 0.0d0
00096 if (dabs(bvxc).le.small) bvxc = 0.0d0
00097 if (dabs(bvyc).le.small) bvyc = 0.0d0
00098 if (dabs(bvzc).le.small) bvzc = 0.0d0
00099
00100 write(14,'(20es14.6)') y(iy), z(iz), &
00101 & rhoc, vxc, vyc, vzc, psic, alphc, bvxc, bvyc, bvzc
00102
00103 end do
00104 end do
00105 close(14)
00106
00107 end subroutine IO_output_cartesian_contour_qeos
