From 1fbfd6ef0334dad81828845c29525cb9da81fb90 Mon Sep 17 00:00:00 2001 From: Brandon Allen Date: Tue, 30 Jun 2020 10:47:42 -0400 Subject: [PATCH] account for ghost elements Co-authored-by: Jeremy E Kozdon --- .../HydrostaticBoussinesq/HydrostaticBoussinesqModel.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/Ocean/HydrostaticBoussinesq/HydrostaticBoussinesqModel.jl b/src/Ocean/HydrostaticBoussinesq/HydrostaticBoussinesqModel.jl index 0c45d81b187..921f0b1a72e 100644 --- a/src/Ocean/HydrostaticBoussinesq/HydrostaticBoussinesqModel.jl +++ b/src/Ocean/HydrostaticBoussinesq/HydrostaticBoussinesqModel.jl @@ -640,10 +640,10 @@ function update_auxiliary_state_gradient!( wz0 = A.data[:, index_wz0, :] # project w(z=0) down the stack - Nq, Nqk, _, _, nelemv, _, nelemh, _ = basic_grid_info(dg) + Nq, Nqk, _, _, nelemv, nelemh, nhorzrealelem, _ = basic_grid_info(dg) data = reshape(A.data, Nq^2, Nqk, number_auxiliary, nelemv, nelemh) - flat_wz0 = @view data[:, end:end, index_w, end:end, :] - boxy_wz0 = @view data[:, :, index_wz0, :, :] + flat_wz0 = @view data[:, end:end, index_w, end:end, 1:nhorzrealelem] + boxy_wz0 = @view data[:, :, index_wz0, :, 1:nhorzrealelem] boxy_wz0 .= flat_wz0 return true