diff --git a/test/IntegrationTests/refs/modal_bernstein_advection_sine/interpolate.h5 b/test/IntegrationTests/refs/modal_bernstein_advection_sine/interpolate.h5
new file mode 100644
index 0000000000000000000000000000000000000000..bf79b76dcc1db77644baf904aa72eeb8aa845e5f
Binary files /dev/null and b/test/IntegrationTests/refs/modal_bernstein_advection_sine/interpolate.h5 differ
diff --git a/test/IntegrationTests/refs/modal_bernstein_advection_sine/modal_bernstein_advection_sine.toml b/test/IntegrationTests/refs/modal_bernstein_advection_sine/modal_bernstein_advection_sine.toml
new file mode 100644
index 0000000000000000000000000000000000000000..6b41c05f43261490bdd4a4db8702a2671847c0d2
--- /dev/null
+++ b/test/IntegrationTests/refs/modal_bernstein_advection_sine/modal_bernstein_advection_sine.toml
@@ -0,0 +1,24 @@
+[Evolution]
+cfl = 0.5
+tend = 2.2
+
+[ScalarEq]
+analyze_error = false
+equation = "advection"
+id = "sine"
+
+[Output]
+substep_every_iteration = 200
+interpolate_variables = ["u"]
+substep_variables = ["u"]
+interpolate_aligned_ts = [0.3, 0.6, 0.8]
+aligned_ts = [0.05, 0.1, 0.15, 0.2]
+interpolate_nodes = "$(collect(range(0.0,1.0,length=10)))"
+variables1d = ["u"]
+
+[Mesh]
+dg_kind = "modal_bernstein"
+range = [0.0, 1.0]
+k = 40
+basis = "lgl"
+n = 4
diff --git a/test/IntegrationTests/refs/modal_bernstein_advection_sine/output1d.h5 b/test/IntegrationTests/refs/modal_bernstein_advection_sine/output1d.h5
new file mode 100644
index 0000000000000000000000000000000000000000..6515693b408270ee7095a7d995b1b9d80881ac4e
Binary files /dev/null and b/test/IntegrationTests/refs/modal_bernstein_advection_sine/output1d.h5 differ
diff --git a/test/IntegrationTests/refs/modal_bernstein_advection_sine/substep_output.h5 b/test/IntegrationTests/refs/modal_bernstein_advection_sine/substep_output.h5
new file mode 100644
index 0000000000000000000000000000000000000000..f67e822f1758cd4bd521039a5362ae5668e170a4
Binary files /dev/null and b/test/IntegrationTests/refs/modal_bernstein_advection_sine/substep_output.h5 differ
diff --git a/test/IntegrationTests/refs/testconfig.toml b/test/IntegrationTests/refs/testconfig.toml
index 5b254ec5ecda84e0d509d04a29d1ce52af4c0721..b5cdbccf82d11f4257b10818dcd97077dd2bfe0b 100644
--- a/test/IntegrationTests/refs/testconfig.toml
+++ b/test/IntegrationTests/refs/testconfig.toml
@@ -7,6 +7,9 @@ default_abstol1d = 1e-8
 [advection_sine]
 variables1d = [ "u", "u_err" ]
 
+[modal_bernstein_advection_sine]
+variables1d = [ "u" ]
+
 [advection_sine_analysis]
 variables0d = [ "u_err_norm", "u_norm2" ]
 variables1d = [ "u" ]