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

write_tp_field modified

Browse source
This commit is contained in:
ignis 2014-07-14 18:19:14 +09:00
parent bb220284a4
commit c6c519555a
1 changed files with 19 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

26
write_stats.f90
View file

@ -19,7 +19,13 @@ subroutine write_tp_field
integer :: i, j, k, ii
if (fdm_sw.eq.0) then
write(tp_field,897) time,nx/spx,ny/spy,nz/spz
if(mod(nx,spx).eq.0) then
write(tp_field,897) time,nx/spx,ny/spy,nz_all/spz
else
write(tp_field,897) time,nx/spx+1,ny/spy+1,nz_all/spz+1
endif
do k=1,nz_all,spz
do j=1,ny,spy
do i=1,nx,spx
@ -27,13 +33,19 @@ subroutine write_tp_field
enddo
enddo
enddo
else
write(tp_field,897) time,nx/svf,ny/svf,nz_all/svf
do ii=1,nz_all,svf
do j=1,ny,svf
do i=1,nx,svf
WRITE(tp_field,'(7e30.20)') REAL(i)*hx,REAL(j)*hx,REAL(ii)*hx, &
else
if(mod(nx,spx).eq.0) then
write(tp_field,897) time,nx/spx,ny/spy,nz_all/spz
else
write(tp_field,897) time,nx/spx+1,ny/spy+1,nz_all/spz+1
endif
do ii=1,nz_all,spz
do j=1,ny,spy
do i=1,nx,spx
write(tp_field,'(7e20.10)') real(i)*hx,real(j)*hx,real(ii)*hx, &
(1.-y1(2,i,j,ii)/y1(1,i,j,ii)),u_(i,j,ii), &
v_(i,j,ii),w_(i,j,ii)
enddo