Fix derzvp

feiyizdx · Jan 20, 2022 · 47867c5 · 47867c5
1 parent d7b5f47
commit 47867c5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/derive.f90 b/src/derive.f90
@@ -5038,7 +5038,7 @@ subroutine derzvp(tz,uz,rz,sz,cfz6,csz6,cwz6,nx,ny,nz,nzm,npaire)
         do j=1,ny
            do i=1,nx
               tz(i,j,nzm-1)=aciz6*(uz(i,j,nzm)-uz(i,j,nzm-1))&
-                   +bciz6*(uz(nz,j,k)-uz(nzm-2,j,k))
+                   +bciz6*(uz(i,j,nz)-uz(i,j,nzm-2))
               tz(i,j,nzm)=aciz6*(uz(i,j,nz)-uz(i,j,nzm))&
                    +bciz6*(uz(i,j,nzm)-uz(i,j,nzm-1))
            enddo