00001 subroutine IO_output_cartesian_contour_CF
00002 use phys_constant, only : long
00003 use def_metric_cartesian
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) :: psic, alphc, bvxdc, bvydc, bvzdc
00009 integer :: ix, iy, iz
00010
00011 open(14,file='rns_contour_xy_CF.dat',status='unknown')
00012 iz = nx_mid
00013 do iy = 1, ny
00014 write(14,*) ' '
00015 do ix = 1, nx
00016 psic = psica(ix,iy,iz)
00017 alphc = alphca(ix,iy,iz)
00018 bvxdc = bvxdca(ix,iy,iz)
00019 bvydc = bvydca(ix,iy,iz)
00020 bvzdc = bvzdca(ix,iy,iz)
00021 if (dabs(psic).le.small) psic = 0.0d0
00022 if (dabs(alphc).le.small) alphc = 0.0d0
00023 if (dabs(bvxdc).le.small) bvxdc = 0.0d0
00024 if (dabs(bvydc).le.small) bvydc = 0.0d0
00025 if (dabs(bvzdc).le.small) bvzdc = 0.0d0
00026 write(14,'(20es14.6)') x(ix), y(iy), psic, alphc, bvxdc, bvydc, bvzdc
00027
00028 end do
00029 end do
00030 close(14)
00031
00032 open(14,file='rns_contour_xz_CF.dat',status='unknown')
00033 iy = nx_mid
00034 do iz = 1, nz
00035 write(14,*) ' '
00036 do ix = 1, nx
00037 psic = psica(ix,iy,iz)
00038 alphc = alphca(ix,iy,iz)
00039 bvxdc = bvxdca(ix,iy,iz)
00040 bvydc = bvydca(ix,iy,iz)
00041 bvzdc = bvzdca(ix,iy,iz)
00042 if (dabs(psic).le.small) psic = 0.0d0
00043 if (dabs(alphc).le.small) alphc = 0.0d0
00044 if (dabs(bvxdc).le.small) bvxdc = 0.0d0
00045 if (dabs(bvydc).le.small) bvydc = 0.0d0
00046 if (dabs(bvzdc).le.small) bvzdc = 0.0d0
00047 write(14,'(20es14.6)') x(ix), z(iz), psic, alphc, bvxdc, bvydc, bvzdc
00048
00049 end do
00050 end do
00051 close(14)
00052
00053 open(14,file='rns_contour_yz_CF.dat',status='unknown')
00054 ix = nx_mid
00055 do iz = 1, nz
00056 write(14,*) ' '
00057 do iy = 1, ny
00058 psic = psica(ix,iy,iz)
00059 alphc = alphca(ix,iy,iz)
00060 bvxdc = bvxdca(ix,iy,iz)
00061 bvydc = bvydca(ix,iy,iz)
00062 bvzdc = bvzdca(ix,iy,iz)
00063 if (dabs(psic).le.small) psic = 0.0d0
00064 if (dabs(alphc).le.small) alphc = 0.0d0
00065 if (dabs(bvxdc).le.small) bvxdc = 0.0d0
00066 if (dabs(bvydc).le.small) bvydc = 0.0d0
00067 if (dabs(bvzdc).le.small) bvzdc = 0.0d0
00068 write(14,'(20es14.6)') y(iy), z(iz), psic, alphc, bvxdc, bvydc, bvzdc
00069
00070 end do
00071 end do
00072 close(14)
00073
00074 end subroutine IO_output_cartesian_contour_CF
