From 434970bfc8cb0ea9c9729e2a93de6eb925188f3f Mon Sep 17 00:00:00 2001 From: Spencer Bryngelson Date: Fri, 9 Oct 2026 23:35:28 -0400 Subject: [PATCH] Use the ghost-cell spacing for the z one-sided viscous gradient s_compute_fd_gradient's one-sided z derivative at idwbuff(3)%beg divided by z_cc(beg + 2) - z_cc(is3_viscous%beg). is3_viscous%beg is idwbuff(3)%beg + 1 there, so the denominator spanned one cell instead of two and the outermost ghost gradient was doubled on a uniform grid. Divide by z_cc(beg + 2) - z_cc(beg), matching the x and y stencils. The 27 3D Viscous tests pass against the existing goldens. --- src/simulation/m_viscous.fpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/simulation/m_viscous.fpp b/src/simulation/m_viscous.fpp index 37afd8988e..2f294a77b9 100644 --- a/src/simulation/m_viscous.fpp +++ b/src/simulation/m_viscous.fpp @@ -1040,7 +1040,7 @@ contains do j = idwbuff(1)%beg, idwbuff(1)%end grad_z%sf(j, k, idwbuff(3)%beg) = (-3._wp*var%sf(j, k, idwbuff(3)%beg) + 4._wp*var%sf(j, k, & & idwbuff(3)%beg + 1) - var%sf(j, k, & - & idwbuff(3)%beg + 2))/(z_cc(idwbuff(3)%beg + 2) - z_cc(is3_viscous%beg)) + & idwbuff(3)%beg + 2))/(z_cc(idwbuff(3)%beg + 2) - z_cc(idwbuff(3)%beg)) grad_z%sf(j, k, idwbuff(3)%end) = (+3._wp*var%sf(j, k, idwbuff(3)%end) - 4._wp*var%sf(j, k, & & idwbuff(3)%end - 1) + var%sf(j, k, & & idwbuff(3)%end - 2))/(z_cc(idwbuff(3)%end) - z_cc(idwbuff(3)%end - 2))