diff --git a/src/GRHD/rhs.jl b/src/GRHD/rhs.jl
index 6f005e7432ed0331499e8dc883377c31b272c310..2992e7a5917007e836007036df22a721862aaf41 100644
--- a/src/GRHD/rhs.jl
+++ b/src/GRHD/rhs.jl
@@ -1825,7 +1825,7 @@ function rhs!(mesh, P::Project{:spherical1d}, hrsc::HRSC.AbstractArtificialVisco
   end
   if :Sr ∉ P.prms.freeze_vars
     compute_rhs_weak_form!(rhs_Sr, flr_Sr, src_Sr, nflr_Sr, mesh)
-    compute_rhs_weak_form!(rhs_Sr, flr_Sr, src_Sr, nflr_Sr, mesh)
+    compute_rhs_weak_form!(rhs_Sr, flr_Sr, src_Sr, nflr_Sr, sdg, freeze_mask)
   end
   if :τ ∉ P.prms.freeze_vars
     compute_rhs_weak_form!(rhs_Ï„,  flr_Ï„,  src_Ï„,  nflr_Ï„,  mesh)