00001 subroutine IO_output_solution_fluid_3D
00002 use def_matter, only : emd, rs, rhof, utf, uxf, uyf, uzf
00003 use def_matter_parameter, only : ome, ber, radi
00004 use grid_parameter, only : nrf, ntf, npf
00005 implicit none
00006 integer :: ir, it, ip
00007
00008
00009 open(12,file='rnsflu_all_3D.las',status='unknown')
00010 write(12,'(5i5)') nrf, ntf, npf
00011 do ip = 0, npf
00012 do it = 0, ntf
00013 do ir = 0, nrf
00014 write(12,'(1p,7e20.12)') emd(ir,it,ip), rs(it,ip), rhof(ir,it,ip), &
00015 & utf(ir,it,ip), uxf(ir,it,ip), uyf(ir,it,ip), uzf(ir,it,ip)
00016 end do
00017 end do
00018 end do
00019 write(12,'(1p,6e20.12)') ome, ber, radi
00020 close(12)
00021
00022 end subroutine IO_output_solution_fluid_3D