fix bug in ddz d2dz

This commit is contained in:
ignis 2019-05-01 17:09:47 +09:00
parent d31328671b
commit 3975534030

View file

@ -93,11 +93,11 @@ contains
integer :: i, j ,k
call tp2(zsrc, src, nyp, nzp*nxp)
call tp2(zsrc, src, nyp*nzp, nxp)
do k = 1,nzp
do i = 1,nxp
call dfp(nzp, hzp, zsrc, zdst, nyp, 3)
call dfp(nzp, hzp, zsrc(:,:,i), zdst(:,:,i), nyp, 3)
end do
@ -149,11 +149,11 @@ contains
integer :: i, j ,k
call tp2(zsrc, src, nyp, nzp*nxp)
call tp2(zsrc, src, nyp*nzp, nxp)
do k = 1,nzp
do i = 1,nxp
call d2fp(nzp, hzp, zsrc, zdst, nyp, 3)
call d2fp(nzp, hzp, zsrc(:,:,i), zdst(:,:,i), nyp, 3)
end do