research/dns-hit3d-fdm
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

gather4 only performs inverse transform

Browse source
This commit is contained in:
ignis 2014-08-22 13:43:31 +09:00
parent 5843ddccf0
commit 344a5b085a
2 changed files with 4 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
gather_4.f90
View file

@ -45,18 +45,15 @@ subroutine gather_4
! velocities
call xFFT3d(-1,1)
fname = 'u.'//file_ext
tmp4(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,1)
call gather_tmp4(u_)
u_(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,1)
call xFFT3d(-1,2)
fname = 'v.'//file_ext
tmp4(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,2)
call gather_tmp4(v_)
v_(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,2)
call xFFT3d(-1,3)
fname = 'w.'//file_ext
tmp4(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,3)
call gather_tmp4(w_)
w_(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,3)
! scalars
if (int_scalars) then

2
main.f90
View file

@ -130,7 +130,7 @@ program x_code
end if
if (fdm_sw.ne.0) then ! Reacting
!call gather_4
call gather_4
call fdm_exe
endif