00001 subroutine IO_output_cartesian_contour_WL
00002 use phys_constant, only : long
00003 use def_metric_hij_cartesian, only : hxxdca, hxydca, hxzdca, &
00004 & hyydca, hyzdca, hzzdca
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) :: hxxdc, hxydc, hxzdc, hyydc, hyzdc, hzzdc
00010 integer :: ix, iy, iz
00011
00012 open(14,file='rns_contour_hij_xy.dat',status='unknown')
00013 iz = nx_mid
00014 do iy = 1, ny
00015 write(14,*) ' '
00016 do ix = 1, nx
00017
00018 hxxdc = hxxdca(ix,iy,iz)
00019 hxydc = hxydca(ix,iy,iz)
00020 hxzdc = hxzdca(ix,iy,iz)
00021 hyydc = hyydca(ix,iy,iz)
00022 hyzdc = hyzdca(ix,iy,iz)
00023 hzzdc = hzzdca(ix,iy,iz)
00024 if (dabs(hxxdc).le.small) hxxdc = 0.0d0
00025 if (dabs(hxydc).le.small) hxydc = 0.0d0
00026 if (dabs(hxzdc).le.small) hxzdc = 0.0d0
00027 if (dabs(hyydc).le.small) hyydc = 0.0d0
00028 if (dabs(hyzdc).le.small) hyzdc = 0.0d0
00029 if (dabs(hzzdc).le.small) hzzdc = 0.0d0
00030
00031 write(14,'(20es14.6)') x(ix), y(iy), &
00032 & hxxdc, hxydc, hxzdc, hyydc, hyzdc, hzzdc
00033
00034 end do
00035 end do
00036 close(14)
00037
00038 open(14,file='rns_contour_hij_xz.dat',status='unknown')
00039 iy = nx_mid
00040 do iz = 1, nz
00041 write(14,*) ' '
00042 do ix = 1, nx
00043
00044 hxxdc = hxxdca(ix,iy,iz)
00045 hxydc = hxydca(ix,iy,iz)
00046 hxzdc = hxzdca(ix,iy,iz)
00047 hyydc = hyydca(ix,iy,iz)
00048 hyzdc = hyzdca(ix,iy,iz)
00049 hzzdc = hzzdca(ix,iy,iz)
00050 if (dabs(hxxdc).le.small) hxxdc = 0.0d0
00051 if (dabs(hxydc).le.small) hxydc = 0.0d0
00052 if (dabs(hxzdc).le.small) hxzdc = 0.0d0
00053 if (dabs(hyydc).le.small) hyydc = 0.0d0
00054 if (dabs(hyzdc).le.small) hyzdc = 0.0d0
00055 if (dabs(hzzdc).le.small) hzzdc = 0.0d0
00056
00057 write(14,'(20es14.6)') x(ix), z(iz), &
00058 & hxxdc, hxydc, hxzdc, hyydc, hyzdc, hzzdc
00059
00060 end do
00061 end do
00062 close(14)
00063
00064 open(14,file='rns_contour_hij_yz.dat',status='unknown')
00065 ix = nx_mid
00066 do iz = 1, nz
00067 write(14,*) ' '
00068 do iy = 1, ny
00069
00070 hxxdc = hxxdca(ix,iy,iz)
00071 hxydc = hxydca(ix,iy,iz)
00072 hxzdc = hxzdca(ix,iy,iz)
00073 hyydc = hyydca(ix,iy,iz)
00074 hyzdc = hyzdca(ix,iy,iz)
00075 hzzdc = hzzdca(ix,iy,iz)
00076 if (dabs(hxxdc).le.small) hxxdc = 0.0d0
00077 if (dabs(hxydc).le.small) hxydc = 0.0d0
00078 if (dabs(hxzdc).le.small) hxzdc = 0.0d0
00079 if (dabs(hyydc).le.small) hyydc = 0.0d0
00080 if (dabs(hyzdc).le.small) hyzdc = 0.0d0
00081 if (dabs(hzzdc).le.small) hzzdc = 0.0d0
00082
00083 write(14,'(20es14.6)') y(iy), z(iz), &
00084 & hxxdc, hxydc, hxzdc, hyydc, hyzdc, hzzdc
00085
00086 end do
00087 end do
00088 close(14)
00089
00090 end subroutine IO_output_cartesian_contour_WL
