diff --git a/src/solvers/dg_3d.jl b/src/solvers/dg_3d.jl
index 77fc3db..a9b28a8 100644
--- a/src/solvers/dg_3d.jl
+++ b/src/solvers/dg_3d.jl
@@ -1285,6 +1285,14 @@ function cuda_interface_flux!(mesh::TreeMesh{3}, nonconservative_terms::False, e
     return nothing
 end
 
+# Pack kernels for calculating interface fluxes
+function cuda_interface_flux!(mesh::TreeMesh{3}, nonconservative_terms::True, equations, dg::DGSEM,
+                              cache)
+    surface_flux, nonconservative_flux = dg.surface_integral.surface_flux
+
+    return nothing
+end
+
 # Dummy function for asserting boundaries
 function cuda_prolong2boundaries!(u, mesh::TreeMesh{3},
                                   boundary_condition::BoundaryConditionPeriodic, equations, cache)
diff --git a/test/test_script.jl b/test/test_script.jl
index 187c043..c20dc37 100644
--- a/test/test_script.jl
+++ b/test/test_script.jl
@@ -62,12 +62,12 @@ TrixiGPU.cuda_prolong2interfaces!(u_gpu, mesh_gpu, equations_gpu, cache_gpu)
 Trixi.prolong2interfaces!(cache, u, mesh, equations, solver.surface_integral, solver)
 @test_approx cache_gpu.interfaces.u ≈ cache.interfaces.u
 
-# # Test `cuda_interface_flux!`
-# TrixiGPU.cuda_interface_flux!(mesh_gpu, Trixi.have_nonconservative_terms(equations_gpu),
-#                               equations_gpu, solver_gpu, cache_gpu)
-# Trixi.calc_interface_flux!(cache.elements.surface_flux_values, mesh,
-#                            Trixi.have_nonconservative_terms(equations), equations,
-#                            solver.surface_integral, solver, cache)
+# Test `cuda_interface_flux!`
+TrixiGPU.cuda_interface_flux!(mesh_gpu, Trixi.have_nonconservative_terms(equations_gpu),
+                              equations_gpu, solver_gpu, cache_gpu)
+Trixi.calc_interface_flux!(cache.elements.surface_flux_values, mesh,
+                           Trixi.have_nonconservative_terms(equations), equations,
+                           solver.surface_integral, solver, cache)
 # @test_approx cache_gpu.elements.surface_flux_values ≈ cache.elements.surface_flux_values
 
 # # Test `cuda_prolong2boundaries!`