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

fix bug rho => rho*D

Browse source
This commit is contained in:
ignis 2019-03-25 04:41:23 +09:00
parent 2934fb0e44
commit 8431f9cac8
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
m_fdm_calc.f90
View file

@ -594,7 +594,7 @@ module m_fdm_calc
CALL pdfp (fbuf1, fzzl, fzzu, fbuf2, hy, xx*yy, zz, yy, 3)
CALL pd2fp(fbuf1, fzzl, fzzu, fbuf3, hy, xx*yy, zz, yy, 3)
fbuf1(:,:,:) = r1_(:,:,:,1) ! rho*D
fbuf1(:,:,:) = r1_(:,:,:,1) * dm(:,:,:) ! rho*D
CALL pdfp (fbuf1, fzzl, fzzu, fbuf4, hy, xx*yy, zz, yy, 3)
fbuf2 = fbuf2 * fbuf4 ! dY/dz * d(rho*D)/dz