00001 subroutine IO_output_solution_IDall_3D
00002 use def_matter, only : rhog, utg, uxg, uyg, uzg
00003 use def_metric, only : psi, alph, bvxd, bvyd, bvzd
00004 use coordinate_grav_r, only : rg
00005 use coordinate_grav_theta, only : thg
00006 use coordinate_grav_phi, only : phig
00007 use grid_parameter, only : nrg, ntg, npg
00008 implicit none
00009 integer :: ir, it, ip
00010
00011
00012 open(12,file='ID/rho.dat',status='unknown')
00013 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00014 write(12,'(1p,7e20.12)') rhog(ir,it,ip)
00015 end do; end do; end do
00016 close(12)
00017
00018 open(12,file='ID/ut.dat',status='unknown')
00019 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00020 write(12,'(1p,7e20.12)') utg(ir,it,ip)
00021 end do; end do; end do
00022 close(12)
00023
00024 open(12,file='ID/ux.dat',status='unknown')
00025 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00026 write(12,'(1p,7e20.12)') uxg(ir,it,ip)
00027 end do; end do; end do
00028 close(12)
00029
00030 open(12,file='ID/uy.dat',status='unknown')
00031 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00032 write(12,'(1p,7e20.12)') uyg(ir,it,ip)
00033 end do; end do; end do
00034 close(12)
00035
00036 open(12,file='ID/uz.dat',status='unknown')
00037 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00038 write(12,'(1p,7e20.12)') uzg(ir,it,ip)
00039 end do; end do; end do
00040 close(12)
00041
00042 open(12,file='ID/psi.dat',status='unknown')
00043 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00044 write(12,'(1p,7e20.12)') psi(ir,it,ip)
00045 end do; end do; end do
00046 close(12)
00047
00048 open(12,file='ID/alpha.dat',status='unknown')
00049 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00050 write(12,'(1p,7e20.12)') alph(ir,it,ip)
00051 end do; end do; end do
00052 close(12)
00053
00054 open(12,file='ID/betax.dat',status='unknown')
00055 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00056 write(12,'(1p,7e20.12)') bvxd(ir,it,ip)
00057 end do; end do; end do
00058 close(12)
00059
00060 open(12,file='ID/betay.dat',status='unknown')
00061 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00062 write(12,'(1p,7e20.12)') bvyd(ir,it,ip)
00063 end do; end do; end do
00064 close(12)
00065
00066 open(12,file='ID/betaz.dat',status='unknown')
00067 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00068 write(12,'(1p,7e20.12)') bvzd(ir,it,ip)
00069 end do; end do; end do
00070 close(12)
00071
00072 open(12,file='ID/SphericalCoordinates.dat',status='unknown')
00073 do ip = 0, npg; do it = 0, ntg; do ir = 0, nrg
00074 write(12,'(1p,7e20.12)') rg(ir), thg(it), phig(ip)
00075 end do; end do; end do
00076 close(12)
00077
00078 end subroutine IO_output_solution_IDall_3D
