diff --git a/model/atmosphere/dycore/tests/dycore/stencil_tests/test_compute_theta_rho_face_values_and_pressure_gradient_and_update_vn.py b/model/atmosphere/dycore/tests/dycore/stencil_tests/test_compute_theta_rho_face_values_and_pressure_gradient_and_update_vn.py
index 350d380792..78dcca681f 100644
--- a/model/atmosphere/dycore/tests/dycore/stencil_tests/test_compute_theta_rho_face_values_and_pressure_gradient_and_update_vn.py
+++ b/model/atmosphere/dycore/tests/dycore/stencil_tests/test_compute_theta_rho_face_values_and_pressure_gradient_and_update_vn.py
@@ -122,7 +122,6 @@ def compute_theta_rho_face_value_by_miura_scheme_numpy(
 
 
 @pytest.mark.embedded_remap_error
-@pytest.mark.skip_value_error
 @pytest.mark.uses_as_offset
 class TestComputeThetaRhoPressureGradientAndUpdateVn(test_helpers.StencilTest):
     PROGRAM = compute_theta_rho_face_values_and_pressure_gradient_and_update_vn
@@ -233,7 +232,7 @@ def reference(
         )
 
         # initialize also nest boundary points with zero
-        if limited_area:
+        if limited_area:  # TODO check this
             (rho_at_edges_on_model_levels, theta_v_at_edges_on_model_levels) = np.where(
                 (horz_idx >= start_edge_lateral_boundary) & (horz_idx < end_edge_halo),
                 (